CCPP Scientific Documentation
v5.0.0
subroutine module_bl_mynn::dmp_mf ( integer, intent(in)  kts,
integer, intent(in)  kte,
real, intent(in)  dt,
real, dimension(kts:kte+1), intent(in)  zw,
real, dimension(kts:kte), intent(in)  dz,
real, dimension(kts:kte), intent(in)  p,
integer, intent(in)  momentum_opt,
integer, intent(in)  tke_opt,
integer, intent(in)  scalar_opt,
real, dimension(kts:kte), intent(in)  u,
real, dimension(kts:kte), intent(in)  v,
real, dimension(kts:kte), intent(in)  w,
real, dimension(kts:kte), intent(in)  th,
real, dimension(kts:kte), intent(in)  thl,
real, dimension(kts:kte), intent(in)  thv,
real, dimension(kts:kte), intent(in)  tk,
real, dimension(kts:kte), intent(in)  qt,
real, dimension(kts:kte), intent(in)  qv,
real, dimension(kts:kte), intent(in)  qc,
real, dimension(kts:kte), intent(in)  qke,
real, dimension(kts:kte), intent(in)  qnc,
real, dimension(kts:kte), intent(in)  qni,
real, dimension(kts:kte), intent(in)  qnwfa,
real, dimension(kts:kte), intent(in)  qnifa,
real, dimension(kts:kte), intent(in)  exner,
real, dimension(kts:kte), intent(inout)  vt,
real, dimension(kts:kte), intent(inout)  vq,
real, dimension(kts:kte), intent(inout)  sgm,
real, intent(in)  ust,
real, intent(in)  flt,
real, intent(in)  flq,
real, intent(in)  flqv,
real, intent(in)  flqc,
real, intent(in)  pblh,
integer, intent(in)  kpbl,
real, intent(in)  DX,
real, intent(in)  landsea,
real, intent(in)  ts 
)

dmp_mf() calculates the nonlocal turbulent transport from the dynamic multiplume mass-flux scheme as well as the shallow-cumulus component of the subgrid clouds. Note that this mass-flux scheme is called when the namelist paramter bl_mynn_edmf is set to 1 (recommended).

Much thanks to Kay Suslj of NASA-JPL for contributing the original version of this mass-flux scheme. Considerable changes have been made from it's original form. Some additions include:

  1. scale-aware tapering as dx -> 0
  2. transport of TKE (extra namelist option)
  3. Chaboureau-Bechtold cloud fraction & coupling to radiation (when icloud_bl > 0)
  4. some extra limits for numerical stability

This scheme remains under development, so consider it experimental code.

References bl_mynn_mixchem, condensation_edmf(), debug_code, env_subs, ep_3, esat_blend(), gtr, onethird, p1000mb, qsat_blend(), svp1, tv0, xl_blend(), and xlvcp.

Referenced by mynn_bl_driver().

Here is the call graph for this function: