icloud()

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) [119] eq. 21 solved implicitly; threshold from wsm6 scheme, Hong et al. (2004) [94], 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.

Definition at line 1544 of file module_gfdl_cloud_microphys.F90.

References acr3d(), linear_prof(), smlt(), and subgrid_z_proc().

Referenced by mpdrv().

Here is the call graph for this function:

Here is the caller graph for this function: