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

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

Location:
LMDZ6/branches/Amaury_dev/libf/dyn3d_common
Files:
29 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/bernoui.F

    r2603 r5086  
    4343c   ---------------------
    4444c
    45       DO 4 ijl = 1,ngrid*nlay
     45      DO ijl = 1,ngrid*nlay
    4646         pbern( ijl ) =  pphi( ijl ) + pecin( ijl )
    47    4  CONTINUE
     47      END DO
    4848c
    4949c-----------------------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.F

    r4593 r5086  
    3232      INCLUDE "comgeom.h"
    3333c
    34       DO 5 l = 1,nbniv
     34      DO l = 1,nbniv
    3535c
    36       DO 2  ij = iip2, ip1jm - 1
     36      DO ij = iip2, ip1jm - 1
    3737      convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   +
    3838     *                      yflu(ij +1,l ) - yflu( ij -iim,l )
    39    2  CONTINUE
     39      END DO
    4040c
    4141c
    … …  
    4545c
    4646CDIR$ IVDEP
    47       DO 3 ij = iip2,ip1jm,iip1
     47      DO ij = iip2,ip1jm,iip1
    4848      convfl( ij,l ) = convfl( ij + iim,l )
    49    3  CONTINUE
     49      END DO
    5050c
    5151c     ......  calcul aux poles  .......
    … …  
    5353      convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
    5454      convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
    55       DO 4 ij = 1,iip1
     55      DO ij = 1,iip1
    5656      convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
    5757      convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
    58    4  CONTINUE
     58      END DO
    5959c
    60    5  CONTINUE
     60      END DO
    6161      RETURN
    6262      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90

    r4257 r5086  
    320320     do l = 1, llm + 1
    321321        read(unit, fmt=*) ap(l), bp(l)
    322      end do
     322     END DO
    323323     close(unit)
    324324     call assert(ap(1) == 0., ap(llm + 1) == 0., bp(1) == 1., &

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F

    r5082 r5086  
    139139         do l=1,llm
    140140            read(99,*) zsig(l)
    141          end do
     141         END DO
    142142         CLOSE(99)
    143143
    … …  
    146146           sig(l) = 0.5 * ( exp(-zsig(l)/scaleheight) +
    147147     &                      exp(-zsig(l-1)/scaleheight) )
    148          end do
     148         END DO
    149149         sig(llm+1) =0
    150150
    … …  
    264264c          zsig(l)= zsig(l-1)-scaleheight*
    265265c    .    log((aps(l) + bps(l)*ps)/(aps(l-1) + bps(l-1)*ps))
    266 c       end do
     266c       END DO
    267267c       write(53,'(I3,50F10.5)') iz, zsig
    268 c      end do
     268c      END DO
    269269c      close(53)
    270270c     --------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.F

    r4593 r5086  
    4141c
    4242c
    43       DO 10 l = 1,klevel
     43      DO l = 1,klevel
    4444c
    4545        DO  ij = iip2, ip1jm - 1
    … …  
    7070         div( ij + ip1jm, l ) =   sumyps
    7171        ENDDO
    72   10  CONTINUE
     72      END DO
    7373c
    7474

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.F

    r4593 r5086  
    4444c
    4545c
    46       DO 10 l = 1,klevel
     46      DO l = 1,klevel
    4747c
    4848        DO  ij = iip2, ip1jm - 1
    … …  
    7474         div( ij + ip1jm, l ) =   sumyps
    7575        ENDDO
    76   10  CONTINUE
     76      END DO
    7777c
    7878

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.F

    r4593 r5086  
    4141c
    4242c
    43       DO 10 l = 1,klevel
     43      DO l = 1,klevel
    4444c
    4545        DO  ij = iip2, ip1jm - 1
    … …  
    7070         div( ij + ip1jm, l ) =   sumyps
    7171        ENDDO
    72   10  CONTINUE
     72      END DO
    7373c
    7474

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F

    r4593 r5086  
    3030c
    3131c
    32       DO 10 l = 1,klevel
     32      DO l = 1,klevel
    3333c
    34       DO 1 ij = iip2, ip1jm - 1
     34      DO ij = iip2, ip1jm - 1
    3535      div( ij + 1, l ) = x(ij+1,l) - x(ij,l)+ y(ij-iim,l)-y(ij+1,l)
    36    1  CONTINUE
     36      END DO
    3737c
    3838c     ....  correction pour  div( 1,j,l)  ......
    … …  
    4040c
    4141CDIR$ IVDEP
    42       DO 3 ij = iip2,ip1jm,iip1
     42      DO ij = iip2,ip1jm,iip1
    4343      div( ij,l ) = div( ij + iim,l )
    44    3  CONTINUE
     44      END DO
    4545c
    4646c     ....  calcul  aux poles  .....
    4747c
    4848c
    49       DO 5 i  = 1,iim
     49      DO i  = 1,iim
    5050      aiy1(i)= y(i,l)
    5151      aiy2(i)= y(i+ip1jmi1,l)
    52    5  CONTINUE
     52      END DO
    5353      sumypn = SSUM ( iim,aiy1,1 )
    5454      sumyps = SSUM ( iim,aiy2,1 )
    55       DO 7 i = 1,iip1
     55      DO i = 1,iip1
    5656      div(     i    , l ) = - sumypn/iim
    5757      div( i + ip1jm, l ) =   sumyps/iim
    58    7  CONTINUE
     58      END DO
    5959c
    60   10  CONTINUE
     60      END DO
    6161      RETURN
    6262      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.F

    r2603 r5086  
    3636      CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
    3737c
    38       DO 10 iter = 1,lh
     38      DO iter = 1,lh
    3939
    4040      CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1  )
    … …  
    4545      CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)
    4646
    47       DO 5 l = 1,klevel
    48       DO 4  ij = 1, ip1jmp1
     47      DO l = 1,klevel
     48      DO ij = 1, ip1jmp1
    4949      divgra( ij,l ) = - cdivh * divgra( ij,l )
    50    4  CONTINUE
    51    5  CONTINUE
     50      END DO
     51      END DO
    5252c
    53   10  CONTINUE
     53      END DO
    5454      RETURN
    5555      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90

    r2598 r5086  
    191191             do while (rlonm025(is2) < - pi .and. is2 < iim)
    192192                is2 = is2 + 1
    193              end do
     193             END DO
    194194
    195195             if (rlonm025(is2) < - pi) then
    … …  
    202202             do while (rlonm025(is2) > pi .and. is2 > 1)
    203203                is2 = is2 - 1
    204              end do
     204             END DO
    205205
    206206             if (rlonm025(is2) > pi) then

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/geopot.F

    r2600 r5086  
    4747c     calcul de phi au niveau 1 pres du sol  .....
    4848
    49       DO   1  ij  = 1, ngrid
     49      DO   ij  = 1, ngrid
    5050      phi( ij,1 ) = phis( ij ) + teta(ij,1) * ( pks(ij) - pk(ij,1) )
    51    1  CONTINUE
     51      END DO
    5252
    5353c     calcul de phi aux niveaux superieurs  .......

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.F

    r4593 r5086  
    2323c
    2424c
    25       DO 6 l = 1,klevel
     25      DO l = 1,klevel
    2626c
    27       DO 2  ij = 1, ip1jmp1 - 1
     27      DO ij = 1, ip1jmp1 - 1
    2828      pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
    29    2  CONTINUE
     29      END DO
    3030c
    3131c    .... correction pour  pgx(ip1,j,l)  ....
    3232c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
    3333CDIR$ IVDEP
    34       DO 3  ij = iip1, ip1jmp1, iip1
     34      DO ij = iip1, ip1jmp1, iip1
    3535      pgx( ij,l ) = pgx( ij -iim,l )
    36    3  CONTINUE
     36      END DO
    3737c
    38       DO 4 ij = 1,ip1jm
     38      DO ij = 1,ip1jm
    3939      pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
    40    4  CONTINUE
     40      END DO
    4141c
    42    6  CONTINUE
     42      END DO
    4343      RETURN
    4444      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.F

    r2603 r5086  
    3737      CALL SCOPY( ip1jm*klevel,  ycov,1,gdy,1 )
    3838c
    39       DO 10 iter = 1,ld
     39      DO iter = 1,ld
    4040c
    4141      CALL  diverg( klevel,  gdx , gdy, div          )
    … …  
    4343      CALL    grad( klevel,  div, gdx, gdy           )
    4444c
    45       DO 5  l = 1, klevel
    46       DO 3 ij = 1, ip1jmp1
     45      DO l = 1, klevel
     46      DO ij = 1, ip1jmp1
    4747      gdx( ij,l ) = - gdx( ij,l ) * cdivu
    48    3  CONTINUE
    49       DO 4 ij = 1, ip1jm
     48      END DO
     49      DO ij = 1, ip1jm
    5050      gdy( ij,l ) = - gdy( ij,l ) * cdivu
    51    4  CONTINUE
    52    5  CONTINUE
     51      END DO
     52      END DO
    5353c
    54   10  CONTINUE
     54      END DO
    5555      RETURN
    5656      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.F

    r5082 r5086  
    224224c
    225225        xo1 = 0.
    226         DO 10 iter = 1, itmax
     226        DO iter = 1, itmax
    227227        x1  = xo1
    228228        f   = x1+ alphax *SIN(x1-pxo)
    … …  
    232232        IF( xdm<=eps )GO TO 11
    233233        xo1 = x1
    234  10     CONTINUE
     234      END DO
    235235 11     CONTINUE
    236236c
    … …  
    241241C
    242242        yo1  = 0.
    243         DO 15 iter = 1,itmay
     243        DO iter = 1,itmay
    244244        y1   = yo1
    245245        f    = y1 + alphay* SIN(y1-pyo)
    … …  
    249249        IF(ydm<=eps) GO TO 17
    250250        yo1  = y1
    251  15     CONTINUE
     251      END DO
    252252c
    253253 17     CONTINUE
    … …  
    346346c
    347347c
    348       DO 35 j = 1, jjp1
     348      DO j = 1, jjp1
    349349c
    350350      IF ( j. eq. 1 )  THEN
    … …  
    356356      radclatm       = 0.5* rad * coslatm
    357357c
    358       DO 30 i = 1, iim
     358      DO i = 1, iim
    359359      xprp           = xprimp025( i )
    360360      xprm           = xprimm025( i )
    … …  
    365365      cvij2  ( i,1 ) = 0.5* rad * yprm
    366366      cvij3  ( i,1 ) = cvij2(i,1)
    367   30  CONTINUE
     367      END DO
    368368c
    369369      DO  i = 1, iim
    … …  
    387387      radclatp            = 0.5* rad * coslatp
    388388c
    389       DO 31 i = 1,iim
     389      DO i = 1,iim
    390390        xprp              = xprimp025( i )
    391391        xprm              = xprimm025( i )
    … …  
    396396        cvij1(i,jjp1)     = 0.5 * rad* yprp
    397397        cvij4(i,jjp1)     = cvij1(i,jjp1)
    398  31   CONTINUE
     398      END DO
    399399c
    400400       DO   i    = 1, iim
    … …  
    428428         ai14            = un4rad2 * coslatp * yprp
    429429         ai23            = un4rad2 * coslatm * yprm
    430          DO 32 i = 1,iim
     430         DO i = 1,iim
    431431         xprp            = xprimp025( i )
    432432         xprm            = xprimm025( i )
    … …  
    444444         cvij3   ( i,j ) = cvij2(i,j)
    445445         cvij4   ( i,j ) = cvij1(i,j)
    446   32     CONTINUE
     446      END DO
    447447c
    448448      END IF
    … …  
    463463         aireij4 (iip1,j) = aireij4 (1,j )
    464464       
    465   35  CONTINUE
     465      END DO
    466466c
    467467c    ..............................................................
    468468c
    469       DO 37 j = 1, jjp1
    470       DO 36 i = 1, iim
     469      DO j = 1, jjp1
     470      DO i = 1, iim
    471471      aire    ( i,j )  = aireij1(i,j) + aireij2(i,j) + aireij3(i,j) +
    472472     *                          aireij4(i,j)
    … …  
    479479      alpha2p3( i,j )  = alpha2 (i,j) + alpha3 (i,j)
    480480      alpha3p4( i,j )  = alpha3 (i,j) + alpha4 (i,j)
    481   36  CONTINUE
     481      END DO
    482482c
    483483c
    … …  
    491491      alpha2p3(iip1,j) = alpha2p3(1,j)
    492492      alpha3p4(iip1,j) = alpha3p4(1,j)
    493   37  CONTINUE
    494 c
    495 
    496       DO 42 j = 1,jjp1
    497       DO 41 i = 1,iim
     493      END DO
     494c
     495
     496      DO j = 1,jjp1
     497      DO i = 1,iim
    498498      aireu       (i,j)= aireij1(i,j) + aireij2(i,j) + aireij4(i+1,j) +
    499499     *                                aireij3(i+1,j)
    … …  
    502502      unsair_gam2( i,j)= unsaire(i,j)** ( - gamdi_h    )
    503503      airesurg   ( i,j)= aire(i,j)/ g
    504   41  CONTINUE
     504      END DO
    505505      aireu     (iip1,j)  = aireu  (1,j)
    506506      unsaire   (iip1,j)  = unsaire(1,j)
    … …  
    508508      unsair_gam2(iip1,j) = unsair_gam2(1,j)
    509509      airesurg   (iip1,j) = airesurg(1,j)
    510   42  CONTINUE
    511 c
    512 c
    513       DO 48 j = 1,jjm
     510      END DO
     511c
     512c
     513      DO j = 1,jjm
    514514c
    515515        DO i=1,iim
    … …  
    529529        unsairz_gam(iip1,j) = unsairz_gam(1,j)
    530530c
    531   48  CONTINUE
     531      END DO
    532532c
    533533c

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90

    r5082 r5086  
    270270          x0   = xxid(idat)
    271271          idat = idat + 1
    272        end do
     272       END DO
    273273       IF (xxim(imod)<xxid(idat)) THEN
    274274          dx   = xxim(imod) - x0
    … …  
    287287          idat = idat + 1
    288288       END IF
    289     end do
     289    END DO
    290290
    291291  END function inter_barx
    … …  
    339339          y0         = yjdat(jdat)
    340340          jdat       = jdat + 1
    341        end do
     341       END DO
    342342       IF (yjmod(jmod) < yjdat(jdat)) THEN
    343343          dy         = yjmod(jmod) - y0
    … …  
    357357          jdat       = jdat + 1
    358358       END IF
    359     end do
     359    END DO
    360360    ! Le test de fin suppose que l'interface 0 est commune aux deux
    361361    ! grilles "yjdat" et "yjmod".

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F

    r2622 r5086  
    4141
    4242c       CALCUL DE LA PRESSION DE SURFACE
    43 c       Les coefficients ap et bp sont passés en common
    44 c       Calcul de la pression au sol en mb optimisée pour
     43c       Les coefficients ap et bp sont pass�s en common
     44c       Calcul de la pression au sol en mb optimis�e pour
    4545c       la vectorialisation
    4646                   
    … …  
    6262             do i=1,iim
    6363                 psppm(i,j)=smass(i,j)/aire(i,j)*g*0.01
    64              end do
    65          end do                       
     64             END DO
     65         END DO
    6666       
    6767c RECONSTRUCTION DES CHAMPS CONTRAVARIANTS
    6868c Le programme ppm3d travaille avec les composantes
    69 c de vitesse et pas les flux, on doit donc passer de l'un à l'autre
    70 c Dans le même temps, on fait le changement d'orientation du vent en v
     69c de vitesse et pas les flux, on doit donc passer de l'un � l'autre
     70c Dans le m�me temps, on fait le changement d'orientation du vent en v
    7171      do l=1,llm
    7272          do j=1,jjm
    … …  
    9999     
    100100c INVERSION DES NIVEAUX
    101 c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
     101c le programme ppm3d travaille avec une 3�me coordonn�e invers�e par rapport
    102102c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
    103 c On passe donc des niveaux du LMDZ à ceux de Lin
     103c On passe donc des niveaux du LMDZ � ceux de Lin
    104104     
    105105      do l=1,llm+1

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90

    r2598 r5086  
    3939       do while (xfi < xf(it) .and. it >= 1)
    4040          it = it - 1
    41        end do
     41       END DO
    4242
    4343       ! Calcul de Xf(xvrai(i))

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.F

    r2597 r5086  
    8282c=======================================================================
    8383
    84       DO   100    l = 1 , llm
     84      DO   l = 1 , llm
    8585c
    8686        DO    ij     = 1, ip1jmp1
    … …  
    103103c       ENDDO
    104104       
    105 100   CONTINUE
     105      END DO
    106106c
    107107      RETURN

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad.F

    r4593 r5086  
    2222c
    2323c
    24       DO 10 l = 1,klevel
     24      DO l = 1,klevel
    2525c
    26       DO 1  ij = 2, ip1jm
     26      DO ij = 2, ip1jm
    2727      y( ij,l ) = (  rot( ij,l ) - rot( ij-1,l )  ) * cvsurcuv( ij )
    28    1  CONTINUE
     28      END DO
    2929c
    3030c    ..... correction pour  y ( 1,j,l )  ......
    … …  
    3232c    ....    y(1,j,l)= y(iip1,j,l) ....
    3333CDIR$ IVDEP
    34       DO 2  ij = 1, ip1jm, iip1
     34      DO ij = 1, ip1jm, iip1
    3535      y( ij,l ) = y( ij +iim,l )
    36    2  CONTINUE
     36      END DO
    3737c
    38       DO 4  ij = iip2,ip1jm
     38      DO ij = iip2,ip1jm
    3939      x( ij,l ) = (  rot( ij,l ) - rot( ij -iip1,l )  ) * cusurcvu( ij )
    40    4  CONTINUE
    41       DO 6 ij = 1,iip1
     40      END DO
     41      DO ij = 1,iip1
    4242      x(    ij    ,l ) = 0.
    4343      x( ij +ip1jm,l ) = 0.
    44    6  CONTINUE
     44      END DO
    4545c
    46   10  CONTINUE
     46      END DO
    4747      RETURN
    4848      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad_gam.F

    r4593 r5086  
    2121      INTEGER   l,ij
    2222c
    23       DO 10 l = 1,klevel
     23      DO l = 1,klevel
    2424c
    25       DO 1  ij = 2, ip1jm
     25      DO ij = 2, ip1jm
    2626      y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij )
    27    1  CONTINUE
     27      END DO
    2828c
    2929c    ..... correction pour  y ( 1,j,l )  ......
    … …  
    3131c    ....    y(1,j,l)= y(iip1,j,l) ....
    3232CDIR$ IVDEP
    33       DO 2  ij = 1, ip1jm, iip1
     33      DO ij = 1, ip1jm, iip1
    3434      y( ij,l ) = y( ij +iim,l )
    35    2  CONTINUE
     35      END DO
    3636c
    37       DO 4  ij = iip2,ip1jm
     37      DO ij = iip2,ip1jm
    3838      x( ij,l ) = (rot( ij,l ) - rot( ij -iip1,l )) * cuscvugam( ij )
    39    4  CONTINUE
    40       DO 6 ij = 1,iip1
     39      END DO
     40      DO ij = 1,iip1
    4141      x(    ij    ,l ) = 0.
    4242      x( ij +ip1jm,l ) = 0.
    43    6  CONTINUE
     43      END DO
    4444c
    45   10  CONTINUE
     45      END DO
    4646      RETURN
    4747      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.F

    r4593 r5086  
    2121      INTEGER l,ij
    2222c
    23       DO 10 l = 1,klevel
     23      DO l = 1,klevel
    2424c
    25       DO 1  ij = 2, ip1jm
     25      DO ij = 2, ip1jm
    2626      y(ij,l)=( rot(ij,l) - rot(ij-1,l))
    27    1  CONTINUE
     27      END DO
    2828c
    2929c    ..... correction pour  y ( 1,j,l )  ......
    … …  
    3131c    ....    y(1,j,l)= y(iip1,j,l) ....
    3232
    33       DO 2  ij = 1, ip1jm, iip1
     33      DO ij = 1, ip1jm, iip1
    3434      y( ij,l ) = y( ij +iim,l )
    35    2  CONTINUE
     35      END DO
    3636c
    37       DO 4  ij = iip2,ip1jm
     37      DO ij = iip2,ip1jm
    3838      x(ij,l)= rot(ij,l)-rot(ij-iip1,l)
    39    4  CONTINUE
    40       DO 6 ij = 1,iip1
     39      END DO
     40      DO ij = 1,iip1
    4141      x(    ij    ,l ) = 0.
    4242      x( ij +ip1jm,l ) = 0.
    43    6  CONTINUE
     43      END DO
    4444c
    45   10  CONTINUE
     45      END DO
    4646      RETURN
    4747      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrarot.F

    r2603 r5086  
    3636      CALL SCOPY (  ip1jm*klevel, ycov, 1, gry, 1 )
    3737c
    38       DO 10 iter = 1,lr
     38      DO iter = 1,lr
    3939      CALL  rotat (klevel,grx, gry, rot )
    4040      CALL filtreg( rot, jjm, klevel, 2,1, .false.,2)
    4141      CALL nxgrad (klevel,rot, grx, gry )
    4242c
    43       DO 5  l = 1, klevel
    44       DO 2 ij = 1, ip1jm
     43      DO l = 1, klevel
     44      DO ij = 1, ip1jm
    4545      gry( ij,l ) = - gry( ij,l ) * crot
    46    2  CONTINUE
    47       DO 3 ij = 1, ip1jmp1
     46      END DO
     47      DO ij = 1, ip1jmp1
    4848      grx( ij,l ) = - grx( ij,l ) * crot
    49    3  CONTINUE
    50    5  CONTINUE
     49      END DO
     50      END DO
    5151c
    52   10  CONTINUE
     52      END DO
    5353      RETURN
    5454      END

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pbar.F

    r4593 r5086  
    8787
    8888
    89       DO 1 ij = 1, ip1jmp1 - 1
     89      DO ij = 1, ip1jmp1 - 1
    9090      pbarx( ij ) = pext(ij) * alpha1p2(ij) + pext(ij+1)*alpha3p4(ij+1)
    91    1  CONTINUE
     91      END DO
    9292
    9393c    .... correction pour pbarx( iip1,j) .....
    … …  
    9595c    ...    pbarx(iip1,j)= pbarx(1,j) ...
    9696CDIR$ IVDEP
    97       DO 2 ij = iip1, ip1jmp1, iip1
     97      DO ij = iip1, ip1jmp1, iip1
    9898      pbarx( ij ) = pbarx( ij - iim )
    99    2  CONTINUE
     99      END DO
    100100
    101101
    102       DO 3 ij = 1,ip1jm
     102      DO ij = 1,ip1jm
    103103      pbary( ij ) = pext(   ij  )   * alpha2p3(   ij   )     +
    104104     *              pext( ij+iip1 ) * alpha1p4( ij+iip1 )
    105    3  CONTINUE
     105      END DO
    106106
    107107
    108       DO 5 ij = 1, ip1jm - 1
     108      DO ij = 1, ip1jm - 1
    109109      pbarxy( ij ) = pext(ij)*alpha2(ij) + pext(ij+1)*alpha3(ij+1) +
    110110     *   pext(ij+iip1)*alpha1(ij+iip1) + pext(ij+iip2)*alpha4(ij+iip2)
    111    5  CONTINUE
     111      END DO
    112112
    113113
    … …  
    116116CDIR$ IVDEP
    117117
    118       DO 7 ij = iip1, ip1jm, iip1
     118      DO ij = iip1, ip1jm, iip1
    119119      pbarxy( ij ) = pbarxy( ij - iim )
    120    7  CONTINUE
     120      END DO
    121121
    122122

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.F

    r5082 r5086  
    149149       call advxp( limit,0.5*dt,pbaru,sm,s0,sx,sy,sz
    150150     .             ,sxx,sxy,sxz,syy,syz,szz,1 )
    151         end do
     151        END DO
    152152        do l=1,llm
    153153        do i=1,iip1
    … …  
    206206       call advxp( limit,0.5*dt,pbaru,sm,s0,sx,sy,sz
    207207     .             ,sxx,sxy,sxz,syy,syz,szz,1 )
    208         end do
     208        END DO
    209209c---------------------------------------------------------
    210210c---------------------------------------------------------

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/psextbar.F

    r4593 r5086  
    8989
    9090
    91       DO     5     ij = 1, ip1jm - 1
     91      DO     ij = 1, ip1jm - 1
    9292      psexbarxy( ij ) = pext(ij)*alpha2(ij) + pext(ij+1)*alpha3(ij+1) +
    9393     *   pext(ij+iip1)*alpha1(ij+iip1) + pext(ij+iip2)*alpha4(ij+iip2)
    94    5  CONTINUE
     94      END DO
    9595
    9696
    … …  
    9999CDIR$ IVDEP
    100100
    101       DO 7 ij = iip1, ip1jm, iip1
     101      DO ij = iip1, ip1jm, iip1
    102102      psexbarxy( ij ) = psexbarxy( ij - iim )
    103    7  CONTINUE
     103      END DO
    104104
    105105

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat.F

    r4593 r5086  
    3030c
    3131c
    32       DO  10 l = 1,klevel
     32      DO  l = 1,klevel
    3333c
    3434        DO   ij = 1, ip1jm - 1
    … …  
    4444        ENDDO
    4545c
    46   10  CONTINUE
     46      END DO
    4747
    4848ccc        CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 )

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat_nfil.F

    r4593 r5086  
    3030c
    3131c
    32       DO  10 l = 1,klevel
     32      DO  l = 1,klevel
    3333c
    3434        DO   ij = 1, ip1jm - 1
    … …  
    4444        ENDDO
    4545c
    46   10  CONTINUE
     46      END DO
    4747
    4848      RETURN

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatf.F

    r4593 r5086  
    3131c
    3232c
    33       DO  10 l = 1,klevel
     33      DO  l = 1,klevel
    3434c
    3535        DO   ij = 1, ip1jm - 1
    … …  
    4545        ENDDO
    4646c
    47   10  CONTINUE
     47      END DO
    4848
    4949        CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 )

  • LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatst.F

    r4593 r5086  
    2424c
    2525c
    26       DO 5 l = 1,klevel
     26      DO l = 1,klevel
    2727c
    28       DO 1 ij = 1, ip1jm - 1
     28      DO ij = 1, ip1jm - 1
    2929      rot( ij,l )  =  (  y( ij+1 , l )  -  y( ij,l )   +
    3030     *                 x(ij +iip1, l )  -  x( ij,l )  )
    31    1  CONTINUE
     31      END DO
    3232c
    3333c    .... correction pour rot( iip1,j,l)  ....
    … …  
    3535c    ....   rot(iip1,j,l)= rot(1,j,l) ...
    3636CDIR$ IVDEP
    37       DO 2 ij = iip1, ip1jm, iip1
     37      DO ij = iip1, ip1jm, iip1
    3838      rot( ij,l ) = rot( ij -iim,l )
    39    2  CONTINUE
     39      END DO
    4040c
    41    5  CONTINUE
     41      END DO
    4242      RETURN
    4343      END
Note: See TracChangeset for help on using the changeset viewer.