5 #ifndef ERF_EULERIANMICROPHYSICS_H
6 #define ERF_EULERIANMICROPHYSICS_H
28 const MoistureType& a_model_type )
32 if (a_model_type == MoistureType::SAM ||
33 a_model_type == MoistureType::SAM_NoIce ||
34 a_model_type == MoistureType::SAM_NoPrecip_NoIce) {
36 }
else if (a_model_type == MoistureType::Kessler ||
37 a_model_type == MoistureType::Kessler_NoRain) {
39 }
else if (a_model_type == MoistureType::Morrison ||
40 a_model_type == MoistureType::Morrison_NoIce) {
42 }
else if (a_model_type == MoistureType::SatAdj) {
44 }
else if (a_model_type == MoistureType::None) {
45 SetModel<NullMoist>();
47 amrex::Abort(
"EulerianMicrophysics: Dont know this moisture_type!") ;
60 const amrex::MultiFab& cons_in,
61 const amrex::BoxArray& grids,
62 const amrex::Geometry& geom,
63 const amrex::Real& dt_advance,
64 std::unique_ptr<amrex::MultiFab>& z_phys_nd,
65 std::unique_ptr<amrex::MultiFab>& detJ_cc )
override
73 const amrex::Real& dt_advance,
77 amrex::Vector<amrex::Vector<amrex::MultiFab>>&,
78 const amrex::Vector<std::unique_ptr<amrex::MultiFab>>&,
79 const amrex::GpuArray<ERF_BC, AMREX_SPACEDIM*2>& )
override
86 amrex::MultiFab& cons_in )
override
93 amrex::MultiFab& cons_in )
override
100 const int& varIdx )
override
127 std::vector<int>& a_idx,
128 std::vector<std::string>& a_names )
const override
130 m_moist_model[a_lev]->Qmoist_Restart_Vars(a_sc, a_idx, a_names);
141 amrex::MultiFab& a_mf,
142 const int a_lev)
const override
150 template<
class NewMoistModel>
Contains the base class for microphysics.
Eulerian microphysics interface.
Definition: ERF_EulerianMicrophysics.H:16
void Get_Qmoist_Restart_Vars(const int a_lev, const SolverChoice &a_sc, std::vector< int > &a_idx, std::vector< std::string > &a_names) const override
get the indices and names of moisture model variables for restart at a given level
Definition: ERF_EulerianMicrophysics.H:125
void Update_Micro_Vars_Lev(const int &lev, amrex::MultiFab &cons_in) override
update microphysics variables from ERF state variables
Definition: ERF_EulerianMicrophysics.H:85
EulerianMicrophysics(const int &nlev, const MoistureType &a_model_type)
Constructor: create the moisture model.
Definition: ERF_EulerianMicrophysics.H:27
~EulerianMicrophysics()=default
default destructor
amrex::Vector< std::unique_ptr< NullMoist > > m_moist_model
Definition: ERF_EulerianMicrophysics.H:159
amrex::MultiFab * Get_Qmoist_Ptr(const int &lev, const int &varIdx) override
get pointer to a moisture variable
Definition: ERF_EulerianMicrophysics.H:99
void Define(const int &lev, SolverChoice &sc) override
Define the moisture model.
Definition: ERF_EulerianMicrophysics.H:52
int Get_Qstate_NonMoist_Size() override
get the number of microphysics conserved non-moist (non-water, i.e., other vapor/condensed species) s...
Definition: ERF_EulerianMicrophysics.H:118
void Update_State_Vars_Lev(const int &lev, amrex::MultiFab &cons_in) override
update ERF state variables from microphysics variables
Definition: ERF_EulerianMicrophysics.H:92
virtual void GetPlotVar(const std::string &a_name, amrex::MultiFab &a_mf, const int a_lev) const override
Fills in a MultiFab for plotting.
Definition: ERF_EulerianMicrophysics.H:140
void SetModel()
Create and set the specified moisture model.
Definition: ERF_EulerianMicrophysics.H:151
virtual void GetPlotVarNames(amrex::Vector< std::string > &a_vec) const override
Returns a list of additional plot variable names.
Definition: ERF_EulerianMicrophysics.H:134
void Advance(const int &lev, const amrex::Real &dt_advance, const int &, const amrex::Real &, const SolverChoice &solverChoice, amrex::Vector< amrex::Vector< amrex::MultiFab >> &, const amrex::Vector< std::unique_ptr< amrex::MultiFab >> &, const amrex::GpuArray< ERF_BC, AMREX_SPACEDIM *2 > &) override
Advance the moisture model for one time step.
Definition: ERF_EulerianMicrophysics.H:72
void Init(const int &lev, const amrex::MultiFab &cons_in, const amrex::BoxArray &grids, const amrex::Geometry &geom, const amrex::Real &dt_advance, std::unique_ptr< amrex::MultiFab > &z_phys_nd, std::unique_ptr< amrex::MultiFab > &detJ_cc) override
Initialize the moisture model.
Definition: ERF_EulerianMicrophysics.H:59
int Get_Qmoist_Size(const int &) override
get the number of moisture model variables
Definition: ERF_EulerianMicrophysics.H:106
int Get_Qstate_Moist_Size() override
get the number of microphysics conserved moist (water-related) state variables
Definition: ERF_EulerianMicrophysics.H:112
EulerianMicrophysics()
Null constructor.
Definition: ERF_EulerianMicrophysics.H:21
Base class for microphysics interface.
Definition: ERF_Microphysics.H:14
static MoistureModelType modelType(const MoistureType a_moisture_type)
query if a specified moisture model is Eulerian or Lagrangian
Definition: ERF_Microphysics.H:86
Definition: ERF_DataStruct.H:91