1251 if (u_star_arr(i,j,k) == 1.E34) { olen_arr(i,j,k) = 1.0e3; }
1256 if (q_surf_arr(i,j,k) > 0.) {
1257 qv_s = q_surf_arr(i,j,k);
1261 qv_s = qvm_arr(i,j,k);
1263 qv_a = qvm_arr(i,j,k);
1273 -ustar *
mdata.
kappa * (qv_a - qv_s) / (C - psi_h);
1274 amrex::Real tflux = -ustar*tstar*(1 + 0.61*qvm_arr(i,j,k)) + 0.61*tm_arr(i,j,k)*qflux;
1275 w_star_arr(i,j,k) =
calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
1277 umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
1278 umm = std::max(umm,
WSMIN);
1282 amrex::Real thv_s = t_surf_arr(i,j,k) * (1.0 + 0.61*qv_s);
1283 amrex::Real thv_a = tm_arr(i,j,k) * (1.0 + 0.61*qv_a);
1285 ( (thv_a - thv_s) / (umm * umm) );
1286 Rib = std::min(std::max(Rib,-4.0),4.0);
1317 }
while ( (std::abs(z0 - z0_old) >
tol_z) && (iter_z <= max_iters) );
1318 AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter_z < max_iters,
1319 "Maximum number of MOST roughness iterations reached.");
1323 num = std::max(C - psi_m, 1.0);
1324 den = std::max(C - psi_h, 1.0);
1330 }
while ( (std::abs(zeta - zeta_old) >
tol) && (iter <= max_iters) );
1331 AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
1332 "Maximum number of MOST iterations reached.");
1337 u_star_arr(i,j,k) =
mdata.
kappa * umm / (C - psi_m);
1338 t_star_arr(i,j,k) =
mdata.
kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) / (C - psi_h);
1341 (u_star_arr(i,j,k) *
mdata.
kappa) + qvm_arr(i,j,k);
1344 q_star_arr(i,j,k) =
mdata.
kappa * (qvm_arr(i,j,k) - q_surf_arr(i,j,k)) / (C - psi_h);
@ num
Definition: ERF_DataStruct.H:22
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real COARE3_roughness(amrex::Real zref, amrex::Real umm, amrex::Real ustar)
Definition: ERF_MOSTRoughness.H:9
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real air_viscosity(amrex::Real T_degK)
Definition: ERF_MOSTStress.H:126
amrex::Real Real
Definition: ERF_ShocInterface.H:16
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real calc_wstar(const amrex::Real &ust, const amrex::Real &tst, const amrex::Real &qst, const amrex::Real &pblh, const amrex::Real &th, const amrex::Real &thv, const amrex::Real &qv=0.0)
Definition: ERF_Wstar.H:13
amrex::Real kappa
von Karman constant
Definition: ERF_MOSTStress.H:17
amrex::Real gravity
Acceleration due to gravity (m/s^2)
Definition: ERF_MOSTStress.H:18
const amrex::Real Bjr_beta
Definition: ERF_MOSTStress.H:29
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_h2(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:64
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_m2(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:44
const amrex::Real tol_z
Definition: ERF_MOSTStress.H:1353
const amrex::Real alpha
Definition: ERF_MOSTStress.H:1354
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:1355
const amrex::Real tol
Definition: ERF_MOSTStress.H:1352
similarity_funs sfuns
Definition: ERF_MOSTStress.H:1351