From f22f0e30494acfb66d2993276313f9f24e40073e Mon Sep 17 00:00:00 2001 From: asalmgren Date: Wed, 27 Nov 2024 02:38:09 +0000 Subject: [PATCH] Deployed from erf-model/ERF --- ERF__ComputeTurbulentViscosity_8cpp.html | 8 +- ERF__DataStruct_8H_source.html | 1137 +++++++++++----------- ERF__DiffusionSrcForState__N_8cpp.html | 8 +- ERF__DiffusionSrcForState__T_8cpp.html | 8 +- ERF__Kessler_8H_source.html | 10 +- ERF__SAM_8H_source.html | 12 +- ERF__TI__substep__fun_8H_source.html | 2 +- ERF__make__buoyancy_8cpp.html | 8 +- ERF__make__mom__sources_8cpp.html | 42 +- ERF__make__sources_8cpp.html | 34 +- ERF__make__tau__terms_8cpp.html | 10 +- ERF__prob__common_8H_source.html | 2 +- ERF__slow__rhs__post_8cpp.html | 22 +- ERF__slow__rhs__pre_8cpp.html | 26 +- classERF.html | 94 +- classKessler.html | 10 +- classProblemBase.html | 2 +- classSAM.html | 12 +- structSolverChoice.html | 867 +++++++++-------- 19 files changed, 1162 insertions(+), 1152 deletions(-) diff --git a/ERF__ComputeTurbulentViscosity_8cpp.html b/ERF__ComputeTurbulentViscosity_8cpp.html index da84f1e9a..b9678b334 100644 --- a/ERF__ComputeTurbulentViscosity_8cpp.html +++ b/ERF__ComputeTurbulentViscosity_8cpp.html @@ -390,10 +390,10 @@

ComputeTurbulentViscosityLES
void ComputeTurbulentViscosityLES(const MultiFab &Tau11, const MultiFab &Tau22, const MultiFab &Tau33, const MultiFab &Tau12, const MultiFab &Tau13, const MultiFab &Tau23, const MultiFab &cons_in, MultiFab &eddyViscosity, MultiFab &Hfx1, MultiFab &Hfx2, MultiFab &Hfx3, MultiFab &Diss, const Geometry &geom, const MultiFab &mapfac_u, const MultiFab &mapfac_v, const std::unique_ptr< MultiFab > &z_phys_nd, const TurbChoice &turbChoice, const Real const_grav, std::unique_ptr< ABLMost > &most, const bool &exp_most)
Definition: ERF_ComputeTurbulentViscosity.cpp:32
Vars::xvel
@ xvel
Definition: ERF_IndexDefines.H:130
Vars::yvel
@ yvel
Definition: ERF_IndexDefines.H:131
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
TurbChoice::les_type
LESType les_type
Definition: ERF_TurbStruct.H:169
diff --git a/ERF__DataStruct_8H_source.html b/ERF__DataStruct_8H_source.html index 4a4e67c3b..4d36f7c3d 100644 --- a/ERF__DataStruct_8H_source.html +++ b/ERF__DataStruct_8H_source.html @@ -231,513 +231,518 @@
143 
144  // Is the terrain none, static or moving?
145  pp.query_enum_case_insensitive("terrain_type",terrain_type);
-
146 
-
147  // Use lagged_delta_rt in the fast integrator?
-
148  pp.query("use_lagged_delta_rt", use_lagged_delta_rt);
-
149 
-
150  // These default to true but are used for unit testing
-
151  pp.query("use_gravity", use_gravity);
-
152  gravity = use_gravity? CONST_GRAV: 0.0;
-
153 
-
154  pp.query("c_p", c_p);
-
155  rdOcp = R_d / c_p;
-
156 
-
157  read_int_string(max_level, "anelastic", anelastic, 0);
+
146  int n_zlevels = pp.countval("terrain_z_levels");
+
147  if (n_zlevels > 0 and terrain_type == TerrainType::None)
+
148  {
+
149  terrain_type = TerrainType::Static;
+
150  }
+
151 
+
152  // Use lagged_delta_rt in the fast integrator?
+
153  pp.query("use_lagged_delta_rt", use_lagged_delta_rt);
+
154 
+
155  // These default to true but are used for unit testing
+
156  pp.query("use_gravity", use_gravity);
+
157  gravity = use_gravity? CONST_GRAV: 0.0;
158 
-
159  // *******************************************************************************
-
160  // Read substepping_type and allow for different values at each level
-
161  // *******************************************************************************
-
162  substepping_type.resize(max_level+1);
+
159  pp.query("c_p", c_p);
+
160  rdOcp = R_d / c_p;
+
161 
+
162  read_int_string(max_level, "anelastic", anelastic, 0);
163 
-
164  for (int i = 0; i <= max_level; i++) {
-
165  substepping_type[i] = SubsteppingType::Implicit;
-
166  }
-
167 
-
168  int nvals = pp.countval("substepping_type");
-
169  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
-
170 
-
171  if (nvals == 1) {
-
172  pp.query_enum_case_insensitive("substepping_type",substepping_type[0]);
-
173  for (int i = 1; i <= max_level; i++) {
-
174  substepping_type[i] = substepping_type[0];
-
175  }
-
176  } else if (nvals > 1) { // in this case we have asserted nvals >= max_level+1
-
177  for (int i = 0; i <= max_level; i++) {
-
178  pp.query_enum_case_insensitive("substepping_type",substepping_type[i],i);
-
179  }
-
180  }
-
181 
-
182  // *******************************************************************************
-
183  // Error check on deprecated input
-
184  // *******************************************************************************
-
185  int nvals_old = pp.countval("no_substepping");
-
186  if (nvals_old > 0) {
-
187  amrex::Abort("The no_substepping flag is deprecated -- set substepping_type instead");
-
188  }
+
164  // *******************************************************************************
+
165  // Read substepping_type and allow for different values at each level
+
166  // *******************************************************************************
+
167  substepping_type.resize(max_level+1);
+
168 
+
169  for (int i = 0; i <= max_level; i++) {
+
170  substepping_type[i] = SubsteppingType::Implicit;
+
171  }
+
172 
+
173  int nvals = pp.countval("substepping_type");
+
174  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
+
175 
+
176  if (nvals == 1) {
+
177  pp.query_enum_case_insensitive("substepping_type",substepping_type[0]);
+
178  for (int i = 1; i <= max_level; i++) {
+
179  substepping_type[i] = substepping_type[0];
+
180  }
+
181  } else if (nvals > 1) { // in this case we have asserted nvals >= max_level+1
+
182  for (int i = 0; i <= max_level; i++) {
+
183  pp.query_enum_case_insensitive("substepping_type",substepping_type[i],i);
+
184  }
+
185  }
+
186 
+
187  // *******************************************************************************
+
188  // Error check on deprecated input
189  // *******************************************************************************
-
190 
-
191  bool any_anelastic = false;
-
192  for (int i = 0; i <= max_level; ++i) {
-
193  if (anelastic[i] == 1) any_anelastic = true;
-
194  }
+
190  int nvals_old = pp.countval("no_substepping");
+
191  if (nvals_old > 0) {
+
192  amrex::Abort("The no_substepping flag is deprecated -- set substepping_type instead");
+
193  }
+
194  // *******************************************************************************
195 
-
196  // If anelastic at all, we do not advect rho -- it is always == rho0
-
197  if (any_anelastic == 1) {
-
198  constant_density = true;
-
199  project_initial_velocity = true;
-
200  buoyancy_type = 3; // (This isn't actually used when anelastic is set)
-
201  } else {
-
202  pp.query("project_initial_velocity", project_initial_velocity);
-
203 
-
204  constant_density = false; // We default to false but allow the user to set it
-
205  pp.query("constant_density", constant_density);
-
206  }
-
207 
-
208  // *******************************************************************************
-
209 
-
210  pp.query("ncorr", ncorr);
-
211  pp.query("poisson_abstol", poisson_abstol);
-
212  pp.query("poisson_reltol", poisson_reltol);
-
213 
-
214  for (int lev = 0; lev <= max_level; lev++) {
-
215  if (anelastic[lev] != 0)
-
216  {
-
217  substepping_type[lev] = SubsteppingType::None;
-
218  }
-
219  }
-
220 
-
221  pp.query("force_stage1_single_substep", force_stage1_single_substep);
-
222 
-
223  // Include Coriolis forcing?
-
224  pp.query("use_coriolis", use_coriolis);
+
196  bool any_anelastic = false;
+
197  for (int i = 0; i <= max_level; ++i) {
+
198  if (anelastic[i] == 1) any_anelastic = true;
+
199  }
+
200 
+
201  // If anelastic at all, we do not advect rho -- it is always == rho0
+
202  if (any_anelastic == 1) {
+
203  constant_density = true;
+
204  project_initial_velocity = true;
+
205  buoyancy_type = 3; // (This isn't actually used when anelastic is set)
+
206  } else {
+
207  pp.query("project_initial_velocity", project_initial_velocity);
+
208 
+
209  constant_density = false; // We default to false but allow the user to set it
+
210  pp.query("constant_density", constant_density);
+
211  }
+
212 
+
213  // *******************************************************************************
+
214 
+
215  pp.query("ncorr", ncorr);
+
216  pp.query("poisson_abstol", poisson_abstol);
+
217  pp.query("poisson_reltol", poisson_reltol);
+
218 
+
219  for (int lev = 0; lev <= max_level; lev++) {
+
220  if (anelastic[lev] != 0)
+
221  {
+
222  substepping_type[lev] = SubsteppingType::None;
+
223  }
+
224  }
225 
-
226  // Include Rayleigh damping (separate flags for each variable)
-
227  pp.query("rayleigh_damp_U", rayleigh_damp_U);
-
228  pp.query("rayleigh_damp_V", rayleigh_damp_V);
-
229  pp.query("rayleigh_damp_W", rayleigh_damp_W);
-
230  pp.query("rayleigh_damp_T", rayleigh_damp_T);
-
231  pp.query("rayleigh_dampcoef", rayleigh_dampcoef);
-
232  pp.query("rayleigh_zdamp", rayleigh_zdamp);
-
233 
-
234  // Flag to do explicit MOST formulation
-
235  pp.query("use_explicit_most",use_explicit_most);
-
236 
-
237  // Flag to do MOST rotations with terrain
-
238  pp.query("use_rotate_most",use_rotate_most);
-
239  if (use_rotate_most) {
-
240  AMREX_ASSERT_WITH_MESSAGE(terrain_type != TerrainType::None,"MOST stress rotations are only valid with terrain!");
-
241  AMREX_ASSERT_WITH_MESSAGE(use_explicit_most, "MOST Stress rotations are only valid with explicit MOST!");
-
242  }
-
243 
-
244  // Which external forcings?
-
245  abl_driver_type = ABLDriverType::None; // Default: no ABL driver for simulating classical fluid dynamics problems
-
246  pp.query_enum_case_insensitive("abl_driver_type",abl_driver_type);
-
247 
-
248  // Which type of inflow turbulent generation
-
249  pert_type = PerturbationType::None; // Default
-
250  pp.query_enum_case_insensitive("perturbation_type",pert_type);
-
251 
-
252  amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
-
253  pp.queryarr("abl_pressure_grad",abl_pressure_grad_in);
-
254  for(int i = 0; i < AMREX_SPACEDIM; ++i) abl_pressure_grad[i] = abl_pressure_grad_in[i];
-
255 
-
256  amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
-
257  if(pp.queryarr("abl_geo_forcing",abl_geo_forcing_in)) {
-
258  amrex::Print() << "Specified abl_geo_forcing: (";
-
259  for (int i = 0; i < AMREX_SPACEDIM; ++i) {
-
260  abl_geo_forcing[i] = abl_geo_forcing_in[i];
-
261  amrex::Print() << abl_geo_forcing[i] << " ";
-
262  }
-
263  amrex::Print() << ")" << std::endl;
-
264  }
-
265 
-
266  if (use_coriolis)
-
267  {
-
268  build_coriolis_forcings();
+
226  pp.query("force_stage1_single_substep", force_stage1_single_substep);
+
227 
+
228  // Include Coriolis forcing?
+
229  pp.query("use_coriolis", use_coriolis);
+
230 
+
231  // Include Rayleigh damping (separate flags for each variable)
+
232  pp.query("rayleigh_damp_U", rayleigh_damp_U);
+
233  pp.query("rayleigh_damp_V", rayleigh_damp_V);
+
234  pp.query("rayleigh_damp_W", rayleigh_damp_W);
+
235  pp.query("rayleigh_damp_T", rayleigh_damp_T);
+
236  pp.query("rayleigh_dampcoef", rayleigh_dampcoef);
+
237  pp.query("rayleigh_zdamp", rayleigh_zdamp);
+
238 
+
239  // Flag to do explicit MOST formulation
+
240  pp.query("use_explicit_most",use_explicit_most);
+
241 
+
242  // Flag to do MOST rotations with terrain
+
243  pp.query("use_rotate_most",use_rotate_most);
+
244  if (use_rotate_most) {
+
245  AMREX_ASSERT_WITH_MESSAGE(terrain_type != TerrainType::None,"MOST stress rotations are only valid with terrain!");
+
246  AMREX_ASSERT_WITH_MESSAGE(use_explicit_most, "MOST Stress rotations are only valid with explicit MOST!");
+
247  }
+
248 
+
249  // Which external forcings?
+
250  abl_driver_type = ABLDriverType::None; // Default: no ABL driver for simulating classical fluid dynamics problems
+
251  pp.query_enum_case_insensitive("abl_driver_type",abl_driver_type);
+
252 
+
253  // Which type of inflow turbulent generation
+
254  pert_type = PerturbationType::None; // Default
+
255  pp.query_enum_case_insensitive("perturbation_type",pert_type);
+
256 
+
257  amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
+
258  pp.queryarr("abl_pressure_grad",abl_pressure_grad_in);
+
259  for(int i = 0; i < AMREX_SPACEDIM; ++i) abl_pressure_grad[i] = abl_pressure_grad_in[i];
+
260 
+
261  amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
+
262  if(pp.queryarr("abl_geo_forcing",abl_geo_forcing_in)) {
+
263  amrex::Print() << "Specified abl_geo_forcing: (";
+
264  for (int i = 0; i < AMREX_SPACEDIM; ++i) {
+
265  abl_geo_forcing[i] = abl_geo_forcing_in[i];
+
266  amrex::Print() << abl_geo_forcing[i] << " ";
+
267  }
+
268  amrex::Print() << ")" << std::endl;
269  }
270 
-
271  pp.query("add_custom_rhotheta_forcing", custom_rhotheta_forcing);
-
272  pp.query("add_custom_moisture_forcing", custom_moisture_forcing);
-
273  pp.query("add_custom_w_subsidence", custom_w_subsidence);
-
274  pp.query("add_custom_geostrophic_profile", custom_geostrophic_profile);
-
275  pp.query("custom_forcing_uses_primitive_vars", custom_forcing_prim_vars);
-
276 
-
277  pp.query("nudging_from_input_sounding", nudging_from_input_sounding);
-
278 
-
279  have_geo_wind_profile = (!abl_geo_wind_table.empty() || custom_geostrophic_profile);
-
280  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(!(!abl_geo_wind_table.empty() && custom_geostrophic_profile),
-
281  "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
-
282 
-
283  pp.query("Ave_Plane", ave_plane);
-
284 
-
285  pp.query("mp_clouds", do_cloud);
-
286  pp.query("mp_precip", do_precip);
-
287  pp.query("use_moist_background", use_moist_background);
-
288 
-
289  // Use numerical diffusion?
-
290  pp.query("use_NumDiff",use_NumDiff);
-
291  if(use_NumDiff) {
-
292  pp.query("NumDiffCoeff",NumDiffCoeff);
-
293  AMREX_ASSERT_WITH_MESSAGE(( (NumDiffCoeff >= 0.) && (NumDiffCoeff <= 1.) ),
-
294  "Numerical diffusion coefficient must be between 0 & 1.");
-
295  NumDiffCoeff *= std::pow(2.0,-6);
-
296  }
-
297 
-
298  // Use monotonic advection?
-
299  pp.query("use_mono_adv",use_mono_adv);
-
300 
-
301  advChoice.init_params();
-
302  diffChoice.init_params();
-
303  spongeChoice.init_params();
-
304 
-
305  turbChoice.resize(max_level+1);
-
306  for (int lev = 0; lev <= max_level; lev++) {
-
307  turbChoice[lev].init_params(lev,max_level);
-
308  }
+
271  if (use_coriolis)
+
272  {
+
273  build_coriolis_forcings();
+
274  }
+
275 
+
276  pp.query("add_custom_rhotheta_forcing", custom_rhotheta_forcing);
+
277  pp.query("add_custom_moisture_forcing", custom_moisture_forcing);
+
278  pp.query("add_custom_w_subsidence", custom_w_subsidence);
+
279  pp.query("add_custom_geostrophic_profile", custom_geostrophic_profile);
+
280  pp.query("custom_forcing_uses_primitive_vars", custom_forcing_prim_vars);
+
281 
+
282  pp.query("nudging_from_input_sounding", nudging_from_input_sounding);
+
283 
+
284  have_geo_wind_profile = (!abl_geo_wind_table.empty() || custom_geostrophic_profile);
+
285  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(!(!abl_geo_wind_table.empty() && custom_geostrophic_profile),
+
286  "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
+
287 
+
288  pp.query("Ave_Plane", ave_plane);
+
289 
+
290  pp.query("mp_clouds", do_cloud);
+
291  pp.query("mp_precip", do_precip);
+
292  pp.query("use_moist_background", use_moist_background);
+
293 
+
294  // Use numerical diffusion?
+
295  pp.query("use_NumDiff",use_NumDiff);
+
296  if(use_NumDiff) {
+
297  pp.query("NumDiffCoeff",NumDiffCoeff);
+
298  AMREX_ASSERT_WITH_MESSAGE(( (NumDiffCoeff >= 0.) && (NumDiffCoeff <= 1.) ),
+
299  "Numerical diffusion coefficient must be between 0 & 1.");
+
300  NumDiffCoeff *= std::pow(2.0,-6);
+
301  }
+
302 
+
303  // Use monotonic advection?
+
304  pp.query("use_mono_adv",use_mono_adv);
+
305 
+
306  advChoice.init_params();
+
307  diffChoice.init_params();
+
308  spongeChoice.init_params();
309 
-
310  // YSU PBL: use consistent coriolis frequency
+
310  turbChoice.resize(max_level+1);
311  for (int lev = 0; lev <= max_level; lev++) {
-
312  if (turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
-
313  if (use_coriolis) {
-
314  turbChoice[lev].pbl_ysu_coriolis_freq = coriolis_factor * sinphi;
-
315  if (lev == 0) {
-
316  amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
-
317  }
-
318  } else {
-
319  amrex::Abort("YSU cannot use ERF coriolis frequency if not using coriolis");
-
320  }
-
321  }
-
322  }
-
323 
-
324 
-
325  // Which type of multilevel coupling
-
326  coupling_type = CouplingType::TwoWay; // Default
-
327  pp.query_enum_case_insensitive("coupling_type",coupling_type);
+
312  turbChoice[lev].init_params(lev,max_level);
+
313  }
+
314 
+
315  // YSU PBL: use consistent coriolis frequency
+
316  for (int lev = 0; lev <= max_level; lev++) {
+
317  if (turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
+
318  if (use_coriolis) {
+
319  turbChoice[lev].pbl_ysu_coriolis_freq = coriolis_factor * sinphi;
+
320  if (lev == 0) {
+
321  amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
+
322  }
+
323  } else {
+
324  amrex::Abort("YSU cannot use ERF coriolis frequency if not using coriolis");
+
325  }
+
326  }
+
327  }
328 
-
329  // Which type of windfarm model
-
330  windfarm_type = WindFarmType::None; // Default
-
331  pp.query_enum_case_insensitive("windfarm_type",windfarm_type);
-
332 
-
333  static std::string windfarm_loc_type_string = "None";
-
334  windfarm_loc_type = WindFarmLocType::None;
-
335  pp.query_enum_case_insensitive("windfarm_loc_type",windfarm_loc_type);
-
336 
-
337  pp.query("windfarm_loc_table", windfarm_loc_table);
-
338  pp.query("windfarm_spec_table", windfarm_spec_table);
-
339  pp.query("windfarm_blade_table", windfarm_blade_table);
-
340  pp.query("windfarm_airfoil_tables", windfarm_airfoil_tables);
-
341  pp.query("windfarm_spec_table_extra", windfarm_spec_table_extra);
-
342 
-
343  // Sampling distance upstream of the turbine to find the
-
344  // incoming free stream velocity as a factor of the diameter of the
-
345  // turbine. ie. the sampling distance will be this number multiplied
-
346  // by the diameter of the turbine
-
347  pp.query("sampling_distance_by_D", sampling_distance_by_D);
-
348  pp.query("turb_disk_angle_from_x", turb_disk_angle);
-
349 
-
350  pp.query("windfarm_x_shift",windfarm_x_shift);
-
351  pp.query("windfarm_y_shift",windfarm_y_shift);
-
352  // Test if time averaged data is to be output
-
353  pp.query("time_avg_vel",time_avg_vel);
+
329 
+
330  // Which type of multilevel coupling
+
331  coupling_type = CouplingType::TwoWay; // Default
+
332  pp.query_enum_case_insensitive("coupling_type",coupling_type);
+
333 
+
334  // Which type of windfarm model
+
335  windfarm_type = WindFarmType::None; // Default
+
336  pp.query_enum_case_insensitive("windfarm_type",windfarm_type);
+
337 
+
338  static std::string windfarm_loc_type_string = "None";
+
339  windfarm_loc_type = WindFarmLocType::None;
+
340  pp.query_enum_case_insensitive("windfarm_loc_type",windfarm_loc_type);
+
341 
+
342  pp.query("windfarm_loc_table", windfarm_loc_table);
+
343  pp.query("windfarm_spec_table", windfarm_spec_table);
+
344  pp.query("windfarm_blade_table", windfarm_blade_table);
+
345  pp.query("windfarm_airfoil_tables", windfarm_airfoil_tables);
+
346  pp.query("windfarm_spec_table_extra", windfarm_spec_table_extra);
+
347 
+
348  // Sampling distance upstream of the turbine to find the
+
349  // incoming free stream velocity as a factor of the diameter of the
+
350  // turbine. ie. the sampling distance will be this number multiplied
+
351  // by the diameter of the turbine
+
352  pp.query("sampling_distance_by_D", sampling_distance_by_D);
+
353  pp.query("turb_disk_angle_from_x", turb_disk_angle);
354 
-
355  check_params(max_level);
-
356  }
-
357 
-
358  void check_params(int max_level)
-
359  {
-
360  // Warn for PBL models and moisture - these may not yet be compatible
-
361  for (int lev = 0; lev <= max_level; lev++) {
-
362  if ((moisture_type != MoistureType::None) && (turbChoice[lev].pbl_type != PBLType::None)) {
-
363  amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
-
364  }
-
365  }
-
366  //
-
367  // Buoyancy type check
-
368  //
-
369  if (buoyancy_type != 1 && buoyancy_type != 2 && buoyancy_type != 3 && buoyancy_type != 4) {
-
370  amrex::Abort("buoyancy_type must be 1, 2, 3 or 4");
-
371  }
-
372 
-
373  if (!use_lagged_delta_rt && !(terrain_type == TerrainType::Moving)) {
-
374  amrex::Error("Can't turn off lagged_delta_rt when terrain not moving");
-
375  }
-
376 
-
377  //
-
378  // Wind farm checks
-
379  //
-
380  if (windfarm_type==WindFarmType::SimpleAD and sampling_distance_by_D < 0.0) {
-
381  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
-
382  " erf.sampling_distance_by_D in the inputs which specifies the upstream"
-
383  " distance as a factor of the turbine diameter at which the incoming free stream"
-
384  " velocity will be computed at.");
-
385  }
-
386  if ( (windfarm_type==WindFarmType::SimpleAD ||
-
387  windfarm_type==WindFarmType::GeneralAD ) && turb_disk_angle < 0.0) {
-
388  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
-
389  " erf.turb_disk_angle_from_x in the inputs which is the angle of the face of the"
-
390  " turbine disk from the x-axis. A turbine facing an oncoming flow in the x-direction"
-
391  " will have turb_disk_angle value of 90 deg.");
-
392  }
-
393  if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
-
394  amrex::Abort("You are using windfarms with latitude-logitude option to position the turbines."
-
395  " For this you should provide the inputs erf.windfarm_x_shift and"
-
396  " erf.windfarm_y_shift which are the values by which the bounding box of the"
-
397  " windfarm is shifted from the x and the y axes.");
-
398  }
-
399  }
-
400 
-
401  void display(int max_level)
-
402  {
-
403  amrex::Print() << "SOLVER CHOICE: " << std::endl;
-
404  amrex::Print() << "force_stage1_single_substep : " << force_stage1_single_substep << std::endl;
-
405  for (int lev = 0; lev <= max_level; lev++) {
-
406  amrex::Print() << "anelastic at level : " << lev << " is " << anelastic[lev] << std::endl;
-
407  if (substepping_type[lev] == SubsteppingType::None) {
-
408  amrex::Print() << "No substepping at level " << lev << std::endl;
-
409  } else if (substepping_type[lev] == SubsteppingType::Explicit) {
-
410  amrex::Print() << "Explicit substepping at level " << lev << std::endl;
-
411  } else if (substepping_type[lev] == SubsteppingType::Implicit) {
-
412  amrex::Print() << "Implicit substepping at level " << lev << std::endl;
-
413  }
-
414  }
-
415  amrex::Print() << "use_coriolis : " << use_coriolis << std::endl;
-
416  amrex::Print() << "use_gravity : " << use_gravity << std::endl;
-
417 
-
418  amrex::Print() << "Terrain Type: " << std::endl;
-
419  if (terrain_type == TerrainType::Static) {
-
420  amrex::Print() << " Static" << std::endl;
-
421  } else if (terrain_type == TerrainType::Moving) {
-
422  amrex::Print() << " Moving" << std::endl;
-
423  } else {
-
424  amrex::Print() << " None" << std::endl;
-
425  }
-
426 
-
427  amrex::Print() << "ABL Driver Type: " << std::endl;
-
428  if (abl_driver_type == ABLDriverType::None) {
+
355  pp.query("windfarm_x_shift",windfarm_x_shift);
+
356  pp.query("windfarm_y_shift",windfarm_y_shift);
+
357  // Test if time averaged data is to be output
+
358  pp.query("time_avg_vel",time_avg_vel);
+
359 
+
360  check_params(max_level);
+
361  }
+
362 
+
363  void check_params(int max_level)
+
364  {
+
365  // Warn for PBL models and moisture - these may not yet be compatible
+
366  for (int lev = 0; lev <= max_level; lev++) {
+
367  if ((moisture_type != MoistureType::None) && (turbChoice[lev].pbl_type != PBLType::None)) {
+
368  amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
+
369  }
+
370  }
+
371  //
+
372  // Buoyancy type check
+
373  //
+
374  if (buoyancy_type != 1 && buoyancy_type != 2 && buoyancy_type != 3 && buoyancy_type != 4) {
+
375  amrex::Abort("buoyancy_type must be 1, 2, 3 or 4");
+
376  }
+
377 
+
378  if (!use_lagged_delta_rt && !(terrain_type == TerrainType::Moving)) {
+
379  amrex::Error("Can't turn off lagged_delta_rt when terrain not moving");
+
380  }
+
381 
+
382  //
+
383  // Wind farm checks
+
384  //
+
385  if (windfarm_type==WindFarmType::SimpleAD and sampling_distance_by_D < 0.0) {
+
386  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
+
387  " erf.sampling_distance_by_D in the inputs which specifies the upstream"
+
388  " distance as a factor of the turbine diameter at which the incoming free stream"
+
389  " velocity will be computed at.");
+
390  }
+
391  if ( (windfarm_type==WindFarmType::SimpleAD ||
+
392  windfarm_type==WindFarmType::GeneralAD ) && turb_disk_angle < 0.0) {
+
393  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
+
394  " erf.turb_disk_angle_from_x in the inputs which is the angle of the face of the"
+
395  " turbine disk from the x-axis. A turbine facing an oncoming flow in the x-direction"
+
396  " will have turb_disk_angle value of 90 deg.");
+
397  }
+
398  if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
+
399  amrex::Abort("You are using windfarms with latitude-logitude option to position the turbines."
+
400  " For this you should provide the inputs erf.windfarm_x_shift and"
+
401  " erf.windfarm_y_shift which are the values by which the bounding box of the"
+
402  " windfarm is shifted from the x and the y axes.");
+
403  }
+
404  }
+
405 
+
406  void display(int max_level)
+
407  {
+
408  amrex::Print() << "SOLVER CHOICE: " << std::endl;
+
409  amrex::Print() << "force_stage1_single_substep : " << force_stage1_single_substep << std::endl;
+
410  for (int lev = 0; lev <= max_level; lev++) {
+
411  amrex::Print() << "anelastic at level : " << lev << " is " << anelastic[lev] << std::endl;
+
412  if (substepping_type[lev] == SubsteppingType::None) {
+
413  amrex::Print() << "No substepping at level " << lev << std::endl;
+
414  } else if (substepping_type[lev] == SubsteppingType::Explicit) {
+
415  amrex::Print() << "Explicit substepping at level " << lev << std::endl;
+
416  } else if (substepping_type[lev] == SubsteppingType::Implicit) {
+
417  amrex::Print() << "Implicit substepping at level " << lev << std::endl;
+
418  }
+
419  }
+
420  amrex::Print() << "use_coriolis : " << use_coriolis << std::endl;
+
421  amrex::Print() << "use_gravity : " << use_gravity << std::endl;
+
422 
+
423  amrex::Print() << "Terrain Type: " << std::endl;
+
424  if (terrain_type == TerrainType::Static) {
+
425  amrex::Print() << " Static" << std::endl;
+
426  } else if (terrain_type == TerrainType::Moving) {
+
427  amrex::Print() << " Moving" << std::endl;
+
428  } else {
429  amrex::Print() << " None" << std::endl;
-
430  } else if (abl_driver_type == ABLDriverType::PressureGradient) {
-
431  amrex::Print() << " Pressure Gradient "
-
432  << amrex::RealVect(abl_pressure_grad[0],abl_pressure_grad[1],abl_pressure_grad[2])
-
433  << std::endl;
-
434  } else if (abl_driver_type == ABLDriverType::GeostrophicWind) {
-
435  amrex::Print() << " Geostrophic Wind "
-
436  << amrex::RealVect(abl_geo_forcing[0],abl_geo_forcing[1],abl_geo_forcing[2])
-
437  << std::endl;
-
438  }
-
439 
-
440  if (max_level > 0) {
-
441  amrex::Print() << "Coupling Type: " << std::endl;
-
442  if (coupling_type == CouplingType::TwoWay) {
-
443  amrex::Print() << " Two-way" << std::endl;
-
444  } else if (coupling_type == CouplingType::OneWay) {
-
445  amrex::Print() << " One-way" << std::endl;
-
446  }
-
447  }
-
448 
-
449  amrex::Print() << "Buoyancy_type : " << buoyancy_type << std::endl;
-
450 
-
451  advChoice.display();
-
452  diffChoice.display();
-
453  spongeChoice.display();
-
454 
-
455  for (int lev = 0; lev <= max_level; lev++) {
-
456  turbChoice[lev].display(lev);
-
457  }
-
458  }
+
430  }
+
431 
+
432  amrex::Print() << "ABL Driver Type: " << std::endl;
+
433  if (abl_driver_type == ABLDriverType::None) {
+
434  amrex::Print() << " None" << std::endl;
+
435  } else if (abl_driver_type == ABLDriverType::PressureGradient) {
+
436  amrex::Print() << " Pressure Gradient "
+
437  << amrex::RealVect(abl_pressure_grad[0],abl_pressure_grad[1],abl_pressure_grad[2])
+
438  << std::endl;
+
439  } else if (abl_driver_type == ABLDriverType::GeostrophicWind) {
+
440  amrex::Print() << " Geostrophic Wind "
+
441  << amrex::RealVect(abl_geo_forcing[0],abl_geo_forcing[1],abl_geo_forcing[2])
+
442  << std::endl;
+
443  }
+
444 
+
445  if (max_level > 0) {
+
446  amrex::Print() << "Coupling Type: " << std::endl;
+
447  if (coupling_type == CouplingType::TwoWay) {
+
448  amrex::Print() << " Two-way" << std::endl;
+
449  } else if (coupling_type == CouplingType::OneWay) {
+
450  amrex::Print() << " One-way" << std::endl;
+
451  }
+
452  }
+
453 
+
454  amrex::Print() << "Buoyancy_type : " << buoyancy_type << std::endl;
+
455 
+
456  advChoice.display();
+
457  diffChoice.display();
+
458  spongeChoice.display();
459 
-
460  void build_coriolis_forcings()
-
461  {
-
462  amrex::ParmParse pp(pp_prefix);
-
463 
-
464  // Read the rotational time period (in seconds)
-
465  amrex::Real rot_time_period = 86400.0;
-
466  pp.query("rotational_time_period", rot_time_period);
-
467 
-
468  coriolis_factor = 2.0 * 2.0 * PI / rot_time_period;
-
469  amrex::Print() << "Coriolis factor = " << coriolis_factor << std::endl;
-
470 
-
471  amrex::Real latitude = 90.0;
-
472  pp.query("latitude", latitude);
-
473 
-
474  pp.query("coriolis_3d", coriolis_3d);
+
460  for (int lev = 0; lev <= max_level; lev++) {
+
461  turbChoice[lev].display(lev);
+
462  }
+
463  }
+
464 
+
465  void build_coriolis_forcings()
+
466  {
+
467  amrex::ParmParse pp(pp_prefix);
+
468 
+
469  // Read the rotational time period (in seconds)
+
470  amrex::Real rot_time_period = 86400.0;
+
471  pp.query("rotational_time_period", rot_time_period);
+
472 
+
473  coriolis_factor = 2.0 * 2.0 * PI / rot_time_period;
+
474  amrex::Print() << "Coriolis factor = " << coriolis_factor << std::endl;
475 
-
476  // Convert to radians
-
477  latitude *= (PI/180.);
-
478  sinphi = std::sin(latitude);
-
479  if (coriolis_3d) {
-
480  cosphi = std::cos(latitude);
-
481  }
-
482 
-
483  amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
-
484 
-
485  if (abl_driver_type == ABLDriverType::GeostrophicWind) {
-
486  // Read in the geostrophic wind -- we only use this to construct
-
487  // the forcing term so no need to keep it
-
488  amrex::Vector<amrex::Real> abl_geo_wind(3);
-
489  pp.queryarr("abl_geo_wind",abl_geo_wind);
-
490 
-
491  if(!pp.query("abl_geo_wind_table",abl_geo_wind_table)) {
-
492  abl_geo_forcing = {
-
493  -coriolis_factor * (abl_geo_wind[1]*sinphi - abl_geo_wind[2]*cosphi),
-
494  coriolis_factor * abl_geo_wind[0]*sinphi,
-
495  -coriolis_factor * abl_geo_wind[0]*cosphi
-
496  };
-
497  } else {
-
498  amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
-
499  }
-
500  }
-
501  }
-
502 
-
503  void read_int_string (int max_level, const char* string_to_read,
-
504  amrex::Vector<int>& vec_to_fill, int default_int)
-
505  {
-
506  amrex::ParmParse pp("erf");
-
507  int nvals = pp.countval(string_to_read);
-
508  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
-
509  amrex::Vector<int> temp; temp.resize(nvals);
-
510  pp.queryarr(string_to_read,temp);
-
511  if (nvals == 0) {
-
512  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(default_int);
-
513  } else if (nvals == 1) {
-
514  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[0]);
-
515  } else {
-
516  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[i]);
-
517  }
-
518  }
-
519 
-
520  inline static
-
521  bool terrain_is_flat = false;
-
522 
-
523  inline static
-
524  TerrainType terrain_type = TerrainType::None;
-
525 
-
526  static
-
527  void set_flat_terrain_flag ()
-
528  {
-
529  terrain_is_flat = true;
-
530  }
-
531 
-
532  AdvChoice advChoice;
-
533  DiffChoice diffChoice;
-
534  SpongeChoice spongeChoice;
-
535  amrex::Vector<TurbChoice> turbChoice;
+
476  amrex::Real latitude = 90.0;
+
477  pp.query("latitude", latitude);
+
478 
+
479  pp.query("coriolis_3d", coriolis_3d);
+
480 
+
481  // Convert to radians
+
482  latitude *= (PI/180.);
+
483  sinphi = std::sin(latitude);
+
484  if (coriolis_3d) {
+
485  cosphi = std::cos(latitude);
+
486  }
+
487 
+
488  amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
+
489 
+
490  if (abl_driver_type == ABLDriverType::GeostrophicWind) {
+
491  // Read in the geostrophic wind -- we only use this to construct
+
492  // the forcing term so no need to keep it
+
493  amrex::Vector<amrex::Real> abl_geo_wind(3);
+
494  pp.queryarr("abl_geo_wind",abl_geo_wind);
+
495 
+
496  if(!pp.query("abl_geo_wind_table",abl_geo_wind_table)) {
+
497  abl_geo_forcing = {
+
498  -coriolis_factor * (abl_geo_wind[1]*sinphi - abl_geo_wind[2]*cosphi),
+
499  coriolis_factor * abl_geo_wind[0]*sinphi,
+
500  -coriolis_factor * abl_geo_wind[0]*cosphi
+
501  };
+
502  } else {
+
503  amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
+
504  }
+
505  }
+
506  }
+
507 
+
508  void read_int_string (int max_level, const char* string_to_read,
+
509  amrex::Vector<int>& vec_to_fill, int default_int)
+
510  {
+
511  amrex::ParmParse pp("erf");
+
512  int nvals = pp.countval(string_to_read);
+
513  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
+
514  amrex::Vector<int> temp; temp.resize(nvals);
+
515  pp.queryarr(string_to_read,temp);
+
516  if (nvals == 0) {
+
517  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(default_int);
+
518  } else if (nvals == 1) {
+
519  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[0]);
+
520  } else {
+
521  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[i]);
+
522  }
+
523  }
+
524 
+
525  inline static
+
526  bool terrain_is_flat = false;
+
527 
+
528  inline static
+
529  TerrainType terrain_type = TerrainType::None;
+
530 
+
531  static
+
532  void set_flat_terrain_flag ()
+
533  {
+
534  terrain_is_flat = true;
+
535  }
536 
-
537  // Default prefix
-
538  std::string pp_prefix {"erf"};
-
539 
-
540  int force_stage1_single_substep = 1;
+
537  AdvChoice advChoice;
+
538  DiffChoice diffChoice;
+
539  SpongeChoice spongeChoice;
+
540  amrex::Vector<TurbChoice> turbChoice;
541 
-
542  amrex::Vector<SubsteppingType> substepping_type;
-
543  amrex::Vector<int> anelastic;
+
542  // Default prefix
+
543  std::string pp_prefix {"erf"};
544 
-
545  int constant_density = 0;
-
546  int ncorr = 1;
-
547  amrex::Real poisson_abstol = 1e-10;
-
548  amrex::Real poisson_reltol = 1e-10;
+
545  int force_stage1_single_substep = 1;
+
546 
+
547  amrex::Vector<SubsteppingType> substepping_type;
+
548  amrex::Vector<int> anelastic;
549 
-
550  bool test_mapfactor = false;
-
551 
-
552  int buoyancy_type = 1; // uses rhoprime directly
-
553 
-
554  // Specify what additional physics/forcing modules we use
-
555  bool use_gravity = false;
-
556  bool use_coriolis = false;
-
557  bool coriolis_3d = true;
+
550  int constant_density = 0;
+
551  int ncorr = 1;
+
552  amrex::Real poisson_abstol = 1e-10;
+
553  amrex::Real poisson_reltol = 1e-10;
+
554 
+
555  bool test_mapfactor = false;
+
556 
+
557  int buoyancy_type = 1; // uses rhoprime directly
558 
-
559  bool rayleigh_damp_U = false;
-
560  bool rayleigh_damp_V = false;
-
561  bool rayleigh_damp_W = false;
-
562  bool rayleigh_damp_T = false;
-
563  amrex::Real rayleigh_dampcoef = 0.2; // inverse time scale [1/s]
-
564  amrex::Real rayleigh_zdamp = 500.0; // damping layer depth [m]
-
565  amrex::Real rayleigh_ztop;
-
566 
-
567  // This defaults to true but can be set to false for moving terrain cases only
-
568  bool use_lagged_delta_rt = true;
-
569 
-
570  // Constants
-
571  amrex::Real gravity;
-
572  amrex::Real c_p = Cp_d; // specific heat at constant pressure for dry air [J/(kg-K)]
-
573  amrex::Real rdOcp;
+
559  // Specify what additional physics/forcing modules we use
+
560  bool use_gravity = false;
+
561  bool use_coriolis = false;
+
562  bool coriolis_3d = true;
+
563 
+
564  bool rayleigh_damp_U = false;
+
565  bool rayleigh_damp_V = false;
+
566  bool rayleigh_damp_W = false;
+
567  bool rayleigh_damp_T = false;
+
568  amrex::Real rayleigh_dampcoef = 0.2; // inverse time scale [1/s]
+
569  amrex::Real rayleigh_zdamp = 500.0; // damping layer depth [m]
+
570  amrex::Real rayleigh_ztop;
+
571 
+
572  // This defaults to true but can be set to false for moving terrain cases only
+
573  bool use_lagged_delta_rt = true;
574 
-
575  // Staggered z levels for vertical grid stretching
-
576  amrex::Real grid_stretching_ratio = 0;
-
577  amrex::Real zsurf = 0.0;
-
578  amrex::Real dz0;
+
575  // Constants
+
576  amrex::Real gravity;
+
577  amrex::Real c_p = Cp_d; // specific heat at constant pressure for dry air [J/(kg-K)]
+
578  amrex::Real rdOcp;
579 
-
580  bool project_initial_velocity = false;
-
581 
-
582  // Coriolis forcing
-
583  amrex::Real coriolis_factor = 0.0;
-
584  amrex::Real cosphi = 0.0;
-
585  amrex::Real sinphi = 0.0;
+
580  // Staggered z levels for vertical grid stretching
+
581  amrex::Real grid_stretching_ratio = 0;
+
582  amrex::Real zsurf = 0.0;
+
583  amrex::Real dz0;
+
584 
+
585  bool project_initial_velocity = false;
586 
-
587  // User-specified forcings in problem definition
-
588  bool custom_rhotheta_forcing = false;
-
589  bool custom_moisture_forcing = false;
-
590  bool custom_w_subsidence = false;
-
591  bool custom_geostrophic_profile = false;
-
592  bool custom_forcing_prim_vars = false;
-
593 
-
594  // Do we use source terms to nudge the solution towards
-
595  // the time-varying data provided in input sounding files?
-
596  bool nudging_from_input_sounding = false;
-
597 
-
598  // User specified MOST BC type
-
599  bool use_explicit_most = false;
-
600 
-
601  // MOST stress rotations
-
602  bool use_rotate_most = false;
-
603 
-
604  // User wishes to output time averaged velocity fields
-
605  bool time_avg_vel = false;
-
606 
-
607  // Type of perturbation
-
608  PerturbationType pert_type;
-
609 
-
610  // Numerical diffusion
-
611  bool use_NumDiff{false};
-
612  amrex::Real NumDiffCoeff{0.};
-
613 
-
614  // Monotonic advection limiter
-
615  bool use_mono_adv{false};
-
616 
-
617  CouplingType coupling_type;
-
618  MoistureType moisture_type;
-
619  WindFarmType windfarm_type;
-
620  WindFarmLocType windfarm_loc_type;
-
621  LandSurfaceType lsm_type;
-
622 
-
623  ABLDriverType abl_driver_type;
-
624  amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> abl_pressure_grad;
-
625  amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> abl_geo_forcing;
-
626  std::string abl_geo_wind_table;
-
627  bool have_geo_wind_profile {false};
-
628 
-
629  int ave_plane {2};
-
630  // Microphysics params
-
631  bool do_cloud {true};
-
632  bool do_precip {true};
-
633  bool use_moist_background {false};
-
634  int RhoQv_comp {-1};
-
635 
-
636  // This component will be model-dependent:
-
637  // if a model with no rain, this will stay -1
-
638  // if Kessler, then it will be set to RhoQ3
-
639  // if SAM, then it will be set to RhoQ4
-
640  int RhoQr_comp {-1};
-
641 
-
642  std::string windfarm_loc_table, windfarm_spec_table, windfarm_spec_table_extra;
-
643  std::string windfarm_blade_table, windfarm_airfoil_tables;
-
644  amrex::Real sampling_distance_by_D = -1.0;
-
645  amrex::Real turb_disk_angle = -1.0;
-
646  amrex::Real windfarm_x_shift = -1.0;
-
647  amrex::Real windfarm_y_shift = -1.0;
-
648 
-
649  // Flag for valid canopy model
-
650  bool do_forest {false};
-
651 };
-
652 #endif
+
587  // Coriolis forcing
+
588  amrex::Real coriolis_factor = 0.0;
+
589  amrex::Real cosphi = 0.0;
+
590  amrex::Real sinphi = 0.0;
+
591 
+
592  // User-specified forcings in problem definition
+
593  bool custom_rhotheta_forcing = false;
+
594  bool custom_moisture_forcing = false;
+
595  bool custom_w_subsidence = false;
+
596  bool custom_geostrophic_profile = false;
+
597  bool custom_forcing_prim_vars = false;
+
598 
+
599  // Do we use source terms to nudge the solution towards
+
600  // the time-varying data provided in input sounding files?
+
601  bool nudging_from_input_sounding = false;
+
602 
+
603  // User specified MOST BC type
+
604  bool use_explicit_most = false;
+
605 
+
606  // MOST stress rotations
+
607  bool use_rotate_most = false;
+
608 
+
609  // User wishes to output time averaged velocity fields
+
610  bool time_avg_vel = false;
+
611 
+
612  // Type of perturbation
+
613  PerturbationType pert_type;
+
614 
+
615  // Numerical diffusion
+
616  bool use_NumDiff{false};
+
617  amrex::Real NumDiffCoeff{0.};
+
618 
+
619  // Monotonic advection limiter
+
620  bool use_mono_adv{false};
+
621 
+
622  CouplingType coupling_type;
+
623  MoistureType moisture_type;
+
624  WindFarmType windfarm_type;
+
625  WindFarmLocType windfarm_loc_type;
+
626  LandSurfaceType lsm_type;
+
627 
+
628  ABLDriverType abl_driver_type;
+
629  amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> abl_pressure_grad;
+
630  amrex::GpuArray<amrex::Real, AMREX_SPACEDIM> abl_geo_forcing;
+
631  std::string abl_geo_wind_table;
+
632  bool have_geo_wind_profile {false};
+
633 
+
634  int ave_plane {2};
+
635  // Microphysics params
+
636  bool do_cloud {true};
+
637  bool do_precip {true};
+
638  bool use_moist_background {false};
+
639  int RhoQv_comp {-1};
+
640 
+
641  // This component will be model-dependent:
+
642  // if a model with no rain, this will stay -1
+
643  // if Kessler, then it will be set to RhoQ3
+
644  // if SAM, then it will be set to RhoQ4
+
645  int RhoQr_comp {-1};
+
646 
+
647  std::string windfarm_loc_table, windfarm_spec_table, windfarm_spec_table_extra;
+
648  std::string windfarm_blade_table, windfarm_airfoil_tables;
+
649  amrex::Real sampling_distance_by_D = -1.0;
+
650  amrex::Real turb_disk_angle = -1.0;
+
651  amrex::Real windfarm_x_shift = -1.0;
+
652  amrex::Real windfarm_y_shift = -1.0;
+
653 
+
654  // Flag for valid canopy model
+
655  bool do_forest {false};
+
656 };
+
657 #endif
ERF_AdvStruct.H
ERF_Constants.H
Cp_d
constexpr amrex::Real Cp_d
Definition: ERF_Constants.H:12
@@ -784,88 +789,88 @@
DiffChoice::init_params
void init_params()
Definition: ERF_DiffStruct.H:21
DiffChoice::display
void display()
Definition: ERF_DiffStruct.H:56
SolverChoice
Definition: ERF_DataStruct.H:78
-
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:562
-
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:578
-
SolverChoice::use_lagged_delta_rt
bool use_lagged_delta_rt
Definition: ERF_DataStruct.H:568
-
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:583
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::windfarm_x_shift
amrex::Real windfarm_x_shift
Definition: ERF_DataStruct.H:646
-
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:560
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:548
-
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:564
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::read_int_string
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:503
-
SolverChoice::windfarm_spec_table
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:642
-
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:615
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:527
-
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:555
-
SolverChoice::build_coriolis_forcings
void build_coriolis_forcings()
Definition: ERF_DataStruct.H:460
-
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:546
-
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:540
-
SolverChoice::windfarm_spec_table_extra
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:642
-
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:584
-
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:621
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::windfarm_loc_table
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:642
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:538
-
SolverChoice::check_params
void check_params(int max_level)
Definition: ERF_DataStruct.H:358
-
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:588
-
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:625
-
SolverChoice::windfarm_loc_type
WindFarmLocType windfarm_loc_type
Definition: ERF_DataStruct.H:620
-
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:590
-
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:596
-
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:559
-
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:591
-
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:565
-
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:576
-
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:585
-
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:627
-
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:650
-
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:624
-
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:612
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:542
-
SolverChoice::coriolis_3d
bool coriolis_3d
Definition: ERF_DataStruct.H:557
+
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:567
+
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:583
+
SolverChoice::use_lagged_delta_rt
bool use_lagged_delta_rt
Definition: ERF_DataStruct.H:573
+
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:588
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::windfarm_x_shift
amrex::Real windfarm_x_shift
Definition: ERF_DataStruct.H:651
+
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:565
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:553
+
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:569
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::read_int_string
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:508
+
SolverChoice::windfarm_spec_table
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:647
+
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:620
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:532
+
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:560
+
SolverChoice::build_coriolis_forcings
void build_coriolis_forcings()
Definition: ERF_DataStruct.H:465
+
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:551
+
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:545
+
SolverChoice::windfarm_spec_table_extra
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:647
+
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:589
+
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:626
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::windfarm_loc_table
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:647
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:543
+
SolverChoice::check_params
void check_params(int max_level)
Definition: ERF_DataStruct.H:363
+
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:593
+
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:630
+
SolverChoice::windfarm_loc_type
WindFarmLocType windfarm_loc_type
Definition: ERF_DataStruct.H:625
+
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:595
+
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:601
+
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:564
+
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:596
+
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:570
+
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:581
+
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:590
+
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:632
+
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:655
+
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:629
+
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:617
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:547
+
SolverChoice::coriolis_3d
bool coriolis_3d
Definition: ERF_DataStruct.H:562
SolverChoice::init_params
void init_params(int max_level)
Definition: ERF_DataStruct.H:80
-
SolverChoice::sampling_distance_by_D
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:644
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:550
-
SolverChoice::display
void display(int max_level)
Definition: ERF_DataStruct.H:401
-
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:556
-
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:611
-
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:589
-
SolverChoice::windfarm_blade_table
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:643
-
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:577
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:580
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:545
-
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:532
-
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:633
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:592
-
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:626
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::abl_driver_type
ABLDriverType abl_driver_type
Definition: ERF_DataStruct.H:623
-
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:561
-
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:608
-
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:534
-
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:619
-
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:552
-
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:547
-
SolverChoice::turb_disk_angle
amrex::Real turb_disk_angle
Definition: ERF_DataStruct.H:645
-
SolverChoice::windfarm_y_shift
amrex::Real windfarm_y_shift
Definition: ERF_DataStruct.H:647
-
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:521
-
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:605
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
-
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:563
-
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:617
-
SolverChoice::windfarm_airfoil_tables
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:643
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::sampling_distance_by_D
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:649
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:555
+
SolverChoice::display
void display(int max_level)
Definition: ERF_DataStruct.H:406
+
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:561
+
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:616
+
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:594
+
SolverChoice::windfarm_blade_table
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:648
+
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:582
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:585
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:550
+
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:537
+
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:638
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:597
+
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:631
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::abl_driver_type
ABLDriverType abl_driver_type
Definition: ERF_DataStruct.H:628
+
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:566
+
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:613
+
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:539
+
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:624
+
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:557
+
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:552
+
SolverChoice::turb_disk_angle
amrex::Real turb_disk_angle
Definition: ERF_DataStruct.H:650
+
SolverChoice::windfarm_y_shift
amrex::Real windfarm_y_shift
Definition: ERF_DataStruct.H:652
+
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:526
+
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:610
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
+
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:568
+
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:622
+
SolverChoice::windfarm_airfoil_tables
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:648
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
SpongeChoice
Definition: ERF_SpongeStruct.H:15
SpongeChoice::display
void display()
Definition: ERF_SpongeStruct.H:45
SpongeChoice::init_params
void init_params()
Definition: ERF_SpongeStruct.H:17
diff --git a/ERF__DiffusionSrcForState__N_8cpp.html b/ERF__DiffusionSrcForState__N_8cpp.html index f2e59c8b6..68197e254 100644 --- a/ERF__DiffusionSrcForState__N_8cpp.html +++ b/ERF__DiffusionSrcForState__N_8cpp.html @@ -1078,10 +1078,10 @@

DiffChoice::alpha_C
amrex::Real alpha_C
Definition: ERF_DiffStruct.H:80
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
MYNNLevel25::B1
amrex::Real B1
Definition: ERF_MYNNStruct.H:43
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_mynn
MYNNLevel25 pbl_mynn
Definition: ERF_TurbStruct.H:194
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
diff --git a/ERF__DiffusionSrcForState__T_8cpp.html b/ERF__DiffusionSrcForState__T_8cpp.html index d0dac52b3..c98038bb6 100644 --- a/ERF__DiffusionSrcForState__T_8cpp.html +++ b/ERF__DiffusionSrcForState__T_8cpp.html @@ -1246,10 +1246,10 @@

DiffChoice::alpha_C
amrex::Real alpha_C
Definition: ERF_DiffStruct.H:80
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
MYNNLevel25::B1
amrex::Real B1
Definition: ERF_MYNNStruct.H:43
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_mynn
MYNNLevel25 pbl_mynn
Definition: ERF_TurbStruct.H:194
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
diff --git a/ERF__Kessler_8H_source.html b/ERF__Kessler_8H_source.html index 124e43a5c..8b3abbe15 100644 --- a/ERF__Kessler_8H_source.html +++ b/ERF__Kessler_8H_source.html @@ -327,11 +327,11 @@
MicVar_Kess::NumVars
@ NumVars
Definition: ERF_Kessler.H:42
MicVar_Kess::qv
@ qv
Definition: ERF_Kessler.H:36
SolverChoice
Definition: ERF_DataStruct.H:78
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
diff --git a/ERF__SAM_8H_source.html b/ERF__SAM_8H_source.html index 0869fe3bc..94c461bf0 100644 --- a/ERF__SAM_8H_source.html +++ b/ERF__SAM_8H_source.html @@ -536,12 +536,12 @@
MicVar::qps
@ qps
Definition: ERF_SAM.H:43
MicVar::qp
@ qp
Definition: ERF_SAM.H:41
SolverChoice
Definition: ERF_DataStruct.H:78
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
diff --git a/ERF__TI__substep__fun_8H_source.html b/ERF__TI__substep__fun_8H_source.html index 34f2cd48c..b7c9d3185 100644 --- a/ERF__TI__substep__fun_8H_source.html +++ b/ERF__TI__substep__fun_8H_source.html @@ -290,7 +290,7 @@
IntVars::cons
@ cons
Definition: ERF_IndexDefines.H:139
IntVars::zmom
@ zmom
Definition: ERF_IndexDefines.H:142
IntVars::xmom
@ xmom
Definition: ERF_IndexDefines.H:140
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
diff --git a/ERF__make__buoyancy_8cpp.html b/ERF__make__buoyancy_8cpp.html index 50b43bd68..3fcb04588 100644 --- a/ERF__make__buoyancy_8cpp.html +++ b/ERF__make__buoyancy_8cpp.html @@ -386,10 +386,10 @@

BaseState::th0_comp
@ th0_comp
Definition: ERF_IndexDefines.H:66
BaseState::r0_comp
@ r0_comp
Definition: ERF_IndexDefines.H:63
IntVars::cons
@ cons
Definition: ERF_IndexDefines.H:139
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:552
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:557

Referenced by make_mom_sources().

diff --git a/ERF__make__mom__sources_8cpp.html b/ERF__make__mom__sources_8cpp.html index f0e4d3d6b..77d98676b 100644 --- a/ERF__make__mom__sources_8cpp.html +++ b/ERF__make__mom__sources_8cpp.html @@ -853,27 +853,27 @@

InputSoundingData::input_sounding_time
amrex::Vector< amrex::Real > input_sounding_time
Definition: ERF_InputSoundingData.H:315

InputSoundingData::tau_nudging
amrex::Real tau_nudging
Definition: ERF_InputSoundingData.H:312
InputSoundingData::U_inp_sound_d
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > U_inp_sound_d
Definition: ERF_InputSoundingData.H:326
-
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:583
-
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:560
-
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:564
-
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:584
-
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:625
-
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:590
-
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:596
-
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:559
-
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:565
-
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:585
-
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:627
-
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:650
-
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:556
-
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:611
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:592
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:561
-
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:534
-
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:563
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:588
+
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:565
+
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:569
+
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:589
+
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:630
+
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:595
+
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:601
+
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:564
+
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:570
+
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:590
+
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:632
+
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:655
+
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:561
+
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:616
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:597
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:566
+
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:539
+
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:568
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
SpongeChoice::sponge_type
std::string sponge_type
Definition: ERF_SpongeStruct.H:61
Here is the call graph for this function:
diff --git a/ERF__make__sources_8cpp.html b/ERF__make__sources_8cpp.html index 7326723bb..7d88bb78b 100644 --- a/ERF__make__sources_8cpp.html +++ b/ERF__make__sources_8cpp.html @@ -689,23 +689,23 @@

InputSoundingData::theta_inp_sound_d
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > theta_inp_sound_d
Definition: ERF_InputSoundingData.H:326
InputSoundingData::input_sounding_time
amrex::Vector< amrex::Real > input_sounding_time
Definition: ERF_InputSoundingData.H:315
InputSoundingData::tau_nudging
amrex::Real tau_nudging
Definition: ERF_InputSoundingData.H:312
-
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:562
-
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:564
-
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:588
-
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:590
-
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:596
-
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:565
-
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:612
-
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:611
-
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:589
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:592
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:608
-
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:534
-
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:563
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:567
+
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:569
+
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:593
+
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:595
+
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:601
+
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:570
+
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:617
+
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:616
+
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:594
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:597
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:613
+
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:539
+
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:568
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
SpongeChoice::sponge_type
std::string sponge_type
Definition: ERF_SpongeStruct.H:61
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
diff --git a/ERF__make__tau__terms_8cpp.html b/ERF__make__tau__terms_8cpp.html index 231d97dac..875709885 100644 --- a/ERF__make__tau__terms_8cpp.html +++ b/ERF__make__tau__terms_8cpp.html @@ -758,11 +758,11 @@

DiffChoice::rho0_trans
amrex::Real rho0_trans
Definition: ERF_DiffStruct.H:83
DiffChoice::dynamicViscosity
amrex::Real dynamicViscosity
Definition: ERF_DiffStruct.H:88
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
TurbChoice::les_type
LESType les_type
Definition: ERF_TurbStruct.H:169
diff --git a/ERF__prob__common_8H_source.html b/ERF__prob__common_8H_source.html index 60d42552f..500e0d417 100644 --- a/ERF__prob__common_8H_source.html +++ b/ERF__prob__common_8H_source.html @@ -611,7 +611,7 @@
ProbParmDefaults::T_0
amrex::Real T_0
Definition: ERF_prob_common.H:16
ProbParmDefaults::rho_0
amrex::Real rho_0
Definition: ERF_prob_common.H:15
SolverChoice
Definition: ERF_DataStruct.H:78
-
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:527
+
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:532
diff --git a/ERF__slow__rhs__post_8cpp.html b/ERF__slow__rhs__post_8cpp.html index fcce66d3c..9e28dcac1 100644 --- a/ERF__slow__rhs__post_8cpp.html +++ b/ERF__slow__rhs__post_8cpp.html @@ -932,17 +932,17 @@

AdvChoice::dryscal_horiz_upw_frac
amrex::Real dryscal_horiz_upw_frac
Definition: ERF_AdvStruct.H:293
DiffChoice
Definition: ERF_DiffStruct.H:19
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:615
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:532
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
-
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:617
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:620
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:537
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
+
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:622
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
TurbChoice::use_KE
bool use_KE
Definition: ERF_TurbStruct.H:206
diff --git a/ERF__slow__rhs__pre_8cpp.html b/ERF__slow__rhs__pre_8cpp.html index 81ba2feac..02d87523f 100644 --- a/ERF__slow__rhs__pre_8cpp.html +++ b/ERF__slow__rhs__pre_8cpp.html @@ -1276,19 +1276,19 @@

AdvChoice::dycore_horiz_upw_frac
amrex::Real dycore_horiz_upw_frac
Definition: ERF_AdvStruct.H:291
DiffChoice
Definition: ERF_DiffStruct.H:19
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:615
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:624
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:545
-
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:532
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
-
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:617
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:620
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:629
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:550
+
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:537
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
+
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:622
TurbChoice
Definition: ERF_TurbStruct.H:29
TurbChoice::pbl_type
PBLType pbl_type
Definition: ERF_TurbStruct.H:192
TurbChoice::les_type
LESType les_type
Definition: ERF_TurbStruct.H:169
diff --git a/classERF.html b/classERF.html index 28ba6ff52..5cd21641a 100644 --- a/classERF.html +++ b/classERF.html @@ -1387,12 +1387,12 @@

Vars::cons
@ cons
Definition: ERF_IndexDefines.H:129
Vars::zvel
@ zvel
Definition: ERF_IndexDefines.H:132
Vars::yvel
@ yvel
Definition: ERF_IndexDefines.H:131
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:608
-
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:619
-
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:605
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:613
+
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:624
+
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:610
TurbulentPerturbation::calc_tpi_update
void calc_tpi_update(const int lev, const amrex::Real dt, amrex::MultiFab &mf_xvel, amrex::MultiFab &mf_yvel, amrex::MultiFab &mf_cons)
Definition: ERF_TurbPertStruct.H:184
TurbulentPerturbation::pb_cell
amrex::MultiFab pb_cell
Definition: ERF_TurbPertStruct.H:536
TurbulentPerturbation::apply_tpi
void apply_tpi(const int &lev, const amrex::Box &vbx, const int &comp, const amrex::IndexType &m_ixtype, const amrex::Array4< amrex::Real > &src_arr, const amrex::Array4< amrex::Real const > &pert_cell)
Definition: ERF_TurbPertStruct.H:245
@@ -1871,22 +1871,22 @@

BaseState::r0_comp
@ r0_comp
Definition: ERF_IndexDefines.H:63
DiffChoice
Definition: ERF_DiffStruct.H:19
DiffChoice::molec_diff_type
MolecDiffType molec_diff_type
Definition: ERF_DiffStruct.H:76
-
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:562
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:560
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:588
-
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:590
-
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:559
-
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:591
-
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:542
-
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:611
-
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:589
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:561
-
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:534
+
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:567
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:565
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:593
+
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:595
+
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:564
+
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:596
+
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:547
+
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:616
+
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:594
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:566
+
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:539
SpongeChoice
Definition: ERF_SpongeStruct.H:15
SpongeChoice::sponge_type
std::string sponge_type
Definition: ERF_SpongeStruct.H:61
TurbChoice
Definition: ERF_TurbStruct.H:29
@@ -1972,7 +1972,7 @@

12 }
ERF::lsm
LandSurface lsm
Definition: ERF.H:777
LandSurface::Advance
void Advance(const int &lev, const amrex::Real &dt_advance)
Definition: ERF_LandSurface.H:51
-
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:621
+
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:626
@@ -2160,7 +2160,7 @@

31  }
32 }
ERF::AverageDownTo
void AverageDownTo(int crse_lev, int scomp, int ncomp)
Definition: ERF_AverageDown.cpp:36
-
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:617
+
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:622
@@ -2662,7 +2662,7 @@

ERF::stretched_dz_d
amrex::Vector< amrex::Gpu::DeviceVector< amrex::Real > > stretched_dz_d
Definition: ERF.H:847
ERF::ax
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ax
Definition: ERF.H:824
ERF::ay
amrex::Vector< std::unique_ptr< amrex::MultiFab > > ay
Definition: ERF.H:825
-
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:521
+
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:526
@@ -4623,8 +4623,8 @@

getRhoThetagivenP
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getRhoThetagivenP(const amrex::Real p, const amrex::Real qv=0.0)
Definition: ERF_EOS.H:175
getExnergivenP
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real getExnergivenP(const amrex::Real P, const amrex::Real rdOcp)
Definition: ERF_EOS.H:144
MicVar_Kess::pres
@ pres
Definition: ERF_Kessler.H:33
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
Here is the call graph for this function:
@@ -4955,9 +4955,9 @@

ERF::initializeMicrophysics
void initializeMicrophysics(const int &)
Definition: ERF.cpp:1214

LandSurface::ReSize
void ReSize(const int &nlev)
Definition: ERF_LandSurface.H:23
amrex::buildInfoGetGitHash
const char * buildInfoGetGitHash(int i)
-
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:578
-
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:576
-
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:577
+
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:583
+
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:581
+
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:582
Here is the call graph for this function:
@@ -5349,7 +5349,7 @@

ERF::fixed_fast_dt
amrex::Vector< amrex::Real > fixed_fast_dt
Definition: ERF.H:915

ERF::cfl
static amrex::Real cfl
Definition: ERF.H:909
MicVar_Kess::rho
@ rho
Definition: ERF_Kessler.H:30
-
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:540
+
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:545
Here is the call graph for this function:
@@ -8063,7 +8063,7 @@

AdvChoice::zero_yflux
amrex::Vector< amrex::IntVect > zero_yflux
Definition: ERF_AdvStruct.H:300

AdvChoice::zero_xflux
amrex::Vector< amrex::IntVect > zero_xflux
Definition: ERF_AdvStruct.H:299
AdvChoice::zero_zflux
amrex::Vector< amrex::IntVect > zero_zflux
Definition: ERF_AdvStruct.H:301
-
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:532
+
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:537
@@ -8211,7 +8211,7 @@

ERF::make_physbcs
void make_physbcs(int lev)
Definition: ERF_make_new_arrays.cpp:550
ERF::init_uniform
void init_uniform(int lev)
Definition: ERF_init_uniform.cpp:17
ERF::turbPert_amplitude
void turbPert_amplitude(const int lev)
Definition: ERF_init_TurbPert.cpp:44
-
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:555
+
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:560
@@ -8615,7 +8615,7 @@

NDRY
#define NDRY
Definition: ERF_IndexDefines.H:13
ERF::ComputeGhostCells
static AMREX_FORCE_INLINE int ComputeGhostCells(const AdvChoice &advChoice, bool use_num_diff)
Definition: ERF.H:1145
BaseState::num_comps
@ num_comps
Definition: ERF_IndexDefines.H:67
-
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:550
+
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:555
TurbulentPerturbation::init_tpi
void init_tpi(const int lev, const amrex::IntVect &nx, const amrex::GpuArray< amrex::Real, 3 > dx, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm, const int ngrow_state)
Definition: ERF_TurbPertStruct.H:29
@@ -9411,10 +9411,10 @@

ERF::samplepoint
amrex::Vector< amrex::IntVect > samplepoint
Definition: ERF.H:1341
ERF::setRecordSamplePointInfo
void setRecordSamplePointInfo(int i, int lev, amrex::IntVect &cell, const std::string &filename)
Definition: ERF.H:1297
Microphysics::modelType
static MoistureModelType modelType(const MoistureType a_moisture_type)
query if a specified moisture model is Eulerian or Lagrangian
Definition: ERF_Microphysics.H:64
-
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:627
-
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:580
-
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:626
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
+
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:632
+
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:585
+
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:631
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
TurbulentPerturbation::debug
void debug(amrex::Real)
Definition: ERF_TurbPertStruct.H:508
Here is the call graph for this function:
@@ -9649,7 +9649,7 @@

127 }
ChopGrids2D
void ChopGrids2D(BoxArray &ba, const Box &domain, int target_size)
Definition: ERF_ChopGrids.cpp:21
ERF::erf_enforce_hse
void erf_enforce_hse(int lev, amrex::MultiFab &dens, amrex::MultiFab &pres, amrex::MultiFab &pi, amrex::MultiFab &th, std::unique_ptr< amrex::MultiFab > &z_cc)
Definition: ERF_init1d.cpp:149
-
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:633
+
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:638
Here is the call graph for this function:
@@ -9906,8 +9906,8 @@

44  }

45 }
ERF::h_rayleigh_ptrs
amrex::Vector< amrex::Vector< amrex::Vector< amrex::Real > > > h_rayleigh_ptrs
Definition: ERF.H:1100
-
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:564
-
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:565
+
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:569
+
SolverChoice::rayleigh_ztop
amrex::Real rayleigh_ztop
Definition: ERF_DataStruct.H:570
@@ -10497,7 +10497,7 @@

ERF::init_zphys
void init_zphys(int lev, amrex::Real time)
Definition: ERF_make_new_arrays.cpp:453
ERF::init_stuff
void init_stuff(int lev, const amrex::BoxArray &ba, const amrex::DistributionMapping &dm, amrex::Vector< amrex::MultiFab > &lev_new, amrex::Vector< amrex::MultiFab > &lev_old, amrex::MultiFab &tmp_base_state, std::unique_ptr< amrex::MultiFab > &tmp_zphys_nd)
Definition: ERF_make_new_arrays.cpp:23
ERF::Define_ERFFillPatchers
void Define_ERFFillPatchers(int lev)
Definition: ERF.cpp:1902
-
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:650
+
SolverChoice::do_forest
bool do_forest
Definition: ERF_DataStruct.H:655
@@ -11892,7 +11892,7 @@

ERF::solve_with_mlmg
void solve_with_mlmg(int lev, amrex::Vector< amrex::MultiFab > &rhs, amrex::Vector< amrex::MultiFab > &p, amrex::Vector< amrex::Array< amrex::MultiFab, AMREX_SPACEDIM >> &fluxes)
Definition: ERF_solve_with_mlmg.cpp:40
ERF::solve_with_gmres
void solve_with_gmres(int lev, amrex::Vector< amrex::MultiFab > &rhs, amrex::Vector< amrex::MultiFab > &p, amrex::Array< amrex::MultiFab, AMREX_SPACEDIM > &fluxes)
Definition: ERF_solve_with_gmres.cpp:12
MicVar::omega
@ omega
Definition: ERF_SAM.H:49
-
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:547
+
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:552
Here is the call graph for this function:
@@ -12161,8 +12161,8 @@

ERF::get_projection_bc
amrex::Array< amrex::LinOpBCType, AMREX_SPACEDIM > get_projection_bc(amrex::Orientation::Side side) const noexcept
Definition: ERF_solve_with_mlmg.cpp:17

ERF::projection_has_dirichlet
bool projection_has_dirichlet(amrex::Array< amrex::LinOpBCType, AMREX_SPACEDIM > bcs) const
Definition: ERF_PoissonSolve_tb.cpp:8
ERF::volWgtSumMF
amrex::Real volWgtSumMF(int lev, const amrex::MultiFab &mf, int comp, const amrex::MultiFab &mapfac, bool local, bool finemask)
Definition: ERF_WriteScalarProfiles.cpp:378
-
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:548
-
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:546
+
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:553
+
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:551
Here is the call graph for this function:
@@ -12921,9 +12921,9 @@

NOAH
Definition: ERF_NOAH.H:20

NullSurf
Definition: ERF_NullSurf.H:8
SLM
Definition: ERF_SLM.H:26
-
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:538
+
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:543
SolverChoice::init_params
void init_params(int max_level)
Definition: ERF_DataStruct.H:80
-
SolverChoice::display
void display(int max_level)
Definition: ERF_DataStruct.H:401
+
SolverChoice::display
void display(int max_level)
Definition: ERF_DataStruct.H:406
Here is the call graph for this function:
diff --git a/classKessler.html b/classKessler.html index f32023356..279b7c73b 100644 --- a/classKessler.html +++ b/classKessler.html @@ -682,7 +682,7 @@

MicVar_Kess::qt
@ qt
Definition: ERF_Kessler.H:35

MicVar_Kess::rain_accum
@ rain_accum
Definition: ERF_Kessler.H:41
MicVar_Kess::qv
@ qv
Definition: ERF_Kessler.H:36
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623

Referenced by Advance().

@@ -927,10 +927,10 @@

Kessler::m_gOcp
amrex::Real m_gOcp
Definition: ERF_Kessler.H:174

Kessler::m_fac_fus
amrex::Real m_fac_fus
Definition: ERF_Kessler.H:172
Kessler::m_axis
int m_axis
Definition: ERF_Kessler.H:165
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
diff --git a/classProblemBase.html b/classProblemBase.html index 7cf263086..ed6fcec8e 100644 --- a/classProblemBase.html +++ b/classProblemBase.html @@ -636,7 +636,7 @@

335 
336  SolverChoice::set_flat_terrain_flag();
337  }
-
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:527
+
SolverChoice::set_flat_terrain_flag
static void set_flat_terrain_flag()
Definition: ERF_DataStruct.H:532
Here is the call graph for this function:
diff --git a/classSAM.html b/classSAM.html index a0d6cdd88..ba5875165 100644 --- a/classSAM.html +++ b/classSAM.html @@ -617,7 +617,7 @@

MicVar::theta
@ theta
Definition: ERF_SAM.H:31

MicVar::qcl
@ qcl
Definition: ERF_SAM.H:38
MicVar::tabs
@ tabs
Definition: ERF_SAM.H:32
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623

Referenced by Advance().

@@ -1116,11 +1116,11 @@

CONST_GRAV
constexpr amrex::Real CONST_GRAV
Definition: ERF_Constants.H:21

SAM::doprecip
bool doprecip
Definition: ERF_SAM.H:311
SAM::docloud
bool docloud
Definition: ERF_SAM.H:311
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
diff --git a/structSolverChoice.html b/structSolverChoice.html index 5ead8f251..f924e7c8b 100644 --- a/structSolverChoice.html +++ b/structSolverChoice.html @@ -327,57 +327,57 @@

-
461  {
-
462  amrex::ParmParse pp(pp_prefix);
-
463 
-
464  // Read the rotational time period (in seconds)
-
465  amrex::Real rot_time_period = 86400.0;
-
466  pp.query("rotational_time_period", rot_time_period);
-
467 
-
468  coriolis_factor = 2.0 * 2.0 * PI / rot_time_period;
-
469  amrex::Print() << "Coriolis factor = " << coriolis_factor << std::endl;
-
470 
-
471  amrex::Real latitude = 90.0;
-
472  pp.query("latitude", latitude);
-
473 
-
474  pp.query("coriolis_3d", coriolis_3d);
+
466  {
+
467  amrex::ParmParse pp(pp_prefix);
+
468 
+
469  // Read the rotational time period (in seconds)
+
470  amrex::Real rot_time_period = 86400.0;
+
471  pp.query("rotational_time_period", rot_time_period);
+
472 
+
473  coriolis_factor = 2.0 * 2.0 * PI / rot_time_period;
+
474  amrex::Print() << "Coriolis factor = " << coriolis_factor << std::endl;
475 
-
476  // Convert to radians
-
477  latitude *= (PI/180.);
-
478  sinphi = std::sin(latitude);
-
479  if (coriolis_3d) {
-
480  cosphi = std::cos(latitude);
-
481  }
-
482 
-
483  amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
-
484 
-
485  if (abl_driver_type == ABLDriverType::GeostrophicWind) {
-
486  // Read in the geostrophic wind -- we only use this to construct
-
487  // the forcing term so no need to keep it
-
488  amrex::Vector<amrex::Real> abl_geo_wind(3);
-
489  pp.queryarr("abl_geo_wind",abl_geo_wind);
-
490 
-
491  if(!pp.query("abl_geo_wind_table",abl_geo_wind_table)) {
-
492  abl_geo_forcing = {
-
493  -coriolis_factor * (abl_geo_wind[1]*sinphi - abl_geo_wind[2]*cosphi),
-
494  coriolis_factor * abl_geo_wind[0]*sinphi,
-
495  -coriolis_factor * abl_geo_wind[0]*cosphi
-
496  };
-
497  } else {
-
498  amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
-
499  }
-
500  }
-
501  }
+
476  amrex::Real latitude = 90.0;
+
477  pp.query("latitude", latitude);
+
478 
+
479  pp.query("coriolis_3d", coriolis_3d);
+
480 
+
481  // Convert to radians
+
482  latitude *= (PI/180.);
+
483  sinphi = std::sin(latitude);
+
484  if (coriolis_3d) {
+
485  cosphi = std::cos(latitude);
+
486  }
+
487 
+
488  amrex::Print() << "Coriolis frequency, f = " << coriolis_factor * sinphi << " 1/s" << std::endl;
+
489 
+
490  if (abl_driver_type == ABLDriverType::GeostrophicWind) {
+
491  // Read in the geostrophic wind -- we only use this to construct
+
492  // the forcing term so no need to keep it
+
493  amrex::Vector<amrex::Real> abl_geo_wind(3);
+
494  pp.queryarr("abl_geo_wind",abl_geo_wind);
+
495 
+
496  if(!pp.query("abl_geo_wind_table",abl_geo_wind_table)) {
+
497  abl_geo_forcing = {
+
498  -coriolis_factor * (abl_geo_wind[1]*sinphi - abl_geo_wind[2]*cosphi),
+
499  coriolis_factor * abl_geo_wind[0]*sinphi,
+
500  -coriolis_factor * abl_geo_wind[0]*cosphi
+
501  };
+
502  } else {
+
503  amrex::Print() << "NOTE: abl_geo_wind_table provided, ignoring input abl_geo_wind" << std::endl;
+
504  }
+
505  }
+
506  }
PI
constexpr amrex::Real PI
Definition: ERF_Constants.H:6
pp
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real pp(amrex::Real y)
Definition: ERF_Microphysics_Utils.H:219
-
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:583
-
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:584
-
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:538
-
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:625
-
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:585
-
SolverChoice::coriolis_3d
bool coriolis_3d
Definition: ERF_DataStruct.H:557
-
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:626
-
SolverChoice::abl_driver_type
ABLDriverType abl_driver_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::coriolis_factor
amrex::Real coriolis_factor
Definition: ERF_DataStruct.H:588
+
SolverChoice::cosphi
amrex::Real cosphi
Definition: ERF_DataStruct.H:589
+
SolverChoice::pp_prefix
std::string pp_prefix
Definition: ERF_DataStruct.H:543
+
SolverChoice::abl_geo_forcing
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_geo_forcing
Definition: ERF_DataStruct.H:630
+
SolverChoice::sinphi
amrex::Real sinphi
Definition: ERF_DataStruct.H:590
+
SolverChoice::coriolis_3d
bool coriolis_3d
Definition: ERF_DataStruct.H:562
+
SolverChoice::abl_geo_wind_table
std::string abl_geo_wind_table
Definition: ERF_DataStruct.H:631
+
SolverChoice::abl_driver_type
ABLDriverType abl_driver_type
Definition: ERF_DataStruct.H:628

Referenced by init_params().

@@ -424,58 +424,58 @@

-
359  {
-
360  // Warn for PBL models and moisture - these may not yet be compatible
-
361  for (int lev = 0; lev <= max_level; lev++) {
-
362  if ((moisture_type != MoistureType::None) && (turbChoice[lev].pbl_type != PBLType::None)) {
-
363  amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
-
364  }
-
365  }
-
366  //
-
367  // Buoyancy type check
-
368  //
-
369  if (buoyancy_type != 1 && buoyancy_type != 2 && buoyancy_type != 3 && buoyancy_type != 4) {
-
370  amrex::Abort("buoyancy_type must be 1, 2, 3 or 4");
-
371  }
-
372 
-
373  if (!use_lagged_delta_rt && !(terrain_type == TerrainType::Moving)) {
-
374  amrex::Error("Can't turn off lagged_delta_rt when terrain not moving");
-
375  }
-
376 
-
377  //
-
378  // Wind farm checks
-
379  //
-
380  if (windfarm_type==WindFarmType::SimpleAD and sampling_distance_by_D < 0.0) {
-
381  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
-
382  " erf.sampling_distance_by_D in the inputs which specifies the upstream"
-
383  " distance as a factor of the turbine diameter at which the incoming free stream"
-
384  " velocity will be computed at.");
-
385  }
-
386  if ( (windfarm_type==WindFarmType::SimpleAD ||
-
387  windfarm_type==WindFarmType::GeneralAD ) && turb_disk_angle < 0.0) {
-
388  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
-
389  " erf.turb_disk_angle_from_x in the inputs which is the angle of the face of the"
-
390  " turbine disk from the x-axis. A turbine facing an oncoming flow in the x-direction"
-
391  " will have turb_disk_angle value of 90 deg.");
-
392  }
-
393  if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
-
394  amrex::Abort("You are using windfarms with latitude-logitude option to position the turbines."
-
395  " For this you should provide the inputs erf.windfarm_x_shift and"
-
396  " erf.windfarm_y_shift which are the values by which the bounding box of the"
-
397  " windfarm is shifted from the x and the y axes.");
-
398  }
-
399  }
-
SolverChoice::use_lagged_delta_rt
bool use_lagged_delta_rt
Definition: ERF_DataStruct.H:568
-
SolverChoice::windfarm_x_shift
amrex::Real windfarm_x_shift
Definition: ERF_DataStruct.H:646
-
SolverChoice::windfarm_loc_type
WindFarmLocType windfarm_loc_type
Definition: ERF_DataStruct.H:620
-
SolverChoice::sampling_distance_by_D
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:644
-
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:535
-
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:618
-
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:524
-
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:619
-
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:552
-
SolverChoice::turb_disk_angle
amrex::Real turb_disk_angle
Definition: ERF_DataStruct.H:645
-
SolverChoice::windfarm_y_shift
amrex::Real windfarm_y_shift
Definition: ERF_DataStruct.H:647
+
364  {
+
365  // Warn for PBL models and moisture - these may not yet be compatible
+
366  for (int lev = 0; lev <= max_level; lev++) {
+
367  if ((moisture_type != MoistureType::None) && (turbChoice[lev].pbl_type != PBLType::None)) {
+
368  amrex::Warning("\n*** WARNING: Moisture may not yet be compatible with PBL models, \n proceed with caution ***");
+
369  }
+
370  }
+
371  //
+
372  // Buoyancy type check
+
373  //
+
374  if (buoyancy_type != 1 && buoyancy_type != 2 && buoyancy_type != 3 && buoyancy_type != 4) {
+
375  amrex::Abort("buoyancy_type must be 1, 2, 3 or 4");
+
376  }
+
377 
+
378  if (!use_lagged_delta_rt && !(terrain_type == TerrainType::Moving)) {
+
379  amrex::Error("Can't turn off lagged_delta_rt when terrain not moving");
+
380  }
+
381 
+
382  //
+
383  // Wind farm checks
+
384  //
+
385  if (windfarm_type==WindFarmType::SimpleAD and sampling_distance_by_D < 0.0) {
+
386  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
+
387  " erf.sampling_distance_by_D in the inputs which specifies the upstream"
+
388  " distance as a factor of the turbine diameter at which the incoming free stream"
+
389  " velocity will be computed at.");
+
390  }
+
391  if ( (windfarm_type==WindFarmType::SimpleAD ||
+
392  windfarm_type==WindFarmType::GeneralAD ) && turb_disk_angle < 0.0) {
+
393  amrex::Abort("To use simplified actuator disks, you need to provide a variable"
+
394  " erf.turb_disk_angle_from_x in the inputs which is the angle of the face of the"
+
395  " turbine disk from the x-axis. A turbine facing an oncoming flow in the x-direction"
+
396  " will have turb_disk_angle value of 90 deg.");
+
397  }
+
398  if (windfarm_loc_type == WindFarmLocType::lat_lon and (windfarm_x_shift < 0.0 or windfarm_y_shift < 0.0)) {
+
399  amrex::Abort("You are using windfarms with latitude-logitude option to position the turbines."
+
400  " For this you should provide the inputs erf.windfarm_x_shift and"
+
401  " erf.windfarm_y_shift which are the values by which the bounding box of the"
+
402  " windfarm is shifted from the x and the y axes.");
+
403  }
+
404  }
+
SolverChoice::use_lagged_delta_rt
bool use_lagged_delta_rt
Definition: ERF_DataStruct.H:573
+
SolverChoice::windfarm_x_shift
amrex::Real windfarm_x_shift
Definition: ERF_DataStruct.H:651
+
SolverChoice::windfarm_loc_type
WindFarmLocType windfarm_loc_type
Definition: ERF_DataStruct.H:625
+
SolverChoice::sampling_distance_by_D
amrex::Real sampling_distance_by_D
Definition: ERF_DataStruct.H:649
+
SolverChoice::turbChoice
amrex::Vector< TurbChoice > turbChoice
Definition: ERF_DataStruct.H:540
+
SolverChoice::moisture_type
MoistureType moisture_type
Definition: ERF_DataStruct.H:623
+
SolverChoice::terrain_type
static TerrainType terrain_type
Definition: ERF_DataStruct.H:529
+
SolverChoice::windfarm_type
WindFarmType windfarm_type
Definition: ERF_DataStruct.H:624
+
SolverChoice::buoyancy_type
int buoyancy_type
Definition: ERF_DataStruct.H:557
+
SolverChoice::turb_disk_angle
amrex::Real turb_disk_angle
Definition: ERF_DataStruct.H:650
+
SolverChoice::windfarm_y_shift
amrex::Real windfarm_y_shift
Definition: ERF_DataStruct.H:652

Referenced by init_params().

@@ -513,75 +513,75 @@

-
402  {
-
403  amrex::Print() << "SOLVER CHOICE: " << std::endl;
-
404  amrex::Print() << "force_stage1_single_substep : " << force_stage1_single_substep << std::endl;
-
405  for (int lev = 0; lev <= max_level; lev++) {
-
406  amrex::Print() << "anelastic at level : " << lev << " is " << anelastic[lev] << std::endl;
-
407  if (substepping_type[lev] == SubsteppingType::None) {
-
408  amrex::Print() << "No substepping at level " << lev << std::endl;
-
409  } else if (substepping_type[lev] == SubsteppingType::Explicit) {
-
410  amrex::Print() << "Explicit substepping at level " << lev << std::endl;
-
411  } else if (substepping_type[lev] == SubsteppingType::Implicit) {
-
412  amrex::Print() << "Implicit substepping at level " << lev << std::endl;
-
413  }
-
414  }
-
415  amrex::Print() << "use_coriolis : " << use_coriolis << std::endl;
-
416  amrex::Print() << "use_gravity : " << use_gravity << std::endl;
-
417 
-
418  amrex::Print() << "Terrain Type: " << std::endl;
-
419  if (terrain_type == TerrainType::Static) {
-
420  amrex::Print() << " Static" << std::endl;
-
421  } else if (terrain_type == TerrainType::Moving) {
-
422  amrex::Print() << " Moving" << std::endl;
-
423  } else {
-
424  amrex::Print() << " None" << std::endl;
-
425  }
-
426 
-
427  amrex::Print() << "ABL Driver Type: " << std::endl;
-
428  if (abl_driver_type == ABLDriverType::None) {
+
407  {
+
408  amrex::Print() << "SOLVER CHOICE: " << std::endl;
+
409  amrex::Print() << "force_stage1_single_substep : " << force_stage1_single_substep << std::endl;
+
410  for (int lev = 0; lev <= max_level; lev++) {
+
411  amrex::Print() << "anelastic at level : " << lev << " is " << anelastic[lev] << std::endl;
+
412  if (substepping_type[lev] == SubsteppingType::None) {
+
413  amrex::Print() << "No substepping at level " << lev << std::endl;
+
414  } else if (substepping_type[lev] == SubsteppingType::Explicit) {
+
415  amrex::Print() << "Explicit substepping at level " << lev << std::endl;
+
416  } else if (substepping_type[lev] == SubsteppingType::Implicit) {
+
417  amrex::Print() << "Implicit substepping at level " << lev << std::endl;
+
418  }
+
419  }
+
420  amrex::Print() << "use_coriolis : " << use_coriolis << std::endl;
+
421  amrex::Print() << "use_gravity : " << use_gravity << std::endl;
+
422 
+
423  amrex::Print() << "Terrain Type: " << std::endl;
+
424  if (terrain_type == TerrainType::Static) {
+
425  amrex::Print() << " Static" << std::endl;
+
426  } else if (terrain_type == TerrainType::Moving) {
+
427  amrex::Print() << " Moving" << std::endl;
+
428  } else {
429  amrex::Print() << " None" << std::endl;
-
430  } else if (abl_driver_type == ABLDriverType::PressureGradient) {
-
431  amrex::Print() << " Pressure Gradient "
-
432  << amrex::RealVect(abl_pressure_grad[0],abl_pressure_grad[1],abl_pressure_grad[2])
-
433  << std::endl;
-
434  } else if (abl_driver_type == ABLDriverType::GeostrophicWind) {
-
435  amrex::Print() << " Geostrophic Wind "
-
436  << amrex::RealVect(abl_geo_forcing[0],abl_geo_forcing[1],abl_geo_forcing[2])
-
437  << std::endl;
-
438  }
-
439 
-
440  if (max_level > 0) {
-
441  amrex::Print() << "Coupling Type: " << std::endl;
-
442  if (coupling_type == CouplingType::TwoWay) {
-
443  amrex::Print() << " Two-way" << std::endl;
-
444  } else if (coupling_type == CouplingType::OneWay) {
-
445  amrex::Print() << " One-way" << std::endl;
-
446  }
-
447  }
-
448 
-
449  amrex::Print() << "Buoyancy_type : " << buoyancy_type << std::endl;
-
450 
-
451  advChoice.display();
-
452  diffChoice.display();
-
453  spongeChoice.display();
-
454 
-
455  for (int lev = 0; lev <= max_level; lev++) {
-
456  turbChoice[lev].display(lev);
-
457  }
-
458  }
+
430  }
+
431 
+
432  amrex::Print() << "ABL Driver Type: " << std::endl;
+
433  if (abl_driver_type == ABLDriverType::None) {
+
434  amrex::Print() << " None" << std::endl;
+
435  } else if (abl_driver_type == ABLDriverType::PressureGradient) {
+
436  amrex::Print() << " Pressure Gradient "
+
437  << amrex::RealVect(abl_pressure_grad[0],abl_pressure_grad[1],abl_pressure_grad[2])
+
438  << std::endl;
+
439  } else if (abl_driver_type == ABLDriverType::GeostrophicWind) {
+
440  amrex::Print() << " Geostrophic Wind "
+
441  << amrex::RealVect(abl_geo_forcing[0],abl_geo_forcing[1],abl_geo_forcing[2])
+
442  << std::endl;
+
443  }
+
444 
+
445  if (max_level > 0) {
+
446  amrex::Print() << "Coupling Type: " << std::endl;
+
447  if (coupling_type == CouplingType::TwoWay) {
+
448  amrex::Print() << " Two-way" << std::endl;
+
449  } else if (coupling_type == CouplingType::OneWay) {
+
450  amrex::Print() << " One-way" << std::endl;
+
451  }
+
452  }
+
453 
+
454  amrex::Print() << "Buoyancy_type : " << buoyancy_type << std::endl;
+
455 
+
456  advChoice.display();
+
457  diffChoice.display();
+
458  spongeChoice.display();
+
459 
+
460  for (int lev = 0; lev <= max_level; lev++) {
+
461  turbChoice[lev].display(lev);
+
462  }
+
463  }
AdvChoice::display
void display()
Definition: ERF_AdvStruct.H:189
DiffChoice::display
void display()
Definition: ERF_DiffStruct.H:56
-
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:533
-
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:555
-
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:540
-
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:624
-
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:542
-
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:556
-
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:543
-
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:532
-
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:534
-
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:617
+
SolverChoice::diffChoice
DiffChoice diffChoice
Definition: ERF_DataStruct.H:538
+
SolverChoice::use_gravity
bool use_gravity
Definition: ERF_DataStruct.H:560
+
SolverChoice::force_stage1_single_substep
int force_stage1_single_substep
Definition: ERF_DataStruct.H:545
+
SolverChoice::abl_pressure_grad
amrex::GpuArray< amrex::Real, AMREX_SPACEDIM > abl_pressure_grad
Definition: ERF_DataStruct.H:629
+
SolverChoice::substepping_type
amrex::Vector< SubsteppingType > substepping_type
Definition: ERF_DataStruct.H:547
+
SolverChoice::use_coriolis
bool use_coriolis
Definition: ERF_DataStruct.H:561
+
SolverChoice::anelastic
amrex::Vector< int > anelastic
Definition: ERF_DataStruct.H:548
+
SolverChoice::advChoice
AdvChoice advChoice
Definition: ERF_DataStruct.H:537
+
SolverChoice::spongeChoice
SpongeChoice spongeChoice
Definition: ERF_DataStruct.H:539
+
SolverChoice::coupling_type
CouplingType coupling_type
Definition: ERF_DataStruct.H:622
SpongeChoice::display
void display()
Definition: ERF_SpongeStruct.H:45
Here is the call graph for this function:
@@ -686,217 +686,222 @@

143 

144  // Is the terrain none, static or moving?
145  pp.query_enum_case_insensitive("terrain_type",terrain_type);
-
146 
-
147  // Use lagged_delta_rt in the fast integrator?
-
148  pp.query("use_lagged_delta_rt", use_lagged_delta_rt);
-
149 
-
150  // These default to true but are used for unit testing
-
151  pp.query("use_gravity", use_gravity);
-
152  gravity = use_gravity? CONST_GRAV: 0.0;
-
153 
-
154  pp.query("c_p", c_p);
-
155  rdOcp = R_d / c_p;
-
156 
-
157  read_int_string(max_level, "anelastic", anelastic, 0);
+
146  int n_zlevels = pp.countval("terrain_z_levels");
+
147  if (n_zlevels > 0 and terrain_type == TerrainType::None)
+
148  {
+
149  terrain_type = TerrainType::Static;
+
150  }
+
151 
+
152  // Use lagged_delta_rt in the fast integrator?
+
153  pp.query("use_lagged_delta_rt", use_lagged_delta_rt);
+
154 
+
155  // These default to true but are used for unit testing
+
156  pp.query("use_gravity", use_gravity);
+
157  gravity = use_gravity? CONST_GRAV: 0.0;
158 
-
159  // *******************************************************************************
-
160  // Read substepping_type and allow for different values at each level
-
161  // *******************************************************************************
-
162  substepping_type.resize(max_level+1);
+
159  pp.query("c_p", c_p);
+
160  rdOcp = R_d / c_p;
+
161 
+
162  read_int_string(max_level, "anelastic", anelastic, 0);
163 
-
164  for (int i = 0; i <= max_level; i++) {
-
165  substepping_type[i] = SubsteppingType::Implicit;
-
166  }
-
167 
-
168  int nvals = pp.countval("substepping_type");
-
169  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
-
170 
-
171  if (nvals == 1) {
-
172  pp.query_enum_case_insensitive("substepping_type",substepping_type[0]);
-
173  for (int i = 1; i <= max_level; i++) {
-
174  substepping_type[i] = substepping_type[0];
-
175  }
-
176  } else if (nvals > 1) { // in this case we have asserted nvals >= max_level+1
-
177  for (int i = 0; i <= max_level; i++) {
-
178  pp.query_enum_case_insensitive("substepping_type",substepping_type[i],i);
-
179  }
-
180  }
-
181 
-
182  // *******************************************************************************
-
183  // Error check on deprecated input
-
184  // *******************************************************************************
-
185  int nvals_old = pp.countval("no_substepping");
-
186  if (nvals_old > 0) {
-
187  amrex::Abort("The no_substepping flag is deprecated -- set substepping_type instead");
-
188  }
+
164  // *******************************************************************************
+
165  // Read substepping_type and allow for different values at each level
+
166  // *******************************************************************************
+
167  substepping_type.resize(max_level+1);
+
168 
+
169  for (int i = 0; i <= max_level; i++) {
+
170  substepping_type[i] = SubsteppingType::Implicit;
+
171  }
+
172 
+
173  int nvals = pp.countval("substepping_type");
+
174  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
+
175 
+
176  if (nvals == 1) {
+
177  pp.query_enum_case_insensitive("substepping_type",substepping_type[0]);
+
178  for (int i = 1; i <= max_level; i++) {
+
179  substepping_type[i] = substepping_type[0];
+
180  }
+
181  } else if (nvals > 1) { // in this case we have asserted nvals >= max_level+1
+
182  for (int i = 0; i <= max_level; i++) {
+
183  pp.query_enum_case_insensitive("substepping_type",substepping_type[i],i);
+
184  }
+
185  }
+
186 
+
187  // *******************************************************************************
+
188  // Error check on deprecated input
189  // *******************************************************************************
-
190 
-
191  bool any_anelastic = false;
-
192  for (int i = 0; i <= max_level; ++i) {
-
193  if (anelastic[i] == 1) any_anelastic = true;
-
194  }
+
190  int nvals_old = pp.countval("no_substepping");
+
191  if (nvals_old > 0) {
+
192  amrex::Abort("The no_substepping flag is deprecated -- set substepping_type instead");
+
193  }
+
194  // *******************************************************************************
195 
-
196  // If anelastic at all, we do not advect rho -- it is always == rho0
-
197  if (any_anelastic == 1) {
-
198  constant_density = true;
-
199  project_initial_velocity = true;
-
200  buoyancy_type = 3; // (This isn't actually used when anelastic is set)
-
201  } else {
-
202  pp.query("project_initial_velocity", project_initial_velocity);
-
203 
-
204  constant_density = false; // We default to false but allow the user to set it
-
205  pp.query("constant_density", constant_density);
-
206  }
-
207 
-
208  // *******************************************************************************
-
209 
-
210  pp.query("ncorr", ncorr);
-
211  pp.query("poisson_abstol", poisson_abstol);
-
212  pp.query("poisson_reltol", poisson_reltol);
-
213 
-
214  for (int lev = 0; lev <= max_level; lev++) {
-
215  if (anelastic[lev] != 0)
-
216  {
-
217  substepping_type[lev] = SubsteppingType::None;
-
218  }
-
219  }
-
220 
-
221  pp.query("force_stage1_single_substep", force_stage1_single_substep);
-
222 
-
223  // Include Coriolis forcing?
-
224  pp.query("use_coriolis", use_coriolis);
+
196  bool any_anelastic = false;
+
197  for (int i = 0; i <= max_level; ++i) {
+
198  if (anelastic[i] == 1) any_anelastic = true;
+
199  }
+
200 
+
201  // If anelastic at all, we do not advect rho -- it is always == rho0
+
202  if (any_anelastic == 1) {
+
203  constant_density = true;
+
204  project_initial_velocity = true;
+
205  buoyancy_type = 3; // (This isn't actually used when anelastic is set)
+
206  } else {
+
207  pp.query("project_initial_velocity", project_initial_velocity);
+
208 
+
209  constant_density = false; // We default to false but allow the user to set it
+
210  pp.query("constant_density", constant_density);
+
211  }
+
212 
+
213  // *******************************************************************************
+
214 
+
215  pp.query("ncorr", ncorr);
+
216  pp.query("poisson_abstol", poisson_abstol);
+
217  pp.query("poisson_reltol", poisson_reltol);
+
218 
+
219  for (int lev = 0; lev <= max_level; lev++) {
+
220  if (anelastic[lev] != 0)
+
221  {
+
222  substepping_type[lev] = SubsteppingType::None;
+
223  }
+
224  }
225 
-
226  // Include Rayleigh damping (separate flags for each variable)
-
227  pp.query("rayleigh_damp_U", rayleigh_damp_U);
-
228  pp.query("rayleigh_damp_V", rayleigh_damp_V);
-
229  pp.query("rayleigh_damp_W", rayleigh_damp_W);
-
230  pp.query("rayleigh_damp_T", rayleigh_damp_T);
-
231  pp.query("rayleigh_dampcoef", rayleigh_dampcoef);
-
232  pp.query("rayleigh_zdamp", rayleigh_zdamp);
-
233 
-
234  // Flag to do explicit MOST formulation
-
235  pp.query("use_explicit_most",use_explicit_most);
-
236 
-
237  // Flag to do MOST rotations with terrain
-
238  pp.query("use_rotate_most",use_rotate_most);
-
239  if (use_rotate_most) {
-
240  AMREX_ASSERT_WITH_MESSAGE(terrain_type != TerrainType::None,"MOST stress rotations are only valid with terrain!");
-
241  AMREX_ASSERT_WITH_MESSAGE(use_explicit_most, "MOST Stress rotations are only valid with explicit MOST!");
-
242  }
-
243 
-
244  // Which external forcings?
-
245  abl_driver_type = ABLDriverType::None; // Default: no ABL driver for simulating classical fluid dynamics problems
-
246  pp.query_enum_case_insensitive("abl_driver_type",abl_driver_type);
-
247 
-
248  // Which type of inflow turbulent generation
-
249  pert_type = PerturbationType::None; // Default
-
250  pp.query_enum_case_insensitive("perturbation_type",pert_type);
-
251 
-
252  amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
-
253  pp.queryarr("abl_pressure_grad",abl_pressure_grad_in);
-
254  for(int i = 0; i < AMREX_SPACEDIM; ++i) abl_pressure_grad[i] = abl_pressure_grad_in[i];
-
255 
-
256  amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
-
257  if(pp.queryarr("abl_geo_forcing",abl_geo_forcing_in)) {
-
258  amrex::Print() << "Specified abl_geo_forcing: (";
-
259  for (int i = 0; i < AMREX_SPACEDIM; ++i) {
-
260  abl_geo_forcing[i] = abl_geo_forcing_in[i];
-
261  amrex::Print() << abl_geo_forcing[i] << " ";
-
262  }
-
263  amrex::Print() << ")" << std::endl;
-
264  }
-
265 
-
266  if (use_coriolis)
-
267  {
-
268  build_coriolis_forcings();
+
226  pp.query("force_stage1_single_substep", force_stage1_single_substep);
+
227 
+
228  // Include Coriolis forcing?
+
229  pp.query("use_coriolis", use_coriolis);
+
230 
+
231  // Include Rayleigh damping (separate flags for each variable)
+
232  pp.query("rayleigh_damp_U", rayleigh_damp_U);
+
233  pp.query("rayleigh_damp_V", rayleigh_damp_V);
+
234  pp.query("rayleigh_damp_W", rayleigh_damp_W);
+
235  pp.query("rayleigh_damp_T", rayleigh_damp_T);
+
236  pp.query("rayleigh_dampcoef", rayleigh_dampcoef);
+
237  pp.query("rayleigh_zdamp", rayleigh_zdamp);
+
238 
+
239  // Flag to do explicit MOST formulation
+
240  pp.query("use_explicit_most",use_explicit_most);
+
241 
+
242  // Flag to do MOST rotations with terrain
+
243  pp.query("use_rotate_most",use_rotate_most);
+
244  if (use_rotate_most) {
+
245  AMREX_ASSERT_WITH_MESSAGE(terrain_type != TerrainType::None,"MOST stress rotations are only valid with terrain!");
+
246  AMREX_ASSERT_WITH_MESSAGE(use_explicit_most, "MOST Stress rotations are only valid with explicit MOST!");
+
247  }
+
248 
+
249  // Which external forcings?
+
250  abl_driver_type = ABLDriverType::None; // Default: no ABL driver for simulating classical fluid dynamics problems
+
251  pp.query_enum_case_insensitive("abl_driver_type",abl_driver_type);
+
252 
+
253  // Which type of inflow turbulent generation
+
254  pert_type = PerturbationType::None; // Default
+
255  pp.query_enum_case_insensitive("perturbation_type",pert_type);
+
256 
+
257  amrex::Vector<amrex::Real> abl_pressure_grad_in = {0.0, 0.0, 0.0};
+
258  pp.queryarr("abl_pressure_grad",abl_pressure_grad_in);
+
259  for(int i = 0; i < AMREX_SPACEDIM; ++i) abl_pressure_grad[i] = abl_pressure_grad_in[i];
+
260 
+
261  amrex::Vector<amrex::Real> abl_geo_forcing_in = {0.0, 0.0, 0.0};
+
262  if(pp.queryarr("abl_geo_forcing",abl_geo_forcing_in)) {
+
263  amrex::Print() << "Specified abl_geo_forcing: (";
+
264  for (int i = 0; i < AMREX_SPACEDIM; ++i) {
+
265  abl_geo_forcing[i] = abl_geo_forcing_in[i];
+
266  amrex::Print() << abl_geo_forcing[i] << " ";
+
267  }
+
268  amrex::Print() << ")" << std::endl;
269  }
270 
-
271  pp.query("add_custom_rhotheta_forcing", custom_rhotheta_forcing);
-
272  pp.query("add_custom_moisture_forcing", custom_moisture_forcing);
-
273  pp.query("add_custom_w_subsidence", custom_w_subsidence);
-
274  pp.query("add_custom_geostrophic_profile", custom_geostrophic_profile);
-
275  pp.query("custom_forcing_uses_primitive_vars", custom_forcing_prim_vars);
-
276 
-
277  pp.query("nudging_from_input_sounding", nudging_from_input_sounding);
-
278 
-
279  have_geo_wind_profile = (!abl_geo_wind_table.empty() || custom_geostrophic_profile);
-
280  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(!(!abl_geo_wind_table.empty() && custom_geostrophic_profile),
-
281  "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
-
282 
-
283  pp.query("Ave_Plane", ave_plane);
-
284 
-
285  pp.query("mp_clouds", do_cloud);
-
286  pp.query("mp_precip", do_precip);
-
287  pp.query("use_moist_background", use_moist_background);
-
288 
-
289  // Use numerical diffusion?
-
290  pp.query("use_NumDiff",use_NumDiff);
-
291  if(use_NumDiff) {
-
292  pp.query("NumDiffCoeff",NumDiffCoeff);
-
293  AMREX_ASSERT_WITH_MESSAGE(( (NumDiffCoeff >= 0.) && (NumDiffCoeff <= 1.) ),
-
294  "Numerical diffusion coefficient must be between 0 & 1.");
-
295  NumDiffCoeff *= std::pow(2.0,-6);
-
296  }
-
297 
-
298  // Use monotonic advection?
-
299  pp.query("use_mono_adv",use_mono_adv);
-
300 
-
301  advChoice.init_params();
-
302  diffChoice.init_params();
-
303  spongeChoice.init_params();
-
304 
-
305  turbChoice.resize(max_level+1);
-
306  for (int lev = 0; lev <= max_level; lev++) {
-
307  turbChoice[lev].init_params(lev,max_level);
-
308  }
+
271  if (use_coriolis)
+
272  {
+
273  build_coriolis_forcings();
+
274  }
+
275 
+
276  pp.query("add_custom_rhotheta_forcing", custom_rhotheta_forcing);
+
277  pp.query("add_custom_moisture_forcing", custom_moisture_forcing);
+
278  pp.query("add_custom_w_subsidence", custom_w_subsidence);
+
279  pp.query("add_custom_geostrophic_profile", custom_geostrophic_profile);
+
280  pp.query("custom_forcing_uses_primitive_vars", custom_forcing_prim_vars);
+
281 
+
282  pp.query("nudging_from_input_sounding", nudging_from_input_sounding);
+
283 
+
284  have_geo_wind_profile = (!abl_geo_wind_table.empty() || custom_geostrophic_profile);
+
285  AMREX_ALWAYS_ASSERT_WITH_MESSAGE(!(!abl_geo_wind_table.empty() && custom_geostrophic_profile),
+
286  "Should not have both abl_geo_wind_table and custom_geostrophic_profile set.");
+
287 
+
288  pp.query("Ave_Plane", ave_plane);
+
289 
+
290  pp.query("mp_clouds", do_cloud);
+
291  pp.query("mp_precip", do_precip);
+
292  pp.query("use_moist_background", use_moist_background);
+
293 
+
294  // Use numerical diffusion?
+
295  pp.query("use_NumDiff",use_NumDiff);
+
296  if(use_NumDiff) {
+
297  pp.query("NumDiffCoeff",NumDiffCoeff);
+
298  AMREX_ASSERT_WITH_MESSAGE(( (NumDiffCoeff >= 0.) && (NumDiffCoeff <= 1.) ),
+
299  "Numerical diffusion coefficient must be between 0 & 1.");
+
300  NumDiffCoeff *= std::pow(2.0,-6);
+
301  }
+
302 
+
303  // Use monotonic advection?
+
304  pp.query("use_mono_adv",use_mono_adv);
+
305 
+
306  advChoice.init_params();
+
307  diffChoice.init_params();
+
308  spongeChoice.init_params();
309 
-
310  // YSU PBL: use consistent coriolis frequency
+
310  turbChoice.resize(max_level+1);
311  for (int lev = 0; lev <= max_level; lev++) {
-
312  if (turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
-
313  if (use_coriolis) {
-
314  turbChoice[lev].pbl_ysu_coriolis_freq = coriolis_factor * sinphi;
-
315  if (lev == 0) {
-
316  amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
-
317  }
-
318  } else {
-
319  amrex::Abort("YSU cannot use ERF coriolis frequency if not using coriolis");
-
320  }
-
321  }
-
322  }
-
323 
-
324 
-
325  // Which type of multilevel coupling
-
326  coupling_type = CouplingType::TwoWay; // Default
-
327  pp.query_enum_case_insensitive("coupling_type",coupling_type);
+
312  turbChoice[lev].init_params(lev,max_level);
+
313  }
+
314 
+
315  // YSU PBL: use consistent coriolis frequency
+
316  for (int lev = 0; lev <= max_level; lev++) {
+
317  if (turbChoice[lev].pbl_ysu_use_consistent_coriolis) {
+
318  if (use_coriolis) {
+
319  turbChoice[lev].pbl_ysu_coriolis_freq = coriolis_factor * sinphi;
+
320  if (lev == 0) {
+
321  amrex::Print() << "YSU PBL using ERF coriolis frequency: " << turbChoice[lev].pbl_ysu_coriolis_freq << std::endl;
+
322  }
+
323  } else {
+
324  amrex::Abort("YSU cannot use ERF coriolis frequency if not using coriolis");
+
325  }
+
326  }
+
327  }
328 
-
329  // Which type of windfarm model
-
330  windfarm_type = WindFarmType::None; // Default
-
331  pp.query_enum_case_insensitive("windfarm_type",windfarm_type);
-
332 
-
333  static std::string windfarm_loc_type_string = "None";
-
334  windfarm_loc_type = WindFarmLocType::None;
-
335  pp.query_enum_case_insensitive("windfarm_loc_type",windfarm_loc_type);
-
336 
-
337  pp.query("windfarm_loc_table", windfarm_loc_table);
-
338  pp.query("windfarm_spec_table", windfarm_spec_table);
-
339  pp.query("windfarm_blade_table", windfarm_blade_table);
-
340  pp.query("windfarm_airfoil_tables", windfarm_airfoil_tables);
-
341  pp.query("windfarm_spec_table_extra", windfarm_spec_table_extra);
-
342 
-
343  // Sampling distance upstream of the turbine to find the
-
344  // incoming free stream velocity as a factor of the diameter of the
-
345  // turbine. ie. the sampling distance will be this number multiplied
-
346  // by the diameter of the turbine
-
347  pp.query("sampling_distance_by_D", sampling_distance_by_D);
-
348  pp.query("turb_disk_angle_from_x", turb_disk_angle);
-
349 
-
350  pp.query("windfarm_x_shift",windfarm_x_shift);
-
351  pp.query("windfarm_y_shift",windfarm_y_shift);
-
352  // Test if time averaged data is to be output
-
353  pp.query("time_avg_vel",time_avg_vel);
+
329 
+
330  // Which type of multilevel coupling
+
331  coupling_type = CouplingType::TwoWay; // Default
+
332  pp.query_enum_case_insensitive("coupling_type",coupling_type);
+
333 
+
334  // Which type of windfarm model
+
335  windfarm_type = WindFarmType::None; // Default
+
336  pp.query_enum_case_insensitive("windfarm_type",windfarm_type);
+
337 
+
338  static std::string windfarm_loc_type_string = "None";
+
339  windfarm_loc_type = WindFarmLocType::None;
+
340  pp.query_enum_case_insensitive("windfarm_loc_type",windfarm_loc_type);
+
341 
+
342  pp.query("windfarm_loc_table", windfarm_loc_table);
+
343  pp.query("windfarm_spec_table", windfarm_spec_table);
+
344  pp.query("windfarm_blade_table", windfarm_blade_table);
+
345  pp.query("windfarm_airfoil_tables", windfarm_airfoil_tables);
+
346  pp.query("windfarm_spec_table_extra", windfarm_spec_table_extra);
+
347 
+
348  // Sampling distance upstream of the turbine to find the
+
349  // incoming free stream velocity as a factor of the diameter of the
+
350  // turbine. ie. the sampling distance will be this number multiplied
+
351  // by the diameter of the turbine
+
352  pp.query("sampling_distance_by_D", sampling_distance_by_D);
+
353  pp.query("turb_disk_angle_from_x", turb_disk_angle);
354 
-
355  check_params(max_level);
-
356  }
+
355  pp.query("windfarm_x_shift",windfarm_x_shift);
+
356  pp.query("windfarm_y_shift",windfarm_y_shift);
+
357  // Test if time averaged data is to be output
+
358  pp.query("time_avg_vel",time_avg_vel);
+
359 
+
360  check_params(max_level);
+
361  }
CONST_GRAV
constexpr amrex::Real CONST_GRAV
Definition: ERF_Constants.H:21
R_d
constexpr amrex::Real R_d
Definition: ERF_Constants.H:10
RhoQ4_comp
#define RhoQ4_comp
Definition: ERF_IndexDefines.H:45
@@ -904,53 +909,53 @@

RhoQ1_comp
#define RhoQ1_comp
Definition: ERF_IndexDefines.H:42

AdvChoice::init_params
void init_params()
Definition: ERF_AdvStruct.H:21
DiffChoice::init_params
void init_params()
Definition: ERF_DiffStruct.H:21
-
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:562
-
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:578
-
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:599
-
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:560
-
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:640
-
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:548
-
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:564
-
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:573
-
SolverChoice::read_int_string
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:503
-
SolverChoice::windfarm_spec_table
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:642
-
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:615
-
SolverChoice::build_coriolis_forcings
void build_coriolis_forcings()
Definition: ERF_DataStruct.H:460
-
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:546
-
SolverChoice::windfarm_spec_table_extra
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:642
-
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:621
-
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:572
-
SolverChoice::windfarm_loc_table
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:642
-
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:571
-
SolverChoice::check_params
void check_params(int max_level)
Definition: ERF_DataStruct.H:358
-
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:588
-
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:590
-
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:596
-
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:559
-
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:591
-
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:576
-
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:627
-
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:612
-
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:632
-
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:634
-
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:550
-
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:611
-
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:589
-
SolverChoice::windfarm_blade_table
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:643
-
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:577
-
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:580
-
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:545
-
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:633
-
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:592
-
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:561
-
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:608
-
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:547
-
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:605
-
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:631
-
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:602
-
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:563
-
SolverChoice::windfarm_airfoil_tables
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:643
-
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:629
+
SolverChoice::rayleigh_damp_T
bool rayleigh_damp_T
Definition: ERF_DataStruct.H:567
+
SolverChoice::dz0
amrex::Real dz0
Definition: ERF_DataStruct.H:583
+
SolverChoice::use_explicit_most
bool use_explicit_most
Definition: ERF_DataStruct.H:604
+
SolverChoice::rayleigh_damp_V
bool rayleigh_damp_V
Definition: ERF_DataStruct.H:565
+
SolverChoice::RhoQr_comp
int RhoQr_comp
Definition: ERF_DataStruct.H:645
+
SolverChoice::poisson_reltol
amrex::Real poisson_reltol
Definition: ERF_DataStruct.H:553
+
SolverChoice::rayleigh_zdamp
amrex::Real rayleigh_zdamp
Definition: ERF_DataStruct.H:569
+
SolverChoice::rdOcp
amrex::Real rdOcp
Definition: ERF_DataStruct.H:578
+
SolverChoice::read_int_string
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:508
+
SolverChoice::windfarm_spec_table
std::string windfarm_spec_table
Definition: ERF_DataStruct.H:647
+
SolverChoice::use_mono_adv
bool use_mono_adv
Definition: ERF_DataStruct.H:620
+
SolverChoice::build_coriolis_forcings
void build_coriolis_forcings()
Definition: ERF_DataStruct.H:465
+
SolverChoice::ncorr
int ncorr
Definition: ERF_DataStruct.H:551
+
SolverChoice::windfarm_spec_table_extra
std::string windfarm_spec_table_extra
Definition: ERF_DataStruct.H:647
+
SolverChoice::lsm_type
LandSurfaceType lsm_type
Definition: ERF_DataStruct.H:626
+
SolverChoice::c_p
amrex::Real c_p
Definition: ERF_DataStruct.H:577
+
SolverChoice::windfarm_loc_table
std::string windfarm_loc_table
Definition: ERF_DataStruct.H:647
+
SolverChoice::gravity
amrex::Real gravity
Definition: ERF_DataStruct.H:576
+
SolverChoice::check_params
void check_params(int max_level)
Definition: ERF_DataStruct.H:363
+
SolverChoice::custom_rhotheta_forcing
bool custom_rhotheta_forcing
Definition: ERF_DataStruct.H:593
+
SolverChoice::custom_w_subsidence
bool custom_w_subsidence
Definition: ERF_DataStruct.H:595
+
SolverChoice::nudging_from_input_sounding
bool nudging_from_input_sounding
Definition: ERF_DataStruct.H:601
+
SolverChoice::rayleigh_damp_U
bool rayleigh_damp_U
Definition: ERF_DataStruct.H:564
+
SolverChoice::custom_geostrophic_profile
bool custom_geostrophic_profile
Definition: ERF_DataStruct.H:596
+
SolverChoice::grid_stretching_ratio
amrex::Real grid_stretching_ratio
Definition: ERF_DataStruct.H:581
+
SolverChoice::have_geo_wind_profile
bool have_geo_wind_profile
Definition: ERF_DataStruct.H:632
+
SolverChoice::NumDiffCoeff
amrex::Real NumDiffCoeff
Definition: ERF_DataStruct.H:617
+
SolverChoice::do_precip
bool do_precip
Definition: ERF_DataStruct.H:637
+
SolverChoice::RhoQv_comp
int RhoQv_comp
Definition: ERF_DataStruct.H:639
+
SolverChoice::test_mapfactor
bool test_mapfactor
Definition: ERF_DataStruct.H:555
+
SolverChoice::use_NumDiff
bool use_NumDiff
Definition: ERF_DataStruct.H:616
+
SolverChoice::custom_moisture_forcing
bool custom_moisture_forcing
Definition: ERF_DataStruct.H:594
+
SolverChoice::windfarm_blade_table
std::string windfarm_blade_table
Definition: ERF_DataStruct.H:648
+
SolverChoice::zsurf
amrex::Real zsurf
Definition: ERF_DataStruct.H:582
+
SolverChoice::project_initial_velocity
bool project_initial_velocity
Definition: ERF_DataStruct.H:585
+
SolverChoice::constant_density
int constant_density
Definition: ERF_DataStruct.H:550
+
SolverChoice::use_moist_background
bool use_moist_background
Definition: ERF_DataStruct.H:638
+
SolverChoice::custom_forcing_prim_vars
bool custom_forcing_prim_vars
Definition: ERF_DataStruct.H:597
+
SolverChoice::rayleigh_damp_W
bool rayleigh_damp_W
Definition: ERF_DataStruct.H:566
+
SolverChoice::pert_type
PerturbationType pert_type
Definition: ERF_DataStruct.H:613
+
SolverChoice::poisson_abstol
amrex::Real poisson_abstol
Definition: ERF_DataStruct.H:552
+
SolverChoice::time_avg_vel
bool time_avg_vel
Definition: ERF_DataStruct.H:610
+
SolverChoice::do_cloud
bool do_cloud
Definition: ERF_DataStruct.H:636
+
SolverChoice::use_rotate_most
bool use_rotate_most
Definition: ERF_DataStruct.H:607
+
SolverChoice::rayleigh_dampcoef
amrex::Real rayleigh_dampcoef
Definition: ERF_DataStruct.H:568
+
SolverChoice::windfarm_airfoil_tables
std::string windfarm_airfoil_tables
Definition: ERF_DataStruct.H:648
+
SolverChoice::ave_plane
int ave_plane
Definition: ERF_DataStruct.H:634
SpongeChoice::init_params
void init_params()
Definition: ERF_SpongeStruct.H:17
Here is the call graph for this function:
@@ -1016,20 +1021,20 @@

-
505  {
-
506  amrex::ParmParse pp("erf");
-
507  int nvals = pp.countval(string_to_read);
-
508  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
-
509  amrex::Vector<int> temp; temp.resize(nvals);
-
510  pp.queryarr(string_to_read,temp);
-
511  if (nvals == 0) {
-
512  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(default_int);
-
513  } else if (nvals == 1) {
-
514  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[0]);
-
515  } else {
-
516  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[i]);
-
517  }
-
518  }
+
510  {
+
511  amrex::ParmParse pp("erf");
+
512  int nvals = pp.countval(string_to_read);
+
513  AMREX_ALWAYS_ASSERT(nvals == 0 || nvals == 1 || nvals >= max_level+1);
+
514  amrex::Vector<int> temp; temp.resize(nvals);
+
515  pp.queryarr(string_to_read,temp);
+
516  if (nvals == 0) {
+
517  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(default_int);
+
518  } else if (nvals == 1) {
+
519  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[0]);
+
520  } else {
+
521  for (int i = 0; i <= max_level; ++i) vec_to_fill.push_back(temp[i]);
+
522  }
+
523  }

Referenced by init_params().

@@ -1075,10 +1080,10 @@

-
528  {
-
529  terrain_is_flat = true;
-
530  }
-
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:521
+
533  {
+
534  terrain_is_flat = true;
+
535  }
+
SolverChoice::terrain_is_flat
static bool terrain_is_flat
Definition: ERF_DataStruct.H:526

Referenced by ProblemBase::init_custom_terrain().