ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Functions | Variables
ERF_UpdateRhoThetaSources_SDMCongestus3D.H File Reference

Go to the source code of this file.

Functions

ParmParse pp_prob ("prob")
 
pp_prob query ("advection_heating_rate", advection_heating_rate)
 
pp_prob query ("advection_heating_rate_base", advection_heating_rate_base)
 
 for (MFIter mfi(*src, TilingIfNotGPU());mfi.isValid();++mfi)
 

Variables

Real advection_heating_rate = amrex::Real(5.4e-5)
 
Real advection_heating_rate_base = amrex::Real(5.4e-5)
 
auto prob_lo = geom.ProbLoArray()
 
auto prob_hi = geom.ProbHiArray()
 
auto dx = geom.CellSizeArray()
 
const Real xc = myhalf * (prob_lo[0] + prob_hi[0])
 
const Real yc = myhalf * (prob_lo[1] + prob_hi[1])
 

Function Documentation

◆ for()

for ( MFIter mfi(*,  TilingIfNotGPU)
16  {
17  const auto &box = mfi.tilebox();
18  const Array4<Real>& src_arr = src->array(mfi);
19  if (box.length(0) != 1)
20  {
21  // if z dimension size is 1, then src is a spatially varying function over x,y at k=0
22  ParallelFor( box, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
23  {
24  auto x = prob_lo[0] + (i + myhalf) * dx[0];
25  auto y = prob_lo[1] + (j + myhalf) * dx[1];
26  auto z = prob_lo[2] + (k + myhalf) * dx[2];
27 
28  auto c = amrex::Real(2890000.0);
29  auto r = std::sqrt((x-xc)*(x-xc) + (y-yc)*(y-yc));
30  auto rsqr = r*r;
31  if (time < 3600) {
32  src_arr(i, j, k) = (advection_heating_rate_base)*exp(-z/100);
33  } else {
34  src_arr(i, j, k) = (advection_heating_rate*exp(-rsqr / c))*exp(-z/100);
35  }
36  });
37  } else {
38  // src is a function over Z
39  ParallelFor(box, [=] AMREX_GPU_DEVICE (int i, int j, int k) {
40  src_arr(i, j, k) = advection_heating_rate;
41  });
42  }
43  }
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:11
Coord::y
@ y
Coord::x
@ x
Coord::z
@ z
ParallelFor
ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const Real *dx=geomdata.CellSize();const Real x=(i+0.5) *dx[0];const Real y=(j+0.5) *dx[1];const Real Omg=erf_vortex_Gaussian(x, y, xc, yc, R, beta, sigma);const Real deltaT=-(gamma - 1.0)/(2.0 *sigma *sigma) *Omg *Omg;const Real rho_norm=std::pow(1.0+deltaT, inv_gm1);const Real T=(1.0+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)=0.25 *(1.0+std::cos(PI *std::min(r2d_xy, R)/R));})
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
prob_lo
auto prob_lo
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:6
advection_heating_rate
Real advection_heating_rate
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:3
yc
const Real yc
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:12
advection_heating_rate_base
Real advection_heating_rate_base
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:4
xc
const Real xc
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:11
dx
auto dx
Definition: ERF_UpdateRhoThetaSources_SDMCongestus3D.H:8
Here is the call graph for this function:

◆ pp_prob()

ParmParse pp_prob ( "prob"  )

◆ query() [1/2]

pp_prob query ( "advection_heating_rate"  ,
advection_heating_rate   
)

◆ query() [2/2]

pp_prob query ( "advection_heating_rate_base"  ,
advection_heating_rate_base   
)

Variable Documentation

◆ advection_heating_rate

Real advection_heating_rate = amrex::Real(5.4e-5)

Referenced by for().

◆ advection_heating_rate_base

Real advection_heating_rate_base = amrex::Real(5.4e-5)

Referenced by for().

◆ dx

auto dx = geom.CellSizeArray()

Referenced by for().

◆ prob_hi

auto prob_hi = geom.ProbHiArray()

◆ prob_lo

auto prob_lo = geom.ProbLoArray()

Referenced by for().

◆ xc

const Real xc = myhalf * (prob_lo[0] + prob_hi[0])

Referenced by for().

◆ yc

const Real yc = myhalf * (prob_lo[1] + prob_hi[1])

Referenced by for().