♻️ DISTP is an integer, not logical

SanderHulst · May 6, 2022 · 91f6696 · 91f6696
1 parent 0155d1a
commit 91f6696
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/model/ftn/w3partmd.ftn b/model/ftn/w3partmd.ftn
@@ -1276,15 +1276,15 @@
         DENOM  = XL*EH - XH*EL
         SIGP   = SIG(IFPMAX(IP))
      !/ --- Parabolic fit around the spectral peak ---
-        IF (.NOT.DISTP .AND. DENOM.NE.0.) THEN
+        IF (DISTP.NE.1 .AND. DENOM.NE.0.) THEN
            SIGP   = SIGP *( 1. + 0.5 * ( XL2*EH - XH2*EL ) &
                          / SIGN ( ABS(DENOM) , DENOM ) )
         END IF
         CALL WAVNU1 ( SIGP, DEPTH, WNP, CGP )
 !
         IK = IFPMAX(IP)
         EFPMAX(IP) = SUMF(IK,IP) * DTH
-        IF (.NOT.DISTP .AND. IK.GT.1 .AND. IK.LT.NK) THEN
+        IF (DISTP.NE.1 .AND. IK.GT.1 .AND. IK.LT.NK) THEN
           EL    = SUMF(IK-1,IP) * DTH
           EH    = SUMF(IK+1,IP) * DTH
           NUMER = 0.125 * ( EL - EH )**2