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 )

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)

Definition at line 2619 of file module_bl_mynn.F90.

References mym_length(), and mym_level2().

Referenced by mynn_bl_driver().

Here is the call graph for this function:

Here is the caller graph for this function: