1196 amrex::Real ustar = 0.0;
1197 amrex::Real wstar = 0.0;
1198 amrex::Real z0 = 0.0;
1199 amrex::Real tflux = 0.0;
1200 amrex::Real zeta = 0.0;
1201 amrex::Real psi_m = 0.0;
1202 amrex::Real psi_h = 0.0;
1203 amrex::Real Olen = 0.0;
1204 amrex::Real Oleno = 0.0;
1205 amrex::Real umm = std::max(umm_arr(i,j,k),
WSMIN);
1206 if (u_star_arr(i,j,k) == 1.E34) {
1209 Olen = olen_arr(i,j,k);
1216 ustar = u_star_arr(i,j,k);
1218 tflux = -(tm_arr(i,j,k) - t_surf_arr(i,j,k)) * ustar *
mdata.
kappa /
1220 tflux *= (1 + 0.61*qvm_arr(i,j,k));
1221 tflux += 0.61*tm_arr(i,j,k) * -ustar*q_star_arr(i,j,k);
1224 w_star_arr(i,j,k) =
calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
1226 umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
1227 umm = std::max(umm,
WSMIN);
1230 if ( ((Olen >= 0.0) && (Oleno <= 0.0)) ||
1231 ((Olen <= 0.0) && (Oleno >= 0.0)) ) {
1232 Olen = 0.5 * (Olen + Oleno);
1240 }
while ((std::abs(u_star_arr(i,j,k) - ustar) >
tol) && iter <= max_iters);
1241 AMREX_ASSERT_WITH_MESSAGE(iter < max_iters,
"Maximum number of MOST iterations reached.");
1243 t_star_arr(i,j,k) =
mdata.
kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) /
1245 olen_arr(i,j,k) = Olen;
AMREX_GPU_DEVICE AMREX_FORCE_INLINE amrex::Real Donelan_roughness(amrex::Real ustar)
Definition: ERF_MOSTRoughness.H:29
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:28
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_m(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:40
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_h(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:54
similarity_funs sfuns
Definition: ERF_MOSTStress.H:1251
const amrex::Real tol
Definition: ERF_MOSTStress.H:1252
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:1253