Changeset 5084 for LMDZ6/trunk/libf/dyn3d_common/advz.F

Timestamp:

Jul 19, 2024, 6:40:44 PM (2 months ago)

Author:

Laurent Fairhead

Message:

Reverting to r4065. Updating fortran standard broke too much stuff. Will do it by smaller chunks
AB, LF

File:

• LMDZ6/trunk/libf/dyn3d_common/advz.F (modified) (9 diffs)

LMDZ6/trunk/libf/dyn3d_common/advz.F

@@ -117 +117 @@
 C  Conversion du flux de masse en kg.s-1
 
-      DO l = 1,llm
-         DO j = 1,jjp1
-            DO i = 1,iip1  
+      DO 500 l = 1,llm
+         DO 500 j = 1,jjp1
+            DO 500 i = 1,iip1  
 c            wgri (i,j,llm+1-l) =  w (i,j,l) / g 
                wgri (i,j,llm+1-l) =  w (i,j,l) 
@@ -125 +125 @@
 c             wgri (i,j,l) = 0.1               !    w (i,j,l) 
 c             wgri (i,j,llm) = 0.              ! a detruire ult.
-      END DO
-      END DO
-      END DO
+  500 CONTINUE
          DO  j = 1,jjp1
             DO i = 1,iip1  
@@ -139 +137 @@
 C  boucle sur les latitudes
 C
-      DO K=1,LAT
+      DO 1 K=1,LAT
 C
 C  place limits on appropriate moments before transport
@@ -146 +144 @@
       IF(.NOT.LIMIT) GO TO 101
 C
-      DO JV=1,NTRA
-      DO L=1,NIV
-         DO I=1,LON
+      DO 10 JV=1,NTRA
+      DO 10 L=1,NIV
+         DO 100 I=1,LON
             sz(I,K,L,JV)=SIGN(AMIN1(AMAX1(S0(I,K,L,JV),0.),
      +                              ABS(sz(I,K,L,JV))),sz(I,K,L,JV))
-      END DO
-      END DO
-      END DO
+ 100     CONTINUE
+ 10   CONTINUE
 C
  101  CONTINUE
@@ -165 +162 @@
 C  2- reajusts moments remaining in the box
 C
-      DO L=1,NIV-1
+      DO 11 L=1,NIV-1
       LP=L+1
 C
-      DO I=1,LON
-C
-         IF(WGRI(I,K,L)<0.) THEN
+      DO 110 I=1,LON
+C
+         IF(WGRI(I,K,L).LT.0.) THEN
            FM(I,L)=-WGRI(I,K,L)*DTZ
            ALF(I)=FM(I,L)/SM(I,K,LP)
@@ -184 +181 @@
          ALF1Q(I)=ALF1(I)*ALF1(I)
 C
-      END DO
-C
-      DO JV=1,NTRA
-      DO I=1,LON
-C
-         IF(WGRI(I,K,L)<0.) THEN
+ 110  CONTINUE
+C
+      DO 111 JV=1,NTRA
+      DO 1110 I=1,LON
+C
+         IF(WGRI(I,K,L).LT.0.) THEN
 C
            F0(I,L,JV)=ALF (I)*( S0(I,K,LP,JV)-ALF1(I)*sz(I,K,LP,JV) )
@@ -215 +212 @@
          ENDIF
 C
-      END DO
-      END DO
-C
-      END DO
+ 1110 CONTINUE
+ 111  CONTINUE
+C
+ 11   CONTINUE
 C
 C  puts the temporary moments Fi into appropriate neighboring boxes
 C
-      DO L=1,NIV-1
+      DO 12 L=1,NIV-1
       LP=L+1
 C
-      DO I=1,LON
-C
-         IF(WGRI(I,K,L)<0.) THEN
+      DO 120 I=1,LON
+C
+         IF(WGRI(I,K,L).LT.0.) THEN
            SM(I,K,L)=SM(I,K,L)+FM(I,L)
            ALF(I)=FM(I,L)/SM(I,K,L)
@@ -239 +236 @@
          ALF1Q(I)=ALF1(I)*ALF1(I)
 C
-      END DO
-C
-      DO JV=1,NTRA
-      DO I=1,LON
-C
-         IF(WGRI(I,K,L)<0.) THEN
+ 120  CONTINUE
+C
+      DO 121 JV=1,NTRA
+      DO 1210 I=1,LON
+C
+         IF(WGRI(I,K,L).LT.0.) THEN
 C
            TEMPTM=-ALF(I)*S0(I,K,L,JV)+ALF1(I)*F0(I,L,JV)
@@ -263 +260 @@
          ENDIF
 C
-      END DO
-      END DO
-C
-      END DO
+ 1210 CONTINUE
+ 121  CONTINUE
+C
+ 12   CONTINUE
 C
 C  fin de la boucle principale sur les latitudes
 C
-      END DO
+ 1    CONTINUE
 C
 C-------------------------------------------------------------