1 #ifndef ERF_SURFACELAYER_H
2 #define ERF_SURFACELAYER_H
4 #include "AMReX_Geometry.H"
5 #include "AMReX_ParmParse.H"
6 #include "AMReX_FArrayBox.H"
7 #include "AMReX_MultiFab.H"
8 #include "AMReX_iMultiFab.H"
9 #include "AMReX_MFInterpolater.H"
36 bool& use_rot_surface_flux,
37 std::string a_pp_prefix,
38 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Qv_prim,
39 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& z_phys_nd,
40 const TerrainType& a_terrain_type,
45 m_ma(geom, (z_phys_nd[0] != nullptr), a_pp_prefix, a_terrain_type)
51 amrex::ParmParse
pp(
"erf");
55 if (use_rot_surface_flux) {
58 std::string flux_string_in;
59 std::string flux_string{
"moeng"};
60 auto read_flux =
pp.query(
"surface_layer.flux_type", flux_string_in);
62 flux_string = amrex::toLower(flux_string_in);
64 if (flux_string ==
"donelan") {
66 }
else if (flux_string ==
"moeng") {
68 }
else if (flux_string ==
"custom") {
71 amrex::Abort(
"Undefined MOST flux type!");
78 std::string pblh_string_in;
79 std::string pblh_string{
"none"};
80 auto read_pblh =
pp.query(
"most.pblh_calc", pblh_string_in);
82 pblh_string = amrex::toLower(pblh_string_in);
84 if (pblh_string ==
"none") {
86 }
else if (pblh_string ==
"mynn25") {
88 }
else if (pblh_string ==
"mynnedmf") {
90 }
else if (pblh_string ==
"ysu") {
92 }
else if (pblh_string ==
"mrf") {
95 amrex::Abort(
"Undefined PBLH calc type!");
99 auto erf_st =
pp.query(
"most.surf_temp",
surf_temp);
117 "Specified custom MOST qv flux without moisture model!");
119 amrex::Print() <<
"Using specified ustar, tstar, qstar for MOST = "
130 amrex::Abort(
"Can only specify one of surf_temp_flux or surf_heating_rate");
136 amrex::Abort(
"Can only specify one of surf_temp_flux or surf_heating_rate");
161 std::string bogus_input;
162 if (
pp.query(
"most.roughness_type", bogus_input) > 0) {
163 amrex::Abort(
"most.roughness_type is deprecated; use "
164 "most.roughness_type_land and/or most.roughness_type_sea");
168 std::string rough_land_string_in;
169 std::string rough_land_string{
"constant"};
170 auto read_rough_land =
171 pp.query(
"most.roughness_type_land", rough_land_string_in);
172 if (read_rough_land) {
173 rough_land_string = amrex::toLower(rough_land_string_in);
175 if (rough_land_string ==
"constant") {
178 amrex::Abort(
"Undefined MOST roughness type for land!");
182 std::string rough_sea_string_in;
183 std::string rough_sea_string{
"charnock"};
184 auto read_rough_sea =
pp.query(
"most.roughness_type_sea", rough_sea_string_in);
185 if (read_rough_sea) {
186 rough_sea_string = amrex::toLower(rough_sea_string_in);
188 if (rough_sea_string ==
"charnock") {
190 pp.query(
"most.charnock_constant",
cnk_a);
191 pp.query(
"most.charnock_viscosity",
cnk_visc);
193 amrex::Print() <<
"If there is water, Charnock relation with C_a="
195 <<
" will be used" << std::endl;
197 amrex::Print() <<
"If there is water, Charnock relation with variable "
198 "Charnock parameter (COARE3.0)"
199 << (
cnk_visc ?
" and viscosity" :
"") <<
" will be used"
202 }
else if (rough_sea_string ==
"coare3.0") {
204 amrex::Print() <<
"If there is water, Charnock relation with variable "
205 "Charnock parameter (COARE3.0)"
206 << (
cnk_visc ?
" and viscosity" :
"") <<
" will be used"
209 }
else if (rough_sea_string ==
"donelan") {
211 }
else if (rough_sea_string ==
"modified_charnock") {
213 pp.query(
"most.modified_charnock_depth",
depth);
214 }
else if (rough_sea_string ==
"wave_coupled") {
216 }
else if (rough_sea_string ==
"constant") {
219 amrex::Abort(
"Undefined MOST roughness type for sea!");
225 const amrex::Vector<amrex::MultiFab*>& mfv,
226 std::unique_ptr<amrex::MultiFab>& Theta_prim,
227 std::unique_ptr<amrex::MultiFab>& Qv_prim,
228 std::unique_ptr<amrex::MultiFab>& Qr_prim,
229 std::unique_ptr<amrex::MultiFab>& z_phys_nd,
230 amrex::MultiFab* Hwave,
231 amrex::MultiFab* Lwave,
232 amrex::MultiFab* eddyDiffs,
233 amrex::Vector<amrex::MultiFab*> lsm_data,
234 amrex::Vector<amrex::MultiFab*> lsm_flux,
235 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& sst_lev,
236 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& tsk_lev,
237 amrex::Vector<std::unique_ptr<amrex::iMultiFab>>& lmask_lev)
241 Theta_prim, Qv_prim, Qr_prim,
245 amrex::MultiFab& mf = *(mfv[0]);
247 amrex::ParmParse
pp(
"erf");
251 const int nghost = 0;
253 amrex::ParallelDescriptor::ReduceIntMin(lmask_min);
257 std::string rough_sea_string{
"charnock"};
258 pp.query(
"most.roughness_type_sea", rough_sea_string);
259 amrex::Print() <<
"Variable sea roughness (type " << rough_sea_string
264 if ((lev == 0) || (lev > nlevs - 1)) {
289 int nt_tot_sst = sst_lev.size();
291 for (
int nt(0); nt < nt_tot_sst; ++nt) {
294 int nt_tot_tsk =
static_cast<int>(tsk_lev.size());
296 for (
int nt(0); nt < nt_tot_tsk; ++nt) {
299 int nt_tot_lmask =
static_cast<int>(lmask_lev.size());
301 for (
int nt(0); nt < nt_tot_lmask; ++nt) {
311 int nvar =
static_cast<int>(lsm_data.size());
314 for (
int n(0); n < nvar; ++n) {
321 bool read_z0 =
false;
324 read_z0 =
pp.query(
"most.roughness_file_name", fname);
330 amrex::BoxArray ba = mf.boxArray();
331 amrex::BoxList bl2d = ba.boxList();
332 for (
auto& b : bl2d) {
335 amrex::BoxArray ba2d(std::move(bl2d));
336 const amrex::DistributionMapping& dm = mf.DistributionMap();
338 amrex::IntVect
ng = mf.nGrowVect();
343 z_0[lev].define(ba2d, dm, ncomp,
ng);
351 u_star[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
352 u_star[lev]->setVal(1.E34);
354 w_star[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
355 w_star[lev]->setVal(1.E34);
357 t_star[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
360 q_star[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
363 olen[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
364 olen[lev]->setVal(1.E34);
366 pblh[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
367 pblh[lev]->setVal(1.E34);
369 t_surf[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
372 q_surf[lev] = std::make_unique<amrex::MultiFab>(ba2d, dm, ncomp,
ng);
390 const amrex::MultiFab& cons_in,
391 const std::unique_ptr<amrex::MultiFab>& z_phys_nd,
392 int max_iters = 100);
394 template <
typename FluxIter>
396 const int& max_iters,
397 const FluxIter& most_flux,
401 amrex::Vector<const amrex::MultiFab*> mfs,
402 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Tau_lev,
403 amrex::MultiFab* xheat_flux,
404 amrex::MultiFab* yheat_flux,
405 amrex::MultiFab* zheat_flux,
406 amrex::MultiFab* xqv_flux,
407 amrex::MultiFab* yqv_flux,
408 amrex::MultiFab* zqv_flux,
409 const amrex::MultiFab* z_phys);
411 template <
typename FluxCalc>
413 amrex::Vector<const amrex::MultiFab*> mfs,
414 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Tau_lev,
415 amrex::MultiFab* xheat_flux,
416 amrex::MultiFab* yheat_flux,
417 amrex::MultiFab* zheat_flux,
418 amrex::MultiFab* xqv_flux,
419 amrex::MultiFab* yqv_flux,
420 amrex::MultiFab* zqv_flux,
421 const amrex::MultiFab* z_phys,
422 const FluxCalc& flux_comp);
428 const amrex::MultiFab& cons_in,
429 const std::unique_ptr<amrex::MultiFab>& z_phys_nd);
435 amrex::Vector<amrex::Vector<amrex::MultiFab>>& vars,
436 amrex::MultiFab* z_phys_cc,
439 template <
typename PBLHeightEstimator>
441 amrex::Vector<amrex::Vector<amrex::MultiFab>>& vars,
442 amrex::MultiFab* z_phys_cc,
443 const PBLHeightEstimator& est,
447 const std::string& fname);
452 int nlevs =
static_cast<int>(
m_geom.size());
453 for (
int lev = 0; lev < nlevs; lev++) {
455 amrex::Print() <<
"Surface temp at t=" << time <<
": "
462 amrex::Vector<amrex::Vector<amrex::MultiFab>>& vars_old,
463 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Theta_prim,
464 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Qv_prim,
465 amrex::Vector<std::unique_ptr<amrex::MultiFab>>& Qr_prim)
493 amrex::MultiFab*
get_z0 (
const int& lev) {
return &
z_0[lev]; }
502 int lmask_min = amrex::ReduceMin(lmask, nghost, [=] AMREX_GPU_HOST_DEVICE(
503 amrex::Box
const& bx, amrex::Array4<int const>
const& lm_arr) ->
int
505 int locmin = std::numeric_limits<int>::max();
506 const auto lo = lbound(bx);
507 const auto hi = ubound(bx);
508 for (
int j = lo.y; j <= hi.y; ++j) {
509 for (
int i = lo.x; i <= hi.x; ++i) {
510 locmin = std::min(locmin, lm_arr(i, j, 0));
578 amrex::Vector<amrex::MultiFab>
z_0;
584 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
u_star;
585 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
w_star;
586 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
t_star;
587 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
q_star;
588 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
olen;
589 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
pblh;
590 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
t_surf;
591 amrex::Vector<std::unique_ptr<amrex::MultiFab>>
q_surf;
593 amrex::Vector<amrex::Vector<amrex::MultiFab*>>
m_sst_lev;
594 amrex::Vector<amrex::Vector<amrex::MultiFab*>>
m_tsk_lev;
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real pp(amrex::Real y)
Definition: ERF_MicrophysicsUtils.H:230
amrex::Real Real
Definition: ERF_ShocInterface.H:16
Definition: ERF_MOSTAverage.H:14
amrex::Real get_zref() const
Definition: ERF_MOSTAverage.H:104
const amrex::MultiFab * get_average(const int &lev, const int &comp) const
Definition: ERF_MOSTAverage.H:101
void update_field_ptrs(const int &lev, amrex::Vector< amrex::Vector< amrex::MultiFab >> &vars_old, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Theta_prim, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Qv_prim, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Qr_prim)
Definition: ERF_MOSTAverage.cpp:246
void make_MOSTAverage_at_level(const int &lev, const amrex::Vector< amrex::MultiFab * > &vars_old, std::unique_ptr< amrex::MultiFab > &Theta_prim, std::unique_ptr< amrex::MultiFab > &Qv_prim, std::unique_ptr< amrex::MultiFab > &Qr_prim, std::unique_ptr< amrex::MultiFab > &z_phys_nd)
Definition: ERF_MOSTAverage.cpp:70
Definition: ERF_SurfaceLayer.H:31
ThetaCalcType theta_type
Definition: ERF_SurfaceLayer.H:549
int lmask_min_reduce(amrex::iMultiFab &lmask, const int &nghost)
Definition: ERF_SurfaceLayer.H:499
bool m_include_wstar
Definition: ERF_SurfaceLayer.H:561
bool specified_rho_surf
Definition: ERF_SurfaceLayer.H:574
bool m_rotate
Definition: ERF_SurfaceLayer.H:558
PBLHeightCalcType pblh_type
Definition: ERF_SurfaceLayer.H:553
amrex::Vector< amrex::Vector< amrex::iMultiFab * > > m_lmask_lev
Definition: ERF_SurfaceLayer.H:595
amrex::iMultiFab * get_lmask(const int &lev)
Definition: ERF_SurfaceLayer.H:497
bool use_moisture
Definition: ERF_SurfaceLayer.H:581
amrex::MultiFab * get_q_surf(const int &lev)
Definition: ERF_SurfaceLayer.H:489
amrex::MultiFab * get_w_star(const int &lev)
Definition: ERF_SurfaceLayer.H:472
void update_surf_temp(const amrex::Real &time)
Definition: ERF_SurfaceLayer.H:449
RoughCalcType rough_type_land
Definition: ERF_SurfaceLayer.H:551
void update_pblh(const int &lev, amrex::Vector< amrex::Vector< amrex::MultiFab >> &vars, amrex::MultiFab *z_phys_cc, const MoistureComponentIndices &moisture_indices)
Definition: ERF_SurfaceLayer.cpp:609
amrex::Vector< std::unique_ptr< amrex::MultiFab > > t_surf
Definition: ERF_SurfaceLayer.H:590
amrex::Real z0_const
Definition: ERF_SurfaceLayer.H:562
amrex::Real cnk_a
Definition: ERF_SurfaceLayer.H:575
amrex::Real surf_temp
Definition: ERF_SurfaceLayer.H:564
void update_mac_ptrs(const int &lev, amrex::Vector< amrex::Vector< amrex::MultiFab >> &vars_old, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Theta_prim, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Qv_prim, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Qr_prim)
Definition: ERF_SurfaceLayer.H:461
void make_SurfaceLayer_at_level(const int &lev, int nlevs, const amrex::Vector< amrex::MultiFab * > &mfv, std::unique_ptr< amrex::MultiFab > &Theta_prim, std::unique_ptr< amrex::MultiFab > &Qv_prim, std::unique_ptr< amrex::MultiFab > &Qr_prim, std::unique_ptr< amrex::MultiFab > &z_phys_nd, amrex::MultiFab *Hwave, amrex::MultiFab *Lwave, amrex::MultiFab *eddyDiffs, amrex::Vector< amrex::MultiFab * > lsm_data, amrex::Vector< amrex::MultiFab * > lsm_flux, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &sst_lev, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &tsk_lev, amrex::Vector< std::unique_ptr< amrex::iMultiFab >> &lmask_lev)
Definition: ERF_SurfaceLayer.H:223
void compute_pblh(const int &lev, amrex::Vector< amrex::Vector< amrex::MultiFab >> &vars, amrex::MultiFab *z_phys_cc, const PBLHeightEstimator &est, const MoistureComponentIndices &moisture_indice)
amrex::Vector< std::unique_ptr< amrex::MultiFab > > q_star
Definition: ERF_SurfaceLayer.H:587
amrex::Vector< amrex::MultiFab * > m_Lwave_lev
Definition: ERF_SurfaceLayer.H:599
void get_lsm_tsurf(const int &lev)
Definition: ERF_SurfaceLayer.cpp:576
void fill_qsurf_with_qsat(const int &lev, const amrex::MultiFab &cons_in, const std::unique_ptr< amrex::MultiFab > &z_phys_nd)
Definition: ERF_SurfaceLayer.cpp:536
amrex::MultiFab * get_olen(const int &lev)
Definition: ERF_SurfaceLayer.H:478
amrex::Vector< amrex::MultiFab > z_0
Definition: ERF_SurfaceLayer.H:578
amrex::Real surf_moist_flux
Definition: ERF_SurfaceLayer.H:569
void compute_SurfaceLayer_bcs(const int &lev, amrex::Vector< const amrex::MultiFab * > mfs, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Tau_lev, amrex::MultiFab *xheat_flux, amrex::MultiFab *yheat_flux, amrex::MultiFab *zheat_flux, amrex::MultiFab *xqv_flux, amrex::MultiFab *yqv_flux, amrex::MultiFab *zqv_flux, const amrex::MultiFab *z_phys, const FluxCalc &flux_comp)
RoughCalcType rough_type_sea
Definition: ERF_SurfaceLayer.H:552
amrex::Real surf_moist
Definition: ERF_SurfaceLayer.H:568
amrex::Vector< std::unique_ptr< amrex::MultiFab > > w_star
Definition: ERF_SurfaceLayer.H:585
amrex::MultiFab * get_u_star(const int &lev)
Definition: ERF_SurfaceLayer.H:470
amrex::Real m_bdy_time_interval
Definition: ERF_SurfaceLayer.H:559
amrex::Vector< amrex::Vector< amrex::MultiFab * > > m_lsm_data_lev
Definition: ERF_SurfaceLayer.H:596
amrex::Real custom_qstar
Definition: ERF_SurfaceLayer.H:572
amrex::Vector< std::unique_ptr< amrex::MultiFab > > u_star
Definition: ERF_SurfaceLayer.H:584
amrex::Real custom_rhosurf
Definition: ERF_SurfaceLayer.H:573
amrex::Vector< std::unique_ptr< amrex::MultiFab > > q_surf
Definition: ERF_SurfaceLayer.H:591
amrex::Vector< amrex::Vector< amrex::MultiFab * > > m_sst_lev
Definition: ERF_SurfaceLayer.H:593
FluxCalcType
Definition: ERF_SurfaceLayer.H:519
@ MOENG
Moeng functional form.
@ CUSTOM
Custom constant flux functional form.
@ ROTATE
Terrain rotation flux functional form.
@ DONELAN
Donelan functional form.
MoistCalcType
Definition: ERF_SurfaceLayer.H:532
@ SURFACE_MOISTURE
Surface Qv specified.
@ MOISTURE_FLUX
Qv-flux specified.
amrex::Real depth
Definition: ERF_SurfaceLayer.H:577
amrex::Vector< amrex::MultiFab * > m_Hwave_lev
Definition: ERF_SurfaceLayer.H:598
amrex::Real default_land_surf_moist
Definition: ERF_SurfaceLayer.H:567
bool m_var_z0
Definition: ERF_SurfaceLayer.H:579
amrex::MultiFab * get_t_star(const int &lev)
Definition: ERF_SurfaceLayer.H:474
bool have_variable_sea_roughness()
Definition: ERF_SurfaceLayer.H:495
void compute_fluxes(const int &lev, const int &max_iters, const FluxIter &most_flux, bool is_land)
Definition: ERF_SurfaceLayer.cpp:190
amrex::MultiFab * get_q_star(const int &lev)
Definition: ERF_SurfaceLayer.H:476
amrex::Vector< amrex::Vector< amrex::MultiFab * > > m_lsm_flux_lev
Definition: ERF_SurfaceLayer.H:597
PBLHeightCalcType
Definition: ERF_SurfaceLayer.H:546
amrex::MultiFab * get_pblh(const int &lev)
Definition: ERF_SurfaceLayer.H:480
void fill_tsurf_with_sst_and_tsk(const int &lev, const amrex::Real &time)
Definition: ERF_SurfaceLayer.cpp:464
amrex::Real surf_temp_flux
Definition: ERF_SurfaceLayer.H:566
amrex::Vector< amrex::Geometry > m_geom
Definition: ERF_SurfaceLayer.H:557
amrex::Vector< std::unique_ptr< amrex::MultiFab > > t_star
Definition: ERF_SurfaceLayer.H:586
amrex::MultiFab * get_z0(const int &lev)
Definition: ERF_SurfaceLayer.H:493
amrex::Vector< amrex::Vector< amrex::MultiFab * > > m_tsk_lev
Definition: ERF_SurfaceLayer.H:594
amrex::Real custom_tstar
Definition: ERF_SurfaceLayer.H:571
bool cnk_visc
Definition: ERF_SurfaceLayer.H:576
amrex::Real surf_heating_rate
Definition: ERF_SurfaceLayer.H:565
RoughCalcType
Definition: ERF_SurfaceLayer.H:538
FluxCalcType flux_type
Definition: ERF_SurfaceLayer.H:548
MoistCalcType moist_type
Definition: ERF_SurfaceLayer.H:550
void impose_SurfaceLayer_bcs(const int &lev, amrex::Vector< const amrex::MultiFab * > mfs, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Tau_lev, amrex::MultiFab *xheat_flux, amrex::MultiFab *yheat_flux, amrex::MultiFab *zheat_flux, amrex::MultiFab *xqv_flux, amrex::MultiFab *yqv_flux, amrex::MultiFab *zqv_flux, const amrex::MultiFab *z_phys)
Definition: ERF_SurfaceLayer.cpp:257
amrex::Vector< amrex::MultiFab * > m_eddyDiffs_lev
Definition: ERF_SurfaceLayer.H:600
const amrex::MultiFab * get_mac_avg(const int &lev, int comp)
Definition: ERF_SurfaceLayer.H:482
amrex::Real custom_ustar
Definition: ERF_SurfaceLayer.H:570
amrex::Vector< std::unique_ptr< amrex::MultiFab > > olen
Definition: ERF_SurfaceLayer.H:588
amrex::MultiFab * get_t_surf(const int &lev)
Definition: ERF_SurfaceLayer.H:487
amrex::Vector< std::unique_ptr< amrex::MultiFab > > pblh
Definition: ERF_SurfaceLayer.H:589
amrex::Real default_land_surf_temp
Definition: ERF_SurfaceLayer.H:563
SurfaceLayer(const amrex::Vector< amrex::Geometry > &geom, bool &use_rot_surface_flux, std::string a_pp_prefix, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &Qv_prim, amrex::Vector< std::unique_ptr< amrex::MultiFab >> &z_phys_nd, const TerrainType &a_terrain_type, amrex::Real bdy_time_interval=0.0)
Definition: ERF_SurfaceLayer.H:35
amrex::Real get_zref()
Definition: ERF_SurfaceLayer.H:491
void update_fluxes(const int &lev, const amrex::Real &time, const amrex::MultiFab &cons_in, const std::unique_ptr< amrex::MultiFab > &z_phys_nd, int max_iters=100)
Definition: ERF_SurfaceLayer.cpp:12
ThetaCalcType
Definition: ERF_SurfaceLayer.H:526
@ SURFACE_TEMPERATURE
Surface temperature specified.
@ HEAT_FLUX
Heat-flux specified.
void read_custom_roughness(const int &lev, const std::string &fname)
Definition: ERF_SurfaceLayer.cpp:636
MOSTAverage m_ma
Definition: ERF_SurfaceLayer.H:583
@ ng
Definition: ERF_Morrison.H:48
real(c_double), parameter epsilon
Definition: ERF_module_model_constants.F90:12
Definition: ERF_DataStruct.H:99