add feature for TAWBUG

umr-lops · Dec 22, 2023 · 1369c6a · 1369c6a
1 parent 51071a7
commit 1369c6a
Showing 1 changed file with 15 additions and 4 deletions.
diff --git a/model/src/w3srcemd.F90 b/model/src/w3srcemd.F90
@@ -1034,10 +1034,14 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD,          &
     TWS = 1./FMEANWS
 #endif
 #ifdef W3_ST4
-    TAUWX=0.
-    TAUWY=0.
-    IF ( IT .eq. 0 ) THEN
+    IF (SINTAILPAR(4).GT.0.5) THEN ! this is designed to keep the bug as an option
+      TAUWX=0.
+      TAUWY=0.
+    END IF
+    IF ( IT .EQ. 0 ) THEN
       LLWS(:) = .TRUE.
+      TAUWX=0.
+      TAUWY=0.
       USTAR=0.
       USTDIR=0.
     ELSE
@@ -1061,7 +1065,7 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD,          &
 #endif
 
 #ifdef W3_ST4
-      CALL W3SIN4 ( SPEC, CG1, WN2, U10ABS, USTAR, DRAT, AS,       &
+      IF (SINTAILPAR(4).GT.0.5) CALL W3SIN4 ( SPEC, CG1, WN2, U10ABS, USTAR, DRAT, AS,       &
            U10DIR, Z0, CD, TAUWX, TAUWY, TAUWAX, TAUWAY,       &
            VSIN, VDIN, LLWS, IX, IY, BRLAMBDA )
     END IF
@@ -1907,6 +1911,13 @@ SUBROUTINE W3SRCE ( srce_call, IT, ISEA, JSEA, IX, IY, IMOD,          &
       CALL W3SIN4 ( SPEC, CG1, WN2, U10ABS, USTAR, DRAT, AS,      &
            U10DIR, Z0, CD, TAUWX, TAUWY, TAUWAX, TAUWAY, &
            VSIN, VDIN, LLWS, IX, IY, BRLAMBDA )
+      IF (SINTAILPAR(4).LT.0.5) CALL W3SPR4 (SPEC, CG1, WN1, EMEAN, FMEAN, FMEAN1, WNMEAN,&
+           AMAX, U10ABS, U10DIR,                          &
+#ifdef W3_FLX5
+           TAUA, TAUADIR, DAIR,                     &
+#endif
+           USTAR, USTDIR,                                 &
+           TAUWX, TAUWY, CD, Z0, CHARN, LLWS, FMEANWS, DLWMEAN)
 #endif
 
       !