From ba5cd68fbb5a5be98e196b51aecda2f4bbb7f7a9 Mon Sep 17 00:00:00 2001
From: Matthew Masarik <86749872+MatthewMasarik-NOAA@users.noreply.github.com>
Date: Mon, 5 Feb 2024 13:45:25 -0500
Subject: [PATCH 1/2] w3fld1md.F90: fix divide by zero in CRIT2 parameter
 (#1184)

---
 model/src/w3fld1md.F90 | 8 ++++++--
 1 file changed, 6 insertions(+), 2 deletions(-)

diff --git a/model/src/w3fld1md.F90 b/model/src/w3fld1md.F90
index 960fd185a..10b2fce08 100644
--- a/model/src/w3fld1md.F90
+++ b/model/src/w3fld1md.F90
@@ -550,9 +550,13 @@ SUBROUTINE W3FLD1( ASPC, FPI, WNDX,WNDY, ZWND,               &
           TAUDIR=atan2(TAUY, TAUX)
           ! Note: add another criterion (stress direction) for iteration.
           CRIT1=(ABS(USTAR-USTRB)*100.0)/((USTAR+USTRB)*0.5) .GT. 0.1
-          CRIT2=(ABS(TAUDIR-TAUDIRB)*100.0/(TAUDIR+TAUDIRB)*0.5) .GT. 0.1
+          IF ((TAUDIR+TAUDIRB).NE.0.) THEN
+              CRIT2=(ABS(TAUDIR-TAUDIRB)*100.0/(TAUDIR+TAUDIRB)*0.5) .GT. 0.1
+          ELSE
+              CRIT2=.TRUE.
+          ENDIF
           IF (CRIT1 .OR. CRIT2) THEN
-            !            IF ((ABS(USTAR-USTRB)*100.0)/((USTAR+USTRB)*0.5) .GT. 0.1) THEN
+            !    IF ((ABS(USTAR-USTRB)*100.0)/((USTAR+USTRB)*0.5) .GT. 0.1) THEN
             USTRB=USTAR
             TAUDIRB=TAUDIR
             CTR=CTR+1

From fd6d559e5dbcea001dec8f138edca90589098c42 Mon Sep 17 00:00:00 2001
From: Matthew Masarik <86749872+MatthewMasarik-NOAA@users.noreply.github.com>
Date: Mon, 5 Feb 2024 13:45:48 -0500
Subject: [PATCH 2/2] ww3_prnc.F90: fix out-of-scope grid index write statement
 (#1185)

---
 model/src/ww3_prnc.F90 | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/model/src/ww3_prnc.F90 b/model/src/ww3_prnc.F90
index e77bbd918..59747d32a 100644
--- a/model/src/ww3_prnc.F90
+++ b/model/src/ww3_prnc.F90
@@ -1059,7 +1059,7 @@ PROGRAM W3PRNC
           ! Manages the simple closure of the grid
           !
           IF (ICLO.EQ.ICLOSE_NONE) THEN
-            IF (IX21(IX,1).LT.1.OR.IX21(IX,1).GT.NXI-1) WRITE(NDSO,1042) IX, IY, X, Y
+            IF (IX21(IX,1).LT.1.OR.IX21(IX,1).GT.NXI-1) WRITE(NDSO,1041) IX, X, Y
             IX21(IX,1) =   MAX ( 1 , MIN(IX21(IX,1),NXI-1) )
             IX22(IX,1) =   IX21(IX,1) + 1
           ELSE
@@ -1067,7 +1067,7 @@ PROGRAM W3PRNC
             IX22(IX,1) =   MOD(IX21(IX,1),NXI)+1
           END IF
           IY21(IX,1) =   1 + INT((Y-Y0I)/SYI)
-          IF (IY21(IX,1).LT.1.OR.IY21(IX,1).GT.NYI-1) WRITE(NDSO,1042) IX, IY, X, Y
+          IF (IY21(IX,1).LT.1.OR.IY21(IX,1).GT.NYI-1) WRITE(NDSO,1041) IX, X, Y
           IY21(IX,1) =   MAX ( 1 , MIN(IY21(IX,1),NYI-1) )
           IY22(IX,1) =   IY21(IX,1) + 1
           !
@@ -2438,6 +2438,9 @@ PROGRAM W3PRNC
        '     2MS2 2MN2 2NK2 MNS2 MSN2 2SM2 3MSN2 '      &
        '     M4 MS4 MN4 M6 2MS6 2MN6'/)
   !
+1041 FORMAT (/' *** WAVEWATCH-III WARNING W3PRNC : '/                &
+       '     GRID POINT ',I6,2F7.2,/                         &
+       ' NOT COVERED BY INPUT GRID.'/)
 1042 FORMAT (/' *** WAVEWATCH-III WARNING W3PRNC : '/                &
        '     GRID POINT ',2I6,2F7.2,/                         &
        ' NOT COVERED BY INPUT GRID.'/)