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
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: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
amrex::Real Rt_min
Definition: ERF_TurbStruct.H:348
amrex::Real pbl_mrf_coriolis_freq
Definition: ERF_TurbStruct.H:400
RANSType rans_type
Definition: ERF_TurbStruct.H:367
amrex::Real pbl_mrf_Ribcr
Definition: ERF_TurbStruct.H:401
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
amrex::Real pbl_ysu_land_Ribcr
Definition: ERF_TurbStruct.H:395
bool mrf_moistvars
Definition: ERF_TurbStruct.H:407
amrex::Real Cs
Definition: ERF_TurbStruct.H:336
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
bool use_smag_stratification
Definition: ERF_TurbStruct.H:364
amrex::Real Pr_t
Definition: ERF_TurbStruct.H:328
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