ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Functions
ERF_TimeAvgVel.cpp File Reference
#include <ERF_Utils.H>
Include dependency graph for ERF_TimeAvgVel.cpp:

Functions

void Time_Avg_Vel_atCC (const Real &dt, Real &t_avg_cnt, MultiFab *vel_t_avg, MultiFab &xvel, MultiFab &yvel, MultiFab &zvel)
 

Function Documentation

◆ Time_Avg_Vel_atCC()

void Time_Avg_Vel_atCC ( const Real &  dt,
Real &  t_avg_cnt,
MultiFab *  vel_t_avg,
MultiFab &  xvel,
MultiFab &  yvel,
MultiFab &  zvel 
)
15 {
16  // Augment the counter
17  t_avg_cnt += dt;
18 
19 #ifdef _OPENMP
20 #pragma omp parallel if (Gpu::notInLaunchRegion())
21 #endif
22  for ( MFIter mfi(*(vel_t_avg),TilingIfNotGPU()); mfi.isValid(); ++mfi)
23  {
24  // CC tilebox
25  Box tbx = mfi.tilebox();
26 
27  // Velocity on faces
28  const Array4<Real>& velx = xvel.array(mfi);
29  const Array4<Real>& vely = yvel.array(mfi);
30  const Array4<Real>& velz = zvel.array(mfi);
31 
32  // Time average at CC
33  Array4<Real> vel_t_avg_arr = vel_t_avg->array(mfi);
34 
35  ParallelFor(tbx, [=] AMREX_GPU_DEVICE (int i, int j, int k)
36  {
37  Real u_cc = myhalf * ( velx(i,j,k) + velx(i+1,j ,k ) );
38  Real v_cc = myhalf * ( vely(i,j,k) + vely(i ,j+1,k ) );
39  Real w_cc = myhalf * ( velz(i,j,k) + velz(i ,j ,k+1) );
40  Real umag_cc = std::sqrt(u_cc*u_cc + v_cc*v_cc + w_cc*w_cc);
41  vel_t_avg_arr(i,j,k,0) += u_cc * dt;
42  vel_t_avg_arr(i,j,k,1) += v_cc * dt;
43  vel_t_avg_arr(i,j,k,2) += w_cc * dt;
44  vel_t_avg_arr(i,j,k,3) += umag_cc * dt;
45  });
46  }
47 }
myhalf
constexpr amrex::Real myhalf
Definition: ERF_Constants.H:11
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
Vars::xvel
@ xvel
Definition: ERF_IndexDefines.H:159
Vars::zvel
@ zvel
Definition: ERF_IndexDefines.H:161
Vars::yvel
@ yvel
Definition: ERF_IndexDefines.H:160

Referenced by ERF::Advance(), and ERF::InitData_post().

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