Changeset 5077 for LMDZ6/trunk


Ignore:
Timestamp:
Jul 19, 2024, 11:01:18 AM (2 months ago)
Author:
abarral
Message:

Fix incorrect alignment of variables in COMMON block comdissip.h
Replace obsolete DO with shared termination in advx.F, advxp.F
(minor) replace obsolete bool operators

Location:
LMDZ6/trunk/libf/dyn3d_common
Files:
3 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/dyn3d_common/advx.F

    r2600 r5077  
    9595
    9696C  -------------------------------------
    97       DO 300 j = 1,jjp1
     97      DO j = 1,jjp1
    9898         NUM(j) = 1
    99   300 CONTINUE
     99      END DO
    100100      sqi = 0.
    101101      sqf = 0.
    … …  
    121121C  ugri est en kg/s
    122122
    123       DO 500 l = 1,llm
    124          DO 500 j = 1,jjm+1
    125             DO 500 i = 1,iip1 
     123      DO l = 1,llm
     124         DO j = 1,jjm+1
     125            DO i = 1,iip1
    126126C            ugri (i,j,llm+1-l) = pbaru (i,j,l) * ( dsig(l) / g )
    127127             ugri (i,j,llm+1-l) = pbaru (i,j,l)
    128   500 CONTINUE
     128      END DO
     129      END DO
     130      END DO
    129131
    130132
    … …  
    137139C  boucle principale sur les niveaux et les latitudes
    138140C
    139       DO 1 L=1,NIV
    140       DO 1 K=lati,latf
     141      DO L=1,NIV
     142      DO K=lati,latf
    141143C
    142144C  initialisation
    … …  
    144146C  program assumes periodic boundaries in X
    145147C
    146       DO 10 I=2,LON
     148      DO I=2,LON
    147149         SMNEW(I)=SM(I,K,L)+(UGRI(I-1,K,L)-UGRI(I,K,L))*DTX
    148  10   CONTINUE
     150      END DO
    149151      SMNEW(1)=SM(1,K,L)+(UGRI(LON,K,L)-UGRI(1,K,L))*DTX
    150152C
    … …  
    154156      LONK=LON/NUMK
    155157C
    156       IF(NUMK.GT.1) THEN
    157 C
    158       DO 111 I=1,LON
     158      IF(NUMK>1) THEN
     159C
     160      DO I=1,LON
    159161         TM(I)=0.
    160  111  CONTINUE
    161       DO 112 JV=1,NTRA
    162       DO 1120 I=1,LON
     162      END DO
     163      DO JV=1,NTRA
     164      DO I=1,LON
    163165         T0(I,JV)=0.
    164166         TX(I,JV)=0.
    165167         TY(I,JV)=0.
    166168         TZ(I,JV)=0.
    167  1120 CONTINUE
    168  112  CONTINUE
    169 C
    170       DO 11 I2=1,NUMK
    171 C
    172          DO 113 I=1,LONK
     169      END DO
     170      END DO
     171C
     172      DO I2=1,NUMK
     173C
     174         DO I=1,LONK
    173175            I3=(I-1)*NUMK+I2
    174176            TM(I)=TM(I)+SM(I3,K,L)
    175177            ALF(I)=SM(I3,K,L)/TM(I)
    176178            ALF1(I)=1.-ALF(I)
    177  113     CONTINUE
     179      END DO
    178180C
    179181         DO  JV=1,NTRA
    … …  
    190192         ENDDO
    191193C
    192  11   CONTINUE
     194      END DO
    193195C
    194196      ELSE
    195197C
    196       DO 115 I=1,LON
     198      DO I=1,LON
    197199         TM(I)=SM(I,K,L)
    198  115  CONTINUE
    199       DO 116 JV=1,NTRA
    200       DO 1160 I=1,LON
     200      END DO
     201      DO JV=1,NTRA
     202      DO I=1,LON
    201203         T0(I,JV)=S0(I,K,L,JV)
    202204         TX(I,JV)=sx(I,K,L,JV)
    203205         TY(I,JV)=sy(I,K,L,JV)
    204206         TZ(I,JV)=sz(I,K,L,JV)
    205  1160 CONTINUE
    206  116  CONTINUE
     207      END DO
     208      END DO
    207209C
    208210      ENDIF
    209211C
    210       DO 117 I=1,LONK
     212      DO I=1,LONK
    211213         UEXT(I)=UGRI(I*NUMK,K,L)
    212  117  CONTINUE
     214      END DO
    213215C
    214216C  place limits on appropriate moments before transport
    … …  
    217219      IF(.NOT.LIMIT) GO TO 13
    218220C
    219       DO 12 JV=1,NTRA
    220       DO 120 I=1,LONK
     221      DO JV=1,NTRA
     222      DO I=1,LONK
    221223        TX(I,JV)=SIGN(AMIN1(AMAX1(T0(I,JV),0.),ABS(TX(I,JV))),TX(I,JV))
    222  120  CONTINUE
    223  12   CONTINUE
     224      END DO
     225      END DO
    224226C
    225227 13   CONTINUE
    … …  
    231233C  flux from IP to I if U(I).lt.0
    232234C
    233       DO 140 I=1,LONK-1
    234          IF(UEXT(I).LT.0.) THEN
     235      DO I=1,LONK-1
     236         IF(UEXT(I)<0.) THEN
    235237           FM(I)=-UEXT(I)*DTX
    236238           ALF(I)=FM(I)/TM(I+1)
    237239           TM(I+1)=TM(I+1)-FM(I)
    238240         ENDIF
    239  140  CONTINUE
     241      END DO
    240242C
    241243      I=LONK
    242       IF(UEXT(I).LT.0.) THEN
     244      IF(UEXT(I)<0.) THEN
    243245        FM(I)=-UEXT(I)*DTX
    244246        ALF(I)=FM(I)/TM(1)
    … …  
    248250C  flux from I to IP if U(I).gt.0
    249251C
    250       DO 141 I=1,LONK
    251          IF(UEXT(I).GE.0.) THEN
     252      DO I=1,LONK
     253         IF(UEXT(I)>=0.) THEN
    252254           FM(I)=UEXT(I)*DTX
    253255           ALF(I)=FM(I)/TM(I)
    254256           TM(I)=TM(I)-FM(I)
    255257         ENDIF
    256  141  CONTINUE
    257 C
    258       DO 142 I=1,LONK
     258      END DO
     259C
     260      DO I=1,LONK
    259261         ALFQ(I)=ALF(I)*ALF(I)
    260262         ALF1(I)=1.-ALF(I)
    261263         ALF1Q(I)=ALF1(I)*ALF1(I)
    262  142  CONTINUE
    263 C
    264       DO 150 JV=1,NTRA
    265       DO 1500 I=1,LONK-1
    266 C
    267          IF(UEXT(I).LT.0.) THEN
     264      END DO
     265C
     266      DO JV=1,NTRA
     267      DO I=1,LONK-1
     268C
     269         IF(UEXT(I)<0.) THEN
    268270C
    269271           F0(I,JV)=ALF (I)* ( T0(I+1,JV)-ALF1(I)*TX(I+1,JV) )
    … …  
    279281         ENDIF
    280282C
    281  1500 CONTINUE
    282  150  CONTINUE
     283      END DO
     284      END DO
    283285C
    284286      I=LONK
    285       IF(UEXT(I).LT.0.) THEN
    286 C
    287         DO 151 JV=1,NTRA
     287      IF(UEXT(I)<0.) THEN
     288C
     289        DO JV=1,NTRA
    288290C
    289291           F0 (I,JV)=ALF (I)* ( T0(1,JV)-ALF1(I)*TX(1,JV) )
    … …  
    297299           TZ(1,JV)=TZ(1,JV)-FZ(I,JV)
    298300C
    299  151    CONTINUE
     301      END DO
    300302C
    301303      ENDIF
    302304C
    303       DO 152 JV=1,NTRA
    304       DO 1520 I=1,LONK
    305 C
    306          IF(UEXT(I).GE.0.) THEN
     305      DO JV=1,NTRA
     306      DO I=1,LONK
     307C
     308         IF(UEXT(I)>=0.) THEN
    307309C
    308310           F0(I,JV)=ALF (I)* ( T0(I,JV)+ALF1(I)*TX(I,JV) )
    … …  
    318320         ENDIF
    319321C
    320  1520 CONTINUE
    321  152  CONTINUE
     322      END DO
     323      END DO
    322324C
    323325C  puts the temporary moments Fi into appropriate neighboring boxes
    324326C
    325       DO 160 I=1,LONK
    326          IF(UEXT(I).LT.0.) THEN
     327      DO I=1,LONK
     328         IF(UEXT(I)<0.) THEN
    327329           TM(I)=TM(I)+FM(I)
    328330           ALF(I)=FM(I)/TM(I)
    329331         ENDIF
    330  160  CONTINUE
    331 C
    332       DO 161 I=1,LONK-1
    333          IF(UEXT(I).GE.0.) THEN
     332      END DO
     333C
     334      DO I=1,LONK-1
     335         IF(UEXT(I)>=0.) THEN
    334336           TM(I+1)=TM(I+1)+FM(I)
    335337           ALF(I)=FM(I)/TM(I+1)
    336338         ENDIF
    337  161  CONTINUE
     339      END DO
    338340C
    339341      I=LONK
    340       IF(UEXT(I).GE.0.) THEN
     342      IF(UEXT(I)>=0.) THEN
    341343        TM(1)=TM(1)+FM(I)
    342344        ALF(I)=FM(I)/TM(1)
    343345      ENDIF
    344346C
    345       DO 162 I=1,LONK
     347      DO I=1,LONK
    346348         ALF1(I)=1.-ALF(I)
    347  162  CONTINUE
    348 C
    349       DO 170 JV=1,NTRA
    350       DO 1700 I=1,LONK
    351 C
    352          IF(UEXT(I).LT.0.) THEN
     349      END DO
     350C
     351      DO JV=1,NTRA
     352      DO I=1,LONK
     353C
     354         IF(UEXT(I)<0.) THEN
    353355C
    354356           TEMPTM=-ALF(I)*T0(I,JV)+ALF1(I)*F0(I,JV)
    … …  
    360362         ENDIF
    361363C
    362  1700 CONTINUE
    363  170  CONTINUE
    364 C
    365       DO 171 JV=1,NTRA
    366       DO 1710 I=1,LONK-1
    367 C
    368          IF(UEXT(I).GE.0.) THEN
     364      END DO
     365      END DO
     366C
     367      DO JV=1,NTRA
     368      DO I=1,LONK-1
     369C
     370         IF(UEXT(I)>=0.) THEN
    369371C
    370372           TEMPTM=ALF(I)*T0(I+1,JV)-ALF1(I)*F0(I,JV)
    … …  
    376378         ENDIF
    377379C
    378  1710 CONTINUE
    379  171  CONTINUE
     380      END DO
     381      END DO
    380382C
    381383      I=LONK
    382       IF(UEXT(I).GE.0.) THEN
    383         DO 172 JV=1,NTRA
     384      IF(UEXT(I)>=0.) THEN
     385        DO JV=1,NTRA
    384386           TEMPTM=ALF(I)*T0(1,JV)-ALF1(I)*F0(I,JV)
    385387           T0(1,JV)=T0(1,JV)+F0(I,JV)
    … …  
    387389           TY(1,JV)=TY(1,JV)+FY(I,JV)
    388390           TZ(1,JV)=TZ(1,JV)+FZ(I,JV)
    389  172    CONTINUE
     391      END DO
    390392      ENDIF
    391393C
    392394C  retour aux mailles d'origine (passage des Tij aux Sij)
    393395C
    394       IF(NUMK.GT.1) THEN
    395 C
    396       DO 180 I2=1,NUMK
    397 C
    398          DO 180 I=1,LONK
     396      IF(NUMK>1) THEN
     397C
     398      DO I2=1,NUMK
     399C
     400         DO I=1,LONK
    399401C
    400402            I3=I2+(I-1)*NUMK
    … …  
    407409            ALF1Q(I)=ALF1(I)*ALF1(I)
    408410C
    409  180     CONTINUE
     411      END DO
     412      END DO
    410413C
    411414         DO  JV=1,NTRA
    … …  
    431434      ELSE
    432435C
    433       DO 190 I=1,LON
     436      DO I=1,LON
    434437         SM(I,K,L)=TM(I)
    435  190  CONTINUE
    436       DO 191 JV=1,NTRA
    437       DO 1910 I=1,LON
     438      END DO
     439      DO JV=1,NTRA
     440      DO I=1,LON
    438441         S0(I,K,L,JV)=T0(I,JV)
    439442         sx(I,K,L,JV)=TX(I,JV)
    440443         sy(I,K,L,JV)=TY(I,JV)
    441444         sz(I,K,L,JV)=TZ(I,JV)
    442  1910 CONTINUE
    443  191  CONTINUE
     445      END DO
     446      END DO
    444447C
    445448      ENDIF
    446449C
    447  1    CONTINUE
     450      END DO
     451      END DO
    448452C
    449453C ----------- AA Test en fin de ADVX ------ Controle des S*

  • LMDZ6/trunk/libf/dyn3d_common/advxp.F

    r2600 r5077  
    126126c test
    127127c  -------------------------------------
    128         DO 300 j =1,jjp1
     128        DO j =1,jjp1
    129129         NUM(j) =1
    130  300  CONTINUE
     130      END DO
    131131c       DO l=1,llm
    132132c      NUM(2,l)=6
    … …  
    150150C  ugri est en kg/s
    151151
    152        DO 500 l = 1,llm
    153        DO 500 j = 1,jjp1
    154        DO 500 i = 1,iip1
     152       DO l = 1,llm
     153       DO j = 1,jjp1
     154       DO i = 1,iip1
    155155       ugri (i,j,llm+1-l) =pbaru (i,j,l)
    156  500   CONTINUE
     156      END DO
     157      END DO
     158      END DO
    157159
    158160C  ---------------------------------------------------------
    … …  
    161163C  boucle principale sur les niveaux et les latitudes
    162164C     
    163       DO 1 L=1,NIV
    164       DO 1 K=lati,latf
     165      DO L=1,NIV
     166      DO K=lati,latf
    165167
    166168C
    … …  
    169171C  program assumes periodic boundaries in X
    170172C
    171       DO 10 I=2,LON
     173      DO I=2,LON
    172174         SMNEW(I)=SM(I,K,L)+(UGRI(I-1,K,L)-UGRI(I,K,L))*DTX
    173  10   CONTINUE
     175      END DO
    174176      SMNEW(1)=SM(1,K,L)+(UGRI(LON,K,L)-UGRI(1,K,L))*DTX
    175177C
    … …  
    179181      LONK=LON/NUMK
    180182C
    181       IF(NUMK.GT.1) THEN
    182 C
    183       DO 111 I=1,LON
     183      IF(NUMK>1) THEN
     184C
     185      DO I=1,LON
    184186         TM(I)=0.
    185  111  CONTINUE
    186       DO 112 JV=1,NTRA
    187       DO 1120 I=1,LON
     187      END DO
     188      DO JV=1,NTRA
     189      DO I=1,LON
    188190         T0 (I,JV)=0.
    189191         TX (I,JV)=0.
    … …  
    196198         TYZ(I,JV)=0.
    197199         TZZ(I,JV)=0.
    198  1120 CONTINUE
    199  112  CONTINUE
    200 C
    201       DO 11 I2=1,NUMK
    202 C
    203          DO 113 I=1,LONK
     200      END DO
     201      END DO
     202C
     203      DO I2=1,NUMK
     204C
     205         DO I=1,LONK
    204206            I3=(I-1)*NUMK+I2
    205207            TM(I)=TM(I)+SM(I3,K,L)
    … …  
    210212            ALF2(I)=ALF1(I)-ALF(I)
    211213            ALF3(I)=ALF(I)*ALF1(I)
    212  113     CONTINUE
    213 C
    214          DO 114 JV=1,NTRA
    215          DO 1140 I=1,LONK
     214      END DO
     215C
     216         DO JV=1,NTRA
     217         DO I=1,LONK
    216218            I3=(I-1)*NUMK+I2
    217219            TEMPTM=-ALF(I)*T0(I,JV)+ALF1(I)*S0(I3,K,L,JV)
    … …  
    229231            TYZ(I,JV)=TYZ(I,JV)+SYZ(I3,K,L,JV)
    230232            TZZ(I,JV)=TZZ(I,JV)+SZZ(I3,K,L,JV)
    231  1140    CONTINUE
    232  114     CONTINUE
    233 C
    234  11   CONTINUE
     233      END DO
     234      END DO
     235C
     236      END DO
    235237C
    236238      ELSE
    237239C
    238       DO 115 I=1,LON
     240      DO I=1,LON
    239241         TM(I)=SM(I,K,L)
    240  115  CONTINUE
    241       DO 116 JV=1,NTRA
    242       DO 1160 I=1,LON
     242      END DO
     243      DO JV=1,NTRA
     244      DO I=1,LON
    243245         T0 (I,JV)=S0 (I,K,L,JV)
    244246         TX (I,JV)=SSX (I,K,L,JV)
    … …  
    251253         TYZ(I,JV)=SYZ(I,K,L,JV)
    252254         TZZ(I,JV)=SZZ(I,K,L,JV)
    253  1160 CONTINUE
    254  116  CONTINUE
     255      END DO
     256      END DO
    255257C
    256258      ENDIF
    257259C
    258       DO 117 I=1,LONK
     260      DO I=1,LONK
    259261         UEXT(I)=UGRI(I*NUMK,K,L)
    260  117  CONTINUE
     262      END DO
    261263C
    262264C  place limits on appropriate moments before transport
    … …  
    265267      IF(.NOT.LIMIT) GO TO 13
    266268C
    267       DO 12 JV=1,NTRA
    268       DO 120 I=1,LONK
    269         IF(T0(I,JV).GT.0.) THEN
     269      DO JV=1,NTRA
     270      DO I=1,LONK
     271        IF(T0(I,JV)>0.) THEN
    270272          SLPMAX=T0(I,JV)
    271273          S1MAX=1.5*SLPMAX
    … …  
    283285          TXZ(I,JV)=0.
    284286        ENDIF
    285  120  CONTINUE
    286  12   CONTINUE
     287      END DO
     288      END DO
    287289C
    288290 13   CONTINUE
    … …  
    294296C  flux from IP to I if U(I).lt.0
    295297C
    296       DO 140 I=1,LONK-1
    297          IF(UEXT(I).LT.0.) THEN
     298      DO I=1,LONK-1
     299         IF(UEXT(I)<0.) THEN
    298300           FM(I)=-UEXT(I)*DTX
    299301           ALF(I)=FM(I)/TM(I+1)
    300302           TM(I+1)=TM(I+1)-FM(I)
    301303         ENDIF
    302  140  CONTINUE
     304      END DO
    303305C
    304306      I=LONK
    305       IF(UEXT(I).LT.0.) THEN
     307      IF(UEXT(I)<0.) THEN
    306308        FM(I)=-UEXT(I)*DTX
    307309        ALF(I)=FM(I)/TM(1)
    … …  
    311313C  flux from I to IP if U(I).gt.0
    312314C
    313       DO 141 I=1,LONK
    314          IF(UEXT(I).GE.0.) THEN
     315      DO I=1,LONK
     316         IF(UEXT(I)>=0.) THEN
    315317           FM(I)=UEXT(I)*DTX
    316318           ALF(I)=FM(I)/TM(I)
    317319           TM(I)=TM(I)-FM(I)
    318320         ENDIF
    319  141  CONTINUE
    320 C
    321       DO 142 I=1,LONK
     321      END DO
     322C
     323      DO I=1,LONK
    322324         ALFQ(I)=ALF(I)*ALF(I)
    323325         ALF1(I)=1.-ALF(I)
    … …  
    326328         ALF3(I)=ALF(I)*ALFQ(I)
    327329         ALF4(I)=ALF1(I)*ALF1Q(I)
    328  142  CONTINUE
    329 C
    330       DO 150 JV=1,NTRA
    331       DO 1500 I=1,LONK-1
    332 C
    333          IF(UEXT(I).LT.0.) THEN
     330      END DO
     331C
     332      DO JV=1,NTRA
     333      DO I=1,LONK-1
     334C
     335         IF(UEXT(I)<0.) THEN
    334336C
    335337           F0 (I,JV)=ALF (I)* ( T0(I+1,JV)-ALF1(I)*
    … …  
    358360         ENDIF
    359361C
    360  1500 CONTINUE
    361  150  CONTINUE
     362      END DO
     363      END DO
    362364C
    363365      I=LONK
    364       IF(UEXT(I).LT.0.) THEN
    365 C
    366         DO 151 JV=1,NTRA
     366      IF(UEXT(I)<0.) THEN
     367C
     368        DO JV=1,NTRA
    367369C
    368370           F0 (I,JV)=ALF (I)* ( T0(1,JV)-ALF1(I)*
    … …  
    389391           TXZ(1,JV)=ALF1Q(I)*TXZ(1,JV)
    390392C
    391  151    CONTINUE
     393      END DO
    392394C
    393395      ENDIF
    394396C
    395       DO 152 JV=1,NTRA
    396       DO 1520 I=1,LONK
    397 C
    398          IF(UEXT(I).GE.0.) THEN
     397      DO JV=1,NTRA
     398      DO I=1,LONK
     399C
     400         IF(UEXT(I)>=0.) THEN
    399401C
    400402           F0 (I,JV)=ALF (I)* ( T0(I,JV)+ALF1(I)*
    … …  
    423425         ENDIF
    424426C
    425  1520 CONTINUE
    426  152  CONTINUE
     427      END DO
     428      END DO
    427429C
    428430C  puts the temporary moments Fi into appropriate neighboring boxes
    429431C
    430       DO 160 I=1,LONK
    431          IF(UEXT(I).LT.0.) THEN
     432      DO I=1,LONK
     433         IF(UEXT(I)<0.) THEN
    432434           TM(I)=TM(I)+FM(I)
    433435           ALF(I)=FM(I)/TM(I)
    434436         ENDIF
    435  160  CONTINUE
    436 C
    437       DO 161 I=1,LONK-1
    438          IF(UEXT(I).GE.0.) THEN
     437      END DO
     438C
     439      DO I=1,LONK-1
     440         IF(UEXT(I)>=0.) THEN
    439441           TM(I+1)=TM(I+1)+FM(I)
    440442           ALF(I)=FM(I)/TM(I+1)
    441443         ENDIF
    442  161  CONTINUE
     444      END DO
    443445C
    444446      I=LONK
    445       IF(UEXT(I).GE.0.) THEN
     447      IF(UEXT(I)>=0.) THEN
    446448        TM(1)=TM(1)+FM(I)
    447449        ALF(I)=FM(I)/TM(1)
    448450      ENDIF
    449451C
    450       DO 162 I=1,LONK
     452      DO I=1,LONK
    451453         ALF1(I)=1.-ALF(I)
    452454         ALFQ(I)=ALF(I)*ALF(I)
    … …  
    454456         ALF2(I)=ALF1(I)-ALF(I)
    455457         ALF3(I)=ALF(I)*ALF1(I)
    456  162  CONTINUE
    457 C
    458       DO 170 JV=1,NTRA
    459       DO 1700 I=1,LONK
    460 C
    461          IF(UEXT(I).LT.0.) THEN
     458      END DO
     459C
     460      DO JV=1,NTRA
     461      DO I=1,LONK
     462C
     463         IF(UEXT(I)<0.) THEN
    462464C
    463465           TEMPTM=-ALF(I)*T0(I,JV)+ALF1(I)*F0(I,JV)
    … …  
    478480         ENDIF
    479481C
    480  1700 CONTINUE
    481  170  CONTINUE
    482 C
    483       DO 171 JV=1,NTRA
    484       DO 1710 I=1,LONK-1
    485 C
    486          IF(UEXT(I).GE.0.) THEN
     482      END DO
     483      END DO
     484C
     485      DO JV=1,NTRA
     486      DO I=1,LONK-1
     487C
     488         IF(UEXT(I)>=0.) THEN
    487489C
    488490           TEMPTM=ALF(I)*T0(I+1,JV)-ALF1(I)*F0(I,JV)
    … …  
    503505         ENDIF
    504506C
    505  1710 CONTINUE
    506  171  CONTINUE
     507      END DO
     508      END DO
    507509C
    508510      I=LONK
    509       IF(UEXT(I).GE.0.) THEN
    510         DO 172 JV=1,NTRA
     511      IF(UEXT(I)>=0.) THEN
     512        DO JV=1,NTRA
    511513           TEMPTM=ALF(I)*T0(1,JV)-ALF1(I)*F0(I,JV)
    512514           T0 (1,JV)=T0(1,JV)+F0(I,JV)
    … …  
    523525           TYZ(1,JV)=TYZ(1,JV)+FYZ(I,JV)
    524526           TZZ(1,JV)=TZZ(1,JV)+FZZ(I,JV)
    525  172    CONTINUE
     527      END DO
    526528      ENDIF
    527529C
    528530C  retour aux mailles d'origine (passage des Tij aux Sij)
    529531C
    530       IF(NUMK.GT.1) THEN
    531 C
    532       DO 18 I2=1,NUMK
    533 C
    534          DO 180 I=1,LONK
     532      IF(NUMK>1) THEN
     533C
     534      DO I2=1,NUMK
     535C
     536         DO I=1,LONK
    535537C
    536538            I3=I2+(I-1)*NUMK
    … …  
    546548            ALF4(I)=ALF1(I)*ALF1Q(I)
    547549C
    548  180     CONTINUE
    549 C
    550          DO 181 JV=1,NTRA
    551          DO 181 I=1,LONK
     550      END DO
     551C
     552         DO JV=1,NTRA
     553         DO I=1,LONK
    552554C
    553555            I3=I2+(I-1)*NUMK
    … …  
    577579            TXZ(I,JV)=ALF1Q(I)*TXZ(I,JV)
    578580C
    579  181     CONTINUE
    580 C
    581  18   CONTINUE
     581      END DO
     582      END DO
     583C
     584      END DO
    582585C
    583586      ELSE
    584587C
    585       DO 190 I=1,LON
     588      DO I=1,LON
    586589         SM(I,K,L)=TM(I)
    587  190  CONTINUE
    588       DO 191 JV=1,NTRA
    589       DO 1910 I=1,LON
     590      END DO
     591      DO JV=1,NTRA
     592      DO I=1,LON
    590593         S0 (I,K,L,JV)=T0 (I,JV)
    591594         SSX (I,K,L,JV)=TX (I,JV)
    … …  
    598601         SYZ(I,K,L,JV)=TYZ(I,JV)
    599602         SZZ(I,K,L,JV)=TZZ(I,JV)
    600  1910 CONTINUE
    601  191  CONTINUE
     603      END DO
     604      END DO
    602605C
    603606      ENDIF
    604607C
    605  1    CONTINUE
     608      END DO
     609      END DO
    606610C
    607611C ----------- AA Test en fin de ADVX ------ Controle des S*

  • LMDZ6/trunk/libf/dyn3d_common/comdissip.h

    r1945 r5077  
    66
    77      COMMON/comdissip/                                                 &
    8      &    niterdis,coefdis,tetavel,tetatemp,gamdissip
     8     &    coefdis,tetavel,tetatemp,gamdissip,niterdis
    99
    1010
Note: See TracChangeset for help on using the changeset viewer.