diff --git a/src/simulation/m_weno.fpp b/src/simulation/m_weno.fpp index 3c787afd8d..b27243d553 100644 --- a/src/simulation/m_weno.fpp +++ b/src/simulation/m_weno.fpp @@ -853,9 +853,10 @@ contains ! Fu, et al. (2016) ! Fu''s code: https://dx.doi.org/10.13140/RG.2.2.36250.34247 tau = abs(beta(2) - beta(0)) - alpha = (1._wp + tau/beta)**6._wp ! Equation 22 (reuse alpha as gamma; pick C=1 & q=6) + alpha = 1._wp + tau/beta ! Equation 22 (reuse alpha as gamma; pick C=1 & q=6) + alpha = (alpha*alpha*alpha)**2._wp ! Equation 22 cont. (some CPU compilers cannot optimize x**6.0) omega = alpha/sum(alpha) ! Equation 25 (reuse omega as xi) - delta = merge(0._wp, 1._wp, omega < teno_CT) ! Equation 26 + delta = merge(0._wp, 1._wp, omega < teno_CT)! Equation 26 alpha = delta*d_cbL_${XYZ}$ (:, j) ! Equation 27 end if @@ -1036,7 +1037,8 @@ contains elseif (teno) then tau = abs(beta(4) - beta(3)) ! Note the reordering of stencils - alpha = (1._wp + tau/beta)**6._wp + alpha = 1._wp + tau/beta + alpha = (alpha*alpha*alpha)**2._wp ! some CPU compilers cannot optimize x**6.0 omega = alpha/sum(alpha) delta = merge(0._wp, 1._wp, omega < teno_CT) alpha = delta*d_cbL_${XYZ}$ (:, j)