CAM5 / vertical diffusion

[andreasb]

Hi,In the calls to qsat in vertical_diffusion.F90, return-values are tem2 (Saturation specific humidity) and ftem2 (Saturation vapor pressure):
call qsat(state%t(:ncol,:), state%pmid(:ncol,:), tem2(:ncol,:), ftem(:ncol,:))
However, in the subroutine qsat (wv_saturation.F90), specific humidity and vapor pressure seem to be the other way round:
subroutine qsat(t, p, es, qs)
  real(r8), intent(out) :: es  ! Saturation vapor pressure
  real(r8), intent(out) :: qs  ! Saturation specific humidity
Am I misunderstanding this or might this be a problem?This is model version CESM1.2.0.thanks!
Andreas

 

[eaton]

Looks like the problem is incorrect comments in the declarations of ftem and tem2.  tem2 is not used, and ftem is used in a context that requires a specific humidity, i.e., as the divisor of the state specific humidity to compute a relative humidity.

 

[andreasb]

excellent. Thanks for your checking.