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::Real Real
Definition: ERF_ShocInterface.H:19
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 B1
Definition: ERF_MYNNStruct.H:43
amrex::Real B2
Definition: ERF_MYNNStruct.H:44
amrex::Real C5
Definition: ERF_MYNNStruct.H:49
amrex::Real A1
Definition: ERF_MYNNStruct.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
amrex::Real Rt_min
Definition: ERF_TurbStruct.H:353
amrex::Real pbl_mrf_coriolis_freq
Definition: ERF_TurbStruct.H:405
RANSType rans_type
Definition: ERF_TurbStruct.H:372
amrex::Real pbl_mrf_Ribcr
Definition: ERF_TurbStruct.H:406
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
amrex::Real pbl_ysu_land_Ribcr
Definition: ERF_TurbStruct.H:400
bool mrf_moistvars
Definition: ERF_TurbStruct.H:412
amrex::Real Cs
Definition: ERF_TurbStruct.H:341
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
bool use_smag_stratification
Definition: ERF_TurbStruct.H:369
amrex::Real Pr_t
Definition: ERF_TurbStruct.H:333
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