docs/ERF__InitCustomPert__TaylorGreenVortex_8H_source.html

     ParmParse pp("prob");

     Real rho_0 =   1.0; pp.query("rho_0", rho_0);

     Real   T_0 = 300.0; pp.query("T_0", T_0);

     Real M_0   =   1.0; pp.query("M_0", M_0);

     Real V_0   =   1.0; pp.query("V_0", V_0);


     ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept

     {

         // Geometry

         const Real* prob_lo = geomdata.ProbLo();

         const Real* 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];


         // Initial potential temperature (actually rho*theta) perturbation

         const Real p = rho_0 * V_0*V_0*

                               ( 1.0 / (Gamma * M_0 * M_0)

                               + (1.0 / 16.0) * (cos(2 * x) + cos(2 * y)) * (cos(2 * z) + 2));

         state_pert(i, j, k, RhoTheta_comp) = getRhoThetagivenP(p) - getRhoThetagivenP(p_hse(i,j,k));


         // Set scalar = 0 everywhere

         state_pert(i, j, k, RhoScalar_comp) = 1.0 * rho_0;

     });

Gamma
constexpr amrex::Real Gamma
Definition: ERF_Constants.H:19

Coord::y
@ y

Coord::x
@ x

Coord::z
@ z

getRhoThetagivenP
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getRhoThetagivenP(const amrex::Real p, const amrex::Real qv=0.0)
Definition: ERF_EOS.H:172

RhoScalar_comp
#define RhoScalar_comp
Definition: ERF_IndexDefines.H:40

RhoTheta_comp
#define RhoTheta_comp
Definition: ERF_IndexDefines.H:37

dx
const Real dx
Definition: ERF_InitCustomPert_ABL.H:23

state_pert
state_pert(i, j, k, RhoTheta_comp)

prob_lo
const amrex::Real * prob_lo
Definition: ERF_InitCustomPert_IsentropicVortex.H:16

p
Real * p
Definition: ERF_InitCustomPert_SquallLine.H:61

M_0
Real M_0
Definition: ERF_InitCustomPert_TaylorGreenVortex.H:5

rho_0
Real rho_0
Definition: ERF_InitCustomPert_TaylorGreenVortex.H:3

pp
ParmParse pp("prob")

T_0
Real T_0
Definition: ERF_InitCustomPert_TaylorGreenVortex.H:4

ParallelFor
ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const Real *prob_lo=geomdata.ProbLo();const Real *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];const Real p=rho_0 *V_0 *V_0 *(1.0/(Gamma *M_0 *M_0)+(1.0/16.0) *(cos(2 *x)+cos(2 *y)) *(cos(2 *z)+2));state_pert(i, j, k, RhoTheta_comp)=getRhoThetagivenP(p) - getRhoThetagivenP(p_hse(i, j, k));state_pert(i, j, k, RhoScalar_comp)=1.0 *rho_0;})

V_0
Real V_0
Definition: ERF_InitCustomPert_TaylorGreenVortex.H:6

Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19