Ignore:
Timestamp:
Jul 19, 2024, 11:28:59 AM (4 months ago)
Author:
abarral
Message:

[coherence with Fortran standards]
Replace obsolete DO with shared termination
(minor) replace obsolete bool operators

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dyn3d_common/advzp.F

    r2600 r5079  
    135135C  Conversion des flux de masses en kg
    136136
    137       DO 500 l = 1,llm
    138          DO 500 j = 1,jjp1
    139             DO 500 i = 1,iip1 
     137      DO l = 1,llm
     138         DO j = 1,jjp1
     139            DO i = 1,iip1 
    140140            wgri (i,j,llm+1-l) = w (i,j,l) 
    141   500 CONTINUE
     141      END DO
     142      END DO
     143      END DO
    142144      do j=1,jjp1
    143145         do i=1,iip1
     
    154156C  boucle sur les latitudes
    155157C
    156       DO 1 K=1,LAT
     158      DO K=1,LAT
    157159C
    158160C  place limits on appropriate moments before transport
     
    161163      IF(.NOT.LIMIT) GO TO 101
    162164C
    163       DO 10 JV=1,NTRA
    164       DO 10 L=1,NIV
    165          DO 100 I=1,LON
    166             IF(S0(I,K,L,JV).GT.0.) THEN
     165      DO JV=1,NTRA
     166      DO L=1,NIV
     167         DO I=1,LON
     168            IF(S0(I,K,L,JV)>0.) THEN
    167169              SLPMAX=S0(I,K,L,JV)
    168170              S1MAX =1.5*SLPMAX
     
    180182              SYZ(I,K,L,JV)=0.
    181183            ENDIF
    182  100     CONTINUE
    183  10   CONTINUE
     184      END DO
     185      END DO
     186      END DO
    184187C
    185188 101  CONTINUE
     
    193196C  2- reajusts moments remaining in the box
    194197C
    195       DO 11 L=1,NIV-1
     198      DO L=1,NIV-1
    196199      LP=L+1
    197200C
    198       DO 110 I=1,LON
    199 C
    200          IF(WGRI(I,K,L).LT.0.) THEN
     201      DO I=1,LON
     202C
     203         IF(WGRI(I,K,L)<0.) THEN
    201204           FM(I,L)=-WGRI(I,K,L)*DTZ
    202205           ALF(I)=FM(I,L)/SM(I,K,LP)
     
    215218         ALF4 (I)=ALF1(I)*ALF1Q(I)
    216219C
    217  110  CONTINUE
    218 C
    219       DO 111 JV=1,NTRA
    220       DO 1110 I=1,LON
    221 C
    222          IF(WGRI(I,K,L).LT.0.) THEN
     220      END DO
     221C
     222      DO JV=1,NTRA
     223      DO I=1,LON
     224C
     225         IF(WGRI(I,K,L)<0.) THEN
    223226C
    224227           F0 (I,L,JV)=ALF (I)* ( S0(I,K,LP,JV)-ALF1(I)*
     
    273276         ENDIF
    274277C
    275  1110 CONTINUE
    276  111  CONTINUE
    277 C
    278  11   CONTINUE
     278      END DO
     279      END DO
     280C
     281      END DO
    279282C
    280283C  puts the temporary moments Fi into appropriate neighboring boxes
    281284C
    282       DO 12 L=1,NIV-1
     285      DO L=1,NIV-1
    283286      LP=L+1
    284287C
    285       DO 120 I=1,LON
    286 C
    287          IF(WGRI(I,K,L).LT.0.) THEN
     288      DO I=1,LON
     289C
     290         IF(WGRI(I,K,L)<0.) THEN
    288291           SM(I,K,L)=SM(I,K,L)+FM(I,L)
    289292           ALF(I)=FM(I,L)/SM(I,K,L)
     
    299302         ALF3(I)=ALF1(I)-ALF(I)
    300303C
    301  120  CONTINUE
    302 C
    303       DO 121 JV=1,NTRA
    304       DO 1210 I=1,LON
    305 C
    306          IF(WGRI(I,K,L).LT.0.) THEN
     304      END DO
     305C
     306      DO JV=1,NTRA
     307      DO I=1,LON
     308C
     309         IF(WGRI(I,K,L)<0.) THEN
    307310C
    308311           TEMPTM=-ALF(I)*S0(I,K,L,JV)+ALF1(I)*F0(I,L,JV)
     
    342345         ENDIF
    343346C
    344  1210 CONTINUE
    345  121  CONTINUE
    346 C
    347  12   CONTINUE
     347      END DO
     348      END DO
     349C
     350      END DO
    348351C
    349352C  fin de la boucle principale sur les latitudes
    350353C
    351  1    CONTINUE
     354      END DO
    352355C
    353356      DO l = 1,llm
Note: See TracChangeset for help on using the changeset viewer.