12 #include <AMReX_FArrayBox.H>
13 #include <AMReX_Geometry.H>
14 #include <AMReX_TableData.H>
15 #include <AMReX_MultiFabUtil.H>
60 using FabPtr = std::shared_ptr<amrex::MultiFab>;
84 Init (
const amrex::MultiFab& cons_in,
85 const amrex::BoxArray& grids,
86 const amrex::Geometry& geom,
88 std::unique_ptr<amrex::MultiFab>& z_phys_nd,
89 std::unique_ptr<amrex::MultiFab>& detJ_cc)
override;
133 AMREX_GPU_HOST_DEVICE
142 const amrex::Array4<amrex::Real>& tabs_array,
143 const amrex::Array4<amrex::Real>& pres_array,
144 const amrex::Array4<amrex::Real>& qv_array,
145 const amrex::Array4<amrex::Real>& qc_array,
146 const amrex::Array4<amrex::Real>& qi_array,
147 const amrex::Array4<amrex::Real>& qn_array,
148 const amrex::Array4<amrex::Real>& qt_array)
203 qsat = omn * qsatw + (1.0-omn) * qsati;
204 dqsat = omn * dqsatw + (1.0-omn) * dqsati
205 + domn * qsatw - domn * qsati;
206 lstar = omn * lstarw + (1.0-omn) * lstari;
207 dlstar = domn * lstarw - domn * lstari;
211 fff = -
tabs + tabs_array(i,j,k) + lstar*(qv_array(i,j,k) - qsat);
214 dfff = -1.0 + dlstar*(qv_array(i,j,k) - qsat) - lstar*dqsat;
222 }
while(std::abs(dtabs) > tol && niter < 20);
228 delta_qv = qv_array(i,j,k) - qsat;
229 delta_qc = std::max(-qc_array(i,j,k), delta_qv * omn);
230 delta_qi = std::max(-qi_array(i,j,k), delta_qv * (1.0-omn));
233 qv_array(i,j,k) = qsat;
234 qc_array(i,j,k) += delta_qc;
235 qi_array(i,j,k) += delta_qi;
236 qn_array(i,j,k) = qc_array(i,j,k) + qi_array(i,j,k);
237 qt_array(i,j,k) = qv_array(i,j,k) + qn_array(i,j,k);
245 std::vector<int>& a_idx,
246 std::vector<std::string>& a_names)
const override
255 a_idx.push_back(0); a_names.push_back(
"Nc");
256 a_idx.push_back(1); a_names.push_back(
"Nr");
257 a_idx.push_back(2); a_names.push_back(
"Ni");
258 a_idx.push_back(3); a_names.push_back(
"Ns");
259 a_idx.push_back(4); a_names.push_back(
"Ng");
260 a_idx.push_back(5); a_names.push_back(
"RainAccum");
261 a_idx.push_back(6); a_names.push_back(
"SnowAccum");
262 a_idx.push_back(7); a_names.push_back(
"GraupAccum");
320 amrex::TableData<amrex::Real, 1>
rho1d;
323 amrex::TableData<amrex::Real, 1>
qt1d;
324 amrex::TableData<amrex::Real, 1>
qv1d;
325 amrex::TableData<amrex::Real, 1>
qn1d;
327 amrex::TableData<amrex::Real, 1>
gamaz;
328 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_Morrison.H:58
virtual ~Morrison()=default
amrex::Geometry m_geom
Definition: ERF_Morrison.H:279
Morrison()
Definition: ERF_Morrison.H:64
void Advance(const amrex::Real &dt_advance, const SolverChoice &sc) override
Definition: ERF_AdvanceMorrison.cpp:503
amrex::TableData< amrex::Real, 1 > evapg2
Definition: ERF_Morrison.H:315
void Qmoist_Restart_Vars(const SolverChoice &, std::vector< int > &a_idx, std::vector< std::string > &a_names) const override
Definition: ERF_Morrison.H:244
amrex::TableData< amrex::Real, 1 > qn1d
Definition: ERF_Morrison.H:325
int Qstate_Moist_Size() override
Definition: ERF_Morrison.H:131
int n_qstate_moist_size
Definition: ERF_Morrison.H:270
void Update_Micro_Vars(amrex::MultiFab &cons_in) override
Definition: ERF_Morrison.H:100
amrex::TableData< amrex::Real, 1 > pres1d
Definition: ERF_Morrison.H:321
amrex::TableData< amrex::Real, 1 > zmid
Definition: ERF_Morrison.H:328
amrex::Real m_gOcp
Definition: ERF_Morrison.H:294
amrex::TableData< amrex::Real, 1 > accrgi
Definition: ERF_Morrison.H:312
static constexpr amrex::Real CFL_MAX
Definition: ERF_Morrison.H:273
amrex::TableData< amrex::Real, 1 > accrgc
Definition: ERF_Morrison.H:313
void Copy_State_to_Micro(const amrex::MultiFab &cons_in) override
Definition: ERF_InitMorrison.cpp:98
amrex::TableData< amrex::Real, 1 > tabs1d
Definition: ERF_Morrison.H:322
amrex::Vector< int > MicVarMap
Definition: ERF_Morrison.H:276
amrex::MultiFab * Qmoist_Ptr(const int &varIdx) override
Definition: ERF_Morrison.H:118
amrex::MultiFab * m_z_phys_nd
Definition: ERF_Morrison.H:299
amrex::TableData< amrex::Real, 1 > qt1d
Definition: ERF_Morrison.H:323
int nlev
Definition: ERF_Morrison.H:285
void Compute_Coefficients()
Definition: ERF_InitMorrison.cpp:152
void Copy_Micro_to_State(amrex::MultiFab &cons_in) override
Definition: ERF_UpdateMorrison.cpp:17
void Define(SolverChoice &sc) override
Definition: ERF_Morrison.H:71
amrex::MultiFab * m_detJ_cc
Definition: ERF_Morrison.H:300
amrex::TableData< amrex::Real, 1 > gamaz
Definition: ERF_Morrison.H:327
amrex::TableData< amrex::Real, 1 > qv1d
Definition: ERF_Morrison.H:324
amrex::Real m_fac_cond
Definition: ERF_Morrison.H:291
amrex::Real m_fac_sub
Definition: ERF_Morrison.H:293
amrex::TableData< amrex::Real, 1 > coefice
Definition: ERF_Morrison.H:309
AMREX_GPU_HOST_DEVICE static AMREX_FORCE_INLINE amrex::Real NewtonIterSat(int &i, int &j, int &k, const amrex::Real &fac_cond, const amrex::Real &fac_fus, 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_Morrison.H:136
int m_axis
Definition: ERF_Morrison.H:288
amrex::TableData< amrex::Real, 1 > accrsi
Definition: ERF_Morrison.H:307
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_InitMorrison.cpp:22
int Qmoist_Size() override
Definition: ERF_Morrison.H:128
amrex::TableData< amrex::Real, 1 > accrrc
Definition: ERF_Morrison.H:306
amrex::TableData< amrex::Real, 1 > evapg1
Definition: ERF_Morrison.H:314
int m_qmoist_size
Definition: ERF_Morrison.H:267
amrex::TableData< amrex::Real, 1 > evapr1
Definition: ERF_Morrison.H:316
int zlo
Definition: ERF_Morrison.H:285
amrex::TableData< amrex::Real, 1 > evaps1
Definition: ERF_Morrison.H:310
amrex::BoxArray m_gtoe
Definition: ERF_Morrison.H:282
amrex::TableData< amrex::Real, 1 > evaps2
Definition: ERF_Morrison.H:311
amrex::Real m_rdOcp
Definition: ERF_Morrison.H:295
amrex::Array< FabPtr, MicVar_Morr::NumVars > mic_fab_vars
Definition: ERF_Morrison.H:303
amrex::TableData< amrex::Real, 1 > rho1d
Definition: ERF_Morrison.H:320
std::shared_ptr< amrex::MultiFab > FabPtr
Definition: ERF_Morrison.H:60
int zhi
Definition: ERF_Morrison.H:285
amrex::TableData< amrex::Real, 1 > evapr2
Definition: ERF_Morrison.H:317
void Update_State_Vars(amrex::MultiFab &cons_in) override
Definition: ERF_Morrison.H:107
amrex::TableData< amrex::Real, 1 > accrsc
Definition: ERF_Morrison.H:308
bool m_do_cond
Definition: ERF_Morrison.H:296
amrex::Real m_fac_fus
Definition: ERF_Morrison.H:292
Definition: ERF_NullMoist.H:8
Definition: ERF_Morrison.H:24
@ qv
Definition: ERF_Morrison.H:34
@ ng
Definition: ERF_Morrison.H:48
@ nc
Definition: ERF_Morrison.H:44
@ qpg
Definition: ERF_Morrison.H:41
@ pres
Definition: ERF_Morrison.H:30
@ nr
Definition: ERF_Morrison.H:45
@ qcl
Definition: ERF_Morrison.H:35
@ tabs
Definition: ERF_Morrison.H:29
@ theta
Definition: ERF_Morrison.H:28
@ qp
Definition: ERF_Morrison.H:38
@ ni
Definition: ERF_Morrison.H:46
@ ns
Definition: ERF_Morrison.H:47
@ omega
Definition: ERF_Morrison.H:53
@ qps
Definition: ERF_Morrison.H:40
@ qn
Definition: ERF_Morrison.H:33
@ graup_accum
Definition: ERF_Morrison.H:52
@ rho
Definition: ERF_Morrison.H:27
@ qpr
Definition: ERF_Morrison.H:39
@ qci
Definition: ERF_Morrison.H:36
@ NumVars
Definition: ERF_Morrison.H:54
@ rain_accum
Definition: ERF_Morrison.H:50
@ snow_accum
Definition: ERF_Morrison.H:51
@ qt
Definition: ERF_Morrison.H:32
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