ERF_InitCustomPert_IsentropicVortex.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp ("prob")
pp	query ("p_inf", p_inf)
pp	query ("T_inf", T_inf)
pp	query ("M_inf", M_inf)
pp	query ("alpha", alpha)
pp	query ("gamma", gamma)
pp	query ("beta", beta)
pp	query ("sigma", sigma)
pp	query ("R", R)
pp	query ("xc", xc_frac)
pp	query ("yc", yc_frac)
	ParallelFor (bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const Real *dx=geomdata.CellSize();const Real x=(i+myhalf) *dx[0];const Real y=(j+myhalf) *dx[1];const Real Omg=erf_vortex_Gaussian(x, y, xc, yc, R, beta, sigma);const Real deltaT=-(gamma - one)/(two *sigma *sigma) *Omg *Omg;const Real rho_norm=std::pow(one+deltaT, inv_gm1);state_pert(i, j, k, Rho_comp)=rho_norm *rho_0 - r_hse(i, j, k);const Real T=(one+deltaT) *T_inf;const Real p=std::pow(rho_norm, Gamma)/Gamma *rho_0 *a_inf *a_inf;const Real rho_theta=rho_0 *rho_norm *(T *std::pow(p_0/p, rdOcp));state_pert(i, j, k, RhoTheta_comp)=rho_theta - getRhoThetagivenP(p_hse(i, j, k));const Real r2d_xy=std::sqrt((x-xc) *(x-xc)+(y-yc) *(y-yc));state_pert(i, j, k, RhoScalar_comp)=fourth *(one+std::cos(PI *std::min(r2d_xy, R)/R));})

Variables
amrex::Real	p_inf = p_0
amrex::Real	T_inf = amrex::Real(300.0)
amrex::Real	M_inf = amrex::Real(0.2)
amrex::Real	alpha = zero
amrex::Real	gamma = Gamma
amrex::Real	beta = amrex::Real(0.01)
amrex::Real	sigma = one
amrex::Real	R = two
amrex::Real	xc_frac = myhalf
amrex::Real	yc_frac = myhalf
const amrex::Real *	prob_lo = geomdata.ProbLo()
const amrex::Real *	prob_hi = geomdata.ProbHi()
amrex::Real	xc = xc_frac * (prob_hi[0] - prob_lo[0])
amrex::Real	yc = yc_frac * (prob_hi[1] - prob_lo[1])
amrex::Real	inv_gm1 = one / (gamma - one)
amrex::Real	rho_0 = p_inf / (R_d * T_inf)
amrex::Real	a_inf = std::sqrt(gamma * R_d * T_inf)
const amrex::Real	rdOcp = sc.rdOcp

Function Documentation

ParallelFor()

ParallelFor	(	bx	,
		[=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const Real *dx=geomdata.CellSize();const Real x=(i+myhalf) *dx[0];const Real y=(j+myhalf) *dx[1];const Real Omg=erf_vortex_Gaussian(x, y, xc, yc, R, beta, sigma);const Real deltaT=-(gamma - one)/(two *sigma *sigma) *Omg *Omg;const Real rho_norm=std::pow(one+deltaT, inv_gm1);state_pert(i, j, k, Rho_comp)=rho_norm *rho_0 - r_hse(i, j, k);const Real T=(one+deltaT) *T_inf;const Real p=std::pow(rho_norm, Gamma)/Gamma *rho_0 *a_inf *a_inf;const Real rho_theta=rho_0 *rho_norm *(T *std::pow(p_0/p, rdOcp));state_pert(i, j, k, RhoTheta_comp)=rho_theta - getRhoThetagivenP(p_hse(i, j, k));const Real r2d_xy=std::sqrt((x-xc) *(x-xc)+(y-yc) *(y-yc));state_pert(i, j, k, RhoScalar_comp)=fourth *(one+std::cos(PI *std::min(r2d_xy, R)/R));}
	)

Here is the call graph for this function:

pp()

ParmParse pp

(

"prob"

)

query() [1/10]

pp query	(	"alpha"	,
		alpha
	)

query() [2/10]

pp query	(	"beta"	,
		beta
	)

query() [3/10]

pp query	(	"gamma"	,
		gamma
	)

query() [4/10]

pp query	(	"M_inf"	,
		M_inf
	)

query() [5/10]

pp query	(	"p_inf"	,
		p_inf
	)

query() [6/10]

pp query	(	"R"	,
		R
	)

query() [7/10]

pp query	(	"sigma"	,
		sigma
	)

query() [8/10]

pp query	(	"T_inf"	,
		T_inf
	)

query() [9/10]

pp query	(	"xc"	,
		xc_frac
	)

query() [10/10]

pp query	(	"yc"	,
		yc_frac
	)

Variable Documentation

a_inf

amrex::Real a_inf = std::sqrt(gamma * R_d * T_inf)

Referenced by ParallelFor().

alpha

amrex::Real alpha = zero

beta

amrex::Real beta = amrex::Real(0.01)

Referenced by TurbulentPerturbation::calc_tpi_amp(), erf_vortex_Gaussian(), interpolate_1d(), orbital_params(), ParallelFor(), ERF::poisson_wall_dist(), and ERF::project_velocity_tb().

gamma

amrex::Real gamma = Gamma

Referenced by HSEutils::compute_dewpoint_temperature(), and ParallelFor().

inv_gm1

amrex::Real inv_gm1 = one / (gamma - one)

Referenced by ParallelFor().

M_inf

amrex::Real M_inf = amrex::Real(0.2)

p_inf

amrex::Real p_inf = p_0

prob_hi

const amrex::Real* prob_hi = geomdata.ProbHi()

Referenced by if(), ParallelFor(), ParallelForRNG(), and writeNCPlotFile().

prob_lo

const amrex::Real* prob_lo = geomdata.ProbLo()

Referenced by ComputeGlobalMinLocation(), ForestDrag::define_drag_field(), ERF::ErrorEst(), ERF::FillForecastStateMultiFabs(), ERF::FillSurfaceStateMultiFabs(), ERF::FindInitialEye(), HurricaneEyeTrackerNotInitial(), HurricaneMaxVelTracker(), HurricaneMinPressureTracker(), if(), init_bx_scalars_from_input_sounding(), init_bx_scalars_from_input_sounding_hse(), init_bx_velocities_from_input_sounding(), ParallelFor(), ParallelForRNG(), ERF::poisson_wall_dist(), ERF::setSpongeRefFromSounding(), tag_on_distance_from_eye(), thinbody_wall_dist(), and writeNCPlotFile().

R

amrex::Real R = two

Referenced by erf_vortex_Gaussian(), and ParallelFor().

rdOcp

const amrex::Real rdOcp = sc.rdOcp

Referenced by ParallelFor().

rho_0

amrex::Real rho_0 = p_inf / (R_d * T_inf)

Referenced by ParallelFor().

sigma

amrex::Real sigma = one

Referenced by ERF::apply_gaussian_smoothing_to_perturbations(), TerminalVelocity< RT >::CloudRainShima(), MaterialPropertiesCore::coeffMolecularDiffusion(), erf_vortex_Gaussian(), SDInitProperties::getDistribution(), and ParallelFor().

T_inf

amrex::Real T_inf = amrex::Real(300.0)

Referenced by ParallelFor().

xc

amrex::Real xc = xc_frac * (prob_hi[0] - prob_lo[0])

Referenced by ParallelFor().

xc_frac

amrex::Real xc_frac = myhalf

yc

amrex::Real yc = yc_frac * (prob_hi[1] - prob_lo[1])

Referenced by ParallelFor().

yc_frac

amrex::Real yc_frac = myhalf