16 #include <AMReX_FArrayBox.H>
17 #include <AMReX_Geometry.H>
18 #include <AMReX_TableData.H>
19 #include <AMReX_MultiFabUtil.H>
55 using FabPtr = std::shared_ptr<amrex::MultiFab>;
91 Init (
const amrex::MultiFab& cons_in,
92 const amrex::BoxArray& grids,
93 const amrex::Geometry& geom,
95 std::unique_ptr<amrex::MultiFab>& z_phys_nd,
96 std::unique_ptr<amrex::MultiFab>& detJ_cc)
override;
148 AMREX_GPU_HOST_DEVICE
152 const int& SAM_moisture_type,
158 const amrex::Array4<amrex::Real>& tabs_array,
159 const amrex::Array4<amrex::Real>& pres_array,
160 const amrex::Array4<amrex::Real>& qv_array,
161 const amrex::Array4<amrex::Real>& qc_array,
162 const amrex::Array4<amrex::Real>& qi_array,
163 const amrex::Array4<amrex::Real>& qn_array,
164 const amrex::Array4<amrex::Real>& qt_array)
203 if (SAM_moisture_type == 1) {
217 }
else if (SAM_moisture_type == 2) {
223 qsat = omn * qsatw + (1.0-omn) * qsati;
224 dqsat = omn * dqsatw + (1.0-omn) * dqsati
225 + domn * qsatw - domn * qsati;
226 lstar = omn * lstarw + (1.0-omn) * lstari;
227 dlstar = domn * lstarw - domn * lstari;
231 fff = -
tabs + tabs_array(i,j,k) + lstar*(qv_array(i,j,k) - qsat);
234 dfff = -1.0 + dlstar*(qv_array(i,j,k) - qsat) - lstar*dqsat;
242 }
while(std::abs(dtabs) > tol && niter < 20);
248 delta_qv = qv_array(i,j,k) - qsat;
249 delta_qc = std::max(-qc_array(i,j,k), delta_qv * omn);
250 delta_qi = std::max(-qi_array(i,j,k), delta_qv * (1.0-omn));
253 qv_array(i,j,k) = qsat;
254 qc_array(i,j,k) += delta_qc;
255 qi_array(i,j,k) += delta_qi;
256 qn_array(i,j,k) = qc_array(i,j,k) + qi_array(i,j,k);
257 qt_array(i,j,k) = qv_array(i,j,k) + qn_array(i,j,k);
265 std::vector<int>& a_idx,
266 std::vector<std::string>& a_names)
const override
272 a_idx.push_back(0); a_names.push_back(
"RainAccum");
273 a_idx.push_back(1); a_names.push_back(
"SnowAccum");
274 a_idx.push_back(2); a_names.push_back(
"GraupAccum");
335 amrex::TableData<amrex::Real, 1>
rho1d;
338 amrex::TableData<amrex::Real, 1>
qt1d;
339 amrex::TableData<amrex::Real, 1>
qv1d;
340 amrex::TableData<amrex::Real, 1>
qn1d;
342 amrex::TableData<amrex::Real, 1>
gamaz;
343 amrex::TableData<amrex::Real, 1>
zmid;
constexpr amrex::Real lfus
Definition: ERF_Constants.H:67
constexpr amrex::Real tbgmax
Definition: ERF_Constants.H:32
constexpr amrex::Real lsub
Definition: ERF_Constants.H:68
constexpr amrex::Real lcond
Definition: ERF_Constants.H:66
constexpr amrex::Real tbgmin
Definition: ERF_Constants.H:31
constexpr amrex::Real CONST_GRAV
Definition: ERF_Constants.H:21
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_dtqsatw(amrex::Real t, amrex::Real p, amrex::Real &dtqsatw)
Definition: ERF_MicrophysicsUtils.H:178
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_dtqsati(amrex::Real t, amrex::Real p, amrex::Real &dtqsati)
Definition: ERF_MicrophysicsUtils.H:170
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_qsatw(amrex::Real t, amrex::Real p, amrex::Real &qsatw)
Definition: ERF_MicrophysicsUtils.H:163
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE void erf_qsati(amrex::Real t, amrex::Real p, amrex::Real &qsati)
Definition: ERF_MicrophysicsUtils.H:152
amrex::Real Real
Definition: ERF_ShocInterface.H:19
Definition: ERF_NullMoist.H:8
void Copy_Micro_to_State(amrex::MultiFab &cons_in) override
Definition: ERF_UpdateSAM.cpp:15
void Advance(const amrex::Real &dt_advance, const SolverChoice &sc) override
Definition: ERF_SAM.H:121
int zlo
Definition: ERF_SAM.H:300
void Precip(const SolverChoice &sc)
Definition: ERF_Precip.cpp:10
amrex::TableData< amrex::Real, 1 > rho1d
Definition: ERF_SAM.H:335
amrex::TableData< amrex::Real, 1 > evapr2
Definition: ERF_SAM.H:332
SAM()
Definition: ERF_SAM.H:59
amrex::MultiFab * m_detJ_cc
Definition: ERF_SAM.H:315
void Copy_State_to_Micro(const amrex::MultiFab &cons_in) override
Definition: ERF_InitSAM.cpp:87
amrex::TableData< amrex::Real, 1 > evapg1
Definition: ERF_SAM.H:329
void Update_Micro_Vars(amrex::MultiFab &cons_in) override
Definition: ERF_SAM.H:107
amrex::TableData< amrex::Real, 1 > accrgi
Definition: ERF_SAM.H:327
std::shared_ptr< amrex::MultiFab > FabPtr
Definition: ERF_SAM.H:55
amrex::TableData< amrex::Real, 1 > evapg2
Definition: ERF_SAM.H:330
void IceFall(const SolverChoice &sc)
Definition: ERF_IceFall.cpp:11
amrex::Array< FabPtr, MicVar::NumVars > mic_fab_vars
Definition: ERF_SAM.H:318
int m_axis
Definition: ERF_SAM.H:303
amrex::Real m_gOcp
Definition: ERF_SAM.H:309
amrex::Real m_rdOcp
Definition: ERF_SAM.H:310
void Cloud(const SolverChoice &sc)
Definition: ERF_CloudSAM.cpp:12
int m_qmoist_size
Definition: ERF_SAM.H:279
amrex::MultiFab * m_z_phys_nd
Definition: ERF_SAM.H:314
amrex::Real m_fac_fus
Definition: ERF_SAM.H:307
amrex::TableData< amrex::Real, 1 > evaps1
Definition: ERF_SAM.H:325
amrex::BoxArray m_gtoe
Definition: ERF_SAM.H:294
amrex::TableData< amrex::Real, 1 > zmid
Definition: ERF_SAM.H:343
amrex::Vector< int > MicVarMap
Definition: ERF_SAM.H:288
bool m_do_cond
Definition: ERF_SAM.H:311
amrex::TableData< amrex::Real, 1 > qv1d
Definition: ERF_SAM.H:339
amrex::Real m_fac_cond
Definition: ERF_SAM.H:306
amrex::TableData< amrex::Real, 1 > accrsi
Definition: ERF_SAM.H:322
amrex::TableData< amrex::Real, 1 > accrrc
Definition: ERF_SAM.H:321
void Compute_Coefficients()
Definition: ERF_InitSAM.cpp:139
amrex::TableData< amrex::Real, 1 > pres1d
Definition: ERF_SAM.H:336
int nlev
Definition: ERF_SAM.H:300
int Qstate_Moist_Size() override
Definition: ERF_SAM.H:146
amrex::Geometry m_geom
Definition: ERF_SAM.H:291
amrex::TableData< amrex::Real, 1 > evapr1
Definition: ERF_SAM.H:331
amrex::Real m_fac_sub
Definition: ERF_SAM.H:308
amrex::TableData< amrex::Real, 1 > evaps2
Definition: ERF_SAM.H:326
void Init(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
Definition: ERF_InitSAM.cpp:22
amrex::TableData< amrex::Real, 1 > qn1d
Definition: ERF_SAM.H:340
amrex::TableData< amrex::Real, 1 > accrgc
Definition: ERF_SAM.H:328
amrex::TableData< amrex::Real, 1 > accrsc
Definition: ERF_SAM.H:323
int n_qstate_moist_size
Definition: ERF_SAM.H:282
amrex::TableData< amrex::Real, 1 > qt1d
Definition: ERF_SAM.H:338
amrex::TableData< amrex::Real, 1 > coefice
Definition: ERF_SAM.H:324
int zhi
Definition: ERF_SAM.H:300
static constexpr amrex::Real CFL_MAX
Definition: ERF_SAM.H:285
void Define(SolverChoice &sc) override
Definition: ERF_SAM.H:78
void Qmoist_Restart_Vars(const SolverChoice &, std::vector< int > &a_idx, std::vector< std::string > &a_names) const override
Definition: ERF_SAM.H:264
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real NewtonIterSat(int &i, int &j, int &k, const int &SAM_moisture_type, const amrex::Real &fac_cond, const amrex::Real &, const amrex::Real &fac_sub, const amrex::Real &an, const amrex::Real &bn, const amrex::Array4< amrex::Real > &tabs_array, const amrex::Array4< amrex::Real > &pres_array, const amrex::Array4< amrex::Real > &qv_array, const amrex::Array4< amrex::Real > &qc_array, const amrex::Array4< amrex::Real > &qi_array, const amrex::Array4< amrex::Real > &qn_array, const amrex::Array4< amrex::Real > &qt_array)
Definition: ERF_SAM.H:151
int Qmoist_Size() override
Definition: ERF_SAM.H:143
void PrecipFall(const SolverChoice &sc)
Definition: ERF_PrecipFall.cpp:16
void Update_State_Vars(amrex::MultiFab &cons_in) override
Definition: ERF_SAM.H:114
amrex::TableData< amrex::Real, 1 > tabs1d
Definition: ERF_SAM.H:337
amrex::TableData< amrex::Real, 1 > gamaz
Definition: ERF_SAM.H:342
amrex::MultiFab * Qmoist_Ptr(const int &varIdx) override
Definition: ERF_SAM.H:133
amrex::Real dt
Definition: ERF_SAM.H:297
@ NumVars
Definition: ERF_SAM.H:49
@ pres
Definition: ERF_SAM.H:33
@ qpr
Definition: ERF_SAM.H:42
@ qci
Definition: ERF_SAM.H:39
@ qv
Definition: ERF_SAM.H:37
@ rho
Definition: ERF_SAM.H:30
@ snow_accum
Definition: ERF_SAM.H:47
@ qpg
Definition: ERF_SAM.H:44
@ qt
Definition: ERF_SAM.H:35
@ qn
Definition: ERF_SAM.H:36
@ rain_accum
Definition: ERF_SAM.H:46
@ theta
Definition: ERF_SAM.H:31
@ qcl
Definition: ERF_SAM.H:38
@ graup_accum
Definition: ERF_SAM.H:48
@ tabs
Definition: ERF_SAM.H:32
@ qps
Definition: ERF_SAM.H:43
@ qp
Definition: ERF_SAM.H:41
Definition: ERF_DataStruct.H:123
amrex::Real rdOcp
Definition: ERF_DataStruct.H:866
amrex::Real c_p
Definition: ERF_DataStruct.H:865
bool use_shoc
Definition: ERF_DataStruct.H:895
int ave_plane
Definition: ERF_DataStruct.H:926