CCPP SciDoc  v6.0.0
Common Community Physics Package Developed at DTC
 
Loading...
Searching...
No Matches

◆ mym_turbulence()

subroutine module_bl_mynn::mym_turbulence ( integer, intent(in)  kts,
integer, intent(in)  kte,
real(kind_phys), intent(in)  xland,
real(kind_phys), intent(in)  closure,
real(kind_phys), dimension(kts:kte), intent(in)  dz,
real(kind_phys), intent(in)  dx,
real(kind_phys), dimension(kts:kte+1), intent(in)  zw,
real(kind_phys), dimension(kts:kte), intent(in)  u,
real(kind_phys), dimension(kts:kte), intent(in)  v,
real(kind_phys), dimension(kts:kte), intent(in)  thl,
real(kind_phys), dimension(kts:kte), intent(in)  thetav,
real(kind_phys), dimension(kts:kte), intent(in)  ql,
real(kind_phys), dimension(kts:kte), intent(in)  qw,
real(kind_phys), dimension(kts:kte), intent(in)  qke,
real(kind_phys), dimension(kts:kte), intent(in)  tsq,
real(kind_phys), dimension(kts:kte), intent(in)  qsq,
real(kind_phys), dimension(kts:kte), intent(in)  cov,
real(kind_phys), dimension(kts:kte), intent(in)  vt,
real(kind_phys), dimension(kts:kte), intent(in)  vq,
real(kind_phys), intent(in)  rmo,
real(kind_phys), intent(in)  flt,
real(kind_phys), intent(in)  fltv,
real(kind_phys), intent(in)  flq,
real(kind_phys), intent(in)  zi,
real(kind_phys), dimension(kts:kte), intent(in)  theta,
real(kind_phys), dimension(kts:kte)  sh,
real(kind_phys), dimension(kts:kte)  sm,
real(kind_phys), dimension(kts:kte), intent(out)  el,
real(kind_phys), dimension(kts:kte), intent(out)  dfm,
real(kind_phys), dimension(kts:kte), intent(out)  dfh,
real(kind_phys), dimension(kts:kte), intent(out)  dfq,
real(kind_phys), dimension(kts:kte), intent(out)  tcd,
real(kind_phys), dimension(kts:kte), intent(out)  qcd,
real(kind_phys), dimension(kts:kte), intent(out)  pdk,
real(kind_phys), dimension(kts:kte), intent(out)  pdt,
real(kind_phys), dimension(kts:kte), intent(out)  pdq,
real(kind_phys), dimension(kts:kte), intent(out)  pdc,
real(kind_phys), dimension(kts:kte), intent(inout)  qwt1d,
real(kind_phys), dimension(kts:kte), intent(inout)  qshear1d,
real(kind_phys), dimension(kts:kte), intent(inout)  qbuoy1d,
real(kind_phys), dimension(kts:kte), intent(inout)  qdiss1d,
integer, intent(in)  tke_budget,
real(kind_phys), intent(in)  psig_bl,
real(kind_phys), intent(in)  psig_shcu,
real(kind_phys), dimension(kts:kte), intent(in)  cldfra_bl1d,
integer, intent(in)  bl_mynn_mixlength,
real(kind_phys), dimension(kts:kte), intent(in)  edmf_w1,
real(kind_phys), dimension(kts:kte), intent(in)  edmf_a1,
real(kind_phys), dimension(kts:kte), intent(in)  tkeprodtd,
integer, intent(in)  spp_pbl,
real(kind_phys), dimension(kts:kte)  rstoch_col 
)

GSD mym_turbulence General Algorithm

Two subroutines mym_level2() and mym_length() are called within this subrouine to collect variable to carry out successive calculations:

  • mym_level2() calculates the level 2 nondimensional wind shear \(G_M\) and vertical temperature gradient \(G_H\) as well as the level 2 stability functions \(S_h\) and \(S_m\).
  • mym_length() calculates the mixing lengths.
  • The stability criteria from Helfand and Labraga (1989) are applied.
  • The stability functions for level 2.5 or level 3.0 are calculated.
  • If level 3.0 is used, counter-gradient terms are calculated.
  • Production terms of TKE, \(\theta^{'2}\), \(q^{'2}\), and \(\theta^{'}q^{'}\) are calculated.
  • Eddy diffusivity \(K_h\) and eddy viscosity \(K_m\) are calculated.
  • TKE budget terms are calculated (if the namelist parameter tke_budget is set to True)