ERF_InitCustomPert_ParticleTests.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp ("prob")
pp	query ("x_c", x_c)
pp	query ("z_c", z_c)
pp	query ("x_r", x_r)
pp	query ("z_r", z_r)
pp	query ("T_pert", T_pert)
	ParallelFor (bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+myhalf) *dx[0];const Real z=z_cc(i, j, k);const Real Tbar_hse=p_hse(i, j, k)/(R_d *r_hse(i, j, k));Real L=std::sqrt(std::pow((x - x_c)/x_r, 2)+std::pow((z - z_c)/z_r, 2));if(L<=one) { Real dT=T_pert *(std::cos(PI *L)+one)/two;Real theta_perturbed=(Tbar_hse+dT) *std::pow(p_0/p_hse(i, j, k), R_d/C_p);state_pert(i, j, k, Rho_comp)=getRhoThetagivenP(p_hse(i, j, k))/theta_perturbed - r_hse(i, j, k);} })

Variables
Real	x_c = zero
Real	z_c = amrex::Real(3000.0)
Real	x_r = amrex::Real(4000.0)
Real	z_r = amrex::Real(2000.0)
Real	T_pert = -amrex::Real(15.0)
Real	C_p = amrex::Real(1004.0)

Function Documentation

ParallelFor()

ParallelFor	(	bx	,
		[=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+myhalf) *dx[0];const Real z=z_cc(i, j, k);const Real Tbar_hse=p_hse(i, j, k)/(R_d *r_hse(i, j, k));Real L=std::sqrt(std::pow((x - x_c)/x_r, 2)+std::pow((z - z_c)/z_r, 2));if(L<=one) { Real dT=T_pert *(std::cos(PI *L)+one)/two;Real theta_perturbed=(Tbar_hse+dT) *std::pow(p_0/p_hse(i, j, k), R_d/C_p);state_pert(i, j, k, Rho_comp)=getRhoThetagivenP(p_hse(i, j, k))/theta_perturbed - r_hse(i, j, k);} }
	)

Here is the call graph for this function:

pp()

ParmParse pp

(

"prob"

)

query() [1/5]

pp query	(	"T_pert"	,
		T_pert
	)

query() [2/5]

pp query	(	"x_c"	,
		x_c
	)

query() [3/5]

pp query	(	"x_r"	,
		x_r
	)

query() [4/5]

pp query	(	"z_c"	,
		z_c
	)

query() [5/5]

pp query	(	"z_r"	,
		z_r
	)

Variable Documentation

C_p

Real C_p = amrex::Real(1004.0)

Referenced by ParallelFor().

T_pert

Real T_pert = -amrex::Real(15.0)

Referenced by ParallelFor().

x_c

Real x_c = zero

Referenced by ParallelFor().

x_r

Real x_r = amrex::Real(4000.0)

Referenced by ParallelFor().

z_c

Real z_c = amrex::Real(3000.0)

Referenced by ParallelFor().

z_r

Real z_r = amrex::Real(2000.0)

Referenced by ParallelFor().