Changeset 5082 for LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F

Timestamp:

Jul 19, 2024, 5:41:58 PM (2 months ago)

Author:

abarral

Message:

(lint) Fix obsolete boolean operators

File:

• LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F (modified) (13 diffs)

LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F

@@ -162 +162 @@
 C
 C***FIRST EXECUTABLE STATEMENT  PCHSP
-      IF ( N.LT.2 )  GO TO 5001
-      IF ( INCFD.LT.1 )  GO TO 5002
+      IF ( N<2 )  GO TO 5001
+      IF ( INCFD<1 )  GO TO 5002
       DO 1  J = 2, N
-         IF ( X(J).LE.X(J-1) )  GO TO 5003
+         IF ( X(J)<=X(J-1) )  GO TO 5003
     1 CONTINUE
 C
@@ -171 +171 @@
       IEND = IC(2)
       IERR = 0
-      IF ( (IBEG.LT.0).OR.(IBEG.GT.4) )  IERR = IERR - 1
-      IF ( (IEND.LT.0).OR.(IEND.GT.4) )  IERR = IERR - 2
-      IF ( IERR.LT.0 )  GO TO 5004
+      IF ( (IBEG<0).OR.(IBEG>4) )  IERR = IERR - 1
+      IF ( (IEND<0).OR.(IEND>4) )  IERR = IERR - 2
+      IF ( IERR<0 )  GO TO 5004
 C
 C  FUNCTION DEFINITION IS OK -- GO ON.
 C
-      IF ( NWK .LT. 2*N )  GO TO 5007
+      IF ( NWK < 2*N )  GO TO 5007
 C
 C  COMPUTE FIRST DIFFERENCES OF X SEQUENCE AND STORE IN WK(1,.). ALSO,
@@ -188 +188 @@
 C  SET TO DEFAULT BOUNDARY CONDITIONS IF N IS TOO SMALL.
 C
-      IF ( IBEG.GT.N )  IBEG = 0
-      IF ( IEND.GT.N )  IEND = 0
+      IF ( IBEG>N )  IBEG = 0
+      IF ( IEND>N )  IEND = 0
 C
 C  SET UP FOR BOUNDARY CONDITIONS.
 C
-      IF ( (IBEG.EQ.1).OR.(IBEG.EQ.2) )  THEN
+      IF ( (IBEG==1).OR.(IBEG==2) )  THEN
          D(1,1) = VC(1)
-      ELSE IF (IBEG .GT. 2)  THEN
+      ELSE IF (IBEG > 2)  THEN
 C        PICK UP FIRST IBEG POINTS, IN REVERSE ORDER.
          DO 10  J = 1, IBEG
@@ -201 +201 @@
 C           INDEX RUNS FROM IBEG DOWN TO 1.
             XTEMP(J) = X(INDEX)
-            IF (J .LT. IBEG)  STEMP(J) = WK(2,INDEX)
+            IF (J < IBEG)  STEMP(J) = WK(2,INDEX)
    10    CONTINUE
 C                 --------------------------------
          D(1,1) = PCHDF (IBEG, XTEMP, STEMP, IERR)
 C                 --------------------------------
-         IF (IERR .NE. 0)  GO TO 5009
+         IF (IERR /= 0)  GO TO 5009
          IBEG = 1
       ENDIF
 C
-      IF ( (IEND.EQ.1).OR.(IEND.EQ.2) )  THEN
+      IF ( (IEND==1).OR.(IEND==2) )  THEN
          D(1,N) = VC(2)
-      ELSE IF (IEND .GT. 2)  THEN
+      ELSE IF (IEND > 2)  THEN
 C        PICK UP LAST IEND POINTS.
          DO 15  J = 1, IEND
@@ -218 +218 @@
 C           INDEX RUNS FROM N+1-IEND UP TO N.
             XTEMP(J) = X(INDEX)
-            IF (J .LT. IEND)  STEMP(J) = WK(2,INDEX+1)
+            IF (J < IEND)  STEMP(J) = WK(2,INDEX+1)
    15    CONTINUE
 C                 --------------------------------
          D(1,N) = PCHDF (IEND, XTEMP, STEMP, IERR)
 C                 --------------------------------
-         IF (IERR .NE. 0)  GO TO 5009
+         IF (IERR /= 0)  GO TO 5009
          IEND = 1
       ENDIF
@@ -237 +237 @@
 C             WK(2,1)*S(1) + WK(1,1)*S(2) = D(1,1)
 C
-      IF (IBEG .EQ. 0)  THEN
-         IF (N .EQ. 2)  THEN
+      IF (IBEG == 0)  THEN
+         IF (N == 2)  THEN
 C           NO CONDITION AT LEFT END AND N = 2.
             WK(2,1) = ONE
@@ -250 +250 @@
      *                        + WK(1,2)**2*WK(2,3)) / WK(1,1)
          ENDIF
-      ELSE IF (IBEG .EQ. 1)  THEN
+      ELSE IF (IBEG == 1)  THEN
 C        SLOPE PRESCRIBED AT LEFT END.
          WK(2,1) = ONE
@@ -266 +266 @@
 C
       NM1 = N-1
-      IF (NM1 .GT. 1)  THEN
+      IF (NM1 > 1)  THEN
          DO 20 J=2,NM1
-            IF (WK(2,J-1) .EQ. ZERO)  GO TO 5008
+            IF (WK(2,J-1) == ZERO)  GO TO 5008
             G = -WK(1,J+1)/WK(2,J-1)
             D(1,J) = G*D(1,J-1)
@@ -282 +282 @@
 C     SUBSTITUTION, SINCE ARRAYS HAPPEN TO BE SET UP JUST RIGHT FOR IT
 C     AT THIS POINT.
-      IF (IEND .EQ. 1)  GO TO 30
-C
-      IF (IEND .EQ. 0)  THEN
-         IF (N.EQ.2 .AND. IBEG.EQ.0)  THEN
+      IF (IEND == 1)  GO TO 30
+C
+      IF (IEND == 0)  THEN
+         IF (N==2 .AND. IBEG==0)  THEN
 C           NOT-A-KNOT AT RIGHT ENDPOINT AND AT LEFT ENDPOINT AND N = 2.
             D(1,2) = WK(2,2)
             GO TO 30
-         ELSE IF ((N.EQ.2) .OR. (N.EQ.3 .AND. IBEG.EQ.0))  THEN
+         ELSE IF ((N==2) .OR. (N==3 .AND. IBEG==0))  THEN
 C           EITHER (N=3 AND NOT-A-KNOT ALSO AT LEFT) OR (N=2 AND *NOT*
 C           NOT-A-KNOT AT LEFT END POINT).
             D(1,N) = TWO*WK(2,N)
             WK(2,N) = ONE
-            IF (WK(2,N-1) .EQ. ZERO)  GO TO 5008
+            IF (WK(2,N-1) == ZERO)  GO TO 5008
             G = -ONE/WK(2,N-1)
          ELSE
@@ -303 +303 @@
             D(1,N) = ((WK(1,N)+TWO*G)*WK(2,N)*WK(1,N-1)
      *                  + WK(1,N)**2*(F(1,N-1)-F(1,N-2))/WK(1,N-1))/G
-            IF (WK(2,N-1) .EQ. ZERO)  GO TO 5008
+            IF (WK(2,N-1) == ZERO)  GO TO 5008
             G = -G/WK(2,N-1)
             WK(2,N) = WK(1,N-1)
@@ -311 +311 @@
          D(1,N) = THREE*WK(2,N) + HALF*WK(1,N)*D(1,N)
          WK(2,N) = TWO
-         IF (WK(2,N-1) .EQ. ZERO)  GO TO 5008
+         IF (WK(2,N-1) == ZERO)  GO TO 5008
          G = -ONE/WK(2,N-1)
       ENDIF
@@ -318 +318 @@
 C
       WK(2,N) = G*WK(1,N-1) + WK(2,N)
-      IF (WK(2,N) .EQ. ZERO)   GO TO 5008
+      IF (WK(2,N) == ZERO)   GO TO 5008
       D(1,N) = (G*D(1,N-1) + D(1,N))/WK(2,N)
 C
@@ -325 +325 @@
    30 CONTINUE
       DO 40 J=NM1,1,-1
-         IF (WK(2,J) .EQ. ZERO)  GO TO 5008
+         IF (WK(2,J) == ZERO)  GO TO 5008
          D(1,J) = (D(1,J) - WK(1,J)*D(1,J+1))/WK(2,J)
    40 CONTINUE