CCPP Scientific Documentation
v4.1.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,
real (kind=kind_phys), dimension(:,:), intent(in)  dz,
real (kind=kind_phys), dimension(:,:), intent(in)  delp,
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,
real (kind=kind_phys), dimension(:), intent(out)  de_lgth 
)
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)
dz(IX,NLAY), layer thickness (km)
delp(IX,NLAY), model layer pressure thickness in mb (100Pa)
ixhorizontal dimention
nlay,nlp1vertical layer/level dimensions
deltaq(ix,nlay), half total water distribution width
supsupersaturation
kdt
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
de_lgth(ix), clouds decorrelation length (km)

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) [80].
  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 climit, climit2, physcons::con_pi, physcons::con_thgni, physcons::con_ttp, gethml(), gfac, gord, iovr, physparam::lcnorm, physparam::lcrick, nf_clds, ptopc, reice_def, reliq_def, rrain_def, and rsnow_def.

Here is the call graph for this function: