ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Functions
derived Namespace Reference

Functions

void erf_derrhodivide (const Box &bx, FArrayBox &derfab, const FArrayBox &datfab, const int scalar_index)
 
void erf_dernull (const Box &, FArrayBox &, int, int, const FArrayBox &, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dersoundspeed (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dertemp (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dermoisttemp (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dertheta (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_derscalar (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_derKE (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dervortx (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &geomdata, Real, const int *, const int)
 
void erf_dervorty (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &geomdata, Real, const int *, const int)
 
void erf_dervortz (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &, const Geometry &geomdata, Real, const int *, const int)
 
void erf_derenstrophysq (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &geomdata, Real, const int *, const int)
 
void erf_dermagvel (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dermagvelsq (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_derreflectivity (const Box &bx, FArrayBox &derfab, int dcomp, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_dermaxreflectivity (const Box &bx, FArrayBox &derfab, int dcomp, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &, Real, const int *, const int)
 
void erf_derlocalhelicity (const Box &bx, FArrayBox &derfab, int dcomp, int, const FArrayBox &datfab, const FArrayBox &, const Geometry &geomdata, Real, const int *, const int)
 
void erf_derhelicity (const Box &bx, FArrayBox &derfab, int dcomp, int, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &geomdata, Real, const int *, const int)
 
void erf_derprecipitable (const Box &bx, FArrayBox &derfab, int dcomp, int, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &, Real, const int *, const int)
 
void erf_dermucape (const Box &bx, FArrayBox &derfab, int dcomp, int ncomp, const FArrayBox &datfab, const FArrayBox &zcc_fab, const Geometry &, Real, const int *, const int)
 
void erf_derrhodivide (const amrex::Box &bx, amrex::FArrayBox &derfab, const amrex::FArrayBox &datfab, const int scalar_index)
 
void erf_dernull (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dersoundspeed (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dertemp (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dermoisttemp (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dertheta (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derscalar (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derKE (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dervortx (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dervorty (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dervortz (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derenstrophysq (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dermagvel (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dermagvelsq (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derreflectivity (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dermaxreflectivity (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derhelicity (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derlocalhelicity (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_derprecipitable (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 
void erf_dermucape (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::FArrayBox &zfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
 

Function Documentation

◆ erf_derenstrophysq() [1/2]

void derived::erf_derenstrophysq ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derenstrophysq() [2/2]

void derived::erf_derenstrophysq ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
453 {
454  AMREX_ALWAYS_ASSERT(dcomp == 0);
455  AMREX_ALWAYS_ASSERT(ncomp == 1);
456 
457  auto const dat = datfab.array(); // cell-centered velocity
458  auto tfab = derfab.array(); // cell-centered vorticity x-component
459  auto z_arr = zcc_fab.array(); // cell-centered height z
460 
461  const Real two_dx = two * geomdata.CellSize(0);
462  const Real two_dy = two * geomdata.CellSize(1);
463 
464  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
465  {
466  Real two_dz = z_arr(i,j,k+1) - z_arr(i,j,k-1);
467 
468  Real vortx = (dat(i,j+1,k,2) - dat(i,j-1,k,2)) / two_dy // dw/dy
469  -(dat(i,j,k+1,1) - dat(i,j,k-1,1)) / two_dz; // dv/dz
470  Real vorty = (dat(i,j,k+1,0) - dat(i,j,k-1,0)) / two_dz // du/dz
471  -(dat(i+1,j,k,2) - dat(i-1,j,k,2)) / two_dx; // dw/dx
472  Real vortz = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / two_dx // dv/dx
473  -(dat(i,j+1,k,0) - dat(i,j-1,k,0)) / two_dy; // du/dy
474 
475  tfab(i,j,k,dcomp) = vortx*vortx + vorty*vorty + vortz*vortz;
476  });
477 }
two
constexpr amrex::Real two
Definition: ERF_Constants.H:10
AMREX_ALWAYS_ASSERT
AMREX_ALWAYS_ASSERT(bx.length()[2]==khi+1)
ParallelFor
ParallelFor(grown_box, [=] AMREX_GPU_DEVICE(int i, int j, int k) { qrcuten_arr(i, j, k)=Real(0);qscuten_arr(i, j, k)=Real(0);qicuten_arr(i, j, k)=Real(0);})
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19

Referenced by ERF::sum_derived_quantities().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derhelicity() [1/2]

void derived::erf_derhelicity ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derhelicity() [2/2]

void derived::erf_derhelicity ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
657 {
658  AMREX_ALWAYS_ASSERT(dcomp == 0);
659  AMREX_ALWAYS_ASSERT(bx.smallEnd(2) == 0);
660 
661  auto const dat = datfab.array(); // cell-centered velocity
662  auto dfab = derfab.array(); // integral of local helicity
663  auto z_arr = zcc_fab.array(); // cell-centered height z
664 
665  const Real dx = geomdata.CellSize(0);
666  const Real dy = geomdata.CellSize(1);
667 
668  // Collapse to i,j box (ignore vertical for now)
669  Box b2d = makeSlab(bx,2,0);
670 
671  ParallelFor(b2d, [=] AMREX_GPU_DEVICE(int i, int j, int ) noexcept
672  {
673  Real int_hel = Real(0.0);
674  for (int k = bx.smallEnd(2); k <= bx.bigEnd(2); ++k)
675  {
676  Real z = z_arr(i,j,k);
677 
678  // Helicity is defined as integral from 2km to 5km in vertical
679  if (z > Real(2000.0) && z < Real(5000.0)) {
680 
681  Real z_hi = myhalf * (z_arr(i,j,k) + z_arr(i,j,k+1));
682  Real z_lo = myhalf * (z_arr(i,j,k) + z_arr(i,j,k-1));
683  Real dz = z_hi - z_lo;
684 
685  Real vortz = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / (two*dx) // dv/dx
686  - (dat(i,j+1,k,0) - dat(i,j-1,k,0)) / (two*dy); // du/dy
687  Real w = dat(i,j,k,2); // vertical velocity
688 
689  int_hel += vortz * w * dz;
690  }
691  }
692 
693  // Store vertical integral into *all* levels for this (i,j)
694  for (int k = bx.smallEnd(2); k <= bx.bigEnd(2); ++k) {
695  dfab(i, j, k, dcomp) = int_hel;
696  }
697  });
698 }
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:13
Coord::z
@ z
dy
const Real dy
Definition: ERF_InitCustomPert_ABL.H:24
dx
const Real dx
Definition: ERF_InitCustomPert_ABL.H:23
WSM6SedCellScratch::dz
@ dz
Definition: ERF_AdvanceWSM6.cpp:104

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derKE() [1/2]

void derived::erf_derKE ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derKE() [2/2]

void derived::erf_derKE ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Function to define the kinetic energy KE by dividing (rho KE) by rho

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds KE @params[in] datfab array of data used to construct derived quantity

352 {
353  erf_derrhodivide(bx, derfab, datfab, RhoKE_comp);
354 }
RhoKE_comp
#define RhoKE_comp
Definition: ERF_IndexDefines.H:38
derived::erf_derrhodivide
void erf_derrhodivide(const Box &bx, FArrayBox &derfab, const FArrayBox &datfab, const int scalar_index)
Definition: ERF_Derive.cpp:168

Referenced by ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derlocalhelicity() [1/2]

void derived::erf_derlocalhelicity ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derlocalhelicity() [2/2]

void derived::erf_derlocalhelicity ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
626 {
627  AMREX_ALWAYS_ASSERT(dcomp == 0);
628 
629  auto const dat = datfab.array(); // cell-centered velocity
630  auto dfab = derfab.array(); // cell-centered local helicity
631 
632  const Real two_dx = Real(2.0)*geomdata.CellSize(0);
633  const Real two_dy = Real(2.0)*geomdata.CellSize(1);
634 
635  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
636  {
637  Real vortz = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / two_dx // dv/dx
638  - (dat(i,j+1,k,0) - dat(i,j-1,k,0)) / two_dy; // du/dy
639  Real w = dat(i,j,k,2);
640 
641  // Helicity
642  dfab(i,j,k,dcomp) = vortz * w;
643  });
644 }

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dermagvel() [1/2]

void derived::erf_dermagvel ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dermagvel() [2/2]

void derived::erf_dermagvel ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
490 {
491  AMREX_ALWAYS_ASSERT(dcomp == 0);
492  AMREX_ALWAYS_ASSERT(ncomp == 1);
493 
494  auto const dat = datfab.array(); // cell-centered velocity
495  auto tfab = derfab.array(); // cell-centered magvel
496 
497  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
498  {
499  Real u = dat(i,j,k,0);
500  Real v = dat(i,j,k,1);
501  Real w = dat(i,j,k,2);
502  tfab(i,j,k,dcomp) = std::sqrt(u*u + v*v + w*w);
503  });
504 }

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dermagvelsq() [1/2]

void derived::erf_dermagvelsq ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dermagvelsq() [2/2]

void derived::erf_dermagvelsq ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
517 {
518  AMREX_ALWAYS_ASSERT(dcomp == 0);
519  AMREX_ALWAYS_ASSERT(ncomp == 1);
520 
521  auto const dat = datfab.array(); // cell-centered velocity
522  auto tfab = derfab.array(); // cell-centered magvel
523 
524  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
525  {
526  Real u = dat(i,j,k,0);
527  Real v = dat(i,j,k,1);
528  Real w = dat(i,j,k,2);
529  tfab(i,j,k,dcomp) = u*u + v*v + w*w;
530  });
531 }

Referenced by ERF::sum_derived_quantities().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dermaxreflectivity() [1/2]

void derived::erf_dermaxreflectivity ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dermaxreflectivity() [2/2]

void derived::erf_dermaxreflectivity ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
576 {
577  AMREX_ALWAYS_ASSERT(dcomp == 0);
578  AMREX_ALWAYS_ASSERT(bx.smallEnd(2) == 0);
579 
580  auto const dat = datfab.array(); // cell-centered state vector
581  auto rfab = derfab.array(); // cell-centered max reflectivity
582 
583  // Collapse to i,j box (ignore vertical for now)
584  Box b2d = makeSlab(bx,2,0);
585 
586  Real l_bogus_large_value = bogus_large_value;
587 
588  ParallelFor(b2d, [=] AMREX_GPU_DEVICE(int i, int j, int) noexcept
589  {
590  Real max_dbz = -l_bogus_large_value;
591 
592  // find max reflectivity over k
593  for (int k = bx.smallEnd(2); k <= bx.bigEnd(2); ++k) {
594 
595  Real rho = dat(i,j,k,Rho_comp);
596  Real qv = std::max(Real(0.0),dat(i,j,k,RhoQ1_comp)/rho);
597  Real qpr = std::max(Real(0.0),dat(i,j,k,RhoQ4_comp)/rho);
598  Real qps = std::max(Real(0.0),dat(i,j,k,RhoQ5_comp)/rho);
599  Real qpg = std::max(Real(0.0),dat(i,j,k,RhoQ6_comp)/rho);
600 
601  Real temp = getTgivenRandRTh(rho, dat(i,j,k,RhoTheta_comp), qv);
602 
603  Real dbz = compute_max_reflectivity_dbz(rho, temp, qpr, qps, qpg,
604  1, 1, 1, 1);
605  max_dbz = amrex::max(max_dbz, dbz);
606  }
607 
608  // Store max_dbz into *all* levels for this (i,j)
609  for (int k = bx.smallEnd(2); k <= bx.bigEnd(2); ++k) {
610  rfab(i, j, k, dcomp) = max_dbz;
611  }
612  });
613 }
bogus_large_value
constexpr amrex::Real bogus_large_value
Definition: ERF_Constants.H:26
for
for(int i=0;i< m_num_species;i++)
Definition: ERF_InitCustomPert_MultiSpeciesBubble.H:48

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dermoisttemp() [1/2]

void derived::erf_dermoisttemp ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dermoisttemp() [2/2]

void derived::erf_dermoisttemp ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
276 {
277  auto const dat = datfab.array();
278  auto tfab = derfab.array();
279 
280  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
281  {
282  const Real rho = dat(i, j, k, Rho_comp);
283  const Real rhotheta = dat(i, j, k, RhoTheta_comp);
284  AMREX_ALWAYS_ASSERT(rhotheta > Real(0.0));
285  const Real qv = dat(i, j, k, RhoQ1_comp) / rho;
286  tfab(i,j,k) = getTgivenRandRTh(rho,rhotheta,qv);
287  });
288 }
getTgivenRandRTh
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getTgivenRandRTh(const amrex::Real rho, const amrex::Real rhotheta, const amrex::Real qv=amrex::Real(0))
Definition: ERF_EOS.H:46
Rho_comp
#define Rho_comp
Definition: ERF_IndexDefines.H:36
RhoTheta_comp
#define RhoTheta_comp
Definition: ERF_IndexDefines.H:37
RhoQ1_comp
#define RhoQ1_comp
Definition: ERF_IndexDefines.H:42
rho
rho
Definition: ERF_InitCustomPert_Bubble.H:106
MicVar_Kess::qv
@ qv
Definition: ERF_Kessler.H:29

Referenced by ERF::Write3DPlotFile(), WriteBndryPlanes::write_planes(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dermucape() [1/2]

void derived::erf_dermucape ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dermucape() [2/2]

void derived::erf_dermucape ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
756 {
757  AMREX_ALWAYS_ASSERT(dcomp == 0);
758  AMREX_ALWAYS_ASSERT(ncomp == 1);
759 
760  auto const dat = datfab.array();
761  auto dfab = derfab.array();
762  auto z_arr = zcc_fab.array();
763  const int ncons = datfab.nComp();
764 
765  Box b2d = bx;
766  b2d.setSmall(2,0);
767  b2d.setBig(2,0);
768 
769  ParallelFor(b2d, [=] AMREX_GPU_DEVICE(int i, int j, int) noexcept
770  {
771  Real mucape = Real(0);
772  int klo = bx.smallEnd(2);
773  int khi = bx.bigEnd(2);
774 
775  if (ncons > RhoQ1_comp) {
776  Real p_sfc = Real(0);
777  for (int k = klo; k <= khi; ++k) {
778  Real rho = dat(i,j,k,Rho_comp);
779  if (rho <= Real(0)) { continue; }
780  Real qv = amrex::max(Real(0), dat(i,j,k,RhoQ1_comp) / rho);
781  p_sfc = amrex::max(p_sfc, getPgivenRTh(dat(i,j,k,RhoTheta_comp), qv));
782  }
783 
784  if (p_sfc > mucape_search_depth_pa()) {
785  Real p_search_min = p_sfc - mucape_search_depth_pa();
786 
787  for (int ks = klo; ks < khi; ++ks) {
788  Real rho_src = dat(i,j,ks,Rho_comp);
789  if (rho_src <= Real(0)) { continue; }
790 
791  Real qv_src = amrex::max(Real(0), dat(i,j,ks,RhoQ1_comp) / rho_src);
792  Real rt_src = dat(i,j,ks,RhoTheta_comp);
793  Real p_src = getPgivenRTh(rt_src, qv_src);
794 
795  if (p_src < p_search_min) { continue; }
796 
797  Real T_src = getTgivenRandRTh(rho_src, rt_src, qv_src);
798  if (T_src <= Real(0)) { continue; }
799 
800  Real Td_src = mucape_dewpoint_temperature(p_src, qv_src, T_src);
801  Real Tlcl = mucape_lcl_temperature(T_src, Td_src);
802  Real plcl = p_src * std::pow(Tlcl / T_src, Cp_d / R_d);
803  plcl = amrex::min(p_src, amrex::max(plcl, mucape_min_pressure_pa()));
804 
805  Real theta_src = getThgivenTandP(T_src, p_src, R_d / Cp_d);
806 
807  Real candidate_cape = Real(0);
808  Real z_prev = z_arr(i,j,ks);
809  Real b_prev = Real(0);
810 
811  bool saturated = false;
812  Real T_sat_prev = Tlcl;
813  Real p_sat_prev = plcl;
814 
815  for (int k = ks + 1; k <= khi; ++k) {
816  Real rho_env = dat(i,j,k,Rho_comp);
817  if (rho_env <= Real(0)) { continue; }
818 
819  Real qv_env = amrex::max(Real(0), dat(i,j,k,RhoQ1_comp) / rho_env);
820  Real rt_env = dat(i,j,k,RhoTheta_comp);
821  Real p_env = getPgivenRTh(rt_env, qv_env);
822  Real T_env = getTgivenRandRTh(rho_env, rt_env, qv_env);
823  Real Tv_env = mucape_virtual_temperature(T_env, qv_env);
824  Real z_env = z_arr(i,j,k);
825 
826  Real T_parcel;
827  Real qv_parcel;
828 
829  if (p_env >= plcl) {
830  T_parcel = getTgivenPandTh(p_env, theta_src, R_d / Cp_d);
831  qv_parcel = qv_src;
832  } else {
833  if (!saturated) {
834  saturated = true;
835  T_sat_prev = Tlcl;
836  p_sat_prev = plcl;
837  }
838 
839  if (p_env < p_sat_prev) {
840  T_sat_prev = mucape_integrate_saturated_temperature(T_sat_prev, p_sat_prev, p_env);
841  p_sat_prev = p_env;
842  }
843 
844  T_parcel = T_sat_prev;
845  qv_parcel = mucape_qsat(T_parcel, p_env);
846  }
847 
848  Real Tv_parcel = mucape_virtual_temperature(T_parcel, qv_parcel);
849  Real buoyancy = CONST_GRAV * (Tv_parcel - Tv_env) /
850  amrex::max(Tv_env, mucape_min_temperature());
851 
852  candidate_cape += mucape_positive_area(b_prev, buoyancy, z_env - z_prev);
853  b_prev = buoyancy;
854  z_prev = z_env;
855  }
856 
857  mucape = amrex::max(mucape, candidate_cape);
858  }
859  }
860  }
861 
862  for (int k = klo; k <= khi; ++k) {
863  dfab(i,j,k,dcomp) = mucape;
864  }
865  });
866 }
Cp_d
constexpr amrex::Real Cp_d
Definition: ERF_Constants.H:44
CONST_GRAV
constexpr amrex::Real CONST_GRAV
Definition: ERF_Constants.H:53
R_d
constexpr amrex::Real R_d
Definition: ERF_Constants.H:42
getThgivenTandP
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getThgivenTandP(const amrex::Real T, const amrex::Real P, const amrex::Real rdOcp)
Definition: ERF_EOS.H:18
getPgivenRTh
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getPgivenRTh(const amrex::Real rhotheta, const amrex::Real qv=amrex::Real(0))
Definition: ERF_EOS.H:81
getTgivenPandTh
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getTgivenPandTh(const amrex::Real P, const amrex::Real th, const amrex::Real rdOcp)
Definition: ERF_EOS.H:32
khi
const int khi
Definition: ERF_InitCustomPert_Bubble.H:21

Referenced by ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dernull() [1/2]

void derived::erf_dernull ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dernull() [2/2]

void derived::erf_dernull ( const Box &  ,
FArrayBox &  ,
int  ,
int  ,
const FArrayBox &  ,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Placeholder function that does nothing

199 { }

Referenced by ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the caller graph for this function:

◆ erf_derprecipitable() [1/2]

void derived::erf_derprecipitable ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derprecipitable() [2/2]

void derived::erf_derprecipitable ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
711 {
712  AMREX_ALWAYS_ASSERT(dcomp == 0);
713  AMREX_ALWAYS_ASSERT(bx.smallEnd(2) == 0);
714 
715  auto const dat = datfab.array(); // cell-centered state vector
716  auto dfab = derfab.array(); // integral of qv to define precipitable water
717 
718  // Collapse to i,j box (ignore vertical for now)
719  Box b2d = makeSlab(bx,2,0);
720 
721  auto z_arr = zcc_fab.array(); // cell-centered height z
722 
723  ParallelFor(b2d, [=] AMREX_GPU_DEVICE(int i, int j, int) noexcept
724  {
725  Real integral_qv = Real(0.0);
726 
727  for (int k = 0; k <= bx.bigEnd(2); ++k)
728  {
729  Real z_hi = myhalf * (z_arr(i,j,k+1) + z_arr(i,j,k));
730  Real z_lo = myhalf * (z_arr(i,j,k-1) + z_arr(i,j,k));
731  Real dz = z_hi - z_lo;
732 
733  Real rhoQ1 = dat(i, j, k, RhoQ1_comp);
734 
735  integral_qv += rhoQ1 * dz;
736  }
737 
738  // Store vertical integral into *all* levels for this (i,j)
739  for (int k = bx.smallEnd(2); k <= bx.bigEnd(2); ++k) {
740  dfab(i, j, k, dcomp) = integral_qv;
741  }
742  });
743 }

Referenced by ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derreflectivity() [1/2]

void derived::erf_derreflectivity ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derreflectivity() [2/2]

void derived::erf_derreflectivity ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)
544 {
545  AMREX_ALWAYS_ASSERT(dcomp == 0);
546 
547  auto const dat = datfab.array(); // cell-centered state vector
548  auto rfab = derfab.array(); // cell-centered reflectivity
549 
550  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
551  {
552  Real rho = dat(i,j,k,Rho_comp);
553  Real qv = std::max(Real(0.0),dat(i,j,k,RhoQ1_comp)/rho);
554  Real qpr = std::max(Real(0.0),dat(i,j,k,RhoQ4_comp)/rho);
555  Real qps = std::max(Real(0.0),dat(i,j,k,RhoQ5_comp)/rho);
556  Real qpg = std::max(Real(0.0),dat(i,j,k,RhoQ6_comp)/rho);
557 
558  Real temp = getTgivenRandRTh(rho, dat(i,j,k,RhoTheta_comp), qv);
559 
560  rfab(i, j, k, dcomp) = compute_max_reflectivity_dbz(rho, temp, qpr, qps, qpg,
561  1, 1, 1, 1);
562  });
563 }
RhoQ4_comp
#define RhoQ4_comp
Definition: ERF_IndexDefines.H:45
RhoQ6_comp
#define RhoQ6_comp
Definition: ERF_IndexDefines.H:47
RhoQ5_comp
#define RhoQ5_comp
Definition: ERF_IndexDefines.H:46
compute_max_reflectivity_dbz
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real compute_max_reflectivity_dbz(amrex::Real rho_air, amrex::Real tmk, amrex::Real qra, amrex::Real qsn, amrex::Real qgr, int in0r, int in0s, int in0g, int iliqskin)
Definition: ERF_StormDiagnostics.H:13
MicVar_Morr::qpg
@ qpg
Definition: ERF_Morrison.H:41
MicVar_Morr::qps
@ qps
Definition: ERF_Morrison.H:40
MicVar_Morr::qpr
@ qpr
Definition: ERF_Morrison.H:39

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derrhodivide() [1/2]

void derived::erf_derrhodivide ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
const amrex::FArrayBox &  datfab,
const int  scalar_index 
)

◆ erf_derrhodivide() [2/2]

void derived::erf_derrhodivide ( const Box &  bx,
FArrayBox &  derfab,
const FArrayBox &  datfab,
const int  scalar_index 
)

Function to define a derived quantity by dividing by density (analogous to cons_to_prim)

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity @params[in] datfab array of data used to construct derived quantity @params[in] scalar_index index of quantity to be divided by density

172 {
173  // This routine divides any cell-centered conserved quantity by density
174  auto const dat = datfab.array();
175  auto primitive = derfab.array();
176 
177  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
178  {
179  const Real rho = dat(i, j, k, Rho_comp);
180  const Real conserved = dat(i, j, k, scalar_index);
181  primitive(i,j,k) = conserved / rho;
182  });
183 }

Referenced by erf_derKE(), erf_derscalar(), erf_dertheta(), and WriteBndryPlanes::write_planes().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_derscalar() [1/2]

void derived::erf_derscalar ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_derscalar() [2/2]

void derived::erf_derscalar ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Function to define a scalar s by dividing (rho s) by rho

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds scalar s @params[in] datfab array of data used to construct derived quantity

330 {
331  erf_derrhodivide(bx, derfab, datfab, RhoScalar_comp);
332 }
RhoScalar_comp
#define RhoScalar_comp
Definition: ERF_IndexDefines.H:40

Referenced by ERF::ErrorEst(), ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dersoundspeed() [1/2]

void derived::erf_dersoundspeed ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dersoundspeed() [2/2]

void derived::erf_dersoundspeed ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Function to define the sound speed by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

219 {
220  auto const dat = datfab.array();
221  auto cfab = derfab.array();
222 
223  // NOTE: we compute the soundspeed of dry air -- we do not account for any moisture effects here
224  Real qv = Real(0.0);
225 
226  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
227  {
228  const Real rhotheta = dat(i, j, k, RhoTheta_comp);
229  const Real rho = dat(i, j, k, Rho_comp);
230  AMREX_ALWAYS_ASSERT(rhotheta > 0);
231  cfab(i,j,k) = std::sqrt(Gamma * getPgivenRTh(rhotheta,qv) / rho);
232  });
233 }
Gamma
constexpr amrex::Real Gamma
Definition: ERF_Constants.H:51

Referenced by ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dertemp() [1/2]

void derived::erf_dertemp ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dertemp() [2/2]

void derived::erf_dertemp ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Function to define the temperature by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

253 {
254  auto const dat = datfab.array();
255  auto tfab = derfab.array();
256 
257  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
258  {
259  const Real rho = dat(i, j, k, Rho_comp);
260  const Real rhotheta = dat(i, j, k, RhoTheta_comp);
261  AMREX_ALWAYS_ASSERT(rhotheta > Real(0.0));
262  tfab(i,j,k) = getTgivenRandRTh(rho,rhotheta);
263  });
264 }

Referenced by ERF::Write3DPlotFile(), WriteBndryPlanes::write_planes(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dertheta() [1/2]

void derived::erf_dertheta ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dertheta() [2/2]

void derived::erf_dertheta ( const Box &  bx,
FArrayBox &  derfab,
int  ,
int  ,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  ,
Real  ,
const int *  ,
const int   
)

Function to define the potential temperature by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

308 {
309  erf_derrhodivide(bx, derfab, datfab, RhoTheta_comp);
310 }

Referenced by ERF::ErrorEst(), ERF::Write3DPlotFile(), and ERF::WriteSubvolume().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dervortx() [1/2]

void derived::erf_dervortx ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dervortx() [2/2]

void derived::erf_dervortx ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
367 {
368  AMREX_ALWAYS_ASSERT(dcomp == 0);
369  AMREX_ALWAYS_ASSERT(ncomp == 1);
370 
371  auto const dat = datfab.array(); // cell-centered velocity
372  auto tfab = derfab.array(); // cell-centered vorticity x-component
373  auto z_arr = zcc_fab.array(); // cell-centered height z
374 
375  const Real two_dy = two * geomdata.CellSize(1);
376 
377  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
378  {
379  Real two_dz = z_arr(i,j,k+1) - z_arr(i,j,k-1);
380  tfab(i,j,k,dcomp) = (dat(i,j+1,k,2) - dat(i,j-1,k,2)) / two_dy // dw/dy
381  - (dat(i,j,k+1,1) - dat(i,j,k-1,1)) / two_dz; // dv/dz
382  });
383 }

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dervorty() [1/2]

void derived::erf_dervorty ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dervorty() [2/2]

void derived::erf_dervorty ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  zcc_fab,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
396 {
397  AMREX_ALWAYS_ASSERT(dcomp == 0);
398  AMREX_ALWAYS_ASSERT(ncomp == 1);
399 
400  auto const dat = datfab.array(); // cell-centered velocity
401  auto tfab = derfab.array(); // cell-centered vorticity y-component
402  auto z_arr = zcc_fab.array(); // cell-centered height z
403 
404  const Real two_dx = two * geomdata.CellSize(0);
405 
406  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
407  {
408  Real two_dz = z_arr(i,j,k+1) - z_arr(i,j,k-1);
409  tfab(i,j,k,dcomp) = (dat(i,j,k+1,0) - dat(i,j,k-1,0)) / two_dz // du/dz
410  - (dat(i+1,j,k,2) - dat(i-1,j,k,2)) / two_dx; // dw/dx
411  });
412 }

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ erf_dervortz() [1/2]

void derived::erf_dervortz ( const amrex::Box &  bx,
amrex::FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const amrex::FArrayBox &  datfab,
const amrex::FArrayBox &  zfab,
const amrex::Geometry &  geomdata,
amrex::Real  time,
const int *  bcrec,
const int  level 
)

◆ erf_dervortz() [2/2]

void derived::erf_dervortz ( const Box &  bx,
FArrayBox &  derfab,
int  dcomp,
int  ncomp,
const FArrayBox &  datfab,
const FArrayBox &  ,
const Geometry &  geomdata,
Real  ,
const int *  ,
const int   
)
425 {
426  AMREX_ALWAYS_ASSERT(dcomp == 0);
427  AMREX_ALWAYS_ASSERT(ncomp == 1);
428 
429  auto const dat = datfab.array(); // cell-centered velocity
430  auto tfab = derfab.array(); // cell-centered vorticity z-component
431 
432  const Real dx = geomdata.CellSize(0);
433  const Real dy = geomdata.CellSize(2);
434 
435  ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
436  {
437  tfab(i,j,k,dcomp) = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / (two*dx) // dv/dx
438  - (dat(i,j+1,k,0) - dat(i,j-1,k,0)) / (two*dy); // du/dy
439  });
440 }

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:
Here is the caller graph for this function: