This module contains the entity of MYNN-EDMF PBL scheme. More...

Functions/Subroutines
subroutine	mym_level2 (kts, kte, dz, u, v, thl, thetav, qw, ql, vt, vq, dtl, dqw, dtv, gm, gh, sm, sh)
	This subroutine calculates the level 2, non-dimensional wind shear \(G_M\) and vertical temperature gradient \(G_H\) as well as the level 2 stability funcitons \(S_h\) and \(S_m\).

subroutine	mym_length (kts, kte, xland, dz, dx, zw, rmo, flt, fltv, flq, vt, vq, u1, v1, qke, dtv, el, zi, theta, qkw, psig_bl, cldfra_bl1d, bl_mynn_mixlength, edmf_w1, edmf_a1)
	This subroutine calculates the mixing lengths.

subroutine	boulac_length0 (k, kts, kte, zw, dz, qtke, theta, lb1, lb2)
	This subroutine was taken from the BouLac scheme in WRF-ARW and modified for integration into the MYNN PBL scheme. WHILE loops were added to reduce the computational expense. This subroutine computes the length scales up and down and then computes the min, average of the up/down length scales, and also considers the distance to the surface.

subroutine	boulac_length (kts, kte, zw, dz, qtke, theta, lb1, lb2)
	This subroutine was taken from the BouLac scheme in WRF-ARW and modified for integration into the MYNN PBL scheme. WHILE loops were added to reduce the computational expense. This subroutine computes the length scales up and down and then computes the min, average of the up/down length scales, and also considers the distance to the surface.

subroutine	mym_turbulence (kts, kte, xland, closure, dz, dx, zw, u, v, thl, thetav, ql, qw, qke, tsq, qsq, cov, vt, vq, rmo, flt, fltv, flq, zi, theta, sh, sm, el, dfm, dfh, dfq, tcd, qcd, pdk, pdt, pdq, pdc, qwt1d, qshear1d, qbuoy1d, qdiss1d, tke_budget, psig_bl, psig_shcu, cldfra_bl1d, bl_mynn_mixlength, edmf_w1, edmf_a1, tkeprodtd, spp_pbl, rstoch_col)
	This subroutine calculates the vertical diffusivity coefficients and the production terms for the turbulent quantities.

subroutine	mym_predict (kts, kte, closure, delt, dz, ust, flt, flq, pmz, phh, el, dfq, rho, pdk, pdt, pdq, pdc, qke, tsq, qsq, cov, s_aw, s_awqke, bl_mynn_edmf_tke, qwt1d, qdiss1d, tke_budget)
	This subroutine predicts the turbulent quantities at the next step.

subroutine	mym_condensation (kts, kte, dx, dz, zw, xland, thl, qw, qv, qc, qi, qs, p, exner, tsq, qsq, cov, sh, el, bl_mynn_cloudpdf, qc_bl1d, qi_bl1d, cldfra_bl1d, pblh1, hfx1, vt, vq, th, sgm, rmo, spp_pbl, rstoch_col)
	This subroutine calculates the nonconvective component of the subgrid cloud fraction and mixing ratio as well as the functions used to calculate the buoyancy flux. Different cloud PDFs can be selected by use of the namelist parameter `bl_mynn_cloudpdf` .

subroutine	mynn_tendencies (kts, kte, i, delt, dz, rho, u, v, th, tk, qv, qc, qi, qs, qnc, qni, psfc, p, exner, thl, sqv, sqc, sqi, sqs, sqw, qnwfa, qnifa, qnbca, ozone, ust, flt, flq, flqv, flqc, wspd, uoce, voce, tsq, qsq, cov, tcd, qcd, dfm, dfh, dfq, du, dv, dth, dqv, dqc, dqi, dqs, dqnc, dqni, dqnwfa, dqnifa, dqnbca, dozone, diss_heat, s_aw, s_awthl, s_awqt, s_awqv, s_awqc, s_awu, s_awv, s_awqnc, s_awqni, s_awqnwfa, s_awqnifa, s_awqnbca, sd_aw, sd_awthl, sd_awqt, sd_awqv, sd_awqc, sd_awu, sd_awv, sub_thl, sub_sqv, sub_u, sub_v, det_thl, det_sqv, det_sqc, det_u, det_v, flag_qc, flag_qi, flag_qnc, flag_qni, flag_qs, flag_qnwfa, flag_qnifa, flag_qnbca, cldfra_bl1d, bl_mynn_cloudmix, bl_mynn_mixqt, bl_mynn_edmf, bl_mynn_edmf_mom, bl_mynn_mixscalars)
	This subroutine solves for tendencies of U, V, \(\theta\), qv, qc, and qi.

subroutine	moisture_check (kte, delt, dp, exner, qv, qc, qi, qs, th, dqv, dqc, dqi, dqs, dth)

subroutine	mynn_mix_chem (kts, kte, i, delt, dz, pblh, nchem, kdvel, ndvel, chem1, vd1, rho, flt, tcd, qcd, dfh, s_aw, s_awchem, emis_ant_no, frp, rrfs_sd, enh_mix, smoke_dbg)

subroutine	retrieve_exchange_coeffs (kts, kte, dfm, dfh, dz, k_m, k_h)

subroutine	tridiag (n, a, b, c, d)

subroutine	tridiag2 (n, a, b, c, d, x)

subroutine	tridiag3 (kte, a, b, c, d, x)

subroutine	dmp_mf (kts, kte, dt, zw, dz, p, rho, momentum_opt, tke_opt, scalar_opt, u, v, w, th, thl, thv, tk, qt, qv, qc, qke, qnc, qni, qnwfa, qnifa, qnbca, exner, vt, vq, sgm, ust, flt, fltv, flq, flqv, pblh, kpbl, dx, landsea, ts, edmf_a, edmf_w, edmf_qt, edmf_thl, edmf_ent, edmf_qc, s_aw, s_awthl, s_awqt, s_awqv, s_awqc, s_awu, s_awv, s_awqke, s_awqnc, s_awqni, s_awqnwfa, s_awqnifa, s_awqnbca, sub_thl, sub_sqv, sub_u, sub_v, det_thl, det_sqv, det_sqc, det_u, det_v, nchem, chem1, s_awchem, mix_chem, qc_bl1d, cldfra_bl1d, qc_bl1d_old, cldfra_bl1d_old, f_qc, f_qi, f_qnc, f_qni, f_qnwfa, f_qnifa, f_qnbca, psig_shcu, maxwidth, ktop, maxmf, ztop, spp_pbl, rstoch_col)
	This subroutine is the Dynamic Multi-Plume (DMP) Mass-Flux Scheme.

subroutine	condensation_edmf (qt, thl, p, zagl, thv, qc)
	This subroutine.

subroutine	condensation_edmf_r (qt, thl, p, zagl, thv, qc)

subroutine	ddmf_jpl (kts, kte, dt, zw, dz, p, u, v, th, thl, thv, tk, qt, qv, qc, rho, exner, ust, wthl, wqt, pblh, kpbl, edmf_a_dd, edmf_w_dd, edmf_qt_dd, edmf_thl_dd, edmf_ent_dd, edmf_qc_dd, sd_aw, sd_awthl, sd_awqt, sd_awqv, sd_awqc, sd_awu, sd_awv, sd_awqke, qc_bl1d, cldfra_bl1d, rthraten)

subroutine	scale_aware (dx, pbl1, psig_bl, psig_shcu)

real(kind_phys) function	esat_blend (t)

real(kind_phys) function	qsat_blend (t, p)
	This function extends function "esat" and returns a "blended" saturation mixing ratio. Tice currently set to 240 K, t0c = 273.15 K.

real(kind_phys) function	xl_blend (t)
	This function interpolates the latent heats of vaporization and sublimation into a single, temperature-dependent, "blended" value, following Chaboureau and Bechtold (2002) [30], Appendix.

real(kind_phys) function	phim (zet)

real(kind_phys) function	phih (zet)

subroutine	topdown_cloudrad (kts, kte, dz1, zw, xland, kpbl, pblh, sqc, sqi, sqw, thl, th1, ex1, p1, rho1, thetav, cldfra_bl1d, rthraten, maxkhtopdown, khtopdown, tkeprodtd)


subroutine	mynn_bl_driver (initflag, restart, cycling, delt, dz, dx, znt, u, v, w, th, sqv3d, sqc3d, sqi3d, sqs3d, qnc, qni, qnwfa, qnifa, qnbca, ozone, p, exner, rho, t3d, xland, ts, qsfc, ps, ust, ch, hfx, qfx, rmol, wspd, uoce, voce, qke, qke_adv, sh3d, sm3d, nchem, kdvel, ndvel, chem3d, vdep, smoke_dbg, frp, emis_ant_no, mix_chem, enh_mix, rrfs_sd, tsq, qsq, cov, rublten, rvblten, rthblten, rqvblten, rqcblten, rqiblten, rqncblten, rqniblten, rqsblten, rqnwfablten, rqnifablten, rqnbcablten, dozone, exch_h, exch_m, pblh, kpbl, el_pbl, dqke, qwt, qshear, qbuoy, qdiss, qc_bl, qi_bl, cldfra_bl, bl_mynn_tkeadvect, tke_budget, bl_mynn_cloudpdf, bl_mynn_mixlength, icloud_bl, closure, bl_mynn_edmf, bl_mynn_edmf_mom, bl_mynn_edmf_tke, bl_mynn_mixscalars, bl_mynn_output, bl_mynn_cloudmix, bl_mynn_mixqt, edmf_a, edmf_w, edmf_qt, edmf_thl, edmf_ent, edmf_qc, sub_thl3d, sub_sqv3d, det_thl3d, det_sqv3d, maxwidth, maxmf, ztop_plume, ktop_plume, spp_pbl, pattern_spp_pbl, rthraten, flag_qc, flag_qi, flag_qnc, flag_qni, flag_qs, flag_qnwfa, flag_qnifa, flag_qnbca, flag_ozone, ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, its, ite, jts, jte, kts, kte)
	This subroutine is the MYNN-EDNF PBL driver routine,which encompassed the majority of the subroutines that comprise the procedures that ultimately solve for tendencies of \(U, V, \theta, q_v, q_c, and q_i\).


subroutine	mym_initialize (kts, kte, xland, dz, dx, zw, u, v, thl, qw, zi, theta, thetav, sh, sm, ust, rmo, el, qke, tsq, qsq, cov, psig_bl, cldfra_bl1d, bl_mynn_mixlength, edmf_w1, edmf_a1, initialize_qke, spp_pbl, rstoch_col)
	This subroutine initializes the mixing length, TKE, \(\theta^{'2}\), \(q^{'2}\), and \(\theta^{'}q^{'}\).


subroutine	get_pblh (kts, kte, zi, thetav1d, qke1d, zw1d, dz1d, landsea, kzi)
	This subroutine calculates hybrid diagnotic boundary-layer height (PBLH).

Variables
real(kind_phys), parameter	cphm_st =5.0

real(kind_phys), parameter	cphm_unst =16.0

real(kind_phys), parameter	cphh_st =5.0

real(kind_phys), parameter	cphh_unst =16.0

real(kind_phys), parameter	pr = 0.74

real(kind_phys), parameter	g1 = 0.235

real(kind_phys), parameter	b1 = 24.0

real(kind_phys), parameter	b2 = 15.0

real(kind_phys), parameter	c2 = 0.729

real(kind_phys), parameter	c3 = 0.340

real(kind_phys), parameter	c4 = 0.0

real(kind_phys), parameter	c5 = 0.2

real(kind_phys), parameter	a1 = b1( 1.0-3.0g1 )/6.0

real(kind_phys), parameter	c1 = g1 -1.0/( 3.0a12.88449914061481660)

real(kind_phys), parameter	a2 = a1( g1-c1 )/( g1pr )

real(kind_phys), parameter	g2 = b2/b1( 1.0-c3 ) +2.0a1/b1( 3.0-2.0c2 )

real(kind_phys), parameter	cc2 = 1.0-c2

real(kind_phys), parameter	cc3 = 1.0-c3

real(kind_phys), parameter	e1c = 3.0a2b2*cc3

real(kind_phys), parameter	e2c = 9.0a1a2*cc2

real(kind_phys), parameter	e3c = 9.0a2a2cc2( 1.0-c5 )

real(kind_phys), parameter	e4c = 12.0a1a2*cc2

real(kind_phys), parameter	e5c = 6.0a1a1

real(kind_phys), parameter	qmin =0.0

real(kind_phys), parameter	zmax =1.0

real(kind_phys), parameter	sqfac =3.0

real(kind_phys), parameter	qkemin =1.e-4

real(kind_phys), parameter	tliq = 269.

real(kind_phys), parameter	rr2 =0.7071068

real(kind_phys), parameter	rrp =0.3989423

real(kind_phys), parameter	ckmod =1.
	Use Canuto/Kitamura mod (remove Ric and negative TKE) (1:yes, 0:no) For more info, see Canuto et al. (2008 JAS) and Kitamura (Journal of the Meteorological Society of Japan, Vol. 88, No. 5, pp. 857-864, 2010). Note that this change required further modification of other parameters above (c2, c3). If you want to remove this option, set c2 and c3 constants (above) back to NN2009 values (see commented out lines next to the parameters above). This only removes the negative TKE problem but does not necessarily improve performance - neutral impact.

real(kind_phys), parameter	scaleaware =1.
	Use Ito et al. (2015, BLM) scale-aware (0: no, 1: yes). Note that this also has impacts on the cloud PDF and mass-flux scheme, using Honnert et al. (2011) similarity function for TKE in the upper PBL/cloud layer.

integer, parameter	bl_mynn_topdown = 0
	Of the following the options, use one OR the other, not both. Adding top-down diffusion driven by cloud-top radiative cooling.

integer, parameter	bl_mynn_edmf_dd = 0
	Option to activate downdrafts, from Elynn Wu (0: deactive, 1: active)

integer, parameter	dheat_opt = 1
	Option to activate heating due to dissipation of TKE (to activate, set to 1.0)

logical, parameter	env_subs = .false.

integer, parameter	bl_mynn_stfunc = 1

logical, parameter	debug_code = .false.

integer, parameter	idbg = 23

integer	mynn_level

Detailed Description

This module contains the entity of MYNN-EDMF PBL scheme.

Function/Subroutine Documentation

◆ condensation_edmf_r()

subroutine module_bl_mynn::condensation_edmf_r	(	real(kind_phys), intent(in)	qt,
		real(kind_phys), intent(out)	thl,
		real(kind_phys), intent(in)	p,
		real(kind_phys), intent(in)	zagl,
		real(kind_phys), intent(in)	thv,
		real(kind_phys), intent(out)	qc
	)

Definition at line 6891 of file module_bl_mynn.F90.

◆ ddmf_jpl()

subroutine module_bl_mynn::ddmf_jpl	(	integer, intent(in)	kts,
		integer, intent(in)	kte,
		real(kind_phys), intent(in)	dt,
		real(kind_phys), dimension(kts:kte+1), intent(in)	zw,
		real(kind_phys), dimension(kts:kte), intent(in)	dz,
		real(kind_phys), dimension(kts:kte), intent(in)	p,
		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)	th,
		real(kind_phys), dimension(kts:kte), intent(in)	thl,
		real(kind_phys), dimension(kts:kte), intent(in)	thv,
		real(kind_phys), dimension(kts:kte), intent(in)	tk,
		real(kind_phys), dimension(kts:kte), intent(in)	qt,
		real(kind_phys), dimension(kts:kte), intent(in)	qv,
		real(kind_phys), dimension(kts:kte), intent(in)	qc,
		real(kind_phys), dimension(kts:kte), intent(in)	rho,
		real(kind_phys), dimension(kts:kte), intent(in)	exner,
		real(kind_phys), intent(in)	ust,
		real(kind_phys), intent(in)	wthl,
		real(kind_phys), intent(in)	wqt,
		real(kind_phys), intent(in)	pblh,
		integer, intent(in)	kpbl,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_a_dd,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_w_dd,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_qt_dd,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_thl_dd,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_ent_dd,
		real(kind_phys), dimension(kts:kte), intent(out)	edmf_qc_dd,
		real(kind_phys), dimension(kts:kte+1)	sd_aw,
		real(kind_phys), dimension(kts:kte+1)	sd_awthl,
		real(kind_phys), dimension(kts:kte+1)	sd_awqt,
		real(kind_phys), dimension(kts:kte+1)	sd_awqv,
		real(kind_phys), dimension(kts:kte+1)	sd_awqc,
		real(kind_phys), dimension(kts:kte+1)	sd_awu,
		real(kind_phys), dimension(kts:kte+1)	sd_awv,
		real(kind_phys), dimension(kts:kte+1)	sd_awqke,
		real(kind_phys), dimension(kts:kte), intent(in)	qc_bl1d,
		real(kind_phys), dimension(kts:kte), intent(in)	cldfra_bl1d,
		real(kind_phys), dimension(kts:kte), intent(in)	rthraten
	)

Definition at line 6933 of file module_bl_mynn.F90.

◆ moisture_check()

subroutine module_bl_mynn::moisture_check	(	integer, intent(in)	kte,
		real(kind_phys), intent(in)	delt,
		real(kind_phys), dimension(kte), intent(in)	dp,
		real(kind_phys), dimension(kte), intent(in)	exner,
		real(kind_phys), dimension(kte), intent(inout)	qv,
		real(kind_phys), dimension(kte), intent(inout)	qc,
		real(kind_phys), dimension(kte), intent(inout)	qi,
		real(kind_phys), dimension(kte), intent(inout)	qs,
		real(kind_phys), dimension(kte), intent(inout)	th,
		real(kind_phys), dimension(kte), intent(inout)	dqv,
		real(kind_phys), dimension(kte), intent(inout)	dqc,
		real(kind_phys), dimension(kte), intent(inout)	dqi,
		real(kind_phys), dimension(kte), intent(inout)	dqs,
		real(kind_phys), dimension(kte), intent(inout)	dth
	)

Definition at line 5133 of file module_bl_mynn.F90.

◆ mynn_mix_chem()

subroutine module_bl_mynn::mynn_mix_chem	(	integer, intent(in)	kts,
		integer, intent(in)	kte,
		integer, intent(in)	i,
		real(kind_phys), intent(in)	delt,
		real(kind_phys), dimension(kts:kte), intent(in)	dz,
		real(kind_phys), intent(in)	pblh,
		integer, intent(in)	nchem,
		integer, intent(in)	kdvel,
		integer, intent(in)	ndvel,
		real(kind_phys), dimension( kts:kte, nchem ), intent(inout)	chem1,
		real(kind_phys), dimension( ndvel ), intent(in)	vd1,
		real(kind_phys), dimension(kts:kte), intent(inout)	rho,
		real(kind_phys), intent(in)	flt,
		real(kind_phys), dimension(kts:kte), intent(in)	tcd,
		real(kind_phys), dimension(kts:kte), intent(in)	qcd,
		real(kind_phys), dimension(kts:kte), intent(in)	dfh,
		real(kind_phys), dimension( kts:kte+1), intent(in)	s_aw,
		real(kind_phys), dimension( kts:kte+1,nchem), intent(in)	s_awchem,
		real(kind_phys), intent(in)	emis_ant_no,
		real(kind_phys), intent(in)	frp,
		logical, intent(in)	rrfs_sd,
		logical, intent(in)	enh_mix,
		logical, intent(in)	smoke_dbg
	)

Definition at line 5224 of file module_bl_mynn.F90.

◆ phih()

real(kind_phys) function module_bl_mynn::phih ( real(kind_phys), intent(in) zet )

Definition at line 7580 of file module_bl_mynn.F90.

◆ phim()

real(kind_phys) function module_bl_mynn::phim ( real(kind_phys), intent(in) zet )

Definition at line 7528 of file module_bl_mynn.F90.

◆ scale_aware()

subroutine module_bl_mynn::scale_aware	(	real(kind_phys), intent(in)	dx,
		real(kind_phys), intent(in)	pbl1,
		real(kind_phys), intent(out)	psig_bl,
		real(kind_phys), intent(out)	psig_shcu
	)

Definition at line 7314 of file module_bl_mynn.F90.

◆ topdown_cloudrad()

subroutine module_bl_mynn::topdown_cloudrad	(	integer, intent(in)	kts,
		integer, intent(in)	kte,
		real(kind_phys), dimension(kts:kte), intent(in)	dz1,
		real(kind_phys), dimension(kts:kte+1), intent(in)	zw,
		real(kind_phys), intent(in)	xland,
		integer, intent(in)	kpbl,
		real(kind_phys), intent(in)	pblh,
		real(kind_phys), dimension(kts:kte), intent(in)	sqc,
		real(kind_phys), dimension(kts:kte), intent(in)	sqi,
		real(kind_phys), dimension(kts:kte), intent(in)	sqw,
		real(kind_phys), dimension(kts:kte), intent(in)	thl,
		real(kind_phys), dimension(kts:kte), intent(in)	th1,
		real(kind_phys), dimension(kts:kte), intent(in)	ex1,
		real(kind_phys), dimension(kts:kte), intent(in)	p1,
		real(kind_phys), dimension(kts:kte), intent(in)	rho1,
		real(kind_phys), dimension(kts:kte), intent(in)	thetav,
		real(kind_phys), dimension(kts:kte), intent(in)	cldfra_bl1d,
		real(kind_phys), dimension(kts:kte), intent(in)	rthraten,
		real(kind_phys), intent(out)	maxkhtopdown,
		real(kind_phys), dimension(kts:kte), intent(out)	khtopdown,
		real(kind_phys), dimension(kts:kte), intent(out)	tkeprodtd
	)

Definition at line 7628 of file module_bl_mynn.F90.

Variable Documentation

◆ a1

real(kind_phys), parameter module_bl_mynn::a1 = b1*( 1.0-3.0*g1 )/6.0

Definition at line 276 of file module_bl_mynn.F90.

◆ a2

real(kind_phys), parameter module_bl_mynn::a2 = a1*( g1-c1 )/( g1*pr )

Definition at line 276 of file module_bl_mynn.F90.

◆ b1

real(kind_phys), parameter module_bl_mynn::b1 = 24.0

Definition at line 276 of file module_bl_mynn.F90.

◆ b2

real(kind_phys), parameter module_bl_mynn::b2 = 15.0

Definition at line 276 of file module_bl_mynn.F90.

◆ bl_mynn_edmf_dd

integer, parameter module_bl_mynn::bl_mynn_edmf_dd = 0

Option to activate downdrafts, from Elynn Wu (0: deactive, 1: active)

Definition at line 331 of file module_bl_mynn.F90.

◆ bl_mynn_stfunc

integer, parameter module_bl_mynn::bl_mynn_stfunc = 1

Definition at line 341 of file module_bl_mynn.F90.

◆ bl_mynn_topdown

integer, parameter module_bl_mynn::bl_mynn_topdown = 0

Of the following the options, use one OR the other, not both. Adding top-down diffusion driven by cloud-top radiative cooling.

Definition at line 329 of file module_bl_mynn.F90.

◆ c1

real(kind_phys), parameter module_bl_mynn::c1 = g1 -1.0/( 3.0*a1*2.88449914061481660)

Definition at line 276 of file module_bl_mynn.F90.

◆ c2

real(kind_phys), parameter module_bl_mynn::c2 = 0.729

Definition at line 276 of file module_bl_mynn.F90.

◆ c3

real(kind_phys), parameter module_bl_mynn::c3 = 0.340

Definition at line 276 of file module_bl_mynn.F90.

◆ c4

real(kind_phys), parameter module_bl_mynn::c4 = 0.0

Definition at line 276 of file module_bl_mynn.F90.

◆ c5

real(kind_phys), parameter module_bl_mynn::c5 = 0.2

Definition at line 276 of file module_bl_mynn.F90.

◆ cc2

real(kind_phys), parameter module_bl_mynn::cc2 = 1.0-c2

Definition at line 291 of file module_bl_mynn.F90.

◆ cc3

real(kind_phys), parameter module_bl_mynn::cc3 = 1.0-c3

Definition at line 291 of file module_bl_mynn.F90.

◆ ckmod

real(kind_phys), parameter module_bl_mynn::ckmod =1.

Use Canuto/Kitamura mod (remove Ric and negative TKE) (1:yes, 0:no) For more info, see Canuto et al. (2008 JAS) and Kitamura (Journal of the Meteorological Society of Japan, Vol. 88, No. 5, pp. 857-864, 2010). Note that this change required further modification of other parameters above (c2, c3). If you want to remove this option, set c2 and c3 constants (above) back to NN2009 values (see commented out lines next to the parameters above). This only removes the negative TKE problem but does not necessarily improve performance - neutral impact.

Definition at line 320 of file module_bl_mynn.F90.

◆ cphh_st

real(kind_phys), parameter module_bl_mynn::cphh_st =5.0

Definition at line 272 of file module_bl_mynn.F90.

◆ cphh_unst

real(kind_phys), parameter module_bl_mynn::cphh_unst =16.0

Definition at line 272 of file module_bl_mynn.F90.

◆ cphm_st

real(kind_phys), parameter module_bl_mynn::cphm_st =5.0

Definition at line 272 of file module_bl_mynn.F90.

◆ cphm_unst

real(kind_phys), parameter module_bl_mynn::cphm_unst =16.0

Definition at line 272 of file module_bl_mynn.F90.

◆ debug_code

logical, parameter module_bl_mynn::debug_code = .false.

Definition at line 344 of file module_bl_mynn.F90.

◆ dheat_opt

integer, parameter module_bl_mynn::dheat_opt = 1

Option to activate heating due to dissipation of TKE (to activate, set to 1.0)

Definition at line 334 of file module_bl_mynn.F90.

◆ e1c

real(kind_phys), parameter module_bl_mynn::e1c = 3.0*a2*b2*cc3

Definition at line 291 of file module_bl_mynn.F90.

◆ e2c

real(kind_phys), parameter module_bl_mynn::e2c = 9.0*a1*a2*cc2

Definition at line 291 of file module_bl_mynn.F90.

◆ e3c

real(kind_phys), parameter module_bl_mynn::e3c = 9.0*a2*a2*cc2*( 1.0-c5 )

Definition at line 291 of file module_bl_mynn.F90.

◆ e4c

real(kind_phys), parameter module_bl_mynn::e4c = 12.0*a1*a2*cc2

Definition at line 291 of file module_bl_mynn.F90.

◆ e5c

real(kind_phys), parameter module_bl_mynn::e5c = 6.0*a1*a1

Definition at line 291 of file module_bl_mynn.F90.

◆ env_subs

logical, parameter module_bl_mynn::env_subs = .false.

Definition at line 337 of file module_bl_mynn.F90.

◆ g1

real(kind_phys), parameter module_bl_mynn::g1 = 0.235

Definition at line 276 of file module_bl_mynn.F90.

◆ g2

real(kind_phys), parameter module_bl_mynn::g2 = b2/b1*( 1.0-c3 ) +2.0*a1/b1*( 3.0-2.0*c2 )

Definition at line 276 of file module_bl_mynn.F90.

◆ idbg

integer, parameter module_bl_mynn::idbg = 23

Definition at line 345 of file module_bl_mynn.F90.

◆ mynn_level

integer module_bl_mynn::mynn_level

Definition at line 348 of file module_bl_mynn.F90.

◆ pr

real(kind_phys), parameter module_bl_mynn::pr = 0.74

Definition at line 276 of file module_bl_mynn.F90.

◆ qkemin

real(kind_phys), parameter module_bl_mynn::qkemin =1.e-4

Definition at line 306 of file module_bl_mynn.F90.

◆ qmin

real(kind_phys), parameter module_bl_mynn::qmin =0.0

Definition at line 302 of file module_bl_mynn.F90.

◆ rr2

real(kind_phys), parameter module_bl_mynn::rr2 =0.7071068

Definition at line 310 of file module_bl_mynn.F90.

◆ rrp

real(kind_phys), parameter module_bl_mynn::rrp =0.3989423

Definition at line 310 of file module_bl_mynn.F90.

◆ scaleaware

real(kind_phys), parameter module_bl_mynn::scaleaware =1.

Use Ito et al. (2015, BLM) scale-aware (0: no, 1: yes). Note that this also has impacts on the cloud PDF and mass-flux scheme, using Honnert et al. (2011) similarity function for TKE in the upper PBL/cloud layer.

Definition at line 325 of file module_bl_mynn.F90.

◆ sqfac

real(kind_phys), parameter module_bl_mynn::sqfac =3.0

Definition at line 302 of file module_bl_mynn.F90.

◆ tliq

real(kind_phys), parameter module_bl_mynn::tliq = 269.

Definition at line 307 of file module_bl_mynn.F90.

◆ zmax

real(kind_phys), parameter module_bl_mynn::zmax =1.0

Definition at line 302 of file module_bl_mynn.F90.

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ condensation_edmf_r()

◆ ddmf_jpl()

◆ moisture_check()

◆ mynn_mix_chem()

◆ phih()

◆ phim()

◆ scale_aware()

◆ topdown_cloudrad()

Variable Documentation

◆ a1

◆ a2

◆ b1

◆ b2

◆ bl_mynn_edmf_dd

◆ bl_mynn_stfunc

◆ bl_mynn_topdown

◆ c1

◆ c2

◆ c3

◆ c4

◆ c5

◆ cc2

◆ cc3

◆ ckmod

◆ cphh_st

◆ cphh_unst

◆ cphm_st

◆ cphm_unst

◆ debug_code

◆ dheat_opt

◆ e1c

◆ e2c

◆ e3c

◆ e4c

◆ e5c

◆ env_subs

◆ g1

◆ g2

◆ idbg

◆ mynn_level

◆ pr

◆ qkemin

◆ qmin

◆ rr2

◆ rrp

◆ scaleaware

◆ sqfac

◆ tliq

◆ zmax