ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
ERF_UpdateRhoThetaSources_RICO.H
Go to the documentation of this file.
1  ParmParse pp_prob("prob");
2 
3  Real advection_heating_rate = zero; pp_prob.query("advection_heating_rate", advection_heating_rate);
4 
5  // If the z coordinate varies in time and or space, then the the height
6  // needs to be calculated at each time step. Here, we assume that only
7  // grid stretching exists.
8 
9  // Only apply temperature source below nominal inversion height
10  for ( amrex::MFIter mfi(*src, amrex::TilingIfNotGPU()); mfi.isValid(); ++mfi )
11  {
12  const auto &box = mfi.tilebox();
13  const Array4<Real>& src_arr = src->array(mfi);
14  if (box.length(0) == 1) {
15  ParallelFor(box, [=] AMREX_GPU_DEVICE (int i, int j, int k) {
16  src_arr(i, j, k) = advection_heating_rate;
17  });
18  }
19  }
zero
constexpr amrex::Real zero
Definition: ERF_Constants.H:6
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
advection_heating_rate
Real advection_heating_rate
Definition: ERF_UpdateRhoThetaSources_RICO.H:3
pp_prob
ParmParse pp_prob("prob")