ERF_InitCustomPert_SDMCongestus3D.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp_prob ("prob")
pp_prob	query ("rho_0", rho_0)
pp_prob	query ("T_0", T_0)
pp_prob	query ("A_0", A_0)
pp_prob	query ("KE_0", KE_0)
pp_prob	query ("T_0_Pert_Mag", T_0_Pert_Mag)
pp_prob	query ("Q_0_Pert_Mag", T_0_Pert_Mag)
pp_prob	query ("pert_ref_height", pert_ref_height)
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)
	ParallelForRNG (bx, [=] AMREX_GPU_DEVICE(int i, int j, int k, const amrex::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 rand_double=amrex::Random(engine);state_pert(i, j, k, RhoTheta_comp)=(rand_double *2.0 - 1.0) *T_0_Pert_Mag;} if((z<=pert_ref_height) &&(Q_0_Pert_Mag !=0.0)) { Real rand_double=amrex::Random(engine);state_pert(i, j, k, RhoQ1_comp)=(rand_double *2.0 - 1.0) *Q_0_Pert_Mag;} state_pert(i, j, k, RhoScalar_comp)=A_0 *exp(-10.*r *r);state_pert(i, j, k, RhoKE_comp)=KE_0;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	rho_0 = 0.0
Real	T_0 = 0.0
Real	A_0 = 1.0
Real	KE_0 = 0.1
Real	T_0_Pert_Mag = 0.0
Real	Q_0_Pert_Mag = 0.0
Real	pert_ref_height = 100.0
Real	pert_deltaT = 0.0
Real	pert_deltaQV = 0.0
Real	pert_periods_T = 0.0
Real	pert_periods_QV = 0.0
Real	prob_lo_x = geomdata.ProbLo(0)
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 amrex::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 rand_double=amrex::Random(engine);state_pert(i, j, k, RhoTheta_comp)=(rand_double *2.0 - 1.0) *T_0_Pert_Mag;} if((z<=pert_ref_height) &&(Q_0_Pert_Mag !=0.0)) { Real rand_double=amrex::Random(engine);state_pert(i, j, k, RhoQ1_comp)=(rand_double *2.0 - 1.0) *Q_0_Pert_Mag;} state_pert(i, j, k, RhoScalar_comp)=A_0 *exp(-10.*r *r);state_pert(i, j, k, RhoKE_comp)=KE_0;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;} }
	)

Here is the call graph for this function:

pp_prob()

ParmParse pp_prob

(

"prob"

)

query() [1/11]

pp_prob query	(	"A_0"	,
		A_0
	)

query() [2/11]

pp_prob query	(	"KE_0"	,
		KE_0
	)

query() [3/11]

pp_prob query	(	"pert_deltaQV"	,
		pert_deltaQV
	)

query() [4/11]

pp_prob query	(	"pert_deltaT"	,
		pert_deltaT
	)

query() [5/11]

pp_prob query	(	"pert_periods_QV"	,
		pert_periods_QV
	)

query() [6/11]

pp_prob query	(	"pert_periods_T"	,
		pert_periods_T
	)

query() [7/11]

pp_prob query	(	"pert_ref_height"	,
		pert_ref_height
	)

query() [8/11]

pp_prob query	(	"Q_0_Pert_Mag"	,
		T_0_Pert_Mag
	)

query() [9/11]

pp_prob query	(	"rho_0"	,
		rho_0
	)

query() [10/11]

pp_prob query	(	"T_0"	,
		T_0
	)

query() [11/11]

pp_prob query	(	"T_0_Pert_Mag"	,
		T_0_Pert_Mag
	)

Variable Documentation

A_0

Real A_0 = 1.0

Referenced by ParallelForRNG().

cval

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

Referenced by ParallelForRNG().

dval

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

Referenced by ParallelForRNG().

KE_0

Real KE_0 = 0.1

Referenced by ParallelForRNG().

pert_deltaQV

Real pert_deltaQV = 0.0

Referenced by ParallelForRNG().

pert_deltaT

Real pert_deltaT = 0.0

Referenced by ParallelForRNG().

pert_periods_QV

Real pert_periods_QV = 0.0

pert_periods_T

Real pert_periods_T = 0.0

pert_ref_height

Real pert_ref_height = 100.0

Referenced by ParallelForRNG().

prob_hi_x

Real prob_hi_x = geomdata.ProbHi(0)

prob_lo_x

Real prob_lo_x = geomdata.ProbLo(0)

Q_0_Pert_Mag

Real Q_0_Pert_Mag = 0.0

Referenced by ParallelForRNG().

qvfac

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

Referenced by ParallelForRNG().

rdOcp

const Real rdOcp = sc.rdOcp

Referenced by ParallelForRNG().

rho_0

Real rho_0 = 0.0

T_0

Real T_0 = 0.0

T_0_Pert_Mag

Real T_0_Pert_Mag = 0.0

Referenced by ParallelForRNG().

tfac

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

Referenced by ParallelForRNG().

use_moisture

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

Referenced by ParallelForRNG().