Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F

Timestamp:

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

Author:

abarral

Message:

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

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F (modified) (4 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F

@@ -57 +57 @@
          deuxjour = 2. * daysec
 
-         DO   1  ij   = 1, ip1jmp1
+         DO   ij   = 1, ip1jmp1
          unsaire2(ij) = unsaire(ij) * unsaire(ij)
-   1     CONTINUE
+      END DO
       END IF
 
@@ -100 +100 @@
 
 c
-      DO 20 l = 1, llmm1
+      DO l = 1, llmm1
 
 
 c       ......   calcul de  - w/2.    au niveau  l+1   .......
 
-      DO 5   ij   = 1, ip1jmp1
+      DO ij   = 1, ip1jmp1
       wsur2( ij ) = - 0.5 * w( ij,l+1 )
-   5  CONTINUE
+      END DO
 
 
 c     .....................     calcul pour  du     ..................
 
-      DO 6 ij = iip2 ,ip1jm-1
+      DO ij = iip2 ,ip1jm-1
       ww        = wsur2 (  ij  )     + wsur2( ij+1 ) 
       uu        = 0.5 * ( ucov(ij,l) + ucov(ij,l+1) )
       du(ij,l)  = du(ij,l)   - ww * ( uu - uav(ij, l ) )/massebx(ij, l )
       du(ij,l+1)= du(ij,l+1) + ww * ( uu - uav(ij,l+1) )/massebx(ij,l+1)
-   6  CONTINUE
+      END DO
 
 c     .....  correction pour  du(iip1,j,l)  ........
@@ -123 +123 @@
 
 CDIR$ IVDEP
-      DO   7  ij   = iip1 +iip1, ip1jm, iip1
+      DO   ij   = iip1 +iip1, ip1jm, iip1
       du( ij, l  ) = du( ij -iim, l  )
       du( ij,l+1 ) = du( ij -iim,l+1 )
-   7  CONTINUE
+      END DO
 
 c     .................    calcul pour   dv      .....................
 
-      DO 8 ij = 1, ip1jm
+      DO ij = 1, ip1jm
       ww        = wsur2( ij+iip1 )   + wsur2( ij )
       vv        = 0.5 * ( vcov(ij,l) + vcov(ij,l+1) )
       dv(ij,l)  = dv(ij, l ) - ww * (vv - vav(ij, l ) )/masseby(ij, l )
       dv(ij,l+1)= dv(ij,l+1) + ww * (vv - vav(ij,l+1) )/masseby(ij,l+1)
-   8  CONTINUE
+      END DO
 
 c
@@ -147 +147 @@
 c                   ...............
 
-        DO 15 ij = 1, ip1jmp1
+        DO ij = 1, ip1jmp1
          ww            = wsur2(ij) * (teta(ij,l) + teta(ij,l+1) )
          dteta(ij, l ) = dteta(ij, l )  -  ww
          dteta(ij,l+1) = dteta(ij,l+1)  +  ww
-  15    CONTINUE
+      END DO
 
       IF( conser)  THEN
-        DO 17 ij = 1,ip1jmp1
+        DO ij = 1,ip1jmp1
         ge(ij)   = wsur2(ij) * wsur2(ij) * unsaire2(ij)
-  17    CONTINUE
+      END DO
         gt       = SSUM( ip1jmp1,ge,1 )
         gtot(l)  = deuxjour * SQRT( gt/ip1jmp1 )
       END IF
 
-  20  CONTINUE
+      END DO
  
       RETURN