Index | Rundeck: P1SoM40 | Created: Thu May 9 03:30:17 EDT 2024 |
CLOUDS | File: CLOUDS2.F90 |
Author : M.S.Yao/A. Del Genio (modifications by Gavin Schmidt)
Version:
Subroutines:
- anvil_optical_thickness
- MC_OPTICAL_THICKNESS calculates anvil cloud optical thickness and droplet radius
- convective_microphysics
- CONVECTIVE_MICROPHYSICS calculates convective microphysics and precip condensate
- lscond
- LSCOND column physics of large scale condensation
- mass_flux
- mc_cloud_fraction
- MC_CLOUD_FRACTION calculates convective cloud fraction
- mc_precip_phase
- MC_PRECIP_PHASE calculates the convective precip phase
- mstcnv
- MSTCNV moist convective processes (precip, convective clouds,...)
- precip_mp
Global Variables:
- rgas : used from constant
- grav : used from constant
- lhe : used from constant
- lhs : used from constant
- lhm : used from constant
- sha : used from constant
- bysha : used from constant
- pi : used from constant
- by6 : used from constant
- lm : used from resolution
- dtsrc : used from model_com
- itime : used from model_com
- seconds_per_hour : used from timeconstants_mod
- ncol : used from clouds_com
- nmom : used from qusdef
- xymoms : used from qusdef
- zmoms : used from qusdef
- zdir : used from qusdef
- ti :
real*8, parameter
- pure ice limit
Initial Value= 233.16d0
- cldmin :
real*8, parameter
- min MC/LSC region
Initial Value= .10d0
- wmu :
real*8, parameter
- critical cloud water content for rapid conversion (g m**-3)
Initial Value= .25
- wmul :
real*8, parameter
- WMU over land
Initial Value= .5
- wmui :
real*8
- WMU for ice clouds
- wmui :
real*8
- WMU for ice clouds
- wmusi :
real*8
- WMU for liquid clouds over sea-ice
- brcld :
real*8, parameter
- for cal. BYBR
Initial Value= .2d0
Used by: | :init_cld |
- fddet :
real*8, parameter
- remainder of downdraft
Initial Value= .25d0
- dtmin1 :
real*8, parameter
- min DT to stop downdraft drop
Initial Value= 1.d0
- slhe :
real*8, parameter
- Initial Value
= LHE*BYSHA
- ccmul :
real*8, parameter
- multiplier for convective cloud cover
Initial Value= 2.
- ccmul1 :
real*8, parameter
- multiplier for deep anvil cloud cover
Initial Value= 5.
- ccmul2 :
real*8, parameter
- multiplier for shallow anvil cloud cover
Initial Value= 3.
- coetau :
real*8, parameter
- multiplier for convective cloud optical thickness
Initial Value= .08d0
- slhs :
real*8, parameter
- Initial Value
= LHS*BYSHA
- rtemp :
real*8
- cmx :
real*8
- wmuix :
real*8
- contce1 :
real*8
- contce2 :
real*8
- tnx :
real*8
- qnx :
real*8
- bybr :
real*8
- factor for converting cloud particle radius to effect. radius
Used by: | :init_cld |
- xmass :
real*8
- dummy variable
Used by: | :init_cld |
- pland :
real*8
- land fraction
- bydtsrc :
real*8
Used by: | :condse | :init_cld |
- rcldlx :
real*8
- multipliers for ice particle/droplet radius
- rcldix :
real*8
- multipliers for ice particle/droplet radius
- lmcm :
integer
- max level for originating MC plumes
Initial Value= -1
Used by: | :init_cld |
- isc :
integer
- integer to turn on computation of stratocumulus clouds
Initial Value= 0
- u00wtrx :
real*8
- multiplies U00ice for critical humidity for water clds
Initial Value= 1.0d0
Used by: | :init_cld |
- u00ice :
real*8
- critical humidity for ice cloud condensation
Initial Value= .7d0
Used by: | :init_cld |
- u00a :
real*8
- tuning knob for U00 above 850 mb without moist convection
Initial Value= 0.55d0
- u00b :
real*8
- tuning knob for U00 below 850 mb and in convective regions
Initial Value= 1.00d0
- maxctop :
real*8
- max cloud top pressure
Initial Value= 50.d0
- funio_denominator :
real*8
- funio denominator
Initial Value= 22.d0
- autoconv_multiplier :
real*8
- autoconversion rate multiplier
Initial Value= 1.d0
- radiusl_multiplier :
real*8
- liquid cloud particle radius multiplier
Initial Value= 1.d0
- radiusi_multiplier :
real*8
- ice cloud particle radius multiplier
Initial Value= 1.d0
- wmu_multiplier :
real*8
- critical liquid cloud water multiplier
Initial Value= 1.d0
- wmui_multiplier :
real*8
- critical ice cloud water multiplier
Initial Value= 1.d0
- entrainment_cont1 :
real*8
- constant for entrainment rate, plume 1
Initial Value= .3d0
- entrainment_cont2 :
real*8
- constant for entrainment rate, plume 2
Initial Value= .6d0
- hrmax :
real*8
- maximum distance an air parcel rises from surface
Initial Value= 1000.d0
- rimax :
real*8
- maximum ice cloud size
Initial Value= 100.d0
- rwmax :
real*8
- maximum water cloud size
Initial Value= 20.d0
- rwcldox :
real*8
- multiplies part.size of water clouds over ocean
Initial Value= 1.d0
- ricldx :
real*8
- multiplies part.size of ice clouds at 1000mb
Initial Value= 1.d0
- do_blu00 :
integer
- =1 if boundary layer U00 is treated differently
Initial Value= 0
- xricld :
real*8
- tautab :
real*8, dimension(0:255)
- look-up table to convert count value to optical thickness
Used by: | :isccp_cloud_types |
- invtau :
integer, dimension(-20:45000)
- look-up table to convert optical thickness to count value
Used by: | :isccp_cloud_types |
- ra :
real*8, dimension(:), allocatable
- ratio of primary grid box to secondary gridbox
- debug :
logical
- v_0 :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- u_0 :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- vm1 :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- um1 :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- vm :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- um :
real*8, dimension(:, :), allocatable
- velocity related variables(UM,VM)=(U,V)*AIRM
- miscellaneous :
- vertical arrays set in driver
- ple :
real*8, dimension(lm+1)
- pressure at layer edge
- lhp :
real*8, dimension(lm+1)
- array of precip phase ! may differ from LHX
- pl :
real*8, dimension(lm)
- layer pressure (mb)
- plk :
real*8, dimension(lm)
- PL**KAPA
- airm :
real*8, dimension(lm)
- the layer's pressure depth (mb)
- byam :
real*8, dimension(lm)
- 1./AIRM
- etal :
real*8, dimension(lm)
- fractional entrainment rate
- tl :
real*8, dimension(lm)
- temperature, specific humidity of the layer
- ql :
real*8, dimension(lm)
- temperature, specific humidity of the layer
- th :
real*8, dimension(lm)
- potential temperature (K)
- rh :
real*8, dimension(lm)
- relative humidity
- qclx :
real*8, dimension(lm)
- liquid cloud water mixing ratio (kg/kg) ! WMX
- qcix :
real*8, dimension(lm)
- ice cloud water mixing ratio (kg/kg)
- dqmtotal :
real*8, dimension(lm)
- dqmshlw :
real*8, dimension(lm)
- dqmdeep :
real*8, dimension(lm)
- pl :
real*8, dimension(lm)
- layer pressure (mb)
- plk :
real*8, dimension(lm)
- PL**KAPA
- airm :
real*8, dimension(lm)
- the layer's pressure depth (mb)
- byam :
real*8, dimension(lm)
- 1./AIRM
- etal :
real*8, dimension(lm)
- fractional entrainment rate
- ql :
real*8, dimension(lm)
- temperature, specific humidity of the layer
- tl :
real*8, dimension(lm)
- temperature, specific humidity of the layer
- th :
real*8, dimension(lm)
- potential temperature (K)
- rh :
real*8, dimension(lm)
- relative humidity
- rh1 :
- relative humidity to compare with the threshold humidity
- qclx :
real*8, dimension(lm)
- liquid cloud water mixing ratio (kg/kg) ! WMX
- qcix :
real*8, dimension(lm)
- ice cloud water mixing ratio (kg/kg)
- vsubl :
- downward vertical velocity due to cumulus subsidence (cm/s)
- dgdqm :
- dummy variables
- dqcond :
- dummy variables
- dphase :
- dummy variables
- dgdsm :
- dummy variables
- mcflx :
- dummy variables
- ddmflx :
- accumulated downdraft mass flux (mb)
- aq :
- time change rate of specific humidity (s**-1)
- dpdt :
- time change rate of pressure (mb/s)
- fssl :
- grid fraction for large-scale clouds
- vlat :
- dummy variable
- wturb :
- turbulent vertical velocity (m)
- precnvl :
real*8, dimension(lm+1)
- convective precip entering the layer top
- sdl :
real*8, dimension(lm)
- vertical velocity in sigma coordinate
- qcll :
real*8, dimension(lm)
- cloud liquid water mixing ratio (kg/kg)
- qcil :
real*8, dimension(lm)
- cloud ice water mixing ratio (kg/kg)
- taumcl :
real*8, dimension(lm)
- convective cloud optical thickness
- svlatl :
real*8, dimension(lm)
- saved LHX for convective cloud
- condpt :
real*8, dimension(lm)
- temporary variable to save CONDP
- cldmcl :
real*8, dimension(lm)
- convective cloud cover
- svlhxl :
real*8, dimension(lm)
- saved LHX for large-scale cloud
- svwmxl :
real*8, dimension(lm)
- saved detrained convective cloud water
- svlat1 :
real*8, dimension(lm)
- csizel :
real*8, dimension(lm)
- cloud particle radius (micron)
- cnvmmrl :
real*8, dimension(lm)
- convective cloud condensate mass mixing ratio
Used by: | :condse |
- cldsal :
real*8, dimension(lm)
- Fractional stratiform cloud area used to derive in-cloud MMR
Used by: | :condse |
- ttoldl :
real*8, dimension(lm)
- previous potential temperature
- cldsavl :
real*8, dimension(lm)
- saved large-scale cloud cover
- cldsv1 :
real*8, dimension(lm)
- dctei :
real*8, dimension(lm)
- height diagnostics of dry and latent heating by MC
- sshr :
real*8, dimension(lm)
- height diagnostics of dry and latent heating by MC
- taussl :
real*8, dimension(lm)
- large-scale cloud optical thickness
- cldssl :
real*8, dimension(lm)
- large-scale cloud cover
- qm :
real*8, dimension(lm)
- Vertical profiles of (T/p**kappa)*AIRM, q*AIRM
- sm :
real*8, dimension(lm)
- Vertical profiles of (T/p**kappa)*AIRM, q*AIRM
- smom :
real*8, dimension(nmom, lm)
- qmom :
real*8, dimension(nmom, lm)
- smommc :
real*8, dimension(nmom, lm)
- qmommc :
real*8, dimension(nmom, lm)
- kmax :
integer
- index for surrounding velocity
- lmcld :
integer
- max cloud top level
- pearth :
real*8
- fraction of land in grid box
- ts :
real*8
- average surface temperture (C)
- ri2 :
real*8
- Richardson numbers
- ri1 :
real*8
- Richardson numbers
- ris :
real*8
- Richardson numbers
- dcl :
integer
- max level of planetary boundary layer
- qs :
real*8
- us :
real*8
- vs :
real*8
- dxypij :
real*8
- roice :
real*8
- prcpmc :
real*8
- precip due to moist convection
- prcpss :
real*8
- precip due to large-scale condensation
- hcndss :
real*8
- heating due to large-scale condensation
- wmsum :
real*8
- cloud liquid water path (kg/m2)
- wisum :
real*8
- cloud ice water path (kg/m2)
- cldslwij :
real*8
- shallow convective cloud cover
- clddepij :
real*8
- deep convective cloud cover
- lmcmax :
integer
- upper-most convective layer
- lmcmin :
integer
- lowerest convective layer
- airxl :
real*8
- is convective mass flux (mb)
- rndssl :
real*8, dimension(3, LM)
- stored random number sequences
- prheat :
real*8
- prebar1 :
real*8, dimension(Lm+1)
- copy of variable prebar
- qiss :
real*8, dimension(Lm)
- stratiform liquid,ice water (cloud+in-cloud precip) available to radiation (kg/kg)
- qlss :
real*8, dimension(Lm)
- stratiform liquid,ice water (cloud+in-cloud precip) available to radiation (kg/kg)
- qimc :
real*8, dimension(Lm)
- convective liquid,ice water (cloud+in-cloud precip) available to radiation (kg/kg)
- qlmc :
real*8, dimension(Lm)
- convective liquid,ice water (cloud+in-cloud precip) available to radiation (kg/kg)