#include <ERF_MYNNStruct.H>

Collaboration diagram for MYNNLevel25:

[legend]

Public Member Functions
AMREX_GPU_DEVICE AMREX_FORCE_INLINE void	calc_stability_funcs (amrex::Real &SM, amrex::Real &SH, amrex::Real &SQ, const amrex::Real GM, const amrex::Real GH, const amrex::Real alphac=1.0) const

Public Attributes
amrex::Real	A1 = 1.18

amrex::Real	A2 = 0.665

amrex::Real	B1 = 24.0

amrex::Real	B2 = 15.0

amrex::Real	C1 = 0.137

amrex::Real	C2 = 0.75

amrex::Real	C3 = 0.352

amrex::Real	C4 = 0.0

amrex::Real	C5 = 0.2

amrex::Real	SMmin = 0.0

amrex::Real	SMmax = 1.0e34

amrex::Real	SHmin = 0.0

amrex::Real	SHmax = 4.0

MYNNConfigType	config

bool	diffuse_moistvars = false

Member Function Documentation

◆ calc_stability_funcs()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE void MYNNLevel25::calc_stability_funcs	(	amrex::Real &	SM,
		amrex::Real &	SH,
		amrex::Real &	SQ,
		const amrex::Real	GM,
		const amrex::Real	GH,
		const amrex::Real	alphac = `1.0`
	)		const

inline

     {
         amrex::Real alphac2 = alphac * alphac;
  
         // Compute non-dimensional parameters (notation follows NN09)
         amrex::Real Phi1 = 1.0  - alphac2*3.0*A2*B2*(1-C3)*GH;
         amrex::Real Phi2 = 1.0  - alphac2*9.0*A1*A2*(1-C2)*GH;
         amrex::Real Phi3 = Phi1 + alphac2*9.0*A2*A2*(1-C2)*(1-C5)*GH;
         amrex::Real Phi4 = Phi1 - alphac2*12.0*A1*A2*(1-C2)*GH;
         amrex::Real Phi5 = 6.0*alphac*A1*A1*GM;
  
         // Compute stability functions
         amrex::Real D = Phi2*Phi4 + Phi5*Phi3;
         SM = alphac * A1 * (Phi3 - 3*C1*Phi4) / D;
         SH = alphac * A2 * (Phi2 + 3*C1*Phi5) / D;
         SQ = 3.0 * SM; // revised in NN09
     }