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)