Bug report on variables QRS, QRL

[simchan_yook@colostate_edu]

Hi all,

I've looked for the temperature budget variables.

I've checked that
PTTEND is not equal to the sum of (QRS+QRL+DTCOND+DTV+TTGWORO)
because
QRS + QRL is not equal to the temperature tendencies calculated by the subroutine radiation_tend

I figured out that this is because
The actual unit of QRS, QRL is not K/s but K*(Pa)/s.

Temperature budget (state%dtdt) is calculated using qrs(i,k)/state%pdel(i,k), qrl(i,k)/state%pdel(i,k)

! convert radiative heating rates from Q*dp to Q for energy conservation
if (conserve_energy) then
!DIR$ CONCURRENT
do k =1 , pver
!DIR$ CONCURRENT
do i = 1, ncol
qrs(i,k) = qrs(i,k)/state%pdel(i,k)
qrl(i,k) = qrl(i,k)/state%pdel(i,k)
end do
end do
end if

While QRS and QRL variables are written without divided by state%pdel(i,k).

As a results
PTTEND budget is not same as the sum of QRS+QRL+DTV+DTCOND+TTGWORO

*** Maybe this is not a sole factor causing the difference, even though I divide the qrs and qrl by pdel some difference still existed.
If I get the QRS and QRL from the radheat_tend subroutine then the difference has been fixed.



Thank you,
Shim

 

[hannay]

See: https://xenforo.cgd.ucar.edu/cesm/threads/cesm1-2-2-pttend-qrs-qrl-dtcond-dtv-ttgworo-tfix.4886/