ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
ERF_SrcHeaders.H
Go to the documentation of this file.
1 #ifndef ERF_INTEGRATION_H_
2 #define ERF_INTEGRATION_H_
3 
4 #include <AMReX_MultiFab.H>
5 #include <AMReX_EBMultiFabUtil.H>
6 
7 #include "ERF_DataStruct.H"
8 #include "ERF_InputSoundingData.H"
9 #include "ERF_TurbPertStruct.H"
10 #include "ERF_EB.H"
11 
12 /**
13  * Function for computing the buoyancy term to be used in the evolution
14  * equation for the z-component of momentum in the slow integrator. There
15  * are three options for how buoyancy is computed (two are the same in the absence of moisture).
16  */
17 void make_buoyancy (const amrex::Vector<amrex::MultiFab>& S_data,
18  const amrex::MultiFab & S_prim,
19  const amrex::MultiFab & qt,
20  amrex::MultiFab& buoyancy,
21  const amrex::Geometry geom,
22  const SolverChoice& solverChoice,
23  const amrex::MultiFab& base_state,
24  const int n_qstate,
25  const eb_& ebfact,
26  const int anelastic);
27 
28 void make_gradp_pert (int level,
29  const SolverChoice& solverChoice,
30  const amrex::Geometry& geom,
31  amrex::Vector<amrex::MultiFab>& S_data,
32  const amrex::MultiFab& p0,
33  const amrex::MultiFab& z_phys_nd,
34  const amrex::MultiFab& z_phys_cc,
35  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
36  const eb_& ebfact,
37  amrex::Vector<amrex::MultiFab>& gradp);
38 
39 void compute_gradp (const amrex::MultiFab& p,
40  const amrex::Geometry& geom,
41  const amrex::MultiFab& z_phys_nd,
42  const amrex::MultiFab& z_phys_cc,
43  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
44  const eb_& ebfact,
45  amrex::Vector<amrex::MultiFab>& gradp,
46  const SolverChoice& solverChoice);
47 
48 void compute_gradp_interpz (const amrex::MultiFab& p,
49  const amrex::Geometry& geom,
50  const amrex::MultiFab& z_phys_nd,
51  const amrex::MultiFab& z_phys_cc,
52  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
53  amrex::Vector<amrex::MultiFab>& gradp,
54  const SolverChoice& solverChoice);
55 
56 void make_sources (int level, int nrk,
57  amrex::Real dt, amrex::Real time,
58  const amrex::Vector<amrex::MultiFab>& S_data,
59  const amrex::MultiFab& S_prim,
60  amrex::MultiFab& cc_source,
61  const amrex::MultiFab& base_state,
62  const amrex::MultiFab* z_phys_cc,
63  const amrex::MultiFab& xvel,
64  const amrex::MultiFab& yvel,
65  const amrex::MultiFab* qheating_rates,
66  amrex::MultiFab* terrain_blank,
67  const amrex::Geometry geom,
68  const SolverChoice& solverChoice,
69  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
70  const amrex::Real* dptr_rhotheta_src,
71  const amrex::Real* dptr_rhoqt_src,
72  const amrex::Real* dptr_wbar_sub,
73  const amrex::Vector<amrex::Real*> d_rayleigh_ptrs_at_lev,
74  InputSoundingData& input_sounding_data,
75  TurbulentPerturbation& turbPert,
76  bool is_slow_step);
77 
78 void make_mom_sources (amrex::Real time,
79  amrex::Real dt,
80  const amrex::Vector<amrex::MultiFab>& S_data,
81  amrex::MultiFab& z_phys_nd,
82  amrex::MultiFab& z_phys_cc,
83  amrex::Vector<amrex::Real>& stretched_dz_h,
84  const amrex::MultiFab& xvel,
85  const amrex::MultiFab& yvel,
86  const amrex::MultiFab& zvel,
87  amrex::MultiFab& xmom_source,
88  amrex::MultiFab& ymom_source,
89  amrex::MultiFab& zmom_source,
90  const amrex::MultiFab& base_state,
91  amrex::MultiFab* forest_drag,
92  amrex::MultiFab* terrain_blank,
93  amrex::MultiFab* cosPhi_mf,
94  amrex::MultiFab* sinPhi_mf,
95  const amrex::Geometry geom,
96  const SolverChoice& solverChoice,
97  amrex::Vector<std::unique_ptr<amrex::MultiFab>>& mapfac,
98  const amrex::Real* dptr_rhotheta_src,
99  const amrex::Real* dptr_rhoqt_src,
100  const amrex::Real* dptr_wbar_sub,
101  const amrex::Vector<amrex::Real*> d_rayleigh_ptrs_at_lev,
102  const amrex::Vector<amrex::Real*> d_sponge_ptrs_at_lev,
103  const amrex::Vector<amrex::MultiFab>* forecast_state_at_lev,
104  InputSoundingData& input_sounding_data,
105  bool is_slow_step);
106 
107 void add_thin_body_sources (amrex::MultiFab& xmom_source,
108  amrex::MultiFab& ymom_source,
109  amrex::MultiFab& zmom_source,
110  std::unique_ptr<amrex::iMultiFab>& xflux_imask_lev,
111  std::unique_ptr<amrex::iMultiFab>& yflux_imask_lev,
112  std::unique_ptr<amrex::iMultiFab>& zflux_imask_lev,
113  std::unique_ptr<amrex::MultiFab>& thin_xforce_lev,
114  std::unique_ptr<amrex::MultiFab>& thin_yforce_lev,
115  std::unique_ptr<amrex::MultiFab>& thin_zforce_lev);
116 
117 #if defined(ERF_USE_NETCDF)
118 void
119 moist_set_rhs (const amrex::Geometry& geom,
120  const amrex::Box& tbx,
121  const amrex::Array4<amrex::Real const>& old_cons,
122  const amrex::Array4<amrex::Real const>& new_cons,
123  const amrex::Array4<amrex::Real >& cell_rhs,
124  const amrex::Real& bdy_time_interval,
125  const amrex::Real& new_stage_time,
126  const amrex::Real& dt,
127  const amrex::Real& stop_time_elapsed,
128  int width, int set_width,
129  const amrex::Box& domain,
130  amrex::Vector<amrex::Vector<amrex::FArrayBox>>& bdy_data_xlo,
131  amrex::Vector<amrex::Vector<amrex::FArrayBox>>& bdy_data_xhi,
132  amrex::Vector<amrex::Vector<amrex::FArrayBox>>& bdy_data_ylo,
133  amrex::Vector<amrex::Vector<amrex::FArrayBox>>& bdy_data_yhi);
134 #endif
135 
136 void ApplySpongeZoneBCsForCC (const SpongeChoice& spongeChoice,
137  const amrex::Geometry geom,
138  const amrex::Box& bx,
139  const amrex::Array4<amrex::Real>& cell_rhs,
140  const amrex::Array4<const amrex::Real>& cell_data,
141  const amrex::Array4<const amrex::Real>& r0,
142  const amrex::Array4<const amrex::Real>& z_phys_cc);
143 
144 void ApplySpongeZoneBCsForMom (const SpongeChoice& spongeChoice,
145  const amrex::Geometry geom,
146  const amrex::Box& tbx,
147  const amrex::Box& tby,
148  const amrex::Box& tbz,
149  const amrex::Array4<amrex::Real>& rho_u_rhs,
150  const amrex::Array4<amrex::Real>& rho_v_rhs,
151  const amrex::Array4<amrex::Real>& rho_w_rhs,
152  const amrex::Array4<const amrex::Real>& rho_u,
153  const amrex::Array4<const amrex::Real>& rho_v,
154  const amrex::Array4<const amrex::Real>& rho_w,
155  const amrex::Array4<const amrex::Real>& r0,
156  const amrex::Array4<const amrex::Real>& z_phys_nd,
157  const amrex::Array4<const amrex::Real>& z_phys_cc);
158 
159 void ApplySpongeZoneBCsForMom_ReadFromFile (const SpongeChoice& spongeChoice,
160  const amrex::Geometry geom,
161  const amrex::Box& tbx,
162  const amrex::Box& tby,
163  const amrex::Array4<const amrex::Real>& cell_data,
164  const amrex::Array4<const amrex::Real>& z_phys_cc,
165  const amrex::Array4<amrex::Real>& rho_u_rhs,
166  const amrex::Array4<amrex::Real>& rho_v_rhs,
167  const amrex::Array4<const amrex::Real>& rho_u,
168  const amrex::Array4<const amrex::Real>& rho_v,
169  const amrex::Vector<amrex::Real*> d_sponge_ptrs_at_lev);
170 
171 void ApplyBndryForcing_Forecast (const SolverChoice& solverChoice,
172  const amrex::Geometry geom,
173  const amrex::Box& tbx,
174  const amrex::Box& tby,
175  const amrex::Box& tbz,
176  const amrex::Array4<const amrex::Real>& z_phys_nd,
177  const amrex::Array4<amrex::Real>& rho_u_rhs,
178  const amrex::Array4<amrex::Real>& rho_v_rhs,
179  const amrex::Array4<amrex::Real>& rho_w_rhs,
180  const amrex::Array4<const amrex::Real>& rho_u,
181  const amrex::Array4<const amrex::Real>& rho_v,
182  const amrex::Array4<const amrex::Real>& rho_w,
183  const amrex::Array4<const amrex::Real>& rho_u_initial_state,
184  const amrex::Array4<const amrex::Real>& rho_v_initial_state,
185  const amrex::Array4<const amrex::Real>& rho_w_initial_state,
186  const amrex::Array4<const amrex::Real>& cons_initial_state);
187 
188 #endif
ERF_DataStruct.H
ERF_EB.H
ERF_InputSoundingData.H
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
make_gradp_pert
void make_gradp_pert(int level, const SolverChoice &solverChoice, const amrex::Geometry &geom, amrex::Vector< amrex::MultiFab > &S_data, const amrex::MultiFab &p0, const amrex::MultiFab &z_phys_nd, const amrex::MultiFab &z_phys_cc, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const eb_ &ebfact, amrex::Vector< amrex::MultiFab > &gradp)
compute_gradp
void compute_gradp(const amrex::MultiFab &p, const amrex::Geometry &geom, const amrex::MultiFab &z_phys_nd, const amrex::MultiFab &z_phys_cc, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const eb_ &ebfact, amrex::Vector< amrex::MultiFab > &gradp, const SolverChoice &solverChoice)
ApplyBndryForcing_Forecast
void ApplyBndryForcing_Forecast(const SolverChoice &solverChoice, const amrex::Geometry geom, const amrex::Box &tbx, const amrex::Box &tby, const amrex::Box &tbz, const amrex::Array4< const amrex::Real > &z_phys_nd, 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_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &rho_w, const amrex::Array4< const amrex::Real > &rho_u_initial_state, const amrex::Array4< const amrex::Real > &rho_v_initial_state, const amrex::Array4< const amrex::Real > &rho_w_initial_state, const amrex::Array4< const amrex::Real > &cons_initial_state)
add_thin_body_sources
void add_thin_body_sources(amrex::MultiFab &xmom_source, amrex::MultiFab &ymom_source, amrex::MultiFab &zmom_source, std::unique_ptr< amrex::iMultiFab > &xflux_imask_lev, std::unique_ptr< amrex::iMultiFab > &yflux_imask_lev, std::unique_ptr< amrex::iMultiFab > &zflux_imask_lev, std::unique_ptr< amrex::MultiFab > &thin_xforce_lev, std::unique_ptr< amrex::MultiFab > &thin_yforce_lev, std::unique_ptr< amrex::MultiFab > &thin_zforce_lev)
make_mom_sources
void make_mom_sources(amrex::Real time, amrex::Real dt, const amrex::Vector< amrex::MultiFab > &S_data, amrex::MultiFab &z_phys_nd, amrex::MultiFab &z_phys_cc, amrex::Vector< amrex::Real > &stretched_dz_h, const amrex::MultiFab &xvel, const amrex::MultiFab &yvel, const amrex::MultiFab &zvel, amrex::MultiFab &xmom_source, amrex::MultiFab &ymom_source, amrex::MultiFab &zmom_source, const amrex::MultiFab &base_state, amrex::MultiFab *forest_drag, amrex::MultiFab *terrain_blank, amrex::MultiFab *cosPhi_mf, amrex::MultiFab *sinPhi_mf, const amrex::Geometry geom, const SolverChoice &solverChoice, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const amrex::Real *dptr_rhotheta_src, const amrex::Real *dptr_rhoqt_src, const amrex::Real *dptr_wbar_sub, const amrex::Vector< amrex::Real * > d_rayleigh_ptrs_at_lev, const amrex::Vector< amrex::Real * > d_sponge_ptrs_at_lev, const amrex::Vector< amrex::MultiFab > *forecast_state_at_lev, InputSoundingData &input_sounding_data, bool is_slow_step)
compute_gradp_interpz
void compute_gradp_interpz(const amrex::MultiFab &p, const amrex::Geometry &geom, const amrex::MultiFab &z_phys_nd, const amrex::MultiFab &z_phys_cc, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, amrex::Vector< amrex::MultiFab > &gradp, const SolverChoice &solverChoice)
make_sources
void make_sources(int level, int nrk, amrex::Real dt, amrex::Real time, const amrex::Vector< amrex::MultiFab > &S_data, const amrex::MultiFab &S_prim, amrex::MultiFab &cc_source, const amrex::MultiFab &base_state, const amrex::MultiFab *z_phys_cc, const amrex::MultiFab &xvel, const amrex::MultiFab &yvel, const amrex::MultiFab *qheating_rates, amrex::MultiFab *terrain_blank, const amrex::Geometry geom, const SolverChoice &solverChoice, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &mapfac, const amrex::Real *dptr_rhotheta_src, const amrex::Real *dptr_rhoqt_src, const amrex::Real *dptr_wbar_sub, const amrex::Vector< amrex::Real * > d_rayleigh_ptrs_at_lev, InputSoundingData &input_sounding_data, TurbulentPerturbation &turbPert, bool is_slow_step)
make_buoyancy
void make_buoyancy(const amrex::Vector< amrex::MultiFab > &S_data, const amrex::MultiFab &S_prim, const amrex::MultiFab &qt, amrex::MultiFab &buoyancy, const amrex::Geometry geom, const SolverChoice &solverChoice, const amrex::MultiFab &base_state, const int n_qstate, const eb_ &ebfact, const int anelastic)
ApplySpongeZoneBCsForMom_ReadFromFile
void ApplySpongeZoneBCsForMom_ReadFromFile(const SpongeChoice &spongeChoice, const amrex::Geometry geom, const amrex::Box &tbx, const amrex::Box &tby, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &z_phys_cc, const amrex::Array4< amrex::Real > &rho_u_rhs, const amrex::Array4< amrex::Real > &rho_v_rhs, const amrex::Array4< const amrex::Real > &rho_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Vector< amrex::Real * > d_sponge_ptrs_at_lev)
ApplySpongeZoneBCsForCC
void ApplySpongeZoneBCsForCC(const SpongeChoice &spongeChoice, const amrex::Geometry geom, const amrex::Box &bx, const amrex::Array4< amrex::Real > &cell_rhs, const amrex::Array4< const amrex::Real > &cell_data, const amrex::Array4< const amrex::Real > &r0, const amrex::Array4< const amrex::Real > &z_phys_cc)
ApplySpongeZoneBCsForMom
void ApplySpongeZoneBCsForMom(const SpongeChoice &spongeChoice, const amrex::Geometry geom, const amrex::Box &tbx, const amrex::Box &tby, const amrex::Box &tbz, 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_u, const amrex::Array4< const amrex::Real > &rho_v, const amrex::Array4< const amrex::Real > &rho_w, const amrex::Array4< const amrex::Real > &r0, const amrex::Array4< const amrex::Real > &z_phys_nd, const amrex::Array4< const amrex::Real > &z_phys_cc)
ERF_TurbPertStruct.H
eb_
Definition: ERF_EB.H:13
MicVar_Kess::qt
@ qt
Definition: ERF_Kessler.H:27
Vars::xvel
@ xvel
Definition: ERF_IndexDefines.H:141
Vars::zvel
@ zvel
Definition: ERF_IndexDefines.H:143
Vars::yvel
@ yvel
Definition: ERF_IndexDefines.H:142
module_model_constants::p0
real(c_double), parameter p0
Definition: ERF_module_model_constants.F90:40
InputSoundingData
Definition: ERF_InputSoundingData.H:22
SolverChoice
Definition: ERF_DataStruct.H:123
SpongeChoice
Definition: ERF_SpongeStruct.H:15
TurbulentPerturbation
Definition: ERF_TurbPertStruct.H:22