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

Go to the source code of this file.

Functions

ParmParse pp_for_pert_vels ("erf")
 
pp_for_pert_vels query_enum_case_insensitive ("init_type", init_type)
 
ParmParse pp_for_pert_vels_prob ("prob")
 
pp_for_pert_vels_prob query ("rho_0", rho_0)
 
 if (init_type==InitType::Uniform)
 

Variables

InitType init_type
 
Real rho_0 = one
 

Function Documentation

◆ if()

if ( init_type  = = InitType::Uniform)
10  {
11  Real rot_time_period;
12  pp_for_pert_vels.get("rotational_time_period", rot_time_period);
13  Real coriolis_factor = amrex::Real(4.0) * PI / rot_time_period;
14 
15  Real Az;
16  pp_for_pert_vels.get("dynamic_viscosity", Az); // dynamic viscosity [kg-m/s]
17  Az = Az / rho_0; // kinematic viscosity [m^2/s]
18 
19  Vector<Real> abl_geo_wind(3);
20  pp_for_pert_vels.queryarr("abl_geo_wind",abl_geo_wind);
21  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(
22  (amrex::Math::abs(abl_geo_wind[1]) < amrex::Real(1.0e-15)) &&
23  (amrex::Math::abs(abl_geo_wind[2]) < amrex::Real(1.0e-15)),
24  "Ekman Spiral uses geostrophic forcing of the form (V_0, 0, 0)");
25  const Real u_0 = abl_geo_wind[0];
26 
27  const Real DE = std::sqrt(two * Az / coriolis_factor);
28  const Real a = one / DE;
29 
30  // Set the x-velocity
31  ParallelFor(xbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
32  const auto dx = geomdata.CellSize();
33  const Real z = (k + myhalf) * dx[2];
34 
35  // Set the x-velocity
36  x_vel_pert(i, j, k) = u_0 * (one - std::exp(-a * z) * std::cos(-a * z));
37  });
38 
39  // Set the y-velocity
40  ParallelFor(ybx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
41  const auto dx = geomdata.CellSize();
42  const Real z = (k + myhalf) * dx[2];
43 
44  // Set the y-velocity
45  y_vel_pert(i, j, k) = -u_0 * std::exp(-a * z) * std::sin(-a * z);
46  });
47 
48  // Set the z-velocity
49  ParallelFor(zbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept {
50  z_vel_pert(i, j, k) = zero;
51  });
52  } // end if uniform
two
constexpr amrex::Real two
Definition: ERF_Constants.H:8
one
constexpr amrex::Real one
Definition: ERF_Constants.H:7
zero
constexpr amrex::Real zero
Definition: ERF_Constants.H:6
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:11
PI
constexpr amrex::Real PI
Definition: ERF_Constants.H:16
Coord::z
@ z
rho_0
Real rho_0
Definition: ERF_InitCustomPertVels_EkmanSpiral.H:7
pp_for_pert_vels
ParmParse pp_for_pert_vels("erf")
dx
const Real dx
Definition: ERF_InitCustomPert_ABL.H:23
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));})
zbx
const Box zbx
Definition: ERF_SetupDiff.H:9
xbx
const Box xbx
Definition: ERF_SetupDiff.H:7
ybx
const Box ybx
Definition: ERF_SetupDiff.H:8
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
Here is the call graph for this function:

◆ pp_for_pert_vels()

ParmParse pp_for_pert_vels ( "erf"  )

Referenced by if().

Here is the caller graph for this function:

◆ pp_for_pert_vels_prob()

ParmParse pp_for_pert_vels_prob ( "prob"  )

◆ query()

pp_for_pert_vels_prob query ( "rho_0"  ,
rho_0   
)

◆ query_enum_case_insensitive()

pp_for_pert_vels query_enum_case_insensitive ( "init_type"  ,
init_type   
)

Variable Documentation

◆ init_type

InitType init_type

◆ rho_0

Real rho_0 = one

Referenced by if().