Functions
ParmParse	pp_prob ("prob")

pp_prob	query ("A_0", A_0)

pp_prob	query ("KE_0", KE_0)

pp_prob	query ("pert_deltaT", pert_deltaT)

pp_prob	query ("pert_deltaQV", pert_deltaQV)

pp_prob	query ("pert_periods_T", pert_periods_T)

pp_prob	query ("pert_periods_QV", pert_periods_QV)

pp_prob	query ("T_0_Pert_Mag", T_0_Pert_Mag)

pp_prob	query ("qv_0_Pert_Mag", qv_0_Pert_Mag)

pp_prob	query ("wbar_sub_max", wbar_sub_max)

pp_prob	query ("pert_ref_height", pert_ref_height)

pp_prob	query ("custom_TKE", custom_TKE)

	ParallelForRNG (bx, [=] AMREX_GPU_DEVICE(int i, int j, int k, const RandomEngine &engine) noexcept { const Real prob_lo=geomdata.ProbLo();const Real prob_hi=geomdata.ProbHi();const Real dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+0.5) dx[0];const Real y=prob_lo[1]+(j+0.5) dx[1];const Real z=prob_lo[2]+(k+0.5) dx[2];const Real xc=0.5 (prob_lo[0]+prob_hi[0]);const Real yc=0.5 (prob_lo[1]+prob_hi[1]);const Real zc=0.5 (prob_lo[2]+prob_hi[2]);const Real r=std::sqrt((x-xc) (x-xc)+(y-yc) (y-yc)+(z-zc) (z-zc));if((z<=pert_ref_height) &&(T_0_Pert_Mag !=0.0)) { Real rhotheta=state(i, j, k, RhoTheta_comp);Real rho=state(i, j, k, Rho_comp);Real qv=state(i, j, k, RhoQ1_comp)/rho;Real Told=getTgivenRandRTh(rho, rhotheta, qv);Real P=getPgivenRTh(rhotheta, qv);Real rand_double=amrex::Random(engine);Real Tpert=(rand_double 2.0 - 1.0) T_0_Pert_Mag;Real Tnew=Told+Tpert;Real theta_new=getThgivenTandP(Tnew, P, rdOcp);Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=0.0;} state_pert(i, j, k, RhoScalar_comp)=A_0 exp(-10.r r);if(custom_TKE) { state_pert(i, j, k, RhoKE_comp)=(1.0 - z/prob_hi[2]) r_hse(i, j, k);} else { state_pert(i, j, k, RhoKE_comp)=KE_0;} if(use_moisture) { if((z<=pert_ref_height) &&(qv_0_Pert_Mag !=0.0)) { Real rhoold=state(i, j, k, Rho_comp);Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;Real rand_double=amrex::Random(engine);Real qvnew=qvold+(rand_double 2.0 - 1.0) qv_0_Pert_Mag;state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} } if(pert_deltaT !=0.0) { const Real zl=z/pert_ref_height;const Real damp=std::exp(-0.5 zl zl);Real rhotheta=state(i, j, k, RhoTheta_comp);Real rho=state(i, j, k, Rho_comp);Real qv=state(i, j, k, RhoQ1_comp)/rho;Real Told=getTgivenRandRTh(rho, rhotheta, qv);Real P=getPgivenRTh(rhotheta, qv);Real Tpert=tfac damp z std::cos(cval (x - xc));Real Tnew=Told+Tpert;Real theta_new=getThgivenTandP(Tnew, P, rdOcp);Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=0.0;} if(use_moisture &&pert_deltaQV !=0.0) { const Real zl=z/pert_ref_height;const Real damp=std::exp(-0.5 zl zl);Real rhoold=state(i, j, k, Rho_comp);Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;Real qvnew=qvold+qvfac damp z std::cos(dval (x - xc));state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} })

Variables
Real	A_0 = 1.0

Real	KE_0 = 0.1

Real	pert_deltaT = 0.0

Real	pert_deltaQV = 0.0

Real	pert_periods_T = 0.0

Real	pert_periods_QV = 0.0

Real	T_0_Pert_Mag = 0.0

Real	qv_0_Pert_Mag = 0.0

Real	wbar_sub_max = -0.65

Real	pert_ref_height = 100.0

bool	custom_TKE = false

const Real	prob_lo_x = geomdata.ProbLo()[0]

const Real	prob_hi_x = geomdata.ProbHi()[0]

Real	cval = pert_periods_T * 2.0 * PI / (prob_hi_x - prob_lo_x)

Real	dval = pert_periods_QV * 2.0 * PI / (prob_hi_x - prob_lo_x)

Real	tfac = pert_deltaT * std::exp(0.5) / pert_ref_height

Real	qvfac = pert_deltaQV * std::exp(0.5) / pert_ref_height

const bool	use_moisture = (sc.moisture_type != MoistureType::None)

const Real	rdOcp = sc.rdOcp

Function Documentation

◆ ParallelForRNG()

ParallelForRNG	(	bx	,
		[=] AMREX_GPU_DEVICE(int i, int j, int k, const RandomEngine &engine) noexcept { const Real prob_lo=geomdata.ProbLo();const Real prob_hi=geomdata.ProbHi();const Real dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+0.5) dx[0];const Real y=prob_lo[1]+(j+0.5) dx[1];const Real z=prob_lo[2]+(k+0.5) dx[2];const Real xc=0.5 (prob_lo[0]+prob_hi[0]);const Real yc=0.5 (prob_lo[1]+prob_hi[1]);const Real zc=0.5 (prob_lo[2]+prob_hi[2]);const Real r=std::sqrt((x-xc) (x-xc)+(y-yc) (y-yc)+(z-zc) (z-zc));if((z<=pert_ref_height) &&(T_0_Pert_Mag !=0.0)) { Real rhotheta=state(i, j, k, RhoTheta_comp);Real rho=state(i, j, k, Rho_comp);Real qv=state(i, j, k, RhoQ1_comp)/rho;Real Told=getTgivenRandRTh(rho, rhotheta, qv);Real P=getPgivenRTh(rhotheta, qv);Real rand_double=amrex::Random(engine);Real Tpert=(rand_double 2.0 - 1.0) T_0_Pert_Mag;Real Tnew=Told+Tpert;Real theta_new=getThgivenTandP(Tnew, P, rdOcp);Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=0.0;} state_pert(i, j, k, RhoScalar_comp)=A_0 exp(-10.r r);if(custom_TKE) { state_pert(i, j, k, RhoKE_comp)=(1.0 - z/prob_hi[2]) r_hse(i, j, k);} else { state_pert(i, j, k, RhoKE_comp)=KE_0;} if(use_moisture) { if((z<=pert_ref_height) &&(qv_0_Pert_Mag !=0.0)) { Real rhoold=state(i, j, k, Rho_comp);Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;Real rand_double=amrex::Random(engine);Real qvnew=qvold+(rand_double 2.0 - 1.0) qv_0_Pert_Mag;state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} } if(pert_deltaT !=0.0) { const Real zl=z/pert_ref_height;const Real damp=std::exp(-0.5 zl zl);Real rhotheta=state(i, j, k, RhoTheta_comp);Real rho=state(i, j, k, Rho_comp);Real qv=state(i, j, k, RhoQ1_comp)/rho;Real Told=getTgivenRandRTh(rho, rhotheta, qv);Real P=getPgivenRTh(rhotheta, qv);Real Tpert=tfac damp z std::cos(cval (x - xc));Real Tnew=Told+Tpert;Real theta_new=getThgivenTandP(Tnew, P, rdOcp);Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=0.0;} if(use_moisture &&pert_deltaQV !=0.0) { const Real zl=z/pert_ref_height;const Real damp=std::exp(-0.5 zl zl);Real rhoold=state(i, j, k, Rho_comp);Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;Real qvnew=qvold+qvfac damp z std::cos(dval (x - xc));state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} }
	)