1111 if (u_star_arr(i,j,k) == 1.E34) { olen_arr(i,j,k) = 1.0e3; }
1116 if (q_surf_arr(i,j,k) > 0.) {
1117 qv_s = q_surf_arr(i,j,k);
1121 qv_s = qvm_arr(i,j,k);
1123 qv_a = qvm_arr(i,j,k);
1133 -ustar *
mdata.
kappa * (qv_a - qv_s) / (C - psi_h);
1134 amrex::Real tflux = -ustar*tstar*(1 + 0.61*qvm_arr(i,j,k)) + 0.61*tm_arr(i,j,k)*qflux;
1135 w_star_arr(i,j,k) =
calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
1137 umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
1138 umm = std::max(umm,
WSMIN);
1142 amrex::Real thv_s = t_surf_arr(i,j,k) * (1.0 + 0.61*qv_s);
1143 amrex::Real thv_a = tm_arr(i,j,k) * (1.0 + 0.61*qv_a);
1145 ( (thv_a - thv_s) / (umm * umm) );
1146 Rib = std::min(std::max(Rib,-4.0),4.0);
1159 num = std::max(C - psi_m, 1.0);
1160 den = std::max(C - psi_h, 1.0);
1166 }
while ( (std::abs(zeta - zeta_old) >
tol) && (iter <= max_iters) );
1167 AMREX_ALWAYS_ASSERT_WITH_MESSAGE(iter < max_iters,
1168 "Maximum number of MOST iterations reached.");
1172 u_star_arr(i,j,k) =
mdata.
kappa * umm / (C - psi_m);
1173 t_star_arr(i,j,k) =
mdata.
kappa * (tm_arr(i,j,k) - t_surf_arr(i,j,k)) / (C - psi_h);
1176 (u_star_arr(i,j,k) *
mdata.
kappa) + qvm_arr(i,j,k);
1179 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::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
similarity_funs sfuns
Definition: ERF_MOSTStress.H:1186
const amrex::Real tol
Definition: ERF_MOSTStress.H:1187
const amrex::Real alpha
Definition: ERF_MOSTStress.H:1188
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:1189