1 #ifndef ERF_TURB_STRUCT_H_
2 #define ERF_TURB_STRUCT_H_
6 AMREX_ENUM(LESType, None, Smagorinsky, Smagorinsky2D, Deardorff);
10 AMREX_ENUM(PBLType, None, MYJ, MYNN25, MYNNEDMF, YSU, MRF);
17 const amrex::ParmParse&
pp,
18 const char* query_string,
23 int count =
pp.countval(query_string);
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);
32 "For parmparse variable " +
pp.prefixedName(query_string) +
33 ": if specified, specify once total or once for each level");
43 void init_params (
int lev,
int max_level, std::string pp_prefix)
45 amrex::ParmParse
pp(pp_prefix);
48 std::string les_type_string =
"None";
52 if (
les_type == LESType::Smagorinsky2D) {
58 std::string rans_type_string =
"None";
62 amrex::Error(
"Hybrid RANS-LES not implemented");
66 static std::string pbl_type_string =
"None";
73 amrex::Print() <<
"Selected a PBL model and an LES model: "
74 <<
"Using PBL for vertical transport, LES for horizontal"
77 if (
les_type == LESType::Smagorinsky) {
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) {
83 "It is not recommended to use Deardorff LES and a PBL model");
112 }
else if (
pbl_type == PBLType::YSU) {
116 pp,
"pbl_ysu_use_consistent_coriolis",
125 }
else if (
pbl_type == PBLType::MRF) {
160 if (
les_type == LESType::Deardorff) {
169 amrex::Print() <<
"Thermal stratification based on gradient of potential temperature" << std::endl;
171 amrex::Print() <<
"Thermal stratification based on gradient of virtual potential temperature" << std::endl;
173 amrex::Print() <<
"Thermal stratification based on gradient of linearized liquid-water potential temperature" << std::endl;
202 if (
les_type == LESType::Smagorinsky) {
204 amrex::Error(
"Need to specify Cs for Smagorsinky LES");
207 amrex::Print() <<
"Turning off mix_isotropic for 2-D Smagorinsky" << std::endl;
215 amrex::Print() <<
"Turbulence Settings at level " << lev << std::endl;
220 amrex::Print() <<
" Using DNS model at level " << lev << std::endl;
221 }
else if (
les_type == LESType::Smagorinsky) {
223 amrex::Print() <<
" Using 2D Smagorinsky LES model at level " << lev << std::endl;
225 amrex::Print() <<
" Using Smagorinsky LES model at level " << lev << std::endl;
228 amrex::Print() <<
" Smagorinsky accounts for moist stratification" << std::endl;
230 }
else if (
les_type == LESType::Deardorff) {
231 amrex::Print() <<
" Using Deardorff LES model at level " << lev << std::endl;
232 }
else if (
rans_type == RANSType::kEqn) {
234 <<
" Using Axell & Liungman one-equation RANS k model at level " << lev << std::endl;
235 }
else if (
pbl_type == PBLType::MYJ) {
236 amrex::Print() <<
" Using MYJ PBL model at level " << lev << std::endl;
237 }
else if (
pbl_type == PBLType::MYNN25) {
238 amrex::Print() <<
" Using MYNN2.5 PBL model at level " << lev << std::endl;
239 }
else if (
pbl_type == PBLType::MYNNEDMF) {
240 amrex::Print() <<
" Using MYNNEDMF PBL model at level " << lev << std::endl;
241 }
else if (
pbl_type == PBLType::YSU) {
242 amrex::Print() <<
" Using YSU PBL model at level " << lev << std::endl;
243 }
else if (
pbl_type == PBLType::MRF) {
244 amrex::Print() <<
" Using MRF PBL model at level " << lev << std::endl;
246 amrex::Error(
"Unknown turbulence model");
250 if (
les_type == LESType::Smagorinsky) {
251 amrex::Print() <<
" Cs : " <<
Cs << std::endl;
253 if (
les_type == LESType::Deardorff) {
254 amrex::Print() <<
" Ce : " <<
Ce << std::endl;
255 amrex::Print() <<
" Ce at wall : " <<
Ce_wall << std::endl;
256 amrex::Print() <<
" Ck : " <<
Ck << std::endl;
257 amrex::Print() <<
" sigma_k : " <<
sigma_k << std::endl;
261 amrex::Print() <<
" equivalent Cs : " << Cs_equiv
270 amrex::Print() <<
"Cmu0 : " <<
Cmu0 << std::endl;
271 amrex::Print() <<
"sigma_k : " <<
sigma_k << std::endl;
272 amrex::Print() <<
"Cb : " <<
Cb << std::endl;
273 amrex::Print() <<
"Rt_crit : " <<
Rt_crit << std::endl;
274 amrex::Print() <<
"Rt_min : " <<
Rt_min << std::endl;
275 amrex::Print() <<
"max_geom_lscale : " <<
l_g_max << std::endl;
279 if ((
les_type == LESType::Deardorff) ||
283 amrex::Print() <<
" reference theta : " <<
theta_ref << std::endl;
285 amrex::Print() <<
" reference theta : n/a" << std::endl;
290 amrex::Print() <<
" Pr_t : " <<
Pr_t << std::endl;
291 amrex::Print() <<
" Sc_t : " <<
Sc_t << std::endl;
295 amrex::Print() <<
" pbl_mynn_A1 : " <<
pbl_mynn.
A1 << std::endl;
296 amrex::Print() <<
" pbl_mynn_A2 : " <<
pbl_mynn.
A2 << std::endl;
297 amrex::Print() <<
" pbl_mynn_B1 : " <<
pbl_mynn.
B1 << std::endl;
298 amrex::Print() <<
" pbl_mynn_B2 : " <<
pbl_mynn.
B2 << std::endl;
299 amrex::Print() <<
" pbl_mynn_C1 : " <<
pbl_mynn.
C1 << std::endl;
300 amrex::Print() <<
" pbl_mynn_C2 : " <<
pbl_mynn.
C2 << std::endl;
301 amrex::Print() <<
" pbl_mynn_C3 : " <<
pbl_mynn.
C3 << std::endl;
302 amrex::Print() <<
" pbl_mynn_C4 : " <<
pbl_mynn.
C4 << std::endl;
303 amrex::Print() <<
" pbl_mynn_C5 : " <<
pbl_mynn.
C5 << std::endl;
304 }
else if (
pbl_type == PBLType::YSU) {
305 amrex::Print() <<
" pbl_ysu_coriolis_freq : "
307 amrex::Print() <<
" pbl_ysu_use_consistent_coriolis : "
309 amrex::Print() <<
" pbl_ysu_force_over_water : "
311 amrex::Print() <<
" pbl_ysu_land_Ribcr : "
313 amrex::Print() <<
" pbl_ysu_unst_Ribcr : "
315 }
else if (
pbl_type == PBLType::MRF) {
322 amrex::Print() <<
" pbl_mrf_sf : " <<
pbl_mrf_sf
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real pp(amrex::Real y)
Definition: ERF_MicrophysicsUtils.H:230
amrex::Real Real
Definition: ERF_ShocInterface.H:19
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:342
amrex::Real sigma_k
Definition: ERF_TurbStruct.H:358
MYNNLevel25 pbl_mynn
Definition: ERF_TurbStruct.H:377
PBLType pbl_type
Definition: ERF_TurbStruct.H:375
amrex::Real pbl_mrf_const_b
Definition: ERF_TurbStruct.H:408
bool use_pbl_tke
Definition: ERF_TurbStruct.H:385
amrex::Real Sc_t_inv
Definition: ERF_TurbStruct.H:338
amrex::Real Rt_min
Definition: ERF_TurbStruct.H:353
bool use_keqn
Definition: ERF_TurbStruct.H:382
amrex::Real pbl_mrf_coriolis_freq
Definition: ERF_TurbStruct.H:405
MYNNLevel2 pbl_mynn_level2
Definition: ERF_TurbStruct.H:378
RANSType rans_type
Definition: ERF_TurbStruct.H:372
StratType strat_type
Definition: ERF_TurbStruct.H:363
amrex::Real pbl_mrf_Ribcr
Definition: ERF_TurbStruct.H:406
bool advect_tke
Definition: ERF_TurbStruct.H:415
amrex::Real Ck
Definition: ERF_TurbStruct.H:347
amrex::Real Cmu0
Definition: ERF_TurbStruct.H:350
bool pbl_ysu_use_consistent_coriolis
Definition: ERF_TurbStruct.H:394
amrex::Real Cb
Definition: ERF_TurbStruct.H:351
amrex::Real pbl_mrf_sf
Definition: ERF_TurbStruct.H:410
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:400
bool mrf_moistvars
Definition: ERF_TurbStruct.H:412
bool use_tke
Definition: ERF_TurbStruct.H:387
bool diffuse_tke_3D
Definition: ERF_TurbStruct.H:417
amrex::Real Cs
Definition: ERF_TurbStruct.H:341
amrex::Real Pr_t_inv
Definition: ERF_TurbStruct.H:334
amrex::Real Ce_wall
Definition: ERF_TurbStruct.H:346
amrex::Real Ce
Definition: ERF_TurbStruct.H:345
LESType les_type
Definition: ERF_TurbStruct.H:330
amrex::Real pbl_ysu_coriolis_freq
Definition: ERF_TurbStruct.H:391
bool pbl_ysu_force_over_water
Definition: ERF_TurbStruct.H:397
void display(int lev)
Definition: ERF_TurbStruct.H:213
bool use_smag_stratification
Definition: ERF_TurbStruct.H:369
amrex::Real Pr_t
Definition: ERF_TurbStruct.H:333
bool use_kturb
Definition: ERF_TurbStruct.H:381
amrex::Real theta_ref
Definition: ERF_TurbStruct.H:361
amrex::Real Sc_t
Definition: ERF_TurbStruct.H:337
amrex::Real Rt_crit
Definition: ERF_TurbStruct.H:352
amrex::Real pbl_ysu_unst_Ribcr
Definition: ERF_TurbStruct.H:402
bool mix_isotropic
Definition: ERF_TurbStruct.H:366
amrex::Real l_g_max
Definition: ERF_TurbStruct.H:354