ww3_gse2: latest work ...

erdc · Jan 22, 2025 · 9c3a6a9 · 9c3a6a9
1 parent 5139344
commit 9c3a6a9
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/model/src/w3profsmd_pdlib.F90 b/model/src/w3profsmd_pdlib.F90
@@ -1767,7 +1767,7 @@ SUBROUTINE PDLIB_W3XYPFSFCT2 ( ISP, C, LCALC, RD10, RD20, DT, ACIN, ACOUT)
           END IF
         END IF
         ST(NI) = ST(NI) + THETA_L(:,IE) ! the 2nd term are the theta values of each node ...
-        THETA_H         = (1./3.+DT/(2.*PDLIB_TRIA(IE)) * KELEM(:,IE) ) * FT ! LAX
+        THETA_H         = (1./3.+DT/ITER(IK,ITH)/(2*PDLIB_TRIA(IE)) * KELEM(:,IE) ) * FT ! LAX
         !        THETA_H = (1./3.+2./3.*KELEM(:,IE)/SUM(MAX(ZERO,KELEM(:,IE))))*FT  ! CENTRAL ... can be tested as well a bit more dispersive then LAX
         THETA_ACE(:,IE) = THETA_H-THETA_L(:,IE)
         PP(NI) =  PP(NI) + MAX(ZERO, -THETA_ACE(:,IE)) * DTSI(NI)