subroutine gfdl_cloud_microphys::cs_profile ( real, dimension (4, km), intent(inout)  a4, real, dimension (km), intent(in)  del, integer, intent(in)  km, logical, intent(in)  do_mono  )

private

GFDL Cloud mpdrv General Algorithm

• Prevent excessive build-up of cloud ice from external sources.
• Convert moist mixing ratios to dry mixing ratios.
• Calculate cloud condensation nuclei (ccn), following klein eq. 15
• Calculate horizontal subgrid variability, which is used in cloud fraction, relative humidity calculation, evaporation and condensation processes. Horizontal sub-grid variability is a function of cell area and land/sea mask:
  Over land:
  $$t_{land}=dw_{land}(\frac{A_{r}}{10^{10}})^{0.25}$$
  
  Over ocean:
  $$t_{ocean}=dw_{ocean}(\frac{A_{r}}{10^{10}})^{0.25}$$
  where $A_{r}$ is cell area. $dw_{land}=0.16$ and $dw_{ocean}=0.10$ are base value for sub-grid variability over land and ocean. The total horizontal sub-grid variability is:
  $$h_{var}=t_{land}\times fr_{land}+t_{ocean}\times (1-fr_{land})$$
  $$h_{var}=min[0.2,max(0.01,h_{var})]$$

• Calculate relative humidity increment.
• If requested, call neg_adj() and fix all negative water species.
• Do loop on cloud microphysics sub time step.
  • Define air density based on hydrostatical property.
  • Call warm_rain() - time-split warm rain processes: 1st pass.
  • Sedimentation of cloud ice, snow, and graupel.
    • Call fall_speed() to calculate the fall velocity of cloud ice, snow and graupel.
    • Call terminal_fall() to calculate the terminal fall speed.
  • Call sedi_heat() to calculate heat transportation during sedimentation.
  • Call warm_rain() to - time-split warm rain processes: 2nd pass
  • Call icloud(): ice-phase microphysics
• Calculate momentum transportation during sedimentation.
• Update moist air mass (actually hydrostatic pressure).
• Update cloud fraction tendency.

References cs_limiters().

Referenced by lagrangian_fall_ppm().

Here is the call graph for this function: