1 #ifndef ERF_TURB_STRUCT_H_
2 #define ERF_TURB_STRUCT_H_
6 AMREX_ENUM(LESType, None, Smagorinsky, Smagorinsky2D, Deardorff);
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) {
158 if (
les_type == LESType::Deardorff) {
167 amrex::Print() <<
"Thermal stratification based on gradient of potential temperature" << std::endl;
169 amrex::Print() <<
"Thermal stratification based on gradient of virtual potential temperature" << std::endl;
171 amrex::Print() <<
"Thermal stratification based on gradient of linearized liquid-water potential temperature" << std::endl;
199 if (
les_type == LESType::Smagorinsky) {
201 amrex::Error(
"Need to specify Cs for Smagorsinky LES");
204 amrex::Print() <<
"Turning off mix_isotropic for 2-D Smagorinsky" << std::endl;
212 amrex::Print() <<
"Turbulence Settings at level " << lev << std::endl;
217 amrex::Print() <<
" Using DNS model at level " << lev << std::endl;
218 }
else if (
les_type == LESType::Smagorinsky) {
220 amrex::Print() <<
" Using 2D Smagorinsky LES model at level " << lev << std::endl;
222 amrex::Print() <<
" Using Smagorinsky LES model at level " << lev << std::endl;
225 amrex::Print() <<
" Smagorinsky accounts for moist stratification" << std::endl;
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) {
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;
241 amrex::Error(
"Unknown turbulence model");
245 if (
les_type == LESType::Smagorinsky) {
246 amrex::Print() <<
" Cs : " <<
Cs << std::endl;
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;
255 amrex::Real Cs_equiv = std::sqrt(
Ck * std::sqrt(
Ck /
Ce));
256 amrex::Print() <<
" equivalent Cs : " << Cs_equiv
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;
274 if ((
les_type == LESType::Deardorff) ||
278 amrex::Print() <<
" reference theta : " <<
theta_ref << std::endl;
280 amrex::Print() <<
" reference theta : n/a" << std::endl;
285 amrex::Print() <<
" Pr_t : " <<
Pr_t << std::endl;
286 amrex::Print() <<
" Sc_t : " <<
Sc_t << std::endl;
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 : "
302 amrex::Print() <<
" pbl_ysu_use_consistent_coriolis : "
304 amrex::Print() <<
" pbl_ysu_force_over_water : "
306 amrex::Print() <<
" pbl_ysu_land_Ribcr : "
308 amrex::Print() <<
" pbl_ysu_unst_Ribcr : "
310 }
else if (
pbl_type == PBLType::MRF) {
317 amrex::Print() <<
" pbl_mrf_sf : " <<
pbl_mrf_sf
328 amrex::Real
Pr_t = amrex::Real(1.0) / amrex::Real(3.0);
336 amrex::Real
Cs = 0.0;
340 amrex::Real
Ce = 0.93;
342 amrex::Real
Ck = 0.1;
346 amrex::Real
Cb = 0.35;
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