ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
ERF_Advection.H
Go to the documentation of this file.
1 #ifndef ERF_ADVECTION_H_
2 #define ERF_ADVECTION_H_
3 
4 #include "AMReX.H"
5 #include "AMReX_MultiFab.H"
6 #include "AMReX_iMultiFab.H"
7 #include "AMReX_BCRec.H"
8 #include "AMReX_YAFluxRegister.H"
9 
10 #include "ERF_DataStruct.H"
11 #include "ERF_IndexDefines.H"
12 #include "ERF_SurfaceLayer.H"
13 #include "ERF_EB.H"
14 
15 /** Compute advection tendency for density and potential temperature */
16 void AdvectionSrcForRho (const amrex::Box& bx,
17  const amrex::Array4<amrex::Real>& src,
18  const amrex::Array4<const amrex::Real>& rho_u, // These are being used
19  const amrex::Array4<const amrex::Real>& rho_v, // to define the fluxes
20  const amrex::Array4<const amrex::Real>& omega,
21  const amrex::Array4< amrex::Real>& avg_xmom, // These are being defined
22  const amrex::Array4< amrex::Real>& avg_ymom, // from the rho fluxes
23  const amrex::Array4< amrex::Real>& avg_zmom,
24  const amrex::Array4<const amrex::Real>& ax_arr,
25  const amrex::Array4<const amrex::Real>& ay_arr,
26  const amrex::Array4<const amrex::Real>& az_arr,
27  const amrex::Array4<const amrex::Real>& detJ,
28  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
29  const amrex::Array4<const amrex::Real>& mf_mx,
30  const amrex::Array4<const amrex::Real>& mf_my,
31  const amrex::Array4<const amrex::Real>& mf_uy,
32  const amrex::Array4<const amrex::Real>& mf_vx,
33  const amrex::GpuArray<const amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_arr,
34  const bool fixed_rho);
35 
36 /** Compute advection tendency for all scalars other than density and potential temperature */
37 void AdvectionSrcForScalars (const amrex::Real& dt,
38  const amrex::Box& bx,
39  const int icomp, const int ncomp,
40  const amrex::Array4<const amrex::Real>& avg_xmom,
41  const amrex::Array4<const amrex::Real>& avg_ymom,
42  const amrex::Array4<const amrex::Real>& avg_zmom,
43  const amrex::Array4<const amrex::Real>& cur_cons,
44  const amrex::Array4<const amrex::Real>& cell_prim,
45  const amrex::Array4<amrex::Real>& src,
46  const bool& use_mono_adv,
47  amrex::Real* max_s_ptr,
48  amrex::Real* min_s_ptr,
49  const amrex::Array4<const amrex::Real>& vf_arr,
50  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
51  const amrex::Array4<const amrex::Real>& mf_mx,
52  const amrex::Array4<const amrex::Real>& mf_my,
53  const AdvType horiz_adv_type, const AdvType vert_adv_type,
54  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
55  const amrex::GpuArray<const amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_arr,
56  const amrex::GpuArray< amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_tmp_arr,
57  const amrex::Box& domain,
58  const amrex::BCRec* bc_ptr_h);
59 
60 /** Compute advection tendencies for all components of momentum */
61 void AdvectionSrcForMom (const amrex::MFIter& mfi,
62  const amrex::Box& bx,
63  const amrex::Box& bxx, const amrex::Box& bxy, const amrex::Box& bxz,
64  const amrex::Vector<amrex::Box>& bxx_grown,
65  const amrex::Vector<amrex::Box>& bxy_grown,
66  const amrex::Vector<amrex::Box>& bxz_grown,
67  const amrex::Array4< amrex::Real>& rho_u_rhs, const amrex::Array4< amrex::Real>& rho_v_rhs,
68  const amrex::Array4< amrex::Real>& rho_w_rhs,
69  const amrex::Array4<const amrex::Real>& rho,
70  const amrex::Array4<const amrex::Real>& u , const amrex::Array4<const amrex::Real>& v,
71  const amrex::Array4<const amrex::Real>& w ,
72  const amrex::Array4<const amrex::Real>& rho_u , const amrex::Array4<const amrex::Real>& rho_v,
73  const amrex::Array4<const amrex::Real>& Omega ,
74  const amrex::Array4<const amrex::Real>& z_nd,
75  const amrex::Array4<const amrex::Real>& ax,
76  const amrex::Array4<const amrex::Real>& ay,
77  const amrex::Array4<const amrex::Real>& az,
78  const amrex::Array4<const amrex::Real>& detJ,
79  amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
80  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
81  const amrex::Array4<const amrex::Real>& mf_mx,
82  const amrex::Array4<const amrex::Real>& mf_ux,
83  const amrex::Array4<const amrex::Real>& mf_vx,
84  const amrex::Array4<const amrex::Real>& mf_my,
85  const amrex::Array4<const amrex::Real>& mf_uy,
86  const amrex::Array4<const amrex::Real>& mf_vy,
87  const AdvType horiz_adv_type, const AdvType vert_adv_type,
88  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
89  MeshType& mesh_type, TerrainType& terrain_type,
90  const eb_& ebfact,
91  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_u_arr,
92  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_v_arr,
93  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_w_arr,
94  const amrex::Vector<amrex::iMultiFab>& physbnd_mask,
95  const bool already_on_centroids,
96  const int lo_z_face, const int hi_z_face,
97  const amrex::Box& domain,
98  const amrex::BCRec* bc_ptr_h);
99 
100 /** Compute advection tendencies for all components of momentum when dz is constant */
101 void AdvectionSrcForMom_ConstantDz (
102  const amrex::Box& bxx, const amrex::Box& bxy, const amrex::Box& bxz,
103  const amrex::Array4< amrex::Real>& rho_u_rhs, const amrex::Array4< amrex::Real>& rho_v_rhs,
104  const amrex::Array4< amrex::Real>& rho_w_rhs,
105  const amrex::Array4<const amrex::Real>& u , const amrex::Array4<const amrex::Real>& v,
106  const amrex::Array4<const amrex::Real>& w ,
107  const amrex::Array4<const amrex::Real>& rho_u , const amrex::Array4<const amrex::Real>& rho_v,
108  const amrex::Array4<const amrex::Real>& Omega ,
109  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
110  const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
111  const amrex::Array4<const amrex::Real>& mf_mx,
112  const amrex::Array4<const amrex::Real>& mf_ux,
113  const amrex::Array4<const amrex::Real>& mf_vx,
114  const amrex::Array4<const amrex::Real>& mf_my,
115  const amrex::Array4<const amrex::Real>& mf_uy,
116  const amrex::Array4<const amrex::Real>& mf_vy,
117  const AdvType horiz_adv_type, const AdvType vert_adv_type,
118  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
119  TerrainType& terrain_type,
120  const int lo_z_face, const int hi_z_face);
121 
122 /** Compute advection tendencies for all components of momentum when dz is stretched */
123 void AdvectionSrcForMom_StretchedDz (
124  const amrex::Box& bxx, const amrex::Box& bxy, const amrex::Box& bxz,
125  const amrex::Array4< amrex::Real>& rho_u_rhs, const amrex::Array4< amrex::Real>& rho_v_rhs,
126  const amrex::Array4< amrex::Real>& rho_w_rhs,
127  const amrex::Array4<const amrex::Real>& u , const amrex::Array4<const amrex::Real>& v,
128  const amrex::Array4<const amrex::Real>& w ,
129  const amrex::Array4<const amrex::Real>& rho_u , const amrex::Array4<const amrex::Real>& rho_v,
130  const amrex::Array4<const amrex::Real>& Omega ,
131  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
132  const amrex::Gpu::DeviceVector<amrex::Real>& stretched_dz_d,
133  const amrex::Array4<const amrex::Real>& mf_mx,
134  const amrex::Array4<const amrex::Real>& mf_ux,
135  const amrex::Array4<const amrex::Real>& mf_vx,
136  const amrex::Array4<const amrex::Real>& mf_my,
137  const amrex::Array4<const amrex::Real>& mf_uy,
138  const amrex::Array4<const amrex::Real>& mf_vy,
139  const AdvType horiz_adv_type, const AdvType vert_adv_type,
140  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
141  const int lo_z_face, const int hi_z_face);
142 
143 /** Compute advection tendencies for all components of momentum for TERRAIN-FITTED COORDS only */
144 void AdvectionSrcForMom_TF (
145  const amrex::Box& bxx, const amrex::Box& bxy, const amrex::Box& bxz,
146  const amrex::Array4< amrex::Real>& rho_u_rhs, const amrex::Array4< amrex::Real>& rho_v_rhs,
147  const amrex::Array4< amrex::Real>& rho_w_rhs,
148  const amrex::Array4<const amrex::Real>& u , const amrex::Array4<const amrex::Real>& v,
149  const amrex::Array4<const amrex::Real>& w ,
150  const amrex::Array4<const amrex::Real>& rho_u , const amrex::Array4<const amrex::Real>& rho_v,
151  const amrex::Array4<const amrex::Real>& Omega ,
152  const amrex::Array4<const amrex::Real>& z_nd,
153  const amrex::Array4<const amrex::Real>& ax,
154  const amrex::Array4<const amrex::Real>& ay,
155  const amrex::Array4<const amrex::Real>& az,
156  const amrex::Array4<const amrex::Real>& detJ,
157  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
158  const amrex::Array4<const amrex::Real>& mf_mx,
159  const amrex::Array4<const amrex::Real>& mf_ux,
160  const amrex::Array4<const amrex::Real>& mf_vx,
161  const amrex::Array4<const amrex::Real>& mf_my,
162  const amrex::Array4<const amrex::Real>& mf_uy,
163  const amrex::Array4<const amrex::Real>& mf_vy,
164  const AdvType horiz_adv_type, const AdvType vert_adv_type,
165  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
166  const int lo_z_face, const int hi_z_face);
167 
168 /** Compute advection tendencies for all components of momentum when using EB for terrain */
169 void AdvectionSrcForMom_EB (const amrex::MFIter& mfi,
170  const amrex::Box& bxx, const amrex::Box& bxy, const amrex::Box& bxz,
171  const amrex::Vector<amrex::Box>& bxx_grown,
172  const amrex::Vector<amrex::Box>& bxy_grown,
173  const amrex::Vector<amrex::Box>& bxz_grown,
174  const amrex::Array4< amrex::Real>& rho_u_rhs, const amrex::Array4< amrex::Real>& rho_v_rhs,
175  const amrex::Array4< amrex::Real>& rho_w_rhs,
176  const amrex::Array4<const amrex::Real>& u , const amrex::Array4<const amrex::Real>& v,
177  const amrex::Array4<const amrex::Real>& w ,
178  const amrex::Array4<const amrex::Real>& rho_u , const amrex::Array4<const amrex::Real>& rho_v,
179  const amrex::Array4<const amrex::Real>& Omega ,
180  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
181  const amrex::Array4<const amrex::Real>& mf_mx,
182  const amrex::Array4<const amrex::Real>& mf_ux,
183  const amrex::Array4<const amrex::Real>& mf_vx,
184  const amrex::Array4<const amrex::Real>& mf_my,
185  const amrex::Array4<const amrex::Real>& mf_uy,
186  const amrex::Array4<const amrex::Real>& mf_vy,
187  const AdvType horiz_adv_type, const AdvType vert_adv_type,
188  const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac,
189  const eb_& ebfact,
190  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_u_arr,
191  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_v_arr,
192  amrex::GpuArray<amrex::Array4<amrex::Real>, AMREX_SPACEDIM>& flx_w_arr,
193  const amrex::Vector<amrex::iMultiFab>& physbnd_mask,
194  const bool already_on_centroids,
195  const int lo_z_face, const int hi_z_face,
196  const amrex::Box& domain);
197 
198 /** Compute advection tendencies for all normal components of momentum */
199 void
200 AdvectionSrcForOpenBC_Normal (const amrex::Box& bx,
201  const int& dir,
202  const amrex::Array4< amrex::Real>& rhs_arr,
203  const amrex::Array4<const amrex::Real>& vel_norm_arr,
204  const amrex::Array4<const amrex::Real>& cell_data_arr,
205  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& dxInv,
206  const bool do_lo=false);
207 
208 /** Compute advection tendencies for u momentum tangent to open bc*/
209 void
210 AdvectionSrcForOpenBC_Tangent_Xmom (const amrex::Box& bxx,
211  const int& dir,
212  const amrex::Array4< amrex::Real>& rho_u_rhs,
213  const amrex::Array4<const amrex::Real>& u,
214  const amrex::Array4<const amrex::Real>& rho_u,
215  const amrex::Array4<const amrex::Real>& rho_v,
216  const amrex::Array4<const amrex::Real>& Omega,
217  const amrex::Array4<const amrex::Real>& ax,
218  const amrex::Array4<const amrex::Real>& az,
219  const amrex::Array4<const amrex::Real>& detJ,
220  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
221  const bool do_lo=false);
222 
223 /** Compute advection tendencies for v momentum tangent to open bc*/
224 void
225 AdvectionSrcForOpenBC_Tangent_Ymom (const amrex::Box& bxy,
226  const int& dir,
227  const amrex::Array4< amrex::Real>& rho_v_rhs,
228  const amrex::Array4<const amrex::Real>& v,
229  const amrex::Array4<const amrex::Real>& rho_u,
230  const amrex::Array4<const amrex::Real>& rho_v,
231  const amrex::Array4<const amrex::Real>& Omega,
232  const amrex::Array4<const amrex::Real>& ay,
233  const amrex::Array4<const amrex::Real>& az,
234  const amrex::Array4<const amrex::Real>& detJ,
235  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
236  const bool do_lo=false);
237 
238 /** Compute advection tendencies for w momentum tangent to open bc*/
239 void
240 AdvectionSrcForOpenBC_Tangent_Zmom (const amrex::Box& bxz,
241  const int& dir,
242  const amrex::Array4< amrex::Real>& rho_w_rhs,
243  const amrex::Array4<const amrex::Real>& w,
244  const amrex::Array4<const amrex::Real>& rho_u,
245  const amrex::Array4<const amrex::Real>& rho_v,
246  const amrex::Array4<const amrex::Real>& Omega,
247  const amrex::Array4<const amrex::Real>& ax,
248  const amrex::Array4<const amrex::Real>& ay,
249  const amrex::Array4<const amrex::Real>& az,
250  const amrex::Array4<const amrex::Real>& detJ,
251  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
252  const int domhi_z,
253  const bool do_lo=false);
254 
255 /** Compute advection tendencies for cons tangential to BC (2nd order)*/
256 void
257 AdvectionSrcForOpenBC_Tangent_Cons (const amrex::Box& bx,
258  const int& dir,
259  const int& icomp,
260  const int& ncomp,
261  const amrex::Array4< amrex::Real>& cell_rhs,
262  const amrex::Array4<const amrex::Real>& cell_prim,
263  const amrex::Array4<const amrex::Real>& avg_xmom,
264  const amrex::Array4<const amrex::Real>& avg_ymom,
265  const amrex::Array4<const amrex::Real>& avg_zmom,
266  const amrex::Array4<const amrex::Real>& detJ,
267  const amrex::GpuArray<amrex::Real, AMREX_SPACEDIM>& cellSizeInv,
268  const bool do_lo=false);
269 
270 /** Compute advection tendencies in normal dir for vars tangent to open bc */
271 AMREX_GPU_HOST_DEVICE
272 AMREX_FORCE_INLINE
273 amrex::Real
274 AdvectionSrcForOpenBC_Tangent (const int& i,
275  const int& j,
276  const int& k,
277  const int& nprim,
278  const int& dir,
279  const amrex::Array4<const amrex::Real>& prim_tang_arr,
280  const amrex::Array4<const amrex::Real>& mom_norm_arr,
281  const amrex::Real& dxInv,
282  const bool do_lo=false);
283 
284 AMREX_GPU_HOST_DEVICE
285 AMREX_FORCE_INLINE
286 AdvType
287 EfficientAdvType (int nrk, AdvType adv_type)
288 {
289  AdvType eff_adv_type;
290 
291  if (nrk == 0) {
292 
293  eff_adv_type = AdvType::Centered_2nd;
294 
295  } else if (nrk == 1) {
296 
297  if ( (adv_type == AdvType::Centered_2nd) ||
298  (adv_type == AdvType::Upwind_3rd) ||
299  (adv_type == AdvType::Centered_4th) ||
300  (adv_type == AdvType::Weno_3) ||
301  (adv_type == AdvType::Weno_3Z) ||
302  (adv_type == AdvType::Weno_3MZQ) )
303  {
304  eff_adv_type = AdvType::Centered_2nd;
305 
306  } else if ( (adv_type == AdvType::Upwind_5th) ||
307  (adv_type == AdvType::Weno_5) ||
308  (adv_type == AdvType::Weno_5Z) )
309  {
310  eff_adv_type = AdvType::Upwind_3rd;
311 
312  } else { // (adv_type == AdvType::Centered_6th)
313 
314  eff_adv_type = AdvType::Centered_4th;
315  }
316 
317  } else {
318 
319  eff_adv_type = adv_type;
320  }
321 
322  return(eff_adv_type);
323 }
324 #endif
AdvectionSrcForRho
void AdvectionSrcForRho(const amrex::Box &bx, const amrex::Array4< amrex::Real > &src, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &omega, const amrex::Array4< amrex::Real > &avg_xmom, const amrex::Array4< amrex::Real > &avg_ymom, const amrex::Array4< amrex::Real > &avg_zmom, const amrex::Array4< const amrex::Real > &ax_arr, const amrex::Array4< const amrex::Real > &ay_arr, const amrex::Array4< const amrex::Real > &az_arr, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::GpuArray< const amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_arr, const bool fixed_rho)
AdvectionSrcForMom_EB
void AdvectionSrcForMom_EB(const amrex::MFIter &mfi, const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Vector< amrex::Box > &bxx_grown, const amrex::Vector< amrex::Box > &bxy_grown, const amrex::Vector< amrex::Box > &bxz_grown, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, const eb_ &ebfact, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_u_arr, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_v_arr, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_w_arr, const amrex::Vector< amrex::iMultiFab > &physbnd_mask, const bool already_on_centroids, const int lo_z_face, const int hi_z_face, const amrex::Box &domain)
AdvectionSrcForMom_StretchedDz
void AdvectionSrcForMom_StretchedDz(const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, const int lo_z_face, const int hi_z_face)
AdvectionSrcForOpenBC_Tangent_Ymom
void AdvectionSrcForOpenBC_Tangent_Ymom(const amrex::Box &bxy, const int &dir, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::Array4< const amrex::Real > &ay, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const bool do_lo=false)
AdvectionSrcForOpenBC_Tangent_Xmom
void AdvectionSrcForOpenBC_Tangent_Xmom(const amrex::Box &bxx, const int &dir, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::Array4< const amrex::Real > &ax, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const bool do_lo=false)
AdvectionSrcForMom_TF
void AdvectionSrcForMom_TF(const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &ax, const amrex::Array4< const amrex::Real > &ay, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, const int lo_z_face, const int hi_z_face)
AdvectionSrcForOpenBC_Normal
void AdvectionSrcForOpenBC_Normal(const amrex::Box &bx, const int &dir, const amrex::Array4< amrex::Real > &rhs_arr, const amrex::Array4< const amrex::Real > &vel_norm_arr, const amrex::Array4< const amrex::Real > &cell_data_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &dxInv, const bool do_lo=false)
AdvectionSrcForOpenBC_Tangent
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real AdvectionSrcForOpenBC_Tangent(const int &i, const int &j, const int &k, const int &nprim, const int &dir, const amrex::Array4< const amrex::Real > &prim_tang_arr, const amrex::Array4< const amrex::Real > &mom_norm_arr, const amrex::Real &dxInv, const bool do_lo=false)
AdvectionSrcForMom
void AdvectionSrcForMom(const amrex::MFIter &mfi, const amrex::Box &bx, const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Vector< amrex::Box > &bxx_grown, const amrex::Vector< amrex::Box > &bxy_grown, const amrex::Vector< amrex::Box > &bxz_grown, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &rho, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::Array4< const amrex::Real > &z_nd, const amrex::Array4< const amrex::Real > &ax, const amrex::Array4< const amrex::Real > &ay, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, MeshType &mesh_type, TerrainType &terrain_type, const eb_ &ebfact, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_u_arr, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_v_arr, amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_w_arr, const amrex::Vector< amrex::iMultiFab > &physbnd_mask, const bool already_on_centroids, const int lo_z_face, const int hi_z_face, const amrex::Box &domain, const amrex::BCRec *bc_ptr_h)
AdvectionSrcForOpenBC_Tangent_Cons
void AdvectionSrcForOpenBC_Tangent_Cons(const amrex::Box &bx, const int &dir, const int &icomp, const int &ncomp, const amrex::Array4< amrex::Real > &cell_rhs, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< const amrex::Real > &avg_xmom, const amrex::Array4< const amrex::Real > &avg_ymom, const amrex::Array4< const amrex::Real > &avg_zmom, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const bool do_lo=false)
AdvectionSrcForOpenBC_Tangent_Zmom
void AdvectionSrcForOpenBC_Tangent_Zmom(const amrex::Box &bxz, const int &dir, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::Array4< const amrex::Real > &ax, const amrex::Array4< const amrex::Real > &ay, const amrex::Array4< const amrex::Real > &az, const amrex::Array4< const amrex::Real > &detJ, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const int domhi_z, const bool do_lo=false)
AdvectionSrcForMom_ConstantDz
void AdvectionSrcForMom_ConstantDz(const amrex::Box &bxx, const amrex::Box &bxy, const amrex::Box &bxz, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< amrex::Real > &rho_w_rhs, const amrex::Array4< const amrex::Real > &u, const amrex::Array4< const amrex::Real > &v, const amrex::Array4< const amrex::Real > &w, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &Omega, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Gpu::DeviceVector< amrex::Real > &stretched_dz_d, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_ux, const amrex::Array4< const amrex::Real > &mf_vx, const amrex::Array4< const amrex::Real > &mf_my, const amrex::Array4< const amrex::Real > &mf_uy, const amrex::Array4< const amrex::Real > &mf_vy, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, TerrainType &terrain_type, const int lo_z_face, const int hi_z_face)
AdvectionSrcForScalars
void AdvectionSrcForScalars(const amrex::Real &dt, const amrex::Box &bx, const int icomp, const int ncomp, const amrex::Array4< const amrex::Real > &avg_xmom, const amrex::Array4< const amrex::Real > &avg_ymom, const amrex::Array4< const amrex::Real > &avg_zmom, const amrex::Array4< const amrex::Real > &cur_cons, const amrex::Array4< const amrex::Real > &cell_prim, const amrex::Array4< amrex::Real > &src, const bool &use_mono_adv, amrex::Real *max_s_ptr, amrex::Real *min_s_ptr, const amrex::Array4< const amrex::Real > &vf_arr, const amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > &cellSizeInv, const amrex::Array4< const amrex::Real > &mf_mx, const amrex::Array4< const amrex::Real > &mf_my, const AdvType horiz_adv_type, const AdvType vert_adv_type, const amrex::Real horiz_upw_frac, const amrex::Real vert_upw_frac, const amrex::GpuArray< const amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_arr, const amrex::GpuArray< amrex::Array4< amrex::Real >, AMREX_SPACEDIM > &flx_tmp_arr, const amrex::Box &domain, const amrex::BCRec *bc_ptr_h)
EfficientAdvType
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE AdvType EfficientAdvType(int nrk, AdvType adv_type)
Definition: ERF_Advection.H:287
ERF_DataStruct.H
ERF_EB.H
ERF_IndexDefines.H
AdvType
AdvType
Definition: ERF_IndexDefines.H:221
AdvType::Upwind_3rd
@ Upwind_3rd
AdvType::Centered_4th
@ Centered_4th
AdvType::Upwind_5th
@ Upwind_5th
AdvType::Weno_3Z
@ Weno_3Z
AdvType::Weno_3
@ Weno_3
AdvType::Centered_2nd
@ Centered_2nd
AdvType::Weno_5Z
@ Weno_5Z
AdvType::Weno_5
@ Weno_5
AdvType::Weno_3MZQ
@ Weno_3MZQ
ERF_SurfaceLayer.H
eb_
Definition: ERF_EB.H:13
MicVar_Kess::rho
@ rho
Definition: ERF_Kessler.H:22
MicVar_Morr::omega
@ omega
Definition: ERF_Morrison.H:53