1 #ifndef ERF_LANDSURFACE_H
2 #define ERF_LANDSURFACE_H
4 #include <AMReX_PlotFileUtil.H>
5 #include <AMReX_MultiFabUtil.H>
26 template<
class NewSurfModel>
31 m_lsm_model[lev] = std::make_unique<NewSurfModel>();
44 const amrex::MultiFab& cons_in,
45 const amrex::Geometry& geom,
46 const amrex::Real& dt_advance)
53 amrex::MultiFab& cons_in,
54 amrex::MultiFab& xvel_in,
55 amrex::MultiFab& yvel_in,
56 amrex::MultiFab* hfx3_out,
57 amrex::MultiFab* qfx3_out,
58 const amrex::Real& dt_advance,
61 m_lsm_model[lev]->Advance(lev, cons_in, xvel_in, yvel_in, hfx3_out, qfx3_out, dt_advance, nstep);
65 Advance (
const int& lev,
const amrex::Real& dt_advance)
99 const amrex::Vector<int> &level_steps,
100 const amrex::Vector<amrex::IntVect> &ref_ratio)
104 amrex::Vector<std::string> varnames;
108 varnames.resize(nvar);
111 std::string plotfilename = amrex::Concatenate(
plot_file_lsm, level_steps[0], 5);
112 for (
int lev(0); lev<nlev; ++lev) {
116 amrex::IntVect ng(0,0,0);
117 amrex::BoxArray ba = mf_lsm->boxArray();
118 amrex::DistributionMapping dm = mf_lsm->DistributionMap();
120 for (
int n(0); n<nvar; ++n) {
126 WriteMultiLevelPlotfile (plotfilename, nlev, GetVecOfConstPtrs(
m_lsm_data_lev),
Definition: ERF_LandSurface.H:15
amrex::Vector< amrex::Geometry > m_lsm_geom_lev
Definition: ERF_LandSurface.H:141
std::string Get_VarName(const int &varIdx)
Definition: ERF_LandSurface.H:95
void Advance(const int &lev, amrex::MultiFab &cons_in, amrex::MultiFab &xvel_in, amrex::MultiFab &yvel_in, amrex::MultiFab *hfx3_out, amrex::MultiFab *qfx3_out, const amrex::Real &dt_advance, const int &nstep)
Definition: ERF_LandSurface.H:52
void Advance(const int &lev, const amrex::Real &dt_advance)
Definition: ERF_LandSurface.H:65
int Get_Data_Size()
Definition: ERF_LandSurface.H:92
amrex::Geometry Get_Lsm_Geom(const int &lev)
Definition: ERF_LandSurface.H:89
void ReSize(const int &nlev)
Definition: ERF_LandSurface.H:24
amrex::Vector< amrex::MultiFab > m_lsm_data_lev
Definition: ERF_LandSurface.H:144
void Plot_Lsm_Data(amrex::Real time, const amrex::Vector< int > &level_steps, const amrex::Vector< amrex::IntVect > &ref_ratio)
Definition: ERF_LandSurface.H:98
amrex::Vector< std::unique_ptr< NullSurf > > m_lsm_model
Definition: ERF_LandSurface.H:135
void SetModel()
Definition: ERF_LandSurface.H:28
amrex::MultiFab * Get_Flux_Ptr(const int &lev, const int &varIdx)
Definition: ERF_LandSurface.H:86
void Init(const int &lev, const amrex::MultiFab &cons_in, const amrex::Geometry &geom, const amrex::Real &dt_advance)
Definition: ERF_LandSurface.H:43
LandSurface()
Definition: ERF_LandSurface.H:19
void Update_State_Vars_Lev(const int &lev, amrex::MultiFab &cons_in)
Definition: ERF_LandSurface.H:77
void Define(const int &lev, SolverChoice &sc)
Definition: ERF_LandSurface.H:36
amrex::MultiFab * Get_Data_Ptr(const int &lev, const int &varIdx)
Definition: ERF_LandSurface.H:83
std::string plot_file_lsm
Definition: ERF_LandSurface.H:138
void Update_Micro_Vars_Lev(const int &lev, amrex::MultiFab &cons_in)
Definition: ERF_LandSurface.H:71
Definition: ERF_DataStruct.H:86