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)