ERF_InitCustomPert_SuperCell.H File Reference

Go to the source code of this file.

Functions
	AMREX_ALWAYS_ASSERT (bx.length()[2]==khi+1)
ParmParse	pp_prob ("prob")
pp_prob	query ("qt_init", q_t)
pp_prob	query ("eq_pot_temp", eq_pot_temp)
pp_prob	query ("use_empirical_psat", use_empirical)
Vector< Real >	h_r (khi+2)
Vector< Real >	h_p (khi+2)
Vector< Real >	h_t (khi+2)
Vector< Real >	h_q_v (khi+2)
amrex::Gpu::DeviceVector< Real >	d_r (khi+2)
amrex::Gpu::DeviceVector< Real >	d_p (khi+2)
amrex::Gpu::DeviceVector< Real >	d_t (khi+2)
amrex::Gpu::DeviceVector< Real >	d_q_v (khi+2)
pp_prob	query ("height", height)
pp_prob	query ("z_tr", z_tr)
pp_prob	query ("x_c", x_c)
pp_prob	query ("z_c", z_c)
pp_prob	query ("x_r", x_r)
pp_prob	query ("y_r", y_r)
pp_prob	query ("z_r", z_r)
pp_prob	query ("T_tr", T_tr)
pp_prob	query ("theta_tr", theta_tr)
pp_prob	query ("theta_c", theta_c)
pp_prob	query ("theta_0", theta_0)
	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+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];Real rad, delta_theta, theta_total, temperature, rho, RH, scalar;rad=std::pow(std::pow((x - x_c)/x_r, 2)+std::pow((y - y_c)/y_r, 2)+std::pow((z - z_c)/z_r, 2), 0.5);if(rad<=r_c) { delta_theta=theta_c *std::pow(cos(PI *rad/2.0), 2);scalar=std::pow(cos(PI *rad/2.0), 2);} else { delta_theta=0.0;scalar=0.0;} Real scaled_height=z/z_tr;int which_zone=-1;if(z<=height) { which_zone=1;} else if(z<=z_tr) { which_zone=2;} else { which_zone=3;} theta_total=t[k]+delta_theta;temperature=getTgivenPandTh(p[k], theta_total,(R_d/Cp_d));Real T_b=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));RH=HSEutils::compute_relative_humidity(p[k], T_b, use_empirical, which_zone, scaled_height);Real q_v_hot;if(z< height) { q_v_hot=0.014;} else { q_v_hot=HSEutils::vapor_mixing_ratio(p[k], T_b, RH, use_empirical, which_zone);} rho=p[k]/(R_d *temperature *(1.0+(R_v/R_d) *q_v_hot));Real temperature_back=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));Real T_back=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));Real RH_back=HSEutils::compute_relative_humidity(p[k], T_back, use_empirical, which_zone, scaled_height);Real q_v_back;if(z< height) { q_v_back=0.014;} else { q_v_back=HSEutils::vapor_mixing_ratio(p[k], T_back, RH_back, use_empirical, which_zone);} Real rho_back=getRhogivenTandPress(temperature_back, p[k], q_v_back);state_pert(i, j, k, RhoTheta_comp)=rho *theta_total - rho_back *t[k] *(1.0+(R_v/R_d) *q_v_back);state_pert(i, j, k, Rho_comp)=rho - rho_back *(1.0+q_v_back);state_pert(i, j, k, RhoScalar_comp)=rho *scalar;if(use_moisture) { state_pert(i, j, k, RhoQ1_comp)=rho *q_v_hot;state_pert(i, j, k, RhoQ2_comp)=0.0;} })

Variables
const bool	use_moisture = (sc.moisture_type != MoistureType::None)
const int	khi = geomdata.Domain().bigEnd()[2]
Real	q_t = 0.02
Real	eq_pot_temp = 320.0
bool	use_empirical = false
const amrex::Real &	dz = geomdata.CellSize()[2]
Real	height = 1200.
Real	z_tr = 12000.
Real	x_c = 0.0
Real	y_c = 0.0
Real	z_c = 1.5e3
Real	r_c = 1.0
Real	x_r = 10.0e3
Real	y_r = 10.0e3
Real	z_r = 1.5e3
Real	T_tr = 213.0
Real	theta_tr = 343.0
Real	theta_c = 3.0
Real	theta_0 = 300.0
Real *	t = d_t.data()
Real *	p = d_p.data()

Function Documentation

AMREX_ALWAYS_ASSERT()

AMREX_ALWAYS_ASSERT

(

bx.

length()[2] = =khi+1

)

d_p()

amrex::Gpu::DeviceVector<Real> d_p

(

khi+

2

)

d_q_v()

amrex::Gpu::DeviceVector<Real> d_q_v

(

khi+

2

)

d_r()

amrex::Gpu::DeviceVector<Real> d_r

(

khi+

2

)

d_t()

amrex::Gpu::DeviceVector<Real> d_t

(

khi+

2

)

h_p()

Vector<Real> h_p

(

khi+

2

)

h_q_v()

Vector<Real> h_q_v

(

khi+

2

)

h_r()

Vector<Real> h_r

(

khi+

2

)

h_t()

Vector<Real> h_t

(

khi+

2

)

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+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];Real rad, delta_theta, theta_total, temperature, rho, RH, scalar;rad=std::pow(std::pow((x - x_c)/x_r, 2)+std::pow((y - y_c)/y_r, 2)+std::pow((z - z_c)/z_r, 2), 0.5);if(rad<=r_c) { delta_theta=theta_c *std::pow(cos(PI *rad/2.0), 2);scalar=std::pow(cos(PI *rad/2.0), 2);} else { delta_theta=0.0;scalar=0.0;} Real scaled_height=z/z_tr;int which_zone=-1;if(z<=height) { which_zone=1;} else if(z<=z_tr) { which_zone=2;} else { which_zone=3;} theta_total=t[k]+delta_theta;temperature=getTgivenPandTh(p[k], theta_total,(R_d/Cp_d));Real T_b=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));RH=HSEutils::compute_relative_humidity(p[k], T_b, use_empirical, which_zone, scaled_height);Real q_v_hot;if(z< height) { q_v_hot=0.014;} else { q_v_hot=HSEutils::vapor_mixing_ratio(p[k], T_b, RH, use_empirical, which_zone);} rho=p[k]/(R_d *temperature *(1.0+(R_v/R_d) *q_v_hot));Real temperature_back=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));Real T_back=getTgivenPandTh(p[k], t[k],(R_d/Cp_d));Real RH_back=HSEutils::compute_relative_humidity(p[k], T_back, use_empirical, which_zone, scaled_height);Real q_v_back;if(z< height) { q_v_back=0.014;} else { q_v_back=HSEutils::vapor_mixing_ratio(p[k], T_back, RH_back, use_empirical, which_zone);} Real rho_back=getRhogivenTandPress(temperature_back, p[k], q_v_back);state_pert(i, j, k, RhoTheta_comp)=rho *theta_total - rho_back *t[k] *(1.0+(R_v/R_d) *q_v_back);state_pert(i, j, k, Rho_comp)=rho - rho_back *(1.0+q_v_back);state_pert(i, j, k, RhoScalar_comp)=rho *scalar;if(use_moisture) { state_pert(i, j, k, RhoQ1_comp)=rho *q_v_hot;state_pert(i, j, k, RhoQ2_comp)=0.0;} }
	)

Here is the call graph for this function:

pp_prob()

ParmParse pp_prob

(

"prob"

)

query() [1/14]

pp_prob query	(	"eq_pot_temp"	,
		eq_pot_temp
	)

query() [2/14]

pp_prob query	(	"height"	,
		height
	)

query() [3/14]

pp_prob query	(	"qt_init"	,
		q_t
	)

query() [4/14]

pp_prob query	(	"T_tr"	,
		T_tr
	)

query() [5/14]

pp_prob query	(	"theta_0"	,
		theta_0
	)

query() [6/14]

pp_prob query	(	"theta_c"	,
		theta_c
	)

query() [7/14]

pp_prob query	(	"theta_tr"	,
		theta_tr
	)

query() [8/14]

pp_prob query	(	"use_empirical_psat"	,
		use_empirical
	)

query() [9/14]

pp_prob query	(	"x_c"	,
		x_c
	)

query() [10/14]

pp_prob query	(	"x_r"	,
		x_r
	)

query() [11/14]

pp_prob query	(	"y_r"	,
		y_r
	)

query() [12/14]

pp_prob query	(	"z_c"	,
		z_c
	)

query() [13/14]

pp_prob query	(	"z_r"	,
		z_r
	)

query() [14/14]

pp_prob query	(	"z_tr"	,
		z_tr
	)

Variable Documentation

dz

const amrex::Real& dz = geomdata.CellSize()[2]

eq_pot_temp

Real eq_pot_temp = 320.0

height

Real height = 1200.

Referenced by ParallelFor().

khi

const int khi = geomdata.Domain().bigEnd()[2]

p

Real* p = d_p.data()

Referenced by ParallelFor().

q_t

Real q_t = 0.02

r_c

Real r_c = 1.0

Referenced by ParallelFor().

t

Real* t = d_t.data()

Referenced by ParallelFor().

T_tr

Real T_tr = 213.0

theta_0

Real theta_0 = 300.0

theta_c

Real theta_c = 3.0

Referenced by ParallelFor().

theta_tr

Real theta_tr = 343.0

use_empirical

bool use_empirical = false

Referenced by ParallelFor().

use_moisture

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

Referenced by ParallelFor().

x_c

Real x_c = 0.0

Referenced by ParallelFor().

x_r

Real x_r = 10.0e3

Referenced by ParallelFor().

y_c

Real y_c = 0.0

Referenced by ParallelFor().

y_r

Real y_r = 10.0e3

Referenced by ParallelFor().

z_c

Real z_c = 1.5e3

Referenced by ParallelFor().

z_r

Real z_r = 1.5e3

Referenced by ParallelFor().

z_tr

Real z_tr = 12000.

Referenced by ParallelFor().