127 amrex::ParmParse
pp(pp_prefix);
130 if (
pp.query(
"use_terrain",
bogus) > 0) {
131 amrex::Abort(
"The input use_terrain is deprecated. Set terrain_type instead.");
137 "The grid stretching ratio must be greater than 1");
144 amrex::Print() <<
"Turning terrain on to enable grid stretching" << std::endl;
149 amrex::Print() <<
"Nominal zsurface height != 0, may result in unexpected behavior"
152 pp.get(
"initial_dz",
dz0);
161 amrex::Abort(
"The input moisture_type is deprecated. Set moisture_model instead.");
174 }
else if ( (
moisture_type == MoistureType::Morrison_NoIce) ||
182 }
else if ( (
moisture_type == MoistureType::SAM_NoPrecip_NoIce) ||
222 pp.query_enum_case_insensitive(
"land_surface_model",
lsm_type);
226 pp.query_enum_case_insensitive(
"radiation_model",
rad_type);
229 #ifndef ERF_USE_RRTMGP
230 if (
rad_type == RadiationType::RRTMGP)
232 amrex::Abort(
"ERF was not compiled with RRTMGP enabled!");
237 std::string terrain_type_temp =
"";
238 pp.query(
"terrain_type", terrain_type_temp);
239 if (terrain_type_temp ==
"Moving") {
240 amrex::Warning(
"erf.terrain_type = Moving is deprecated; please replace Moving by MovingFittedMesh");
242 }
else if (terrain_type_temp ==
"Static") {
243 amrex::Warning(
"erf.terrain_type = Static is deprecated; please replace Static by StaticFittedMesh");
252 std::string init_type_temp_string;
253 pp.query(
"init_type",init_type_temp_string);
254 if ( (init_type_temp_string ==
"Real") || (init_type_temp_string ==
"real") ) {
255 amrex::Warning(
"erf.init_type = Real is deprecated; please replace Real by WRFInput");
259 }
else if ( (init_type_temp_string ==
"Ideal") || (init_type_temp_string ==
"ideal") ) {
260 amrex::Warning(
"erf.init_type = Ideal is deprecated; please replace Ideal by WRFInput");
263 }
else if (init_type_temp_string ==
"hindcast") {
268 pp.query_enum_case_insensitive(
"init_type",
init_type);
277 amrex::Abort(
"Only terrain_type = StaticFittedMesh are allowed with init_type = WRFInput or Metgrid");
285 amrex::Abort(
"Can't have moisture_type = None with init_type = WRFInput");
289 std::string nc_bdy_file_temp_string;
290 bool has_bdy =
pp.query(
"nc_bdy_file", nc_bdy_file_temp_string);
294 pp.query(
"use_real_bcs", use_real_bcs_temp);
309 if (
init_type == InitType::Input_Sounding) {
318 int n_zlevels =
pp.countval(
"terrain_z_levels");
336 pp.query(
"c_p",
c_p);
346 for (
int i = 0; i <= max_level; i++) {
350 int nvals =
pp.countval(
"substepping_type");
351 AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
355 for (
int i = 1; i <= max_level; i++) {
358 }
else if (nvals > 1) {
359 for (
int i = 0; i <= max_level; i++) {
367 int nvals_old =
pp.countval(
"no_substepping");
369 amrex::Abort(
"The no_substepping flag is deprecated -- set substepping_type instead");
373 bool any_anelastic =
false;
374 for (
int i = 0; i <= max_level; ++i) {
375 if (
anelastic[i] == 1) any_anelastic =
true;
378 if (any_anelastic == 1) {
394 for (
int lev = 0; lev <= max_level; lev++) {
424 AMREX_ASSERT_WITH_MESSAGE(
terrain_type != TerrainType::None,
"MOST stress rotations are only valid with terrain!");
438 pp.query_enum_case_insensitive(
"perturbation_type",
pert_type);
440 amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
441 pp.queryarr(
"abl_pressure_grad",abl_pressure_grad_in);
442 for(
int i = 0; i < AMREX_SPACEDIM; ++i)
abl_pressure_grad[i] = abl_pressure_grad_in[i];
444 amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
445 if(
pp.queryarr(
"abl_geo_forcing",abl_geo_forcing_in)) {
446 amrex::Print() <<
"Specified abl_geo_forcing: (";
447 for (
int i = 0; i < AMREX_SPACEDIM; ++i) {
451 amrex::Print() <<
")" << std::endl;
469 "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
478 "Numerical diffusion coefficient must be between 0 & 1.");
481 amrex::Print() <<
"6th-order numerical diffusion turned on with coefficient = "
491 for (
int lev = 0; lev <= max_level; lev++) {
492 turbChoice[lev].init_params(lev,max_level,pp_prefix);
496 for (
int lev = 0; lev <= max_level; lev++) {
497 if (
turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
501 amrex::Print() <<
"YSU PBL using ERF coriolis frequency: " <<
turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
504 amrex::Abort(
"YSU cannot use ERF coriolis frequency if not using coriolis");
509 for (
int lev = 0; lev <= max_level; lev++) {
510 if (
turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
514 amrex::Print() <<
"MRF PBL using ERF coriolis frequency: " <<
turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
517 amrex::Abort(
"MRF cannot use ERF coriolis frequency if not using coriolis");
526 amrex::Abort(
"You set use_shoc to true but didn't build with SHOC; you must rebuild the executable");
538 static std::string windfarm_loc_type_string =
"None";
565 amrex::Abort(
"ERROR: Missing input parameter 'erf.hindcast_boundary_data_dir' for boundary data for lateral forcing");
577 amrex::Abort(
"ERROR: Missing input parameter 'erf.hindcast_lateral_sponge_strength' or it is specified to be less than zero");
581 amrex::Abort(
"ERROR: Missing input parameter 'erf.hindcast_lateral_sponge_length' or it is specified to be less than zero");
585 amrex::Abort(
"ERROR: Missing input parameter 'erf.hindcast_zhi_sponge_strength' or it is specified to be less than zero");
589 amrex::Abort(
"ERROR: Missing input parameter 'erf.hindcast_zhi_sponge_strength' or it is specified to be less than zero");
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
#define RhoQ6_comp
Definition: ERF_IndexDefines.H:47
#define RhoQ5_comp
Definition: ERF_IndexDefines.H:46
@ 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
Definition: ERF_DataStruct.H:99
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:850
amrex::Real dz0
Definition: ERF_DataStruct.H:871
amrex::Real const_massflux_layer_lo
Definition: ERF_DataStruct.H:948
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:848
bool rebalance_wrfinput
Definition: ERF_DataStruct.H:908
amrex::Real hindcast_lateral_sponge_strength
Definition: ERF_DataStruct.H:955
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:836
bool rayleigh_damp_substep
Definition: ERF_DataStruct.H:854
void build_coriolis_forcings_const_lat(std::string pp_prefix)
Definition: ERF_DataStruct.H:740
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:852
amrex::Real rdOcp
Definition: ERF_DataStruct.H:866
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:782
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:934
amrex::Real hindcast_zhi_sponge_length
Definition: ERF_DataStruct.H:956
amrex::Real const_massflux_v
Definition: ERF_DataStruct.H:946
int ncorr
Definition: ERF_DataStruct.H:834
bool hindcast_zhi_sponge_damping
Definition: ERF_DataStruct.H:957
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:934
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:914
amrex::Real c_p
Definition: ERF_DataStruct.H:865
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:934
amrex::Real gravity
Definition: ERF_DataStruct.H:864
void check_params(int max_level)
Definition: ERF_DataStruct.H:596
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:881
amrex::Real hindcast_lateral_sponge_length
Definition: ERF_DataStruct.H:955
bool use_shoc
Definition: ERF_DataStruct.H:895
bool hindcast_lateral_forcing
Definition: ERF_DataStruct.H:954
bool moisture_tight_coupling
Definition: ERF_DataStruct.H:932
bool custom_w_subsidence
Definition: ERF_DataStruct.H:883
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:889
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:847
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:884
bool immersed_forcing_substep
Definition: ERF_DataStruct.H:857
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:869
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:921
amrex::Real const_massflux_u
Definition: ERF_DataStruct.H:945
bool use_num_diff
Definition: ERF_DataStruct.H:904
amrex::Real hindcast_zhi_sponge_strength
Definition: ERF_DataStruct.H:956
bool test_mapfactor
Definition: ERF_DataStruct.H:838
static SoundingType sounding_type
Definition: ERF_DataStruct.H:803
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:882
amrex::Real num_diff_coeff
Definition: ERF_DataStruct.H:905
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:935
amrex::Real zsurf
Definition: ERF_DataStruct.H:870
bool variable_coriolis
Definition: ERF_DataStruct.H:924
bool project_initial_velocity
Definition: ERF_DataStruct.H:873
static bool upwind_real_bcs
Definition: ERF_DataStruct.H:812
bool use_moist_background
Definition: ERF_DataStruct.H:929
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:885
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:849
PerturbationType pert_type
Definition: ERF_DataStruct.H:901
static InitType init_type
Definition: ERF_DataStruct.H:800
amrex::Real const_massflux_layer_hi
Definition: ERF_DataStruct.H:949
static bool use_real_bcs
Definition: ERF_DataStruct.H:809
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:835
MoistureComponentIndices moisture_indices
Definition: ERF_DataStruct.H:930
bool has_lat_lon
Definition: ERF_DataStruct.H:923
bool use_rotate_surface_flux
Definition: ERF_DataStruct.H:892
amrex::Real const_massflux_tau
Definition: ERF_DataStruct.H:947
bool time_avg_vel
Definition: ERF_DataStruct.H:898
bool forest_substep
Definition: ERF_DataStruct.H:858
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:851
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:935
int ave_plane
Definition: ERF_DataStruct.H:926
std::string hindcast_boundary_data_dir
Definition: ERF_DataStruct.H:953
void init_params(std::string pp_prefix)
Definition: ERF_SpongeStruct.H:17