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

Timestamp:

Jul 19, 2024, 7:54:50 PM (2 months ago)

Author:

abarral

Message:

convert labeled do (f77) to do .. end do

File:

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

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

@@ -164 +164 @@
       IF ( N<2 )  GO TO 5001
       IF ( INCFD<1 )  GO TO 5002
-      DO 1  J = 2, N
+      DO J = 2, N
          IF ( X(J)<=X(J-1) )  GO TO 5003
-    1 CONTINUE
+      END DO
 C
       IBEG = IC(1)
@@ -181 +181 @@
 C  COMPUTE FIRST DIFFERENCES OF X SEQUENCE AND STORE IN WK(1,.). ALSO,
 C  COMPUTE FIRST DIVIDED DIFFERENCE OF DATA AND STORE IN WK(2,.).
-      DO 5  J=2,N
+      DO J=2,N
          WK(1,J) = X(J) - X(J-1)
          WK(2,J) = (F(1,J) - F(1,J-1))/WK(1,J)
-    5 CONTINUE
+      END DO
 C
 C  SET TO DEFAULT BOUNDARY CONDITIONS IF N IS TOO SMALL.
@@ -197 +197 @@
       ELSE IF (IBEG > 2)  THEN
 C        PICK UP FIRST IBEG POINTS, IN REVERSE ORDER.
-         DO 10  J = 1, IBEG
+         DO J = 1, IBEG
             INDEX = IBEG-J+1
 C           INDEX RUNS FROM IBEG DOWN TO 1.
             XTEMP(J) = X(INDEX)
             IF (J < IBEG)  STEMP(J) = WK(2,INDEX)
-   10    CONTINUE
+      END DO
 C                 --------------------------------
          D(1,1) = PCHDF (IBEG, XTEMP, STEMP, IERR)
@@ -214 +214 @@
       ELSE IF (IEND > 2)  THEN
 C        PICK UP LAST IEND POINTS.
-         DO 15  J = 1, IEND
+         DO J = 1, IEND
             INDEX = N-IEND+J
 C           INDEX RUNS FROM N+1-IEND UP TO N.
             XTEMP(J) = X(INDEX)
             IF (J < IEND)  STEMP(J) = WK(2,INDEX+1)
-   15    CONTINUE
+      END DO
 C                 --------------------------------
          D(1,N) = PCHDF (IEND, XTEMP, STEMP, IERR)
@@ -267 +267 @@
       NM1 = N-1
       IF (NM1 > 1)  THEN
-         DO 20 J=2,NM1
+         DO J=2,NM1
             IF (WK(2,J-1) == ZERO)  GO TO 5008
             G = -WK(1,J+1)/WK(2,J-1)
@@ -273 +273 @@
      *                  + THREE*(WK(1,J)*WK(2,J+1) + WK(1,J+1)*WK(2,J))
             WK(2,J) = G*WK(1,J-1) + TWO*(WK(1,J) + WK(1,J+1))
-   20    CONTINUE
+      END DO
       ENDIF
 C
@@ -324 +324 @@
 C
    30 CONTINUE
-      DO 40 J=NM1,1,-1
+      DO J=NM1,1,-1
          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
+      END DO
 C --------------------(  END  CODING FROM CUBSPL )--------------------
 C