Ignore:
Timestamp:
Jul 19, 2024, 7:54:50 PM (2 months ago)
Author:
abarral
Message:

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

Location:
LMDZ6/branches/Amaury_dev/libf/dyn3d
Files:
8 edited

Legend:

Unmodified
Added
Removed

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

    r2622 r5086  
    5757         deuxjour = 2. * daysec
    5858
    59          DO   1  ij   = 1, ip1jmp1
     59         DO   ij   = 1, ip1jmp1
    6060         unsaire2(ij) = unsaire(ij) * unsaire(ij)
    61    1     CONTINUE
     61      END DO
    6262      END IF
    6363
    … …  
    100100
    101101c
    102       DO 20 l = 1, llmm1
     102      DO l = 1, llmm1
    103103
    104104
    105105c       ......   calcul de  - w/2.    au niveau  l+1   .......
    106106
    107       DO 5   ij   = 1, ip1jmp1
     107      DO ij   = 1, ip1jmp1
    108108      wsur2( ij ) = - 0.5 * w( ij,l+1 )
    109    5  CONTINUE
     109      END DO
    110110
    111111
    112112c     .....................     calcul pour  du     ..................
    113113
    114       DO 6 ij = iip2 ,ip1jm-1
     114      DO ij = iip2 ,ip1jm-1
    115115      ww        = wsur2 (  ij  )     + wsur2( ij+1 )
    116116      uu        = 0.5 * ( ucov(ij,l) + ucov(ij,l+1) )
    117117      du(ij,l)  = du(ij,l)   - ww * ( uu - uav(ij, l ) )/massebx(ij, l )
    118118      du(ij,l+1)= du(ij,l+1) + ww * ( uu - uav(ij,l+1) )/massebx(ij,l+1)
    119    6  CONTINUE
     119      END DO
    120120
    121121c     .....  correction pour  du(iip1,j,l)  ........
    … …  
    123123
    124124CDIR$ IVDEP
    125       DO   7  ij   = iip1 +iip1, ip1jm, iip1
     125      DO   ij   = iip1 +iip1, ip1jm, iip1
    126126      du( ij, l  ) = du( ij -iim, l  )
    127127      du( ij,l+1 ) = du( ij -iim,l+1 )
    128    7  CONTINUE
     128      END DO
    129129
    130130c     .................    calcul pour   dv      .....................
    131131
    132       DO 8 ij = 1, ip1jm
     132      DO ij = 1, ip1jm
    133133      ww        = wsur2( ij+iip1 )   + wsur2( ij )
    134134      vv        = 0.5 * ( vcov(ij,l) + vcov(ij,l+1) )
    135135      dv(ij,l)  = dv(ij, l ) - ww * (vv - vav(ij, l ) )/masseby(ij, l )
    136136      dv(ij,l+1)= dv(ij,l+1) + ww * (vv - vav(ij,l+1) )/masseby(ij,l+1)
    137    8  CONTINUE
     137      END DO
    138138
    139139c
    … …  
    147147c                   ...............
    148148
    149         DO 15 ij = 1, ip1jmp1
     149        DO ij = 1, ip1jmp1
    150150         ww            = wsur2(ij) * (teta(ij,l) + teta(ij,l+1) )
    151151         dteta(ij, l ) = dteta(ij, l )  -  ww
    152152         dteta(ij,l+1) = dteta(ij,l+1)  +  ww
    153   15    CONTINUE
     153      END DO
    154154
    155155      IF( conser)  THEN
    156         DO 17 ij = 1,ip1jmp1
     156        DO ij = 1,ip1jmp1
    157157        ge(ij)   = wsur2(ij) * wsur2(ij) * unsaire2(ij)
    158   17    CONTINUE
     158      END DO
    159159        gt       = SSUM( ip1jmp1,ge,1 )
    160160        gtot(l)  = deuxjour * SQRT( gt/ip1jmp1 )
    161161      END IF
    162162
    163   20  CONTINUE
     163      END DO
    164164 
    165165      RETURN

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dteta1.F

    r2603 r5086  
    3232c
    3333
    34       DO 5 l = 1,llm
     34      DO l = 1,llm
    3535
    36       DO 1  ij = iip2, ip1jm - 1
     36      DO ij = iip2, ip1jm - 1
    3737      hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( teta(ij,l) + teta(ij+1,l) )
    38    1  CONTINUE
     38      END DO
    3939
    4040c    .... correction pour  hbxu(iip1,j,l)  .....
    … …  
    4242
    4343CDIR$ IVDEP
    44       DO 2 ij = iip1+ iip1, ip1jm, iip1
     44      DO ij = iip1+ iip1, ip1jm, iip1
    4545      hbxu( ij, l ) = hbxu( ij - iim, l )
    46    2  CONTINUE
     46      END DO
    4747
    4848
    49       DO 3 ij = 1,ip1jm
     49      DO ij = 1,ip1jm
    5050      hbyv(ij,l)= pbarv(ij, l)* 0.5 * ( teta(ij, l)+ teta(ij +iip1,l) )
    51    3  CONTINUE
     51      END DO
    5252
    53    5  CONTINUE
     53      END DO
    5454
    5555

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv1.F

    r1907 r5086  
    2727c
    2828c
    29       DO 10 l = 1,llm
     29      DO l = 1,llm
    3030c
    31       DO 2  ij = iip2, ip1jm - 1
     31      DO ij = iip2, ip1jm - 1
    3232      du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) *
    3333     *                    (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) +
    3434     *                        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
    35    2  CONTINUE
     35      END DO
    3636c
    37       DO 3 ij = 1, ip1jm - 1
     37      DO ij = 1, ip1jm - 1
    3838      dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) *
    3939     *                        (   pbaru(ij, l)  +  pbaru(ij+1   , l) +
    4040     *                       pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
    41    3  CONTINUE
     41      END DO
    4242c
    4343c    .... correction  pour  dv( 1,j,l )  .....
    … …  
    4545c
    4646CDIR$ IVDEP
    47       DO 4 ij = 1, ip1jm, iip1
     47      DO ij = 1, ip1jm, iip1
    4848      dv( ij,l ) = dv( ij + iim, l )
    49    4  CONTINUE
     49      END DO
    5050c
    51   10  CONTINUE
     51      END DO
    5252      RETURN
    5353      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F

    r2600 r5086  
    3434c
    3535c
    36       DO 5 l = 1,llm
     36      DO l = 1,llm
    3737c
    38       DO 2  ij  = iip2, ip1jm - 1
     38      DO ij  = iip2, ip1jm - 1
    3939       du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
    4040     * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
    41    2  CONTINUE
     41      END DO
    4242c
    4343c
    … …  
    4646c
    4747CDIR$ IVDEP
    48       DO 3 ij = iip1+ iip1, ip1jm, iip1
     48      DO ij = iip1+ iip1, ip1jm, iip1
    4949      du( ij,l ) = du( ij - iim,l )
    50    3  CONTINUE
     50      END DO
    5151c
    5252c
    53       DO 4 ij  = 1,ip1jm
     53      DO ij  = 1,ip1jm
    5454      dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
    5555     *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
    5656     *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
    57    4  CONTINUE
     57      END DO
    5858c
    59    5  CONTINUE
     59      END DO
    6060c
    6161      RETURN

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/iniinterp_horiz.F

    r5082 r5086  
    7474         a(i) = rlonuo(i-1)
    7575         b(i) =  rlonuo(i)
    76       end do
     76      END DO
    7777
    7878      d(1) = pi/2
    … …  
    8080         c(j) = rlatvo(j)
    8181         d(j+1) = rlatvo(j)
    82       end do
     82      END DO
    8383      c(jmo+1) = -pi/2
    8484     
    … …  
    9191         an(i) = rlonun(i-1)
    9292         bn(i) =  rlonun(i)
    93       end do
     93      END DO
    9494
    9595      dn(1) = pi/2
    … …  
    9797         cn(j) = rlatvn(j)
    9898         dn(j+1) = rlatvn(j)
    99       end do
     99      END DO
    100100      cn(jmn+1) = -pi/2
    101101
    … …  
    105105        do jj = 1,jmn+1
    106106               airen(ii,jj) = (bn(ii)-an(ii))*(sin(dn(jj))-sin(cn(jj)))
    107         end do
    108       end do
     107        END DO
     108      END DO
    109109
    110110c Calcul de la surface des intersections
    … …  
    151151                      intersec(ktotal)=(bb-aa)*(sin(dd)-sin(cc))
    152152                     end if
    153                 end do
    154                end do
     153                END DO
     154               END DO
    155155             end if
    156          end do
    157        end do       
     156         END DO
     157       END DO
    158158
    159159

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/interp_horiz.F

    r1907 r5086  
    7373        do ii=1, imn+1
    7474          varn(ii,jj,l) =0.
    75         end do
    76        end do
    77       end do
     75        END DO
     76       END DO
     77      END DO
    7878       
    7979c Interpolation horizontale
    … …  
    8888         varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l)
    8989     &        + varo(ik(k), jk(k),l)*intersec(k)/airen(iik(k),jjk(k))
    90         end do
    91       end do
     90        END DO
     91      END DO
    9292
    9393c Une seule valeur au pole pour les variables ! :
    … …  
    9999             totn = totn + varn(ii,1,l)
    100100             tots = tots + varn (ii,jmn+1,l)
    101            end do
     101           END DO
    102102           do ii =1, imn+1
    103103             varn(ii,1,l) = totn/REAL(imn+1)
    104104             varn(ii,jmn+1,l) = tots/REAL(imn+1)
    105            end do
    106        end do
     105           END DO
     106       END DO
    107107           
    108108
    … …  
    115115!!         do ii=1, imn+1
    116116!!           airetest(ii,jj) =0.
    117 !!         end do
    118 !!       end do
     117!!         END DO
     118!!       END DO
    119119!!       PRINT *, 'ktotal = ', ktotal
    120120!!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
    … …  
    136136!!              aire_ok = .false.
    137137!!          end if
    138 !!        end do
    139 !!       end do
     138!!        END DO
     139!!       END DO
    140140!! !      if (aire_ok) write(*,*) 'INTERP. HORIZ. : AIRES OK'
    141141!!  99   continue

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F

    r5082 r5086  
    6565c Eventuellement, faire l'extrapolation a partir des deux couches
    6666c les plus basses ou les deux couches les plus hautes:
    67       DO 130 i = 1, ilon
     67      DO i = 1, ilon
    6868cIM      IF ( ABS(pres-pgcm(i,ilev) ) .LT.
    6969         IF ( ABS(pres-pgcm(i,ilev) ) >
    … …  
    7777cIM   PRINT*,'i, ABS(pres-pgcm),ABS(pres-pgcm)',
    7878cIM  .i, ABS(pres-pgcm(i,ilev)),ABS(pres-pgcm(i,1))
    79   130 CONTINUE
    80       DO 150 k = 1, ilev-1
    81          DO 140 i = 1, ilon
     79      END DO
     80      DO k = 1, ilev-1
     81         DO i = 1, ilon
    8282            pbot = pgcm(i,k)
    8383            ptop = pgcm(i,k+1)
    … …  
    8787               lb(i) = k
    8888            ENDIF
    89   140    CONTINUE
    90   150 CONTINUE
     89      END DO
     90      END DO
    9191c
    9292c Interpolation lineaire:

  • LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j1.F

    r5082 r5086  
    6565c Eventuellement, faire l'extrapolation a partir des deux couches
    6666c les plus basses ou les deux couches les plus hautes:
    67       DO 130 i = 1, ilon
     67      DO i = 1, ilon
    6868cIM      IF ( ABS(pres-pgcm(i,ilev) ) .LT.
    6969         IF ( ABS(pres-pgcm(i,ilev) ) >
    … …  
    7777cIM   PRINT*,'i, ABS(pres-pgcm),ABS(pres-pgcm)',
    7878cIM  .i, ABS(pres-pgcm(i,ilev)),ABS(pres-pgcm(i,1))
    79   130 CONTINUE
    80       DO 150 k = 1, ilev-1
    81          DO 140 i = 1, ilon
     79      END DO
     80      DO k = 1, ilev-1
     81         DO i = 1, ilon
    8282            pbot = pgcm(i,k)
    8383            ptop = pgcm(i,k+1)
    … …  
    8787               lb(i) = k
    8888            ENDIF
    89   140    CONTINUE
    90   150 CONTINUE
     89      END DO
     90      END DO
    9191c
    9292c Interpolation lineaire:
Note: See TracChangeset for help on using the changeset viewer.