199 amrex::Print() <<
"Turbulence Settings at level " << lev << std::endl;
204 amrex::Print() <<
" Using DNS model at level " << lev << std::endl;
205 }
else if (
les_type == LESType::Smagorinsky) {
207 amrex::Print() <<
" Using 2D Smagorinsky LES model at level " << lev
210 amrex::Print() <<
" Using Smagorinsky LES model at level " << lev
213 }
else if (
les_type == LESType::Deardorff) {
214 amrex::Print() <<
" Using Deardorff LES model at level " << lev
216 }
else if (
rans_type == RANSType::kEqn) {
218 <<
" Using Axell & Liungman one-equation RANS k model at level "
220 }
else if (
pbl_type == PBLType::MYNN25) {
221 amrex::Print() <<
" Using MYNN2.5 PBL model at level " << lev
223 }
else if (
pbl_type == PBLType::MYNNEDMF) {
224 amrex::Print() <<
" Using MYNNEDMF PBL model at level " << lev
226 }
else if (
pbl_type == PBLType::YSU) {
227 amrex::Print() <<
" Using YSU PBL model at level " << lev << std::endl;
228 }
else if (
pbl_type == PBLType::MRF) {
229 amrex::Print() <<
" Using MRF PBL model at level " << lev << std::endl;
231 amrex::Error(
"Unknown turbulence model");
235 if (
les_type == LESType::Smagorinsky) {
236 amrex::Print() <<
"Cs : " <<
Cs << std::endl;
238 if (
les_type == LESType::Deardorff) {
239 amrex::Print() <<
"Ce : " <<
Ce << std::endl;
240 amrex::Print() <<
"Ce at wall : " <<
Ce_wall
242 amrex::Print() <<
"Ck : " <<
Ck << std::endl;
243 amrex::Print() <<
"sigma_k : " <<
sigma_k
245 amrex::Print() <<
"reference theta : " <<
theta_ref
248 amrex::Real Cs_equiv = std::sqrt(
Ck * std::sqrt(
Ck /
Ce));
249 amrex::Print() <<
"equivalent Cs : " << Cs_equiv
258 amrex::Print() <<
"Cmu0 : " <<
Cmu0 << std::endl;
259 amrex::Print() <<
"sigma_k : " <<
sigma_k
261 amrex::Print() <<
"Cb : " <<
Cb << std::endl;
262 amrex::Print() <<
"Rt_crit : " <<
Rt_crit
264 amrex::Print() <<
"Rt_min : " <<
Rt_min
266 amrex::Print() <<
"max_geom_lscale : " <<
l_g_max
268 amrex::Print() <<
"reference theta : " <<
theta_ref
274 amrex::Print() <<
"Pr_t : " <<
Pr_t << std::endl;
275 amrex::Print() <<
"Sc_t : " <<
Sc_t << std::endl;
279 amrex::Print() <<
" pbl_mynn_A1 : " <<
pbl_mynn.
A1 << std::endl;
280 amrex::Print() <<
" pbl_mynn_A2 : " <<
pbl_mynn.
A2 << std::endl;
281 amrex::Print() <<
" pbl_mynn_B1 : " <<
pbl_mynn.
B1 << std::endl;
282 amrex::Print() <<
" pbl_mynn_B2 : " <<
pbl_mynn.
B2 << std::endl;
283 amrex::Print() <<
" pbl_mynn_C1 : " <<
pbl_mynn.
C1 << std::endl;
284 amrex::Print() <<
" pbl_mynn_C2 : " <<
pbl_mynn.
C2 << std::endl;
285 amrex::Print() <<
" pbl_mynn_C3 : " <<
pbl_mynn.
C3 << std::endl;
286 amrex::Print() <<
" pbl_mynn_C4 : " <<
pbl_mynn.
C4 << std::endl;
287 amrex::Print() <<
" pbl_mynn_C5 : " <<
pbl_mynn.
C5 << std::endl;
288 }
else if (
pbl_type == PBLType::YSU) {
289 amrex::Print() <<
" pbl_ysu_coriolis_freq : "
291 amrex::Print() <<
" pbl_ysu_use_consistent_coriolis : "
293 amrex::Print() <<
" pbl_ysu_force_over_water : "
295 amrex::Print() <<
" pbl_ysu_land_Ribcr : "
297 amrex::Print() <<
" pbl_ysu_unst_Ribcr : "
299 }
else if (
pbl_type == PBLType::MRF) {
306 amrex::Print() <<
" pbl_mrf_sf : " <<
pbl_mrf_sf
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:326
amrex::Real sigma_k
Definition: ERF_TurbStruct.H:342
MYNNLevel25 pbl_mynn
Definition: ERF_TurbStruct.H:356
PBLType pbl_type
Definition: ERF_TurbStruct.H:354
amrex::Real pbl_mrf_const_b
Definition: ERF_TurbStruct.H:387
amrex::Real Rt_min
Definition: ERF_TurbStruct.H:337
amrex::Real pbl_mrf_coriolis_freq
Definition: ERF_TurbStruct.H:384
RANSType rans_type
Definition: ERF_TurbStruct.H:351
amrex::Real pbl_mrf_Ribcr
Definition: ERF_TurbStruct.H:385
amrex::Real Ck
Definition: ERF_TurbStruct.H:331
amrex::Real Cmu0
Definition: ERF_TurbStruct.H:334
bool pbl_ysu_use_consistent_coriolis
Definition: ERF_TurbStruct.H:373
amrex::Real Cb
Definition: ERF_TurbStruct.H:335
amrex::Real pbl_mrf_sf
Definition: ERF_TurbStruct.H:389
amrex::Real pbl_ysu_land_Ribcr
Definition: ERF_TurbStruct.H:379
bool mrf_moistvars
Definition: ERF_TurbStruct.H:391
amrex::Real Cs
Definition: ERF_TurbStruct.H:325
amrex::Real Ce_wall
Definition: ERF_TurbStruct.H:330
amrex::Real Ce
Definition: ERF_TurbStruct.H:329
LESType les_type
Definition: ERF_TurbStruct.H:314
amrex::Real pbl_ysu_coriolis_freq
Definition: ERF_TurbStruct.H:370
bool pbl_ysu_force_over_water
Definition: ERF_TurbStruct.H:376
amrex::Real Pr_t
Definition: ERF_TurbStruct.H:317
amrex::Real theta_ref
Definition: ERF_TurbStruct.H:345
amrex::Real Sc_t
Definition: ERF_TurbStruct.H:321
amrex::Real Rt_crit
Definition: ERF_TurbStruct.H:336
amrex::Real pbl_ysu_unst_Ribcr
Definition: ERF_TurbStruct.H:381
bool mix_isotropic
Definition: ERF_TurbStruct.H:348
amrex::Real l_g_max
Definition: ERF_TurbStruct.H:338