subroutine gfdl_cloud_microphys_mod::icloud ( integer, intent(in)  ktop, integer, intent(in)  kbot, real, dimension (ktop:kbot), intent(inout)  tzk, real, dimension (ktop:kbot), intent(in)  p1, real, dimension (ktop:kbot), intent(inout)  qvk, real, dimension (ktop:kbot), intent(inout)  qlk, real, dimension (ktop:kbot), intent(inout)  qrk, real, dimension (ktop:kbot), intent(inout)  qik, real, dimension (ktop:kbot), intent(inout)  qsk, real, dimension (ktop:kbot), intent(inout)  qgk, real, dimension (ktop:kbot), intent(in)  dp1, real, dimension (ktop:kbot), intent(in)  den, real, dimension (ktop:kbot), intent(in)  denfac, real, dimension (ktop:kbot), intent(in)  vts, real, dimension (ktop:kbot), intent(in)  vtg, real, dimension (ktop:kbot), intent(in)  vtr, real, dimension (ktop:kbot), intent(inout)  qak, real, intent(in)  rh_adj, real, intent(in)  rh_rain, real, intent(in)  dts, real, intent(in)  h_var  )

private

Author

Shian-Jiann Lin, GFDL

This scheme is featured with:

• bulk cloud microphysics
• processes splitting with some un-split sub-grouping
• time implicit (when possible) accretion and autoconversion

GFDL icloud Detailed Algorithm

• Define conversion scalar/factor.
• Define heat capacity and latend heat coefficient.
• Calculate $P_{imlt}$: instant melting of cloud ice.
• Calculate $P_{ihom}$: homogeneous freezing of cloud water into cloud ice. This is the 1st occurance of liquid water freezing in the split MP process.
• Call linear_prof() to calculate vertical subgrid variability of cloud ice.
• Update capacity heat and latend heat coefficient.
• Melting of snow:
  • $P_{sacw}$: accretion of cloud water by snow
  • $P_{sacr}$: accretion of rain by melted snow
  • $P_{racs}$: accretion of snow by rain
  • Total snow sink: $P_{smlt}$: snow melt (smlt(); due to rain accretion)
• Update capacity heat and latend heat coefficient.
• Melting of graupel:
  • $P_{gacr}$: accretion of rain by graupel
  • $P_{gacw}$: accretion of cloud water by graupel
  • $P_{gmlt}$: graupel melt (gmlt())
• Cloud ice processes:
  • $P_{saci}$: accretion of cloud ice by snow
  • $P_{saut}$: autoconversion: cloud ice $\rightarrow$ snow.
    similar to Lin et al.(1983) [107] eq. 21 solved implicitly; threshold from wsm6 scheme, Hong et al. (2004) [88], eq (13) : qi0_crt ~0.8e-4.
  • $P_{gaci}$: accretion of cloud ice by graupel
• Cold-rain processes:
  • $P_{sacr}$: accretion of rain by snow (acr3d())
  • $P_{gfr}$: rain freezing $\rightarrow$ graupel
  • Calculate total sink to $q_r$ :
    Sink terms to $q_r$: $P_{sacr}+P_{gfr}$
    source term to $q_s$: $P_{sacr}$
    source term to $q_g$: $P_{gfr}$
• Update capacity heat and latend heat coefficient.
• Graupel production terms:
  • accretion: snow $\rightarrow$ graupel (acr3d())
  • autoconversion: snow $\rightarrow$ graupel
  • $P_{gacw}$: accretion of cloud water by graupel
  • $P_{gacr}$: accretion of rain by graupel (acr3d())
• Call subgrid_z_proc() for subgrid cloud microphysics.

References acco, acr3d(), c_air, c_ice, c_liq, c_vap, ces0, cgaci, cgacr, cgacs, cgacw, cgfr, cgmlt, const_vi, cracs, csaci, csacr, csacw, csmlt, d0_vap, dc_ice, dt_fr, gmlt(), li00, linear_prof(), lv00, p_min, qcmin, qi0_crt, qi_gen, qi_lim, ql_mlt, qrmin, qs0_crt, qs_mlt, smlt(), subgrid_z_proc(), t_wfr, tau_g2v, tau_i2s, tau_imlt, tau_v2g, tice, tice0, and z_slope_ice.

Referenced by mpdrv().

Here is the call graph for this function: