Common Community Physics Package (CCPP) Scientific Documentation  Version 2.0
subroutine, public module_radiation_clouds::progcld3 ( real (kind=kind_phys), dimension(:,:), intent(in)  plyr,
real (kind=kind_phys), dimension(:,:), intent(in)  plvl,
real (kind=kind_phys), dimension(:,:), intent(in)  tlyr,
real (kind=kind_phys), dimension(:,:), intent(in)  tvly,
real (kind=kind_phys), dimension(:,:), intent(in)  qlyr,
real (kind=kind_phys), dimension(:,:), intent(in)  qstl,
real (kind=kind_phys), dimension(:,:), intent(in)  rhly,
real (kind=kind_phys), dimension(:,:), intent(in)  clw,
real (kind=kind_phys), dimension(:,:)  cnvw,
real (kind=kind_phys), dimension(:,:)  cnvc,
real (kind=kind_phys), dimension(:), intent(in)  xlat,
real (kind=kind_phys), dimension(:), intent(in)  xlon,
real (kind=kind_phys), dimension(:), intent(in)  slmsk,
integer, intent(in)  ix,
integer, intent(in)  nlay,
integer, intent(in)  nlp1,
real (kind=kind_phys), dimension(:,:)  deltaq,
real (kind=kind_phys), intent(in)  sup,
integer, intent(in)  kdt,
integer  me,
real (kind=kind_phys), dimension(:,:,:), intent(out)  clouds,
real (kind=kind_phys), dimension(:,:), intent(out)  clds,
integer, dimension(:,:), intent(out)  mtop,
integer, dimension(:,:), intent(out)  mbot 
)
Parameters
plyr(ix,nlay), model layer mean pressure in mb (100pa)
plvl(ix,nlp1), model level pressure in mb (100pa)
tlyr(ix,nlay), model layer mean temperature in K
tvly(ix,nlay), model layer virtual temperature in K
qlyr(ix,nlay), layer specific humidity in gm/gm
qstl(ix,nlay), layer saturate humidity in gm/gm
rhly(ix,nlay), layer relative humidity (=qlyr/qstl)
clw(ix,nlay), layer cloud condensate amount
cnvw(ix,nlay), layer convective cloud condensate
cnvc(ix,nlay), layer convective cloud cover
xlat(ix), grid latitude in radians, default to pi/2 -> -pi/2 range, otherwise see in-line comment
xlon(ix), grid longitude in radians (not used)
slmsk(ix), sea/land mask array (sea:0,land:1,sea-ice:2)
ixhorizontal dimention
nlayvertical layer dimension
nlp1vertical level dimension
deltaq(ix,nlay), half total water distribution width
supsupersaturation
kdtcurrent forecast iteration
meprint control flag
clouds(ix,nlay,nf_clds), cloud profiles
(:,:,1) - layer total cloud fraction
(:,:,2) - layer cloud liq water path (g/m**2)
(:,:,3) - mean eff radius for liq cloud (micron)
(:,:,4) - layer cloud ice water path (g/m**2)
(:,:,5) - mean eff radius for ice cloud (micron)
(:,:,6) - layer rain drop water path not assigned
(:,:,7) - mean eff radius for rain drop (micron)
(:,:,8) - layer snow flake water path not assigned
(:,:,9) - mean eff radius for snow flake(micron)
clds(ix,5), fraction of clouds for low, mid, hi, tot, bl
mtop(ix,3), vertical indices for low, mid, hi cloud tops
mbot(ix,3), vertical indices for low, mid, hi cloud bases

progcld3 General Algorithm

  1. Find top pressure (ptopc) for each cloud domain for given latitude.
  2. Calculate liquid/ice condensate path in \( g/m^2 \)
  3. Calculate effective liquid cloud droplet radius over land.
  4. Calculate layer cloud fraction.
  5. Calculate effective ice cloud droplet radius following Heymsfield and McFarquhar (1996) [44].
  6. Call gethml() to compute low,mid,high,total, and boundary layer cloud fractions and clouds top/bottom layer indices for low, mid, and high clouds.

References physcons::con_pi, physcons::con_thgni, physcons::con_ttp, gethml(), physparam::ivflip, physparam::lcnorm, and physparam::lcrick.

Here is the call graph for this function: