103 amrex::ParmParse
pp(pp_prefix);
106 if (
pp.query(
"use_terrain",
bogus) > 0) {
107 amrex::Abort(
"The input use_terrain is deprecated. Set terrain_type instead.");
113 "The grid stretching ratio must be greater than 1");
120 amrex::Print() <<
"Turning terrain on to enable grid stretching" << std::endl;
125 amrex::Print() <<
"Nominal zsurface height != 0, may result in unexpected behavior"
128 pp.get(
"initial_dz",
dz0);
137 amrex::Abort(
"The input moisture_type is deprecated. Set moisture_model instead.");
147 }
else if ( (
moisture_type == MoistureType::SAM_NoPrecip_NoIce) ||
180 pp.query_enum_case_insensitive(
"land_surface_model",
lsm_type);
184 pp.query_enum_case_insensitive(
"radiation_model",
rad_type);
187 #ifndef ERF_USE_RRTMGP
188 if (
rad_type == RadiationType::RRTMGP)
190 amrex::Abort(
"ERF was not compiled with RRTMGP enabled!");
195 std::string terrain_type_temp =
"";
196 pp.query(
"terrain_type", terrain_type_temp);
197 if (terrain_type_temp ==
"Moving") {
198 amrex::Warning(
"erf.terrain_type = Moving is deprecated; please replace Moving by MovingFittedMesh");
200 }
else if (terrain_type_temp ==
"Static") {
201 amrex::Warning(
"erf.terrain_type = Static is deprecated; please replace Static by StaticFittedMesh");
210 std::string init_type_temp_string;
211 pp.query(
"init_type",init_type_temp_string);
212 if ( (init_type_temp_string ==
"Real") || (init_type_temp_string ==
"real") ) {
213 amrex::Warning(
"erf.init_type = Real is deprecated; please replace Real by WRFInput");
216 }
else if ( (init_type_temp_string ==
"Ideal") || (init_type_temp_string ==
"ideal") ) {
217 amrex::Warning(
"erf.init_type = Ideal is deprecated; please replace Ideal by WRFInput");
221 pp.query_enum_case_insensitive(
"init_type",
init_type);
227 amrex::Abort(
"Only terrain_type = StaticFittedMesh are allowed with init_type = WRFInput or Metgrid");
233 amrex::Abort(
"Can't have moisture_type = None with init_type = WRFInput");
237 std::string nc_bdy_file_temp_string;
238 bool has_bdy =
pp.query(
"nc_bdy_file", nc_bdy_file_temp_string);
242 pp.query(
"use_real_bcs", use_real_bcs_temp);
259 int n_zlevels =
pp.countval(
"terrain_z_levels");
277 pp.query(
"c_p",
c_p);
287 for (
int i = 0; i <= max_level; i++) {
291 int nvals =
pp.countval(
"substepping_type");
292 AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
296 for (
int i = 1; i <= max_level; i++) {
299 }
else if (nvals > 1) {
300 for (
int i = 0; i <= max_level; i++) {
308 int nvals_old =
pp.countval(
"no_substepping");
310 amrex::Abort(
"The no_substepping flag is deprecated -- set substepping_type instead");
314 bool any_anelastic =
false;
315 for (
int i = 0; i <= max_level; ++i) {
316 if (
anelastic[i] == 1) any_anelastic =
true;
319 if (any_anelastic == 1) {
335 for (
int lev = 0; lev <= max_level; lev++) {
365 AMREX_ASSERT_WITH_MESSAGE(
terrain_type != TerrainType::None,
"MOST stress rotations are only valid with terrain!");
379 pp.query_enum_case_insensitive(
"perturbation_type",
pert_type);
381 amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
382 pp.queryarr(
"abl_pressure_grad",abl_pressure_grad_in);
383 for(
int i = 0; i < AMREX_SPACEDIM; ++i)
abl_pressure_grad[i] = abl_pressure_grad_in[i];
385 amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
386 if(
pp.queryarr(
"abl_geo_forcing",abl_geo_forcing_in)) {
387 amrex::Print() <<
"Specified abl_geo_forcing: (";
388 for (
int i = 0; i < AMREX_SPACEDIM; ++i) {
392 amrex::Print() <<
")" << std::endl;
410 "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
419 "Numerical diffusion coefficient must be between 0 & 1.");
422 amrex::Print() <<
"6th-order numerical diffusion turned on with coefficient = "
435 for (
int lev = 0; lev <= max_level; lev++) {
436 turbChoice[lev].init_params(lev,max_level,pp_prefix);
440 for (
int lev = 0; lev <= max_level; lev++) {
441 if (
turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
445 amrex::Print() <<
"YSU PBL using ERF coriolis frequency: " <<
turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
448 amrex::Abort(
"YSU cannot use ERF coriolis frequency if not using coriolis");
453 for (
int lev = 0; lev <= max_level; lev++) {
454 if (
turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
458 amrex::Print() <<
"MRF PBL using ERF coriolis frequency: " <<
turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
461 amrex::Abort(
"MRF cannot use ERF coriolis frequency if not using coriolis");
474 static std::string windfarm_loc_type_string =
"None";
constexpr amrex::Real CONST_GRAV
Definition: ERF_Constants.H:21
constexpr amrex::Real R_d
Definition: ERF_Constants.H:10
#define RhoQ4_comp
Definition: ERF_IndexDefines.H:45
#define RhoQ2_comp
Definition: ERF_IndexDefines.H:43
#define RhoQ3_comp
Definition: ERF_IndexDefines.H:44
#define RhoQ1_comp
Definition: ERF_IndexDefines.H:42
@ bogus
Definition: ERF_IndexDefines.H:204
void init_params(std::string pp_prefix)
Definition: ERF_AdvStruct.H:21
void init_params(std::string pp_prefix)
Definition: ERF_DiffStruct.H:21
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:743
amrex::Real dz0
Definition: ERF_DataStruct.H:764
amrex::Real const_massflux_layer_lo
Definition: ERF_DataStruct.H:846
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:741
int RhoQr_comp
Definition: ERF_DataStruct.H:830
bool rebalance_wrfinput
Definition: ERF_DataStruct.H:801
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:729
bool rayleigh_damp_substep
Definition: ERF_DataStruct.H:747
void build_coriolis_forcings_const_lat(std::string pp_prefix)
Definition: ERF_DataStruct.H:639
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:745
amrex::Real rdOcp
Definition: ERF_DataStruct.H:759
void read_int_string(int max_level, const char *string_to_read, amrex::Vector< int > &vec_to_fill, int default_int)
Definition: ERF_DataStruct.H:681
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:832
bool use_mono_adv
Definition: ERF_DataStruct.H:798
amrex::Real const_massflux_v
Definition: ERF_DataStruct.H:844
int ncorr
Definition: ERF_DataStruct.H:727
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:832
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:807
amrex::Real c_p
Definition: ERF_DataStruct.H:758
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:832
amrex::Real gravity
Definition: ERF_DataStruct.H:757
void check_params(int max_level)
Definition: ERF_DataStruct.H:499
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:774
bool custom_w_subsidence
Definition: ERF_DataStruct.H:776
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:782
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:740
int RhoQc_comp
Definition: ERF_DataStruct.H:824
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:777
bool immersed_forcing_substep
Definition: ERF_DataStruct.H:750
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:762
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:814
amrex::Real const_massflux_u
Definition: ERF_DataStruct.H:843
bool use_num_diff
Definition: ERF_DataStruct.H:794
int RhoQv_comp
Definition: ERF_DataStruct.H:823
bool test_mapfactor
Definition: ERF_DataStruct.H:731
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:775
amrex::Real num_diff_coeff
Definition: ERF_DataStruct.H:795
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:833
amrex::Real zsurf
Definition: ERF_DataStruct.H:763
bool variable_coriolis
Definition: ERF_DataStruct.H:817
bool project_initial_velocity
Definition: ERF_DataStruct.H:766
bool use_moist_background
Definition: ERF_DataStruct.H:822
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:778
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:742
PerturbationType pert_type
Definition: ERF_DataStruct.H:791
static InitType init_type
Definition: ERF_DataStruct.H:699
amrex::Real const_massflux_layer_hi
Definition: ERF_DataStruct.H:847
static bool use_real_bcs
Definition: ERF_DataStruct.H:705
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:728
bool has_lat_lon
Definition: ERF_DataStruct.H:816
bool use_rotate_surface_flux
Definition: ERF_DataStruct.H:785
amrex::Real const_massflux_tau
Definition: ERF_DataStruct.H:845
bool time_avg_vel
Definition: ERF_DataStruct.H:788
bool forest_substep
Definition: ERF_DataStruct.H:751
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:744
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:833
int ave_plane
Definition: ERF_DataStruct.H:819
void init_params(std::string pp_prefix)
Definition: ERF_SpongeStruct.H:17