#include <ERF_MYNNStruct.H>

Collaboration diagram for MYNNLevel2:

Public Member Functions
void	init_coeffs (amrex::Real A1_lvl25, amrex::Real A2_lvl25, amrex::Real B1, amrex::Real B2, amrex::Real C1, amrex::Real C2, amrex::Real C3, amrex::Real, amrex::Real C5)

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_Rf (const amrex::Real GM, const amrex::Real GH) const

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_SH (const amrex::Real Rf) const

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real	calc_SM (const amrex::Real Rf) const

Public Attributes
bool	initialized {false}

amrex::Real	A1

amrex::Real	A2

amrex::Real	F1

amrex::Real	F2

const amrex::Real	gam1 = 0.235

amrex::Real	gam2

amrex::Real	Rf1

amrex::Real	Rf2

amrex::Real	Rfc

amrex::Real	Ri1

amrex::Real	Ri2

amrex::Real	Ri3

const amrex::Real	eps = std::numeric_limits<amrex::Real>::epsilon()

Member Function Documentation

◆ calc_Rf()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_Rf	(	const amrex::Real	GM,
		const amrex::Real	GH
	)		const

inline

     {
         AMREX_ASSERT(initialized);
         amrex::Real Ri = -GH/(GM + eps);
         return Ri1*(Ri + Ri2 - std::sqrt(Ri*Ri - Ri3*Ri + Ri2*Ri2));
     }

◆ calc_SH()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_SH ( const amrex::Real Rf ) const

inline

     {
         AMREX_ASSERT(initialized);
         return 3.0*A2*(gam1+gam2)*(Rfc-Rf)/(1.0-Rf);
     }

Referenced by calc_SM().

Here is the caller graph for this function:

◆ calc_SM()

AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real MYNNLevel2::calc_SM ( const amrex::Real Rf ) const

inline

     {
         AMREX_ASSERT(initialized);
         return A1*F1/(A2*F2) * (Rf1-Rf) / (Rf2-Rf) * calc_SH(Rf);
     }

Here is the call graph for this function:

◆ init_coeffs()

void MYNNLevel2::init_coeffs	(	amrex::Real	A1_lvl25,
		amrex::Real	A2_lvl25,
		amrex::Real	B1,
		amrex::Real	B2,
		amrex::Real	C1,
		amrex::Real	C2,
		amrex::Real	C3,
		amrex::Real	,
		amrex::Real	C5
	)

inline

     {
         A1 = A1_lvl25;
         A2 = A2_lvl25;
         gam2 = (2.0*A1*(3.0-2.0*C2) + B2*(1.0-C3)) / B1;
         F1  = B1*(gam1-C1) + 2.0*A1*(3.0-2.0*C2) + 3.0*A2*(1.0-C2)*(1.0-C5);
         F2  = B1*(gam1+gam2) - 3.0*A1*(1.0-C2);
         Rf1 = B1*(gam1-C1) / F1;
         Rf2 = B1*gam1 / F2;
         Rfc = gam1 / (gam1 + gam2);
         Ri1 = 0.5*A2*F2/(A1*F1);
         Ri2 = 0.5*Rf1/Ri1;
         Ri3 = (2.0*Rf2-Rf1)/Ri1;
  
         initialized = true;
     }