Functions
ParmParse	pp ("prob")

pp	query ("A_0", A_0)

pp	query ("KE_0", KE_0)

pp	query ("T_0_Pert_Mag", T_0_Pert_Mag)

pp	query ("qv_0_Pert_Mag", qv_0_Pert_Mag)

pp	query ("pert_ref_height", pert_ref_height)

pp	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+myhalf) dx[0];const Real y=prob_lo[1]+(j+myhalf) dx[1];const Real z=use_terrain ? z_cc(i, j, k) :prob_lo[2]+(k+myhalf) dx[2];const Real z_agl=z - erf_prob_heights::z_ground_cc(i, j, z_sfc_nd);const Real xc=myhalf (prob_lo[0]+prob_hi[0]);const Real yc=myhalf (prob_lo[1]+prob_hi[1]);const Real zc=myhalf (phys_z_lo+phys_z_hi);const Real r=std::sqrt((x-xc) (x-xc)+(y-yc) (y-yc)+(z-zc) (z-zc));if(pert_ref_height > zero &&z_agl<=pert_ref_height &&T_0_Pert_Mag !=zero) { const Real rhotheta=state(i, j, k, RhoTheta_comp);const Real rho=state(i, j, k, Rho_comp);const Real qv=(use_moisture) ? state(i, j, k, RhoQ1_comp)/rho :zero;const Real Told=getTgivenRandRTh(rho, rhotheta, qv);const Real P=getPgivenRTh(rhotheta, qv);const Real rand_double=amrex::Random(engine);const Real Tnew=Told+T_0_Pert_Mag (one - two rand_double) (one - z_agl/pert_ref_height);const Real theta_new=getThgivenTandP(Tnew, P, rdOcp);const Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);if(const_rho) { state_pert(i, j, k, Rho_comp)=zero;state_pert(i, j, k, RhoTheta_comp)=r_hse(i, j, k) (theta_new -(rhotheta/rho));} else { state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=zero;} } state_pert(i, j, k, RhoScalar_comp)=A_0 exp(-amrex::Real(10.0) r r);if(custom_TKE) { state_pert(i, j, k, RhoKE_comp)=(one - z/phys_z_hi) r_hse(i, j, k);} else { state_pert(i, j, k, RhoKE_comp)=KE_0;} if(use_moisture) { state_pert(i, j, k, RhoQ1_comp)=zero;state_pert(i, j, k, RhoQ2_comp)=zero;if(z_agl<=pert_ref_height &&qv_0_Pert_Mag !=zero) { const Real rhoold=state(i, j, k, Rho_comp);const Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);const Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;const Real rand_double=amrex::Random(engine);const Real qvnew=qvold+(two rand_double - one) qv_0_Pert_Mag;state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} } })

Variables
const Real	phys_z_lo = phys_zlo(geomdata, lev)

const Real	phys_z_hi = phys_zhi(geomdata, lev)

Real	A_0 = amrex::Real(1.0)

Real	KE_0 = amrex::Real(0.1)

Real	T_0_Pert_Mag = amrex::Real(0.0)

Real	qv_0_Pert_Mag = amrex::Real(0.0)

Real	pert_ref_height = amrex::Real(200.0)

bool	custom_TKE = false

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

const bool	use_terrain = (sc.terrain_type != TerrainType::None)

const bool	const_rho = (sc.fixed_density[lev] == 1)

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+myhalf) dx[0];const Real y=prob_lo[1]+(j+myhalf) dx[1];const Real z=use_terrain ? z_cc(i, j, k) :prob_lo[2]+(k+myhalf) dx[2];const Real z_agl=z - erf_prob_heights::z_ground_cc(i, j, z_sfc_nd);const Real xc=myhalf (prob_lo[0]+prob_hi[0]);const Real yc=myhalf (prob_lo[1]+prob_hi[1]);const Real zc=myhalf (phys_z_lo+phys_z_hi);const Real r=std::sqrt((x-xc) (x-xc)+(y-yc) (y-yc)+(z-zc) (z-zc));if(pert_ref_height > zero &&z_agl<=pert_ref_height &&T_0_Pert_Mag !=zero) { const Real rhotheta=state(i, j, k, RhoTheta_comp);const Real rho=state(i, j, k, Rho_comp);const Real qv=(use_moisture) ? state(i, j, k, RhoQ1_comp)/rho :zero;const Real Told=getTgivenRandRTh(rho, rhotheta, qv);const Real P=getPgivenRTh(rhotheta, qv);const Real rand_double=amrex::Random(engine);const Real Tnew=Told+T_0_Pert_Mag (one - two rand_double) (one - z_agl/pert_ref_height);const Real theta_new=getThgivenTandP(Tnew, P, rdOcp);const Real rhonew=getRhogivenThetaPress(theta_new, P, rdOcp, qv);if(const_rho) { state_pert(i, j, k, Rho_comp)=zero;state_pert(i, j, k, RhoTheta_comp)=r_hse(i, j, k) (theta_new -(rhotheta/rho));} else { state_pert(i, j, k, Rho_comp)=rhonew - rho;state_pert(i, j, k, RhoTheta_comp)=zero;} } state_pert(i, j, k, RhoScalar_comp)=A_0 exp(-amrex::Real(10.0) r r);if(custom_TKE) { state_pert(i, j, k, RhoKE_comp)=(one - z/phys_z_hi) r_hse(i, j, k);} else { state_pert(i, j, k, RhoKE_comp)=KE_0;} if(use_moisture) { state_pert(i, j, k, RhoQ1_comp)=zero;state_pert(i, j, k, RhoQ2_comp)=zero;if(z_agl<=pert_ref_height &&qv_0_Pert_Mag !=zero) { const Real rhoold=state(i, j, k, Rho_comp);const Real rhonew=rhoold+state_pert(i, j, k, Rho_comp);const Real qvold=state(i, j, k, RhoQ1_comp)/rhoold;const Real rand_double=amrex::Random(engine);const Real qvnew=qvold+(two rand_double - one) qv_0_Pert_Mag;state_pert(i, j, k, RhoQ1_comp)=rhonew qvnew - rhoold qvold;} } }
	)