ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Macros Pages
ERF_TurbStruct.H
Go to the documentation of this file.
1 #ifndef ERF_TURB_STRUCT_H_
2 #define ERF_TURB_STRUCT_H_
3 
4 #include <ERF_MYNNStruct.H>
5 
6 AMREX_ENUM(LESType, None, Smagorinsky, Smagorinsky2D, Deardorff);
7 
8 AMREX_ENUM(RANSType, None, kEqn);
9 
10 AMREX_ENUM(PBLType, None, MYNN25, MYNNEDMF, YSU, MRF);
11 
12 AMREX_ENUM(StratType, theta, thetav, thetal);
13 
14 template <typename T>
15 void
17  const amrex::ParmParse& pp,
18  const char* query_string,
19  T& query_var,
20  const int lev,
21  const int maxlev)
22 {
23  int count = pp.countval(query_string);
24  if (count == 0) {
25  return; // nothing to do
26  } else if (count == 1) {
27  pp.query(query_string, query_var);
28  } else if (count == maxlev + 1) {
29  pp.query(query_string, query_var, lev);
30  } else {
31  amrex::Error(
32  "For parmparse variable " + pp.prefixedName(query_string) +
33  ": if specified, specify once total or once for each level");
34  }
35 }
36 
37 /**
38  * Container holding quantities related to turbulence parametrizations
39  */
40 struct TurbChoice
41 {
42 public:
43  void init_params(int lev, int max_level, std::string pp_prefix)
44  {
45  amrex::ParmParse pp(pp_prefix);
46 
47  // Which LES closure?
48  std::string les_type_string = "None";
49  query_one_or_per_level(pp, "les_type", les_type, lev, max_level);
50 
51  // Handle 2-D Smag
52  if (les_type == LESType::Smagorinsky2D) {
53  les_type = LESType::Smagorinsky;
54  smag2d = true;
55  }
56 
57  // Which RANS closure?
58  std::string rans_type_string = "None";
59  query_one_or_per_level(pp, "rans_type", rans_type, lev, max_level);
60 
61  if ((rans_type != RANSType::None) && (les_type != LESType::None)) {
62  amrex::Error("Hybrid RANS-LES not implemented");
63  }
64 
65  // Which PBL Closure
66  static std::string pbl_type_string = "None";
67  query_one_or_per_level(pp, "pbl_type", pbl_type, lev, max_level);
68 
69  // Do some more stuff for PBL Modeling
70  if (pbl_type != PBLType::None) {
71  // Check for compatibility between PBL, LES, Molec Transport
72  if (les_type != LESType::None) {
73  amrex::Print() << "Selected a PBL model and an LES model: "
74  << "Using PBL for vertical transport, LES for horizontal"
75  << std::endl;
76  }
77  if (les_type == LESType::Smagorinsky) {
78  if (!smag2d)
79  amrex::Error("If using Smagorinsky with a PBL model, the 2-D "
80  "formulation should be used");
81  } else if (les_type == LESType::Deardorff) {
82  amrex::Error(
83  "It is not recommended to use Deardorff LES and a PBL model");
84  }
85 
86  if (pbl_type == PBLType::MYNN25 || pbl_type == PBLType::MYNNEDMF) {
87  query_one_or_per_level(pp, "pbl_mynn_A1", pbl_mynn.A1, lev, max_level);
88  query_one_or_per_level(pp, "pbl_mynn_A2", pbl_mynn.A2, lev, max_level);
89  query_one_or_per_level(pp, "pbl_mynn_B1", pbl_mynn.B1, lev, max_level);
90  query_one_or_per_level(pp, "pbl_mynn_B2", pbl_mynn.B2, lev, max_level);
91  query_one_or_per_level(pp, "pbl_mynn_C1", pbl_mynn.C1, lev, max_level);
92  query_one_or_per_level(pp, "pbl_mynn_C2", pbl_mynn.C2, lev, max_level);
93  query_one_or_per_level(pp, "pbl_mynn_C3", pbl_mynn.C3, lev, max_level);
94  query_one_or_per_level(pp, "pbl_mynn_C4", pbl_mynn.C4, lev, max_level);
95  query_one_or_per_level(pp, "pbl_mynn_C5", pbl_mynn.C5, lev, max_level);
100  pp, "pbl_mynn_diffuse_moistvars", pbl_mynn.diffuse_moistvars, lev,
101  max_level);
103  pp, "pbl_mynn_SMmin", pbl_mynn.SMmin, lev, max_level);
105  pp, "pbl_mynn_SMmax", pbl_mynn.SMmax, lev, max_level);
107  pp, "pbl_mynn_SHmin", pbl_mynn.SHmin, lev, max_level);
109  pp, "pbl_mynn_SHmax", pbl_mynn.SHmax, lev, max_level);
111  pp, "pbl_mynn_SQfactor", pbl_mynn.SQfac, lev, max_level);
112  } else if (pbl_type == PBLType::YSU) {
114  pp, "pbl_ysu_coriolis_freq", pbl_ysu_coriolis_freq, lev, max_level);
116  pp, "pbl_ysu_use_consistent_coriolis",
117  pbl_ysu_use_consistent_coriolis, lev, max_level);
119  pp, "pbl_ysu_force_over_water", pbl_ysu_force_over_water, lev,
120  max_level);
122  pp, "pbl_ysu_land_Ribcr", pbl_ysu_land_Ribcr, lev, max_level);
124  pp, "pbl_ysu_unst_Ribcr", pbl_ysu_unst_Ribcr, lev, max_level);
125  } else if (pbl_type == PBLType::MRF) {
127  pp, "pbl_mrf_coriolis_freq", pbl_mrf_coriolis_freq, lev, max_level);
129  pp, "pbl_mrf_Ribcr", pbl_mrf_Ribcr, lev, max_level);
131  pp, "pbl_mrf_const_b", pbl_mrf_const_b, lev, max_level);
132  query_one_or_per_level(pp, "pbl_mrf_sf", pbl_mrf_sf, lev, max_level);
134  pp, "mrf_moistvars", mrf_moistvars, lev, max_level);
135  }
136  }
137 
138  // Right now, solving the QKE equation is only supported when MYNN PBL is
139  // turned on
140  if (pbl_type == PBLType::MYNN25 || pbl_type == PBLType::MYNNEDMF) {
141  use_pbl_tke = true;
142 
143  query_one_or_per_level(pp, "advect_tke", advect_tke, lev, max_level);
145  pp, "diffuse_tke_3D", diffuse_tke_3D, lev, max_level);
146  }
147 
148  // LES constants...
149  query_one_or_per_level(pp, "Cs", Cs, lev, max_level);
150 
151  query_one_or_per_level(pp, "Pr_t", Pr_t, lev, max_level);
152  query_one_or_per_level(pp, "Sc_t", Sc_t, lev, max_level);
153 
154  // Compute relevant forms of diffusion parameters
155  Pr_t_inv = amrex::Real(1.0) / Pr_t;
156  Sc_t_inv = amrex::Real(1.0) / Sc_t;
157 
158  if (les_type == LESType::Deardorff) {
159  query_one_or_per_level(pp, "Ck", Ck, lev, max_level);
160  query_one_or_per_level(pp, "Ce", Ce, lev, max_level);
161  query_one_or_per_level(pp, "Ce_wall", Ce_wall, lev, max_level);
162  }
163 
164  // To quantify atmospheric stability for subgrid modeling
165  query_one_or_per_level(pp, "thermal_stratification", strat_type, lev, max_level);
166  if (strat_type == StratType::theta) {
167  amrex::Print() << "Thermal stratification based on gradient of potential temperature" << std::endl;
168  } else if (strat_type == StratType::thetav) {
169  amrex::Print() << "Thermal stratification based on gradient of virtual potential temperature" << std::endl;
170  } else if (strat_type == StratType::thetal) {
171  amrex::Print() << "Thermal stratification based on gradient of linearized liquid-water potential temperature" << std::endl;
172  }
173 
174  // k-eqn constants
175  query_one_or_per_level(pp, "Cmu0", Cmu0, lev, max_level);
176  query_one_or_per_level(pp, "Cb", Cb, lev, max_level);
177  query_one_or_per_level(pp, "Rt_crit", Rt_crit, lev, max_level);
178  query_one_or_per_level(pp, "Rt_min", Rt_min, lev, max_level);
179  query_one_or_per_level(pp, "max_geom_lscale", l_g_max, lev, max_level);
180 
181  // Common inputs (LES or RANS)
182  query_one_or_per_level(pp, "sigma_k", sigma_k, lev, max_level);
183  query_one_or_per_level(pp, "theta_ref", theta_ref, lev, max_level);
184 
185  query_one_or_per_level(pp, "mix_isotropic", mix_isotropic, lev, max_level);
186  query_one_or_per_level(pp, "use_smag_stratification", use_smag_stratification, lev, max_level);
187 
188  // Set common flags
189  use_kturb =
190  ((les_type != LESType::None) || (rans_type != RANSType::None) ||
191  (pbl_type != PBLType::None));
192  use_keqn =
193  ((les_type == LESType::Deardorff) || (rans_type == RANSType::kEqn));
194  use_tke =
195  ((les_type == LESType::Deardorff) || (rans_type == RANSType::kEqn) ||
196  (pbl_type == PBLType::MYNN25) || (pbl_type == PBLType::MYNNEDMF));
197 
198  // Validate inputs
199  if (les_type == LESType::Smagorinsky) {
200  if (Cs == 0) {
201  amrex::Error("Need to specify Cs for Smagorsinky LES");
202  }
203  if (smag2d && mix_isotropic) {
204  amrex::Print() << "Turning off mix_isotropic for 2-D Smagorinsky" << std::endl;
205  mix_isotropic = false;
206  }
207  }
208  }
209 
210  void display(int lev)
211  {
212  amrex::Print() << "Turbulence Settings at level " << lev << std::endl;
213 
214  if (
215  les_type == LESType::None && rans_type == RANSType::None &&
216  pbl_type == PBLType::None) {
217  amrex::Print() << " Using DNS model at level " << lev << std::endl;
218  } else if (les_type == LESType::Smagorinsky) {
219  if (smag2d) {
220  amrex::Print() << " Using 2D Smagorinsky LES model at level " << lev << std::endl;
221  } else {
222  amrex::Print() << " Using Smagorinsky LES model at level " << lev << std::endl;
223  }
225  amrex::Print() << " Smagorinsky accounts for moist stratification" << std::endl;
226  }
227  } else if (les_type == LESType::Deardorff) {
228  amrex::Print() << " Using Deardorff LES model at level " << lev << std::endl;
229  } else if (rans_type == RANSType::kEqn) {
230  amrex::Print()
231  << " Using Axell & Liungman one-equation RANS k model at level " << lev << std::endl;
232  } else if (pbl_type == PBLType::MYNN25) {
233  amrex::Print() << " Using MYNN2.5 PBL model at level " << lev << std::endl;
234  } else if (pbl_type == PBLType::MYNNEDMF) {
235  amrex::Print() << " Using MYNNEDMF PBL model at level " << lev << std::endl;
236  } else if (pbl_type == PBLType::YSU) {
237  amrex::Print() << " Using YSU PBL model at level " << lev << std::endl;
238  } else if (pbl_type == PBLType::MRF) {
239  amrex::Print() << " Using MRF PBL model at level " << lev << std::endl;
240  } else {
241  amrex::Error("Unknown turbulence model");
242  }
243 
244  if (les_type != LESType::None) {
245  if (les_type == LESType::Smagorinsky) {
246  amrex::Print() << " Cs : " << Cs << std::endl;
247  }
248  if (les_type == LESType::Deardorff) {
249  amrex::Print() << " Ce : " << Ce << std::endl;
250  amrex::Print() << " Ce at wall : " << Ce_wall << std::endl;
251  amrex::Print() << " Ck : " << Ck << std::endl;
252  amrex::Print() << " sigma_k : " << sigma_k << std::endl;
253 
254  // Sullivan et al 1994, Eqn 14
255  amrex::Real Cs_equiv = std::sqrt(Ck * std::sqrt(Ck / Ce));
256  amrex::Print() << " equivalent Cs : " << Cs_equiv
257  << std::endl;
258  }
259  amrex::Print() << " isotropic mixing : " << mix_isotropic
260  << std::endl;
261  }
262 
263  if (rans_type != RANSType::None) {
264  if (rans_type == RANSType::kEqn) {
265  amrex::Print() << "Cmu0 : " << Cmu0 << std::endl;
266  amrex::Print() << "sigma_k : " << sigma_k << std::endl;
267  amrex::Print() << "Cb : " << Cb << std::endl;
268  amrex::Print() << "Rt_crit : " << Rt_crit << std::endl;
269  amrex::Print() << "Rt_min : " << Rt_min << std::endl;
270  amrex::Print() << "max_geom_lscale : " << l_g_max << std::endl;
271  }
272  }
273 
274  if ((les_type == LESType::Deardorff) ||
275  (les_type == LESType::Smagorinsky && use_smag_stratification) ||
276  (rans_type == RANSType::kEqn)) {
277  if (theta_ref > 0) {
278  amrex::Print() << " reference theta : " << theta_ref << std::endl;
279  } else {
280  amrex::Print() << " reference theta : n/a" << std::endl;
281  }
282  }
283 
284  if ((les_type != LESType::None) || (rans_type != RANSType::None)) {
285  amrex::Print() << " Pr_t : " << Pr_t << std::endl;
286  amrex::Print() << " Sc_t : " << Sc_t << std::endl;
287  }
288 
289  if (pbl_type == PBLType::MYNN25 || pbl_type == PBLType::MYNNEDMF) {
290  amrex::Print() << " pbl_mynn_A1 : " << pbl_mynn.A1 << std::endl;
291  amrex::Print() << " pbl_mynn_A2 : " << pbl_mynn.A2 << std::endl;
292  amrex::Print() << " pbl_mynn_B1 : " << pbl_mynn.B1 << std::endl;
293  amrex::Print() << " pbl_mynn_B2 : " << pbl_mynn.B2 << std::endl;
294  amrex::Print() << " pbl_mynn_C1 : " << pbl_mynn.C1 << std::endl;
295  amrex::Print() << " pbl_mynn_C2 : " << pbl_mynn.C2 << std::endl;
296  amrex::Print() << " pbl_mynn_C3 : " << pbl_mynn.C3 << std::endl;
297  amrex::Print() << " pbl_mynn_C4 : " << pbl_mynn.C4 << std::endl;
298  amrex::Print() << " pbl_mynn_C5 : " << pbl_mynn.C5 << std::endl;
299  } else if (pbl_type == PBLType::YSU) {
300  amrex::Print() << " pbl_ysu_coriolis_freq : "
301  << pbl_ysu_coriolis_freq << std::endl;
302  amrex::Print() << " pbl_ysu_use_consistent_coriolis : "
303  << pbl_ysu_use_consistent_coriolis << std::endl;
304  amrex::Print() << " pbl_ysu_force_over_water : "
305  << pbl_ysu_force_over_water << std::endl;
306  amrex::Print() << " pbl_ysu_land_Ribcr : "
307  << pbl_ysu_land_Ribcr << std::endl;
308  amrex::Print() << " pbl_ysu_unst_Ribcr : "
309  << pbl_ysu_unst_Ribcr << std::endl;
310  } else if (pbl_type == PBLType::MRF) {
311  amrex::Print() << " pbl_mrf_coriolis_freq : " << pbl_mrf_coriolis_freq
312  << std::endl;
313  amrex::Print() << " pbl_mrf_Ribcr : " << pbl_mrf_Ribcr
314  << std::endl;
315  amrex::Print() << " pbl_mrf_const_b : " << pbl_mrf_const_b
316  << std::endl;
317  amrex::Print() << " pbl_mrf_sf : " << pbl_mrf_sf
318  << std::endl;
319  amrex::Print() << " mrf_moistvars : " << mrf_moistvars
320  << std::endl;
321  }
322  }
323 
324  // LES model
325  LESType les_type;
326 
327  // Turbulent Prandtl number
328  amrex::Real Pr_t = amrex::Real(1.0) / amrex::Real(3.0);
329  amrex::Real Pr_t_inv = amrex::Real(3.0);
330 
331  // Turbulent Schmidt number
332  amrex::Real Sc_t = 1.0;
333  amrex::Real Sc_t_inv = 1.0;
334 
335  // Smagorinsky
336  amrex::Real Cs = 0.0;
337  bool smag2d = false;
338 
339  // Deardorff
340  amrex::Real Ce = 0.93;
341  amrex::Real Ce_wall = 0.0; // if > 0, then set Ce to this at k=0
342  amrex::Real Ck = 0.1;
343 
344  // k-eqn RANS coefficients (Axell & Liungman 2001)
345  amrex::Real Cmu0 = 0.5562;
346  amrex::Real Cb = 0.35;
347  amrex::Real Rt_crit = -1.0;
348  amrex::Real Rt_min = -3.0;
349  amrex::Real l_g_max = 1e99;
350 
351  // Deardorff or k-eqn RANS
352  // - diffusivity of tke is 1/sigma_k times the eddy viscosity
353  amrex::Real sigma_k = 0.5;
354  // - reference potential temperature used to quantify the stratification in
355  // a stable region
356  amrex::Real theta_ref = 0.0;
357  // - how to quantify stratification effects
359 
360  // Anistropic length scales
361  bool mix_isotropic = true;
362 
363  // Use the stratification adjustment for the Smagorinsky Scheme
365 
366  // RANS type
367  RANSType rans_type;
368 
369  // PBL model
370  PBLType pbl_type;
371 
373  MYNNLevel2 pbl_mynn_level2; // for filtering
374 
375  // Common Flags
376  bool use_kturb = false; // Any turbulence modeling?
377  bool use_keqn =
378  false; // Any microscale turbulence modeling (LES, RANS) with TKE closure?
379  // Then need to populate SmnSmn_lev for production term.
380  bool use_pbl_tke =
381  false; // Any mesoscale turbulence modeling (PBL) with TKE closure?
382  bool use_tke = false; // Any TKE closure (meso or microscale)?
383 
384  // Model coefficients - YSU
385  // TODO: Add parmparse for all of these above
386  amrex::Real pbl_ysu_coriolis_freq =
387  1.0e-4; // 1e-4 is hardcoded in WRF, we let the user specify or tske the
388  // value from ERF coriolis forcing
390  false; // ignore input pbl_ysu_coriolis_freq, take value from ERF coriolis
391  // forcing instead
393  false; // Force YSU to act as if it is over water regardless of other inputs
394  // (for testing)
395  amrex::Real pbl_ysu_land_Ribcr =
396  0.25; // Critical Bulk Richardson number of Land for stable conditions
397  amrex::Real pbl_ysu_unst_Ribcr =
398  0.0; // Critical Bulk Richardson number for unstable conditions
399  // Model coefficients - MRF
400  amrex::Real pbl_mrf_coriolis_freq = 1.0e-4;
401  amrex::Real pbl_mrf_Ribcr =
402  0.5; // Critical Bulk Richardson number for MRF PBL model
403  amrex::Real pbl_mrf_const_b =
404  7.8; // Constant b in MRF PBL model, used to compute the PBL height
405  amrex::Real pbl_mrf_sf =
406  0.1; // MRF surface flux, used to compute the PBL height
407  bool mrf_moistvars = false; // if true, adds turbulence to moisture
408  // MYNN2.5 PBL model
409  // TKE/QKE stuff
410  bool advect_tke = true; // if MYNN2.5 PBL is used default is turb transport in
411  // Z-direction only
412  bool diffuse_tke_3D = true; // if numerical diffusion is turned on
413 };
414 #endif
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real pp(amrex::Real y)
Definition: ERF_MicrophysicsUtils.H:230
AMREX_ENUM(LESType, None, Smagorinsky, Smagorinsky2D, Deardorff)
void query_one_or_per_level(const amrex::ParmParse &pp, const char *query_string, T &query_var, const int lev, const int maxlev)
Definition: ERF_TurbStruct.H:16
@ theta
Definition: ERF_MM5.H:20
@ T
Definition: ERF_IndexDefines.H:110
Definition: ERF_MYNNStruct.H:9
amrex::Real SMmax
Definition: ERF_MYNNStruct.H:56
amrex::Real SHmax
Definition: ERF_MYNNStruct.H:58
amrex::Real SQfac
Definition: ERF_MYNNStruct.H:52
amrex::Real C4
Definition: ERF_MYNNStruct.H:48
amrex::Real C1
Definition: ERF_MYNNStruct.H:45
amrex::Real C3
Definition: ERF_MYNNStruct.H:47
amrex::Real C2
Definition: ERF_MYNNStruct.H:46
amrex::Real A2
Definition: ERF_MYNNStruct.H:42
amrex::Real SHmin
Definition: ERF_MYNNStruct.H:57
amrex::Real B1
Definition: ERF_MYNNStruct.H:43
amrex::Real B2
Definition: ERF_MYNNStruct.H:44
amrex::Real C5
Definition: ERF_MYNNStruct.H:49
amrex::Real SMmin
Definition: ERF_MYNNStruct.H:55
amrex::Real A1
Definition: ERF_MYNNStruct.H:41
bool diffuse_moistvars
Definition: ERF_MYNNStruct.H:63
Definition: ERF_MYNNStruct.H:66
void init_coeffs(amrex::Real A1_lvl25, amrex::Real A2_lvl25, amrex::Real B1, amrex::Real B2, amrex::Real C1, amrex::Real C2, amrex::Real C3, amrex::Real, amrex::Real C5)
Definition: ERF_MYNNStruct.H:67
Definition: ERF_TurbStruct.H:41
bool smag2d
Definition: ERF_TurbStruct.H:337
amrex::Real sigma_k
Definition: ERF_TurbStruct.H:353
MYNNLevel25 pbl_mynn
Definition: ERF_TurbStruct.H:372
PBLType pbl_type
Definition: ERF_TurbStruct.H:370
amrex::Real pbl_mrf_const_b
Definition: ERF_TurbStruct.H:403
bool use_pbl_tke
Definition: ERF_TurbStruct.H:380
amrex::Real Sc_t_inv
Definition: ERF_TurbStruct.H:333
amrex::Real Rt_min
Definition: ERF_TurbStruct.H:348
bool use_keqn
Definition: ERF_TurbStruct.H:377
amrex::Real pbl_mrf_coriolis_freq
Definition: ERF_TurbStruct.H:400
MYNNLevel2 pbl_mynn_level2
Definition: ERF_TurbStruct.H:373
RANSType rans_type
Definition: ERF_TurbStruct.H:367
StratType strat_type
Definition: ERF_TurbStruct.H:358
amrex::Real pbl_mrf_Ribcr
Definition: ERF_TurbStruct.H:401
bool advect_tke
Definition: ERF_TurbStruct.H:410
amrex::Real Ck
Definition: ERF_TurbStruct.H:342
amrex::Real Cmu0
Definition: ERF_TurbStruct.H:345
bool pbl_ysu_use_consistent_coriolis
Definition: ERF_TurbStruct.H:389
amrex::Real Cb
Definition: ERF_TurbStruct.H:346
amrex::Real pbl_mrf_sf
Definition: ERF_TurbStruct.H:405
void init_params(int lev, int max_level, std::string pp_prefix)
Definition: ERF_TurbStruct.H:43
amrex::Real pbl_ysu_land_Ribcr
Definition: ERF_TurbStruct.H:395
bool mrf_moistvars
Definition: ERF_TurbStruct.H:407
bool use_tke
Definition: ERF_TurbStruct.H:382
bool diffuse_tke_3D
Definition: ERF_TurbStruct.H:412
amrex::Real Cs
Definition: ERF_TurbStruct.H:336
amrex::Real Pr_t_inv
Definition: ERF_TurbStruct.H:329
amrex::Real Ce_wall
Definition: ERF_TurbStruct.H:341
amrex::Real Ce
Definition: ERF_TurbStruct.H:340
LESType les_type
Definition: ERF_TurbStruct.H:325
amrex::Real pbl_ysu_coriolis_freq
Definition: ERF_TurbStruct.H:386
bool pbl_ysu_force_over_water
Definition: ERF_TurbStruct.H:392
void display(int lev)
Definition: ERF_TurbStruct.H:210
bool use_smag_stratification
Definition: ERF_TurbStruct.H:364
amrex::Real Pr_t
Definition: ERF_TurbStruct.H:328
bool use_kturb
Definition: ERF_TurbStruct.H:376
amrex::Real theta_ref
Definition: ERF_TurbStruct.H:356
amrex::Real Sc_t
Definition: ERF_TurbStruct.H:332
amrex::Real Rt_crit
Definition: ERF_TurbStruct.H:347
amrex::Real pbl_ysu_unst_Ribcr
Definition: ERF_TurbStruct.H:397
bool mix_isotropic
Definition: ERF_TurbStruct.H:361
amrex::Real l_g_max
Definition: ERF_TurbStruct.H:349