895 amrex::Real ustar = 0.0;
896 amrex::Real wstar = 0.0;
897 amrex::Real tflux = 0.0;
898 amrex::Real qflux = 0.0;
899 amrex::Real z0 = 0.0;
900 amrex::Real zeta = 0.0;
901 amrex::Real psi_m = 0.0;
902 amrex::Real psi_h = 0.0;
903 amrex::Real Olen = 0.0;
905 ie = i < lbound(eta_arr).x ? lbound(eta_arr).x : i;
906 je = j < lbound(eta_arr).y ? lbound(eta_arr).y : j;
907 ie = ie > ubound(eta_arr).x ? ubound(eta_arr).x : ie;
908 je = je > ubound(eta_arr).y ? ubound(eta_arr).y : je;
909 amrex::Real umm = std::max(umm_arr(i,j,k),
WSMIN);
910 if (u_star_arr(i,j,k) == 1.E34) {
913 Olen = olen_arr(i,j,k);
919 ustar = u_star_arr(i,j,k);
920 z0 = std::min( std::max(1200.0 * Hwave_arr(i,j,k) * std::pow( Hwave_arr(i,j,k)/(Lwave_arr(i,j,k)+
eps), 4.5 )
923 -(qvm_arr(i,j,k) - q_surf_arr(i,j,k)) * ustar *
mdata.
kappa /
924 (std::log(
mdata.
zref / z0_arr(i,j,k)) - psi_h);
928 w_star_arr(i,j,k) =
calc_wstar(tflux, pblh_arr(i,j,k), tvm_arr(i,j,k));
930 umm = std::sqrt(umm_arr(i,j,k)*umm_arr(i,j,k) + wstar*wstar);
931 umm = std::max(umm,
WSMIN);
939 }
while ((std::abs(u_star_arr(i,j,k) - ustar) >
tol) && iter <= max_iters);
940 AMREX_ASSERT_WITH_MESSAGE(iter < max_iters,
"Maximum number of MOST iterations reached.");
943 (u_star_arr(i,j,k) *
mdata.
kappa) + tm_arr(i,j,k);
945 olen_arr(i,j,k) = Olen;
949 (u_star_arr(i,j,k) *
mdata.
kappa) + qvm_arr(i,j,k);
952 q_star_arr(i,j,k) =
mdata.
kappa * (qvm_arr(i,j,k) - q_surf_arr(i,j,k)) /
953 (std::log(
mdata.
zref / z0_arr(i,j,k)) - psi_h);
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
@ Mom_v
Definition: ERF_IndexDefines.H:175
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_m(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:41
AMREX_GPU_HOST_DEVICE AMREX_FORCE_INLINE amrex::Real calc_psi_h(amrex::Real zeta) const
Definition: ERF_MOSTStress.H:55
const amrex::Real WSMIN
Definition: ERF_MOSTStress.H:965
const amrex::Real z0_max
Definition: ERF_MOSTStress.H:964
const amrex::Real z0_eps
Definition: ERF_MOSTStress.H:963
similarity_funs sfuns
Definition: ERF_MOSTStress.H:960
const amrex::Real tol
Definition: ERF_MOSTStress.H:961
const amrex::Real eps
Definition: ERF_MOSTStress.H:962