Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/advect.F	(revision 5086)
@@ -57,7 +57,7 @@
          deuxjour = 2. * daysec
 
-         DO   1  ij   = 1, ip1jmp1
+         DO   ij   = 1, ip1jmp1
          unsaire2(ij) = unsaire(ij) * unsaire(ij)
-   1     CONTINUE
+      END DO
       END IF
 
@@ -100,22 +100,22 @@
 
 c
-      DO 20 l = 1, llmm1
+      DO l = 1, llmm1
 
 
 c       ......   calcul de  - w/2.    au niveau  l+1   .......
 
-      DO 5   ij   = 1, ip1jmp1
+      DO ij   = 1, ip1jmp1
       wsur2( ij ) = - 0.5 * w( ij,l+1 )
-   5  CONTINUE
+      END DO
 
 
 c     .....................     calcul pour  du     ..................
 
-      DO 6 ij = iip2 ,ip1jm-1
+      DO ij = iip2 ,ip1jm-1
       ww        = wsur2 (  ij  )     + wsur2( ij+1 ) 
       uu        = 0.5 * ( ucov(ij,l) + ucov(ij,l+1) )
       du(ij,l)  = du(ij,l)   - ww * ( uu - uav(ij, l ) )/massebx(ij, l )
       du(ij,l+1)= du(ij,l+1) + ww * ( uu - uav(ij,l+1) )/massebx(ij,l+1)
-   6  CONTINUE
+      END DO
 
 c     .....  correction pour  du(iip1,j,l)  ........
@@ -123,17 +123,17 @@
 
 CDIR$ IVDEP
-      DO   7  ij   = iip1 +iip1, ip1jm, iip1
+      DO   ij   = iip1 +iip1, ip1jm, iip1
       du( ij, l  ) = du( ij -iim, l  )
       du( ij,l+1 ) = du( ij -iim,l+1 )
-   7  CONTINUE
+      END DO
 
 c     .................    calcul pour   dv      .....................
 
-      DO 8 ij = 1, ip1jm
+      DO ij = 1, ip1jm
       ww        = wsur2( ij+iip1 )   + wsur2( ij )
       vv        = 0.5 * ( vcov(ij,l) + vcov(ij,l+1) )
       dv(ij,l)  = dv(ij, l ) - ww * (vv - vav(ij, l ) )/masseby(ij, l )
       dv(ij,l+1)= dv(ij,l+1) + ww * (vv - vav(ij,l+1) )/masseby(ij,l+1)
-   8  CONTINUE
+      END DO
 
 c
@@ -147,19 +147,19 @@
 c                   ...............
 
-        DO 15 ij = 1, ip1jmp1
+        DO ij = 1, ip1jmp1
          ww            = wsur2(ij) * (teta(ij,l) + teta(ij,l+1) )
          dteta(ij, l ) = dteta(ij, l )  -  ww
          dteta(ij,l+1) = dteta(ij,l+1)  +  ww
-  15    CONTINUE
+      END DO
 
       IF( conser)  THEN
-        DO 17 ij = 1,ip1jmp1
+        DO ij = 1,ip1jmp1
         ge(ij)   = wsur2(ij) * wsur2(ij) * unsaire2(ij)
-  17    CONTINUE
+      END DO
         gt       = SSUM( ip1jmp1,ge,1 )
         gtot(l)  = deuxjour * SQRT( gt/ip1jmp1 )
       END IF
 
-  20  CONTINUE
+      END DO
  
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/dteta1.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/dteta1.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/dteta1.F	(revision 5086)
@@ -32,9 +32,9 @@
 c
 
-      DO 5 l = 1,llm
+      DO l = 1,llm
 
-      DO 1  ij = iip2, ip1jm - 1
+      DO ij = iip2, ip1jm - 1
       hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( teta(ij,l) + teta(ij+1,l) )
-   1  CONTINUE
+      END DO
 
 c    .... correction pour  hbxu(iip1,j,l)  .....
@@ -42,14 +42,14 @@
 
 CDIR$ IVDEP
-      DO 2 ij = iip1+ iip1, ip1jm, iip1
+      DO ij = iip1+ iip1, ip1jm, iip1
       hbxu( ij, l ) = hbxu( ij - iim, l )
-   2  CONTINUE
+      END DO
 
 
-      DO 3 ij = 1,ip1jm
+      DO ij = 1,ip1jm
       hbyv(ij,l)= pbarv(ij, l)* 0.5 * ( teta(ij, l)+ teta(ij +iip1,l) )
-   3  CONTINUE
+      END DO
 
-   5  CONTINUE
+      END DO
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv1.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv1.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv1.F	(revision 5086)
@@ -27,17 +27,17 @@
 c
 c
-      DO 10 l = 1,llm
+      DO l = 1,llm
 c
-      DO 2  ij = iip2, ip1jm - 1
+      DO ij = iip2, ip1jm - 1
       du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) *
      *                    (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) +
      *                        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
-   2  CONTINUE
+      END DO
 c
-      DO 3 ij = 1, ip1jm - 1
+      DO ij = 1, ip1jm - 1
       dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) *
      *                        (   pbaru(ij, l)  +  pbaru(ij+1   , l) +
      *                       pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
-   3  CONTINUE
+      END DO
 c
 c    .... correction  pour  dv( 1,j,l )  .....
@@ -45,9 +45,9 @@
 c
 CDIR$ IVDEP
-      DO 4 ij = 1, ip1jm, iip1
+      DO ij = 1, ip1jm, iip1
       dv( ij,l ) = dv( ij + iim, l )
-   4  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/dudv2.F	(revision 5086)
@@ -34,10 +34,10 @@
 c
 c
-      DO 5 l = 1,llm
+      DO l = 1,llm
 c
-      DO 2  ij  = iip2, ip1jm - 1
+      DO ij  = iip2, ip1jm - 1
        du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
      * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
-   2  CONTINUE
+      END DO
 c
 c
@@ -46,16 +46,16 @@
 c
 CDIR$ IVDEP
-      DO 3 ij = iip1+ iip1, ip1jm, iip1
+      DO ij = iip1+ iip1, ip1jm, iip1
       du( ij,l ) = du( ij - iim,l )
-   3  CONTINUE
+      END DO
 c
 c
-      DO 4 ij  = 1,ip1jm
+      DO ij  = 1,ip1jm
       dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
      *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
      *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
-   4  CONTINUE
+      END DO
 c
-   5  CONTINUE
+      END DO
 c
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/iniinterp_horiz.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/iniinterp_horiz.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/iniinterp_horiz.F	(revision 5086)
@@ -74,5 +74,5 @@
          a(i) = rlonuo(i-1)
          b(i) =  rlonuo(i)
-      end do
+      END DO
 
       d(1) = pi/2 
@@ -80,5 +80,5 @@
          c(j) = rlatvo(j) 
          d(j+1) = rlatvo(j)
-      end do
+      END DO
       c(jmo+1) = -pi/2 
       
@@ -91,5 +91,5 @@
          an(i) = rlonun(i-1)
          bn(i) =  rlonun(i)
-      end do
+      END DO
 
       dn(1) = pi/2 
@@ -97,5 +97,5 @@
          cn(j) = rlatvn(j)
          dn(j+1) = rlatvn(j)
-      end do
+      END DO
       cn(jmn+1) = -pi/2 
 
@@ -105,6 +105,6 @@
         do jj = 1,jmn+1
                airen(ii,jj) = (bn(ii)-an(ii))*(sin(dn(jj))-sin(cn(jj)))
-        end do
-      end do
+        END DO
+      END DO
 
 c Calcul de la surface des intersections
@@ -151,9 +151,9 @@
                       intersec(ktotal)=(bb-aa)*(sin(dd)-sin(cc))
                      end if
-                end do
-               end do
+                END DO
+               END DO
              end if
-         end do
-       end do       
+         END DO
+       END DO
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/interp_horiz.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/interp_horiz.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/interp_horiz.F	(revision 5086)
@@ -73,7 +73,7 @@
         do ii=1, imn+1
           varn(ii,jj,l) =0.
-        end do
-       end do
-      end do 
+        END DO
+       END DO
+      END DO
        
 c Interpolation horizontale
@@ -88,6 +88,6 @@
          varn(iik(k),jjk(k),l) = varn(iik(k),jjk(k),l) 
      &        + varo(ik(k), jk(k),l)*intersec(k)/airen(iik(k),jjk(k))
-        end do
-      end do
+        END DO
+      END DO
 
 c Une seule valeur au pole pour les variables ! :
@@ -99,10 +99,10 @@
              totn = totn + varn(ii,1,l)
              tots = tots + varn (ii,jmn+1,l)
-           end do 
+           END DO
            do ii =1, imn+1
              varn(ii,1,l) = totn/REAL(imn+1)
              varn(ii,jmn+1,l) = tots/REAL(imn+1)
-           end do 
-       end do
+           END DO
+       END DO
            
 
@@ -115,6 +115,6 @@
 !!         do ii=1, imn+1
 !!           airetest(ii,jj) =0.
-!!         end do
-!!       end do 
+!!         END DO
+!!       END DO
 !!       PRINT *, 'ktotal = ', ktotal
 !!       PRINT *, 'jmn+1 =', jmn+1, 'imn+1', imn+1
@@ -136,6 +136,6 @@
 !!              aire_ok = .false.
 !!          end if
-!!        end do
-!!       end do
+!!        END DO
+!!       END DO
 !! !      if (aire_ok) write(*,*) 'INTERP. HORIZ. : AIRES OK'
 !!  99   continue
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F	(revision 5086)
@@ -65,5 +65,5 @@
 c Eventuellement, faire l'extrapolation a partir des deux couches
 c les plus basses ou les deux couches les plus hautes:
-      DO 130 i = 1, ilon
+      DO i = 1, ilon
 cIM      IF ( ABS(pres-pgcm(i,ilev) ) .LT.
          IF ( ABS(pres-pgcm(i,ilev) ) >
@@ -77,7 +77,7 @@
 cIM   PRINT*,'i, ABS(pres-pgcm),ABS(pres-pgcm)',
 cIM  .i, ABS(pres-pgcm(i,ilev)),ABS(pres-pgcm(i,1))
-  130 CONTINUE
-      DO 150 k = 1, ilev-1
-         DO 140 i = 1, ilon
+      END DO
+      DO k = 1, ilev-1
+         DO i = 1, ilon
             pbot = pgcm(i,k)
             ptop = pgcm(i,k+1)
@@ -87,6 +87,6 @@
                lb(i) = k
             ENDIF
-  140    CONTINUE
-  150 CONTINUE
+      END DO
+      END DO
 c
 c Interpolation lineaire:
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j1.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j1.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j1.F	(revision 5086)
@@ -65,5 +65,5 @@
 c Eventuellement, faire l'extrapolation a partir des deux couches
 c les plus basses ou les deux couches les plus hautes:
-      DO 130 i = 1, ilon
+      DO i = 1, ilon
 cIM      IF ( ABS(pres-pgcm(i,ilev) ) .LT.
          IF ( ABS(pres-pgcm(i,ilev) ) >
@@ -77,7 +77,7 @@
 cIM   PRINT*,'i, ABS(pres-pgcm),ABS(pres-pgcm)',
 cIM  .i, ABS(pres-pgcm(i,ilev)),ABS(pres-pgcm(i,1))
-  130 CONTINUE
-      DO 150 k = 1, ilev-1
-         DO 140 i = 1, ilon
+      END DO
+      DO k = 1, ilev-1
+         DO i = 1, ilon
             pbot = pgcm(i,k)
             ptop = pgcm(i,k+1)
@@ -87,6 +87,6 @@
                lb(i) = k
             ENDIF
-  140    CONTINUE
-  150 CONTINUE
+      END DO
+      END DO
 c
 c Interpolation lineaire:
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/bernoui.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/bernoui.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/bernoui.F	(revision 5086)
@@ -43,7 +43,7 @@
 c   ---------------------
 c
-      DO 4 ijl = 1,ngrid*nlay
+      DO ijl = 1,ngrid*nlay
          pbern( ijl ) =  pphi( ijl ) + pecin( ijl )
-   4  CONTINUE
+      END DO
 c
 c-----------------------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/convflu.F	(revision 5086)
@@ -32,10 +32,10 @@
       INCLUDE "comgeom.h"
 c
-      DO 5 l = 1,nbniv
+      DO l = 1,nbniv
 c
-      DO 2  ij = iip2, ip1jm - 1
+      DO ij = iip2, ip1jm - 1
       convfl( ij + 1,l ) =  xflu(   ij,l ) - xflu( ij +  1,l )   +
      *                      yflu(ij +1,l ) - yflu( ij -iim,l )
-   2  CONTINUE
+      END DO
 c
 c
@@ -45,7 +45,7 @@
 c
 CDIR$ IVDEP
-      DO 3 ij = iip2,ip1jm,iip1
+      DO ij = iip2,ip1jm,iip1
       convfl( ij,l ) = convfl( ij + iim,l )
-   3  CONTINUE
+      END DO
 c
 c     ......  calcul aux poles  .......
@@ -53,10 +53,10 @@
       convpn =   SSUM( iim, yflu(     1    ,l ),  1 )
       convps = - SSUM( iim, yflu( ip1jm-iim,l ),  1 )
-      DO 4 ij = 1,iip1
+      DO ij = 1,iip1
       convfl(     ij   ,l ) = convpn * aire(   ij     ) / apoln
       convfl( ij+ ip1jm,l ) = convps * aire( ij+ ip1jm) / apols
-   4  CONTINUE
+      END DO
 c
-   5  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert.F90	(revision 5086)
@@ -320,5 +320,5 @@
      do l = 1, llm + 1
         read(unit, fmt=*) ap(l), bp(l)
-     end do
+     END DO
      close(unit)
      call assert(ap(1) == 0., ap(llm + 1) == 0., bp(1) == 1., &
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/disvert_noterre.F	(revision 5086)
@@ -139,5 +139,5 @@
          do l=1,llm
             read(99,*) zsig(l)
-         end do
+         END DO
          CLOSE(99)
 
@@ -146,5 +146,5 @@
            sig(l) = 0.5 * ( exp(-zsig(l)/scaleheight) + 
      &                      exp(-zsig(l-1)/scaleheight) )
-         end do
+         END DO
          sig(llm+1) =0
 
@@ -264,7 +264,7 @@
 c          zsig(l)= zsig(l-1)-scaleheight*
 c    .    log((aps(l) + bps(l)*ps)/(aps(l-1) + bps(l-1)*ps))
-c       end do
+c       END DO
 c       write(53,'(I3,50F10.5)') iz, zsig
-c      end do
+c      END DO
 c      close(53)
 c     --------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg.F	(revision 5086)
@@ -41,5 +41,5 @@
 c
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = iip2, ip1jm - 1
@@ -70,5 +70,5 @@
          div( ij + ip1jm, l ) =   sumyps
         ENDDO
-  10  CONTINUE
+      END DO
 c
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/diverg_gam.F	(revision 5086)
@@ -44,5 +44,5 @@
 c
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = iip2, ip1jm - 1
@@ -74,5 +74,5 @@
          div( ij + ip1jm, l ) =   sumyps 
         ENDDO
-  10  CONTINUE
+      END DO
 c
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergf.F	(revision 5086)
@@ -41,5 +41,5 @@
 c
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = iip2, ip1jm - 1
@@ -70,5 +70,5 @@
          div( ij + ip1jm, l ) =   sumyps
         ENDDO
-  10  CONTINUE
+      END DO
 c
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divergst.F	(revision 5086)
@@ -30,9 +30,9 @@
 c
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 1 ij = iip2, ip1jm - 1
+      DO ij = iip2, ip1jm - 1
       div( ij + 1, l ) = x(ij+1,l) - x(ij,l)+ y(ij-iim,l)-y(ij+1,l)
-   1  CONTINUE
+      END DO
 c
 c     ....  correction pour  div( 1,j,l)  ......
@@ -40,23 +40,23 @@
 c
 CDIR$ IVDEP
-      DO 3 ij = iip2,ip1jm,iip1
+      DO ij = iip2,ip1jm,iip1
       div( ij,l ) = div( ij + iim,l )
-   3  CONTINUE
+      END DO
 c
 c     ....  calcul  aux poles  .....
 c
 c
-      DO 5 i  = 1,iim
+      DO i  = 1,iim
       aiy1(i)= y(i,l)
       aiy2(i)= y(i+ip1jmi1,l)
-   5  CONTINUE
+      END DO
       sumypn = SSUM ( iim,aiy1,1 )
       sumyps = SSUM ( iim,aiy2,1 )
-      DO 7 i = 1,iip1
+      DO i = 1,iip1
       div(     i    , l ) = - sumypn/iim
       div( i + ip1jm, l ) =   sumyps/iim
-   7  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/divgrad.F	(revision 5086)
@@ -36,5 +36,5 @@
       CALL SCOPY ( ip1jmp1*klevel,h,1,divgra,1 )
 c
-      DO 10 iter = 1,lh
+      DO iter = 1,lh
 
       CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1  )
@@ -45,11 +45,11 @@
       CALL filtreg ( divgra,jjp1,klevel,2,1,.true.,1)
 
-      DO 5 l = 1,klevel
-      DO 4  ij = 1, ip1jmp1
+      DO l = 1,klevel
+      DO ij = 1, ip1jmp1
       divgra( ij,l ) = - cdivh * divgra( ij,l )
-   4  CONTINUE
-   5  CONTINUE
+      END DO
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/fxhyp_m.F90	(revision 5086)
@@ -191,5 +191,5 @@
              do while (rlonm025(is2) < - pi .and. is2 < iim)
                 is2 = is2 + 1
-             end do
+             END DO
 
              if (rlonm025(is2) < - pi) then
@@ -202,5 +202,5 @@
              do while (rlonm025(is2) > pi .and. is2 > 1)
                 is2 = is2 - 1
-             end do
+             END DO
 
              if (rlonm025(is2) > pi) then
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/geopot.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/geopot.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/geopot.F	(revision 5086)
@@ -47,7 +47,7 @@
 c     calcul de phi au niveau 1 pres du sol  .....
 
-      DO   1  ij  = 1, ngrid
+      DO   ij  = 1, ngrid
       phi( ij,1 ) = phis( ij ) + teta(ij,1) * ( pks(ij) - pk(ij,1) )
-   1  CONTINUE
+      END DO
 
 c     calcul de phi aux niveaux superieurs  .......
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/grad.F	(revision 5086)
@@ -23,22 +23,22 @@
 c
 c
-      DO 6 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 2  ij = 1, ip1jmp1 - 1
+      DO ij = 1, ip1jmp1 - 1
       pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
-   2  CONTINUE
+      END DO
 c
 c    .... correction pour  pgx(ip1,j,l)  ....
 c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
 CDIR$ IVDEP
-      DO 3  ij = iip1, ip1jmp1, iip1
+      DO ij = iip1, ip1jmp1, iip1
       pgx( ij,l ) = pgx( ij -iim,l )
-   3  CONTINUE
+      END DO
 c
-      DO 4 ij = 1,ip1jm
+      DO ij = 1,ip1jm
       pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
-   4  CONTINUE
+      END DO
 c
-   6  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/gradiv.F	(revision 5086)
@@ -37,5 +37,5 @@
       CALL SCOPY( ip1jm*klevel,  ycov,1,gdy,1 )
 c
-      DO 10 iter = 1,ld
+      DO iter = 1,ld
 c
       CALL  diverg( klevel,  gdx , gdy, div          )
@@ -43,14 +43,14 @@
       CALL    grad( klevel,  div, gdx, gdy           )
 c
-      DO 5  l = 1, klevel
-      DO 3 ij = 1, ip1jmp1
+      DO l = 1, klevel
+      DO ij = 1, ip1jmp1
       gdx( ij,l ) = - gdx( ij,l ) * cdivu
-   3  CONTINUE
-      DO 4 ij = 1, ip1jm
+      END DO
+      DO ij = 1, ip1jm
       gdy( ij,l ) = - gdy( ij,l ) * cdivu
-   4  CONTINUE
-   5  CONTINUE
+      END DO
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigeom.F	(revision 5086)
@@ -224,5 +224,5 @@
 c
         xo1 = 0.
-        DO 10 iter = 1, itmax
+        DO iter = 1, itmax
         x1  = xo1
         f   = x1+ alphax *SIN(x1-pxo)
@@ -232,5 +232,5 @@
         IF( xdm<=eps )GO TO 11
         xo1 = x1
- 10     CONTINUE
+      END DO
  11     CONTINUE
 c
@@ -241,5 +241,5 @@
 C
         yo1  = 0.
-        DO 15 iter = 1,itmay
+        DO iter = 1,itmay
         y1   = yo1
         f    = y1 + alphay* SIN(y1-pyo)
@@ -249,5 +249,5 @@
         IF(ydm<=eps) GO TO 17
         yo1  = y1
- 15     CONTINUE
+      END DO
 c
  17     CONTINUE
@@ -346,5 +346,5 @@
 c
 c
-      DO 35 j = 1, jjp1
+      DO j = 1, jjp1
 c
       IF ( j. eq. 1 )  THEN
@@ -356,5 +356,5 @@
       radclatm       = 0.5* rad * coslatm
 c
-      DO 30 i = 1, iim
+      DO i = 1, iim
       xprp           = xprimp025( i )
       xprm           = xprimm025( i )
@@ -365,5 +365,5 @@
       cvij2  ( i,1 ) = 0.5* rad * yprm
       cvij3  ( i,1 ) = cvij2(i,1)
-  30  CONTINUE
+      END DO
 c
       DO  i = 1, iim
@@ -387,5 +387,5 @@
       radclatp            = 0.5* rad * coslatp
 c
-      DO 31 i = 1,iim
+      DO i = 1,iim
         xprp              = xprimp025( i )
         xprm              = xprimm025( i )
@@ -396,5 +396,5 @@
         cvij1(i,jjp1)     = 0.5 * rad* yprp
         cvij4(i,jjp1)     = cvij1(i,jjp1)
- 31   CONTINUE
+      END DO
 c
        DO   i    = 1, iim
@@ -428,5 +428,5 @@
          ai14            = un4rad2 * coslatp * yprp
          ai23            = un4rad2 * coslatm * yprm
-         DO 32 i = 1,iim
+         DO i = 1,iim
          xprp            = xprimp025( i )
          xprm            = xprimm025( i )
@@ -444,5 +444,5 @@
          cvij3   ( i,j ) = cvij2(i,j)
          cvij4   ( i,j ) = cvij1(i,j)
-  32     CONTINUE
+      END DO
 c
       END IF
@@ -463,10 +463,10 @@
          aireij4 (iip1,j) = aireij4 (1,j )
         
-  35  CONTINUE
+      END DO
 c
 c    ..............................................................
 c
-      DO 37 j = 1, jjp1
-      DO 36 i = 1, iim
+      DO j = 1, jjp1
+      DO i = 1, iim
       aire    ( i,j )  = aireij1(i,j) + aireij2(i,j) + aireij3(i,j) +
      *                          aireij4(i,j)
@@ -479,5 +479,5 @@
       alpha2p3( i,j )  = alpha2 (i,j) + alpha3 (i,j)
       alpha3p4( i,j )  = alpha3 (i,j) + alpha4 (i,j)
-  36  CONTINUE
+      END DO
 c
 c
@@ -491,9 +491,9 @@
       alpha2p3(iip1,j) = alpha2p3(1,j)
       alpha3p4(iip1,j) = alpha3p4(1,j)
-  37  CONTINUE
-c
-
-      DO 42 j = 1,jjp1
-      DO 41 i = 1,iim
+      END DO
+c
+
+      DO j = 1,jjp1
+      DO i = 1,iim
       aireu       (i,j)= aireij1(i,j) + aireij2(i,j) + aireij4(i+1,j) +
      *                                aireij3(i+1,j)
@@ -502,5 +502,5 @@
       unsair_gam2( i,j)= unsaire(i,j)** ( - gamdi_h    )
       airesurg   ( i,j)= aire(i,j)/ g
-  41  CONTINUE
+      END DO
       aireu     (iip1,j)  = aireu  (1,j)
       unsaire   (iip1,j)  = unsaire(1,j)
@@ -508,8 +508,8 @@
       unsair_gam2(iip1,j) = unsair_gam2(1,j)
       airesurg   (iip1,j) = airesurg(1,j)
-  42  CONTINUE
-c
-c
-      DO 48 j = 1,jjm
+      END DO
+c
+c
+      DO j = 1,jjm
 c
         DO i=1,iim
@@ -529,5 +529,5 @@
         unsairz_gam(iip1,j) = unsairz_gam(1,j)
 c
-  48  CONTINUE
+      END DO
 c
 c
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inter_barxy_m.F90	(revision 5086)
@@ -270,5 +270,5 @@
           x0   = xxid(idat)
           idat = idat + 1
-       end do
+       END DO
        IF (xxim(imod)<xxid(idat)) THEN
           dx   = xxim(imod) - x0
@@ -287,5 +287,5 @@
           idat = idat + 1
        END IF
-    end do
+    END DO
 
   END function inter_barx
@@ -339,5 +339,5 @@
           y0         = yjdat(jdat)
           jdat       = jdat + 1
-       end do
+       END DO
        IF (yjmod(jmod) < yjdat(jdat)) THEN
           dy         = yjmod(jmod) - y0
@@ -357,5 +357,5 @@
           jdat       = jdat + 1
        END IF
-    end do
+    END DO
     ! Le test de fin suppose que l'interface 0 est commune aux deux
     ! grilles "yjdat" et "yjmod".
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/interpre.F	(revision 5086)
@@ -41,6 +41,6 @@
 
 c       CALCUL DE LA PRESSION DE SURFACE
-c       Les coefficients ap et bp sont passés en common 
-c       Calcul de la pression au sol en mb optimisée pour 
+c       Les coefficients ap et bp sont pass�s en common 
+c       Calcul de la pression au sol en mb optimis�e pour 
 c       la vectorialisation
                    
@@ -62,11 +62,11 @@
              do i=1,iim
                  psppm(i,j)=smass(i,j)/aire(i,j)*g*0.01
-             end do
-         end do                        
+             END DO
+         END DO
        
 c RECONSTRUCTION DES CHAMPS CONTRAVARIANTS
 c Le programme ppm3d travaille avec les composantes
-c de vitesse et pas les flux, on doit donc passer de l'un à l'autre
-c Dans le même temps, on fait le changement d'orientation du vent en v
+c de vitesse et pas les flux, on doit donc passer de l'un � l'autre
+c Dans le m�me temps, on fait le changement d'orientation du vent en v
       do l=1,llm
           do j=1,jjm
@@ -99,7 +99,7 @@
       
 c INVERSION DES NIVEAUX
-c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
+c le programme ppm3d travaille avec une 3�me coordonn�e invers�e par rapport
 c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
-c On passe donc des niveaux du LMDZ à ceux de Lin
+c On passe donc des niveaux du LMDZ � ceux de Lin
      
       do l=1,llm+1
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/invert_zoom_x_m.F90	(revision 5086)
@@ -39,5 +39,5 @@
        do while (xfi < xf(it) .and. it >= 1)
           it = it - 1
-       end do
+       END DO
 
        ! Calcul de Xf(xvrai(i)) 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/massdair.F	(revision 5086)
@@ -82,5 +82,5 @@
 c=======================================================================
 
-      DO   100    l = 1 , llm
+      DO   l = 1 , llm
 c
         DO    ij     = 1, ip1jmp1
@@ -103,5 +103,5 @@
 c       ENDDO
        
-100   CONTINUE
+      END DO
 c
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad.F	(revision 5086)
@@ -22,9 +22,9 @@
 c
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 1  ij = 2, ip1jm
+      DO ij = 2, ip1jm
       y( ij,l ) = (  rot( ij,l ) - rot( ij-1,l )  ) * cvsurcuv( ij )
-   1  CONTINUE
+      END DO
 c
 c    ..... correction pour  y ( 1,j,l )  ......
@@ -32,17 +32,17 @@
 c    ....    y(1,j,l)= y(iip1,j,l) ....
 CDIR$ IVDEP
-      DO 2  ij = 1, ip1jm, iip1
+      DO ij = 1, ip1jm, iip1
       y( ij,l ) = y( ij +iim,l )
-   2  CONTINUE
+      END DO
 c
-      DO 4  ij = iip2,ip1jm
+      DO ij = iip2,ip1jm
       x( ij,l ) = (  rot( ij,l ) - rot( ij -iip1,l )  ) * cusurcvu( ij )
-   4  CONTINUE
-      DO 6 ij = 1,iip1
+      END DO
+      DO ij = 1,iip1
       x(    ij    ,l ) = 0.
       x( ij +ip1jm,l ) = 0.
-   6  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad_gam.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad_gam.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrad_gam.F	(revision 5086)
@@ -21,9 +21,9 @@
       INTEGER   l,ij
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 1  ij = 2, ip1jm
+      DO ij = 2, ip1jm
       y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij )
-   1  CONTINUE
+      END DO
 c
 c    ..... correction pour  y ( 1,j,l )  ......
@@ -31,17 +31,17 @@
 c    ....    y(1,j,l)= y(iip1,j,l) ....
 CDIR$ IVDEP
-      DO 2  ij = 1, ip1jm, iip1
+      DO ij = 1, ip1jm, iip1
       y( ij,l ) = y( ij +iim,l )
-   2  CONTINUE
+      END DO
 c
-      DO 4  ij = iip2,ip1jm
+      DO ij = iip2,ip1jm
       x( ij,l ) = (rot( ij,l ) - rot( ij -iip1,l )) * cuscvugam( ij )
-   4  CONTINUE
-      DO 6 ij = 1,iip1
+      END DO
+      DO ij = 1,iip1
       x(    ij    ,l ) = 0.
       x( ij +ip1jm,l ) = 0.
-   6  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgradst.F	(revision 5086)
@@ -21,9 +21,9 @@
       INTEGER l,ij
 c
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 1  ij = 2, ip1jm
+      DO ij = 2, ip1jm
       y(ij,l)=( rot(ij,l) - rot(ij-1,l))
-   1  CONTINUE
+      END DO
 c
 c    ..... correction pour  y ( 1,j,l )  ......
@@ -31,17 +31,17 @@
 c    ....    y(1,j,l)= y(iip1,j,l) ....
 
-      DO 2  ij = 1, ip1jm, iip1
+      DO ij = 1, ip1jm, iip1
       y( ij,l ) = y( ij +iim,l )
-   2  CONTINUE
+      END DO
 c
-      DO 4  ij = iip2,ip1jm
+      DO ij = iip2,ip1jm
       x(ij,l)= rot(ij,l)-rot(ij-iip1,l)
-   4  CONTINUE
-      DO 6 ij = 1,iip1
+      END DO
+      DO ij = 1,iip1
       x(    ij    ,l ) = 0.
       x( ij +ip1jm,l ) = 0.
-   6  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrarot.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrarot.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/nxgrarot.F	(revision 5086)
@@ -36,19 +36,19 @@
       CALL SCOPY (  ip1jm*klevel, ycov, 1, gry, 1 )
 c
-      DO 10 iter = 1,lr
+      DO iter = 1,lr
       CALL  rotat (klevel,grx, gry, rot )
       CALL filtreg( rot, jjm, klevel, 2,1, .false.,2)
       CALL nxgrad (klevel,rot, grx, gry )
 c
-      DO 5  l = 1, klevel
-      DO 2 ij = 1, ip1jm
+      DO l = 1, klevel
+      DO ij = 1, ip1jm
       gry( ij,l ) = - gry( ij,l ) * crot
-   2  CONTINUE
-      DO 3 ij = 1, ip1jmp1
+      END DO
+      DO ij = 1, ip1jmp1
       grx( ij,l ) = - grx( ij,l ) * crot
-   3  CONTINUE
-   5  CONTINUE
+      END DO
+      END DO
 c
-  10  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pbar.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pbar.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/pbar.F	(revision 5086)
@@ -87,7 +87,7 @@
 
 
-      DO 1 ij = 1, ip1jmp1 - 1
+      DO ij = 1, ip1jmp1 - 1
       pbarx( ij ) = pext(ij) * alpha1p2(ij) + pext(ij+1)*alpha3p4(ij+1)
-   1  CONTINUE
+      END DO
 
 c    .... correction pour pbarx( iip1,j) .....
@@ -95,19 +95,19 @@
 c    ...    pbarx(iip1,j)= pbarx(1,j) ...
 CDIR$ IVDEP
-      DO 2 ij = iip1, ip1jmp1, iip1
+      DO ij = iip1, ip1jmp1, iip1
       pbarx( ij ) = pbarx( ij - iim )
-   2  CONTINUE
+      END DO
 
 
-      DO 3 ij = 1,ip1jm
+      DO ij = 1,ip1jm
       pbary( ij ) = pext(   ij  )   * alpha2p3(   ij   )     +
      *              pext( ij+iip1 ) * alpha1p4( ij+iip1 )
-   3  CONTINUE
+      END DO
 
 
-      DO 5 ij = 1, ip1jm - 1
+      DO ij = 1, ip1jm - 1
       pbarxy( ij ) = pext(ij)*alpha2(ij) + pext(ij+1)*alpha3(ij+1) +
      *   pext(ij+iip1)*alpha1(ij+iip1) + pext(ij+iip2)*alpha4(ij+iip2)
-   5  CONTINUE
+      END DO
 
 
@@ -116,7 +116,7 @@
 CDIR$ IVDEP
 
-      DO 7 ij = iip1, ip1jm, iip1
+      DO ij = iip1, ip1jm, iip1
       pbarxy( ij ) = pbarxy( ij - iim )
-   7  CONTINUE
+      END DO
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/prather.F	(revision 5086)
@@ -149,5 +149,5 @@
        call advxp( limit,0.5*dt,pbaru,sm,s0,sx,sy,sz
      .             ,sxx,sxy,sxz,syy,syz,szz,1 )
-        end do
+        END DO
         do l=1,llm
         do i=1,iip1
@@ -206,5 +206,5 @@
        call advxp( limit,0.5*dt,pbaru,sm,s0,sx,sy,sz
      .             ,sxx,sxy,sxz,syy,syz,szz,1 )
-        end do
+        END DO
 c---------------------------------------------------------
 c---------------------------------------------------------
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/psextbar.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/psextbar.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/psextbar.F	(revision 5086)
@@ -89,8 +89,8 @@
 
 
-      DO     5     ij = 1, ip1jm - 1
+      DO     ij = 1, ip1jm - 1
       psexbarxy( ij ) = pext(ij)*alpha2(ij) + pext(ij+1)*alpha3(ij+1) +
      *   pext(ij+iip1)*alpha1(ij+iip1) + pext(ij+iip2)*alpha4(ij+iip2)
-   5  CONTINUE
+      END DO
 
 
@@ -99,7 +99,7 @@
 CDIR$ IVDEP
 
-      DO 7 ij = iip1, ip1jm, iip1
+      DO ij = iip1, ip1jm, iip1
       psexbarxy( ij ) = psexbarxy( ij - iim )
-   7  CONTINUE
+      END DO
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat.F	(revision 5086)
@@ -30,5 +30,5 @@
 c
 c
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = 1, ip1jm - 1
@@ -44,5 +44,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 
 ccc        CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 )
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat_nfil.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat_nfil.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotat_nfil.F	(revision 5086)
@@ -30,5 +30,5 @@
 c
 c
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = 1, ip1jm - 1
@@ -44,5 +44,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatf.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatf.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatf.F	(revision 5086)
@@ -31,5 +31,5 @@
 c
 c
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = 1, ip1jm - 1
@@ -45,5 +45,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 
         CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 )
Index: LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatst.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatst.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3d_common/rotatst.F	(revision 5086)
@@ -24,10 +24,10 @@
 c
 c
-      DO 5 l = 1,klevel
+      DO l = 1,klevel
 c
-      DO 1 ij = 1, ip1jm - 1
+      DO ij = 1, ip1jm - 1
       rot( ij,l )  =  (  y( ij+1 , l )  -  y( ij,l )   +
      *                 x(ij +iip1, l )  -  x( ij,l )  )
-   1  CONTINUE
+      END DO
 c
 c    .... correction pour rot( iip1,j,l)  ....
@@ -35,9 +35,9 @@
 c    ....   rot(iip1,j,l)= rot(1,j,l) ...
 CDIR$ IVDEP
-      DO 2 ij = iip1, ip1jm, iip1
+      DO ij = iip1, ip1jm, iip1
       rot( ij,l ) = rot( ij -iim,l )
-   2  CONTINUE
+      END DO
 c
-   5  CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/advect_new_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/advect_new_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/advect_new_loc.F	(revision 5086)
@@ -63,7 +63,7 @@
          deuxjour = 2. * daysec
 
-         DO   1  ij   = 1, ip1jmp1
+         DO   ij   = 1, ip1jmp1
          unsaire2(ij) = unsaire(ij) * unsaire(ij)
-   1     CONTINUE
+      END DO
       END IF
 
@@ -175,5 +175,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO 20 l = 1, llmm1
+      DO l = 1, llmm1
 
 
@@ -183,7 +183,7 @@
       if (pole_sud)  ije=ij_end
       
-      DO 5   ij   = ijb, ije
+      DO ij   = ijb, ije
       wsur2( ij ) = - 0.5 * w( ij,l+1 )
-   5  CONTINUE
+      END DO
 
 
@@ -195,10 +195,10 @@
       if (pole_sud)  ije=ije-iip1
          
-      DO 6 ij = ijb ,ije-1
+      DO ij = ijb ,ije-1
       ww        = wsur2 (  ij  )     + wsur2( ij+1 ) 
       uu        = 0.5 * ( ucov(ij,l) + ucov(ij,l+1) )
       du1(ij,l)  =  ww * ( uu - uav(ij, l ) )/massebx(ij, l )
       du2(ij,l+1)=  ww * ( uu - uav(ij,l+1) )/massebx(ij,l+1)
-   6  CONTINUE
+      END DO
 
 c     .................    calcul pour   dv      .....................
@@ -207,10 +207,10 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 8 ij = ijb, ije
+      DO ij = ijb, ije
       ww        = wsur2( ij+iip1 )   + wsur2( ij )
       vv        = 0.5 * ( vcov(ij,l) + vcov(ij,l+1) )
       dv1(ij,l)  =  ww * (vv - vav(ij, l ) )/masseby(ij, l )
       dv2(ij,l+1)=  ww * (vv - vav(ij,l+1) )/masseby(ij,l+1)
-   8  CONTINUE
+      END DO
 
 c
@@ -226,9 +226,9 @@
         ije=ij_end
         
-        DO 15 ij = ijb, ije
+        DO ij = ijb, ije
          ww            = wsur2(ij) * (teta(ij,l) + teta(ij,l+1) )
          dteta1(ij, l ) =   ww
          dteta2(ij,l+1) =   ww
-  15    CONTINUE
+      END DO
 
 c ym ---> conser a voir plus tard
@@ -243,5 +243,5 @@
 c      END IF
 
-  20  CONTINUE
+      END DO
 c$OMP END DO
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bernoui_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bernoui_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/bernoui_loc.F	(revision 5086)
@@ -54,7 +54,7 @@
       DO l=1,llm
     
-        DO 4 ij = ijb,ije
+        DO ij = ijb,ije
           pbern( ij,l ) =  pphi( ij,l ) + pecin( ij,l )
-   4    CONTINUE
+      END DO
        
        ENDDO
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/convflu_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/convflu_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/convflu_loc.F	(revision 5086)
@@ -34,5 +34,5 @@
      
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)          
-      DO 5 l = 1,nbniv
+      DO l = 1,nbniv
 c
         ijb=ij_begin
@@ -42,8 +42,8 @@
         IF (pole_sud)  ije=ij_end-iip1
         
-        DO 2  ij = ijb , ije - 1
+        DO ij = ijb , ije - 1
           convfl(ij+1,l) = xflu(ij,l) - xflu(ij+ 1,l)   +
      *                     yflu(ij +1,l ) - yflu( ij -iim,l )
-   2    CONTINUE
+      END DO
 c
 c
@@ -53,7 +53,7 @@
 c
 CDIR$ IVDEP
-        DO 3 ij = ijb,ije,iip1
+        DO ij = ijb,ije,iip1
           convfl( ij,l ) = convfl( ij + iim,l )
-   3    CONTINUE
+      END DO
 c
 c     ......  calcul aux poles  .......
@@ -79,5 +79,5 @@
         ENDIF
       
-   5  CONTINUE
+      END DO
 c$OMP END DO NOWAIT   
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/covcont_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/covcont_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/covcont_loc.F	(revision 5086)
@@ -44,15 +44,15 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)    
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 
-      DO 2  ij = ijb_u,ije_u
+      DO ij = ijb_u,ije_u
       ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
-   2  CONTINUE
+      END DO
 
-      DO 4 ij = ijb_v,ije_v
+      DO ij = ijb_v,ije_v
       vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
-   4  CONTINUE
+      END DO
 
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_gam_loc.F	(revision 5086)
@@ -51,5 +51,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = ijb, ije - 1
@@ -91,5 +91,5 @@
           ENDDO
        endif
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 c
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/diverg_p.F	(revision 5086)
@@ -47,5 +47,5 @@
       
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = ijb, ije - 1
@@ -88,5 +88,5 @@
 
 
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 c
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/divergf_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/divergf_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/divergf_loc.F	(revision 5086)
@@ -49,5 +49,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
         DO  ij = ijb, ije - 1
@@ -95,5 +95,5 @@
         endif
         
-  10    CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dteta1_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dteta1_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dteta1_loc.F	(revision 5086)
@@ -39,5 +39,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)    
-      DO 5 l = 1,llm
+      DO l = 1,llm
       
       ijb=ij_begin
@@ -47,7 +47,7 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 1  ij = ijb, ije - 1
+      DO ij = ijb, ije - 1
         hbxu(ij,l) = pbaru(ij,l) * 0.5 * ( teta(ij,l)+teta(ij+1,l) )
-   1  CONTINUE
+      END DO
 
 c    .... correction pour  hbxu(iip1,j,l)  .....
@@ -55,14 +55,14 @@
 
 CDIR$ IVDEP
-      DO 2 ij = ijb+iip1-1, ije, iip1
+      DO ij = ijb+iip1-1, ije, iip1
         hbxu( ij, l ) = hbxu( ij - iim, l )
-   2  CONTINUE
+      END DO
 
       ijb=ij_begin-iip1
       if (pole_nord) ijb=ij_begin
       
-      DO 3 ij = ijb,ije
+      DO ij = ijb,ije
         hbyv(ij,l)= pbarv(ij, l)* 0.5 * ( teta(ij, l)+teta(ij+iip1,l) )
-   3  CONTINUE
+      END DO
 
        if (.not. pole_sud) then
@@ -71,5 +71,5 @@
 	endif
 	
-   5  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
        
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv1_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv1_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv1_loc.F	(revision 5086)
@@ -28,5 +28,5 @@
       
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
-      DO 10 l = 1,llm
+      DO l = 1,llm
 c
       ijb=ij_begin
@@ -36,9 +36,9 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 2  ij = ijb, ije-1 
+      DO ij = ijb, ije-1
       du( ij,l ) = 0.125 *(  vorpot(ij-iip1, l) + vorpot( ij, l)  ) *
      *                    (   pbarv(ij-iip1, l) + pbarv(ij-iim,  l) +
      *                        pbarv(   ij  , l) + pbarv(ij+ 1 ,  l)   )
-   2  CONTINUE
+      END DO
    
  
@@ -46,9 +46,9 @@
       if (pole_nord) ijb=ij_begin
       
-      DO 3 ij = ijb, ije-1 
+      DO ij = ijb, ije-1
       dv( ij+1,l ) = - 0.125 *(  vorpot(ij, l)  + vorpot(ij+1, l)  ) *
      *                        (   pbaru(ij, l)  +  pbaru(ij+1   , l) +
      *                       pbaru(ij+iip1, l)  +  pbaru(ij+iip2, l)  )
-   3  CONTINUE
+      END DO
 c
 c    .... correction  pour  dv( 1,j,l )  .....
@@ -56,9 +56,9 @@
 c
 CDIR$ IVDEP
-      DO 4 ij = ijb, ije, iip1
+      DO ij = ijb, ije, iip1
       dv( ij,l ) = dv( ij + iim, l )
-   4  CONTINUE
+      END DO
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv2_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv2_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/dudv2_loc.F	(revision 5086)
@@ -33,5 +33,5 @@
 c
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO 5 l = 1,llm
+      DO l = 1,llm
 c
       ijb=ij_begin
@@ -40,8 +40,8 @@
       if (pole_sud)  ije=ije-iip1
 
-      DO 2  ij  = ijb, ije - 1
+      DO ij  = ijb, ije - 1
        du(ij,l) = du(ij,l) + 0.5* ( teta( ij,l ) + teta( ij+1,l ) ) *
      * ( pkf( ij,l ) - pkf(ij+1,l) )  + bern(ij,l) - bern(ij+1,l)
-   2  CONTINUE
+      END DO
 c
 c
@@ -50,18 +50,18 @@
 c
 CDIR$ IVDEP
-      DO 3 ij = ijb+iip1-1, ije, iip1
+      DO ij = ijb+iip1-1, ije, iip1
       du( ij,l ) = du( ij - iim,l )
-   3  CONTINUE
+      END DO
 c
 c
       if (pole_nord) ijb=ijb-iip1
 
-      DO 4 ij  = ijb,ije
+      DO ij  = ijb,ije
       dv( ij,l) = dv(ij,l) + 0.5 * ( teta(ij,l) + teta( ij+iip1,l ) ) *
      *                             ( pkf(ij+iip1,l) - pkf(  ij,l  ) )
      *                           +   bern( ij+iip1,l ) - bern( ij  ,l )
-   4  CONTINUE
+      END DO
 c
-   5  CONTINUE
+      END DO
 c$OMP END DO NOWAIT 
 c
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_loc.F	(revision 5086)
@@ -23,18 +23,18 @@
 c
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO 6 l = 1,klevel
+      DO l = 1,klevel
 c
       ijb=ij_begin
       ije=ij_end
-      DO 2  ij = ijb, ije - 1
+      DO ij = ijb, ije - 1
         pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
-   2  CONTINUE
+      END DO
 c
 c    .... correction pour  pgx(ip1,j,l)  ....
 c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
 CDIR$ IVDEP
-      DO 3  ij = ijb+iip1-1, ije, iip1
+      DO ij = ijb+iip1-1, ije, iip1
         pgx( ij,l ) = pgx( ij -iim,l )
-   3  CONTINUE
+      END DO
 c
       ijb=ij_begin-iip1
@@ -43,9 +43,9 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 4 ij = ijb,ije
+      DO ij = ijb,ije
         pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
-   4  CONTINUE
+      END DO
 c
-   6  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_p.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_p.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/grad_p.F	(revision 5086)
@@ -23,18 +23,18 @@
 c
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO 6 l = 1,klevel
+      DO l = 1,klevel
 c
       ijb=ij_begin
       ije=ij_end
-      DO 2  ij = ijb, ije - 1
+      DO ij = ijb, ije - 1
         pgx( ij,l ) = pg( ij +1,l ) - pg( ij,l )
-   2  CONTINUE
+      END DO
 c
 c    .... correction pour  pgx(ip1,j,l)  ....
 c    ...    pgx(iip1,j,l)= pgx(1,j,l)  ....
 CDIR$ IVDEP
-      DO 3  ij = ijb+iip1-1, ije, iip1
+      DO ij = ijb+iip1-1, ije, iip1
         pgx( ij,l ) = pgx( ij -iim,l )
-   3  CONTINUE
+      END DO
 c
       ijb=ij_begin-iip1
@@ -43,9 +43,9 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 4 ij = ijb,ije
+      DO ij = ijb,ije
         pgy( ij,l ) = pg( ij,l ) - pg( ij +iip1,l )
-   4  CONTINUE
+      END DO
 c
-   6  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/integrd_loc.F	(revision 5086)
@@ -124,9 +124,9 @@
 
 c$OMP DO SCHEDULE(STATIC)
-      DO 2 ij = ijb,ije
+      DO ij = ijb,ije
        pscr (ij)    = ps0(ij)
        ps (ij)      = psm1(ij) + dt * dp(ij)     
 
-   2  CONTINUE
+      END DO
 
 c$OMP END DO  
@@ -250,5 +250,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 
-      DO 10 l = 1,llm
+      DO l = 1,llm
       
       ijb=ij_begin
@@ -257,8 +257,8 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 4 ij = ijb,ije
+      DO ij = ijb,ije
       uscr( ij )   =  ucov( ij,l )
       ucov( ij,l ) = ucovm1( ij,l ) + dt * du( ij,l )
-   4  CONTINUE
+      END DO
 
       ijb=ij_begin
@@ -266,17 +266,17 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 5 ij = ijb,ije
+      DO ij = ijb,ije
       vscr( ij )   =  vcov( ij,l )
       vcov( ij,l ) = vcovm1( ij,l ) + dt * dv( ij,l )
-   5  CONTINUE
+      END DO
       
       ijb=ij_begin
       ije=ij_end
       
-      DO 6 ij = ijb,ije
+      DO ij = ijb,ije
       hscr( ij )    =  teta(ij,l)
       teta ( ij,l ) = tetam1(ij,l) *  massem1(ij,l) / masse(ij,l) 
      $                + dt * dteta(ij,l) / masse(ij,l)
-   6  CONTINUE
+      END DO
 
 c   ....  Calcul de la valeur moyenne, unique  aux poles pour  teta    ......
@@ -324,5 +324,5 @@
       END IF
 
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/massdair_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/massdair_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/massdair_loc.F	(revision 5086)
@@ -92,5 +92,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)    
-      DO   100    l = 1 , llm
+      DO   l = 1 , llm
 c
         DO    ij     = ijb, ije
@@ -113,5 +113,5 @@
 c       ENDDO
        
-100   CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 c
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/mod_hallo.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/mod_hallo.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/mod_hallo.F90	(revision 5086)
@@ -169,5 +169,5 @@
     do while (Buffer_Pos(Index_Pos)==-1 .and. Index_Pos>1)
       Index_Pos=Index_Pos-1
-    end do
+    END DO
 
   end subroutine deallocate_buffer  
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_gam_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_gam_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_gam_loc.F	(revision 5086)
@@ -25,5 +25,5 @@
 c
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
       ijb=ij_begin
@@ -31,7 +31,7 @@
       if(pole_sud) ije=ij_end-iip1
       
-      DO 1  ij = ijb+1, ije
+      DO ij = ijb+1, ije
       y( ij,l ) = (rot( ij,l ) - rot( ij-1,l )) * cvscuvgam( ij )
-   1  CONTINUE
+      END DO
 c
 c    ..... correction pour  y ( 1,j,l )  ......
@@ -39,7 +39,7 @@
 c    ....    y(1,j,l)= y(iip1,j,l) ....
 CDIR$ IVDEP
-      DO 2  ij = ijb, ije, iip1
+      DO ij = ijb, ije, iip1
       y( ij,l ) = y( ij +iim,l )
-   2  CONTINUE
+      END DO
 c
       ijb=ij_begin
@@ -48,7 +48,7 @@
       if(pole_sud) ije=ij_end-iip1
       
-      DO 4  ij = ijb,ije
+      DO ij = ijb,ije
       x( ij,l ) = (rot( ij,l ) - rot( ij -iip1,l )) * cuscvugam( ij )
-   4  CONTINUE
+      END DO
     
       if (pole_nord) then
@@ -64,5 +64,5 @@
       endif
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/nxgrad_loc.F	(revision 5086)
@@ -23,5 +23,5 @@
 c
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      DO 10 l = 1,klevel
+      DO l = 1,klevel
 c
       ijb=ij_begin
@@ -29,7 +29,7 @@
       if (pole_sud)  ije=ij_end-iip1
        
-      DO 1  ij = ijb+1, ije
+      DO ij = ijb+1, ije
       y( ij,l ) = (  rot( ij,l ) - rot( ij-1,l )  ) * cvsurcuv( ij )
-   1  CONTINUE
+      END DO
 c
 c    ..... correction pour  y ( 1,j,l )  ......
@@ -37,7 +37,7 @@
 c    ....    y(1,j,l)= y(iip1,j,l) ....
 CDIR$ IVDEP
-      DO 2  ij = ijb, ije, iip1
+      DO ij = ijb, ije, iip1
       y( ij,l ) = y( ij +iim,l )
-   2  CONTINUE
+      END DO
 c
       ijb=ij_begin
@@ -47,7 +47,7 @@
       if (pole_sud)  ije=ij_end-iip1
       
-      DO 4  ij = ijb,ije
+      DO ij = ijb,ije
       x( ij,l ) = (  rot( ij,l ) - rot( ij -iip1,l )  ) * cusurcvu( ij )
-   4  CONTINUE
+      END DO
    
       if (pole_nord) then 
@@ -63,5 +63,5 @@
       endif
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_nfil_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_nfil_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_nfil_loc.F	(revision 5086)
@@ -33,5 +33,5 @@
       if(pole_sud) ije=ij_end-iip1
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = ijb, ije - 1
@@ -47,5 +47,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_p.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_p.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotat_p.F	(revision 5086)
@@ -34,5 +34,5 @@
       
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = ijb, ije - 1
@@ -48,5 +48,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
 ccc        CALL filtreg( rot, jjm, klevel, 2, 2, .FALSE., 1 )
Index: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotatf_loc.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotatf_loc.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dyn3dmem/rotatf_loc.F	(revision 5086)
@@ -35,5 +35,5 @@
 
 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)      
-      DO  10 l = 1,klevel
+      DO  l = 1,klevel
 c
         DO   ij = ijb, ije - 1
@@ -49,5 +49,5 @@
         ENDDO
 c
-  10  CONTINUE
+      END DO
 c$OMP END DO NOWAIT
         jjb=jj_begin
Index: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/calfis.F	(revision 5086)
@@ -347,7 +347,7 @@
 c   ------------
 
-      DO 50 l=1,llm
-
-         DO 25 j=2,jjm
+      DO l=1,llm
+
+         DO j=2,jjm
             ig0 = 1+(j-2)*iim
             zufi(ig0+1,l)= 0.5 * 
@@ -355,13 +355,13 @@
             pcvgu(ig0+1,l)= 0.5 * 
      $      ( pducov(iim,j,l)/cu(iim,j) + pducov(1,j,l)/cu(1,j) )
-            DO 10 i=2,iim
+            DO i=2,iim
                zufi(ig0+i,l)= 0.5 *
      $         ( pucov(i-1,j,l)/cu(i-1,j) + pucov(i,j,l)/cu(i,j) )
                pcvgu(ig0+i,l)= 0.5 *
      $         ( pducov(i-1,j,l)/cu(i-1,j) + pducov(i,j,l)/cu(i,j) )
-10         CONTINUE
-25      CONTINUE
-
-50    CONTINUE
+      END DO
+      END DO
+
+      END DO
 
 
Index: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/test_disvert_m.F90	(revision 5086)
@@ -55,6 +55,6 @@
                      p(i, l + 1) / 100., " hPa"
              end if
-          end do
-       end do
+          END DO
+       END DO
        call abort_physic("test_disvert", "bad order of pressure values", 1)
     end if
Index: LMDZ6/branches/Amaury_dev/libf/filtrez/eigen.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/filtrez/eigen.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/filtrez/eigen.F	(revision 5086)
@@ -10,10 +10,10 @@
       im=iim
 c
-      DO 48 i = 1,im
+      DO i = 1,im
 	 asm( i ) = d( im-i+1 )
- 48   CONTINUE
-      DO 49 i = 1,iim
+      END DO
+      DO i = 1,iim
 	 d( i ) = asm( i )
- 49   CONTINUE
+      END DO
 c
 c     PRINT 70,d
@@ -21,12 +21,12 @@
 		print *
 c
-      DO 51 i = 1,im
-	 DO 52 j = 1,im
+      DO i = 1,im
+	 DO j = 1,im
             asm( j ) = e( i , im-j+1 )
- 52      CONTINUE
-	 DO 50 j = 1,im
+      END DO
+	 DO j = 1,im
 	    e( i,j ) = asm( j )
- 50      CONTINUE
- 51   CONTINUE
+      END DO
+      END DO
 
       RETURN
Index: LMDZ6/branches/Amaury_dev/libf/misc/arth_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/arth_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/arth_m.F90	(revision 5086)
@@ -34,9 +34,9 @@
        do k=2,n
           arth_r(k)=arth_r(k-1)+increment
-       end do
+       END DO
     else
        do k=2,NPAR2_ARTH
           arth_r(k)=arth_r(k-1)+increment
-       end do
+       END DO
        temp=increment*NPAR2_ARTH
        k=NPAR2_ARTH
@@ -47,5 +47,5 @@
           temp=temp+temp
           k=k2
-       end do
+       END DO
     end if
 
@@ -68,9 +68,9 @@
        do k=2,n
           arth_i(k)=arth_i(k-1)+increment
-       end do
+       END DO
     else
        do k=2,NPAR2_ARTH
           arth_i(k)=arth_i(k-1)+increment
-       end do
+       END DO
        temp=increment*NPAR2_ARTH
        k=NPAR2_ARTH
@@ -81,5 +81,5 @@
           temp=temp+temp
           k=k2
-       end do
+       END DO
     end if
 
Index: LMDZ6/branches/Amaury_dev/libf/misc/chfev.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/chfev.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/chfev.F	(revision 5086)
@@ -124,5 +124,5 @@
 C  EVALUATION LOOP.
 C
-      DO 500  I = 1, NE
+      DO I = 1, NE
          X = XE(I) - X1
          FE(I) = F1 + X*(D1 + X*(C2 + X*C3))
@@ -131,5 +131,5 @@
          IF ( X>XMA )  NEXT(2) = NEXT(2) + 1
 C        (NOTE REDUNDANCY--IF EITHER CONDITION IS TRUE, OTHER IS FALSE.)
-  500 CONTINUE
+      END DO
 C
 C  NORMAL RETURN.
Index: LMDZ6/branches/Amaury_dev/libf/misc/cray.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/cray.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/cray.F	(revision 5086)
@@ -36,5 +36,5 @@
          ssum=ssum+sx(ix)
          ix=ix+incx
-      end do
+      END DO
 c
       return
Index: LMDZ6/branches/Amaury_dev/libf/misc/interpolation.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/interpolation.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/interpolation.F90	(revision 5086)
@@ -40,5 +40,5 @@
           ju=jm ! or the upper limit, as appropriate.
        end if
-    end do
+    END DO
     ! {ju == jl + 1}
 
@@ -102,5 +102,5 @@
                 inc=inc+inc ! so double the increment
              end if
-          end do ! and try again.
+          END DO ! and try again.
        else ! Hunt down:
           jhi=jlo
@@ -115,5 +115,5 @@
                 inc=inc+inc ! so double the increment
              end if
-          end do ! and try again.
+          END DO ! and try again.
        end if
     end if ! Done hunting, value bracketed.
@@ -132,5 +132,5 @@
           end if
        end if
-    end do
+    END DO
 
   END SUBROUTINE hunt
Index: LMDZ6/branches/Amaury_dev/libf/misc/ismax.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/ismax.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/ismax.F	(revision 5086)
@@ -12,5 +12,5 @@
       ismax=1
       sxmax=sx(1)
-      do 10 i=1,n-1
+      do i=1,n-1
        ix=ix+incx
        if(sx(ix)>sxmax) then
@@ -18,5 +18,5 @@
          ismax=i+1
        endif
-10    continue
+      END DO
 c
       return
Index: LMDZ6/branches/Amaury_dev/libf/misc/new_unit_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/new_unit_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/new_unit_m.F90	(revision 5086)
@@ -19,5 +19,5 @@
        if (exist .and. .not. opened) exit
        unit = unit + 1
-    end do
+    END DO
 
   end subroutine new_unit
Index: LMDZ6/branches/Amaury_dev/libf/misc/pchdf.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/pchdf.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/pchdf.F	(revision 5086)
@@ -76,16 +76,16 @@
 C  COMPUTE COEFFICIENTS OF INTERPOLATING POLYNOMIAL.
 C
-      DO 10  J = 2, K-1
-         DO 9  I = 1, K-J
+      DO J = 2, K-1
+         DO I = 1, K-J
             S(I) = (S(I+1)-S(I))/(X(I+J)-X(I))
-    9    CONTINUE
-   10 CONTINUE
+      END DO
+      END DO
 C
 C  EVALUATE DERIVATIVE AT X(K).
 C
       VALUE = S(1)
-      DO 20  I = 2, K-1
+      DO I = 2, K-1
          VALUE = S(I) + VALUE*(X(K)-X(I))
-   20 CONTINUE
+      END DO
 C
 C  NORMAL RETURN.
Index: LMDZ6/branches/Amaury_dev/libf/misc/pchfe.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/pchfe.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/pchfe.F	(revision 5086)
@@ -145,7 +145,7 @@
       IF ( N<2 )  GO TO 5001
       IF ( INCFD<1 )  GO TO 5002
-      DO 1  I = 2, N
+      DO I = 2, N
          IF ( X(I)<=X(I-1) )  GO TO 5003
-    1 CONTINUE
+      END DO
 C
 C  FUNCTION DEFINITION IS OK, GO ON.
@@ -168,7 +168,7 @@
 C     LOCATE ALL POINTS IN INTERVAL.
 C
-         DO 20  J = JFIRST, NE
+         DO J = JFIRST, NE
             IF (XE(J) >= X(IR))  GO TO 30
-   20    CONTINUE
+      END DO
          J = NE + 1
          GO TO 40
@@ -228,7 +228,7 @@
 C
 C              FIRST, LOCATE FIRST POINT TO LEFT OF X(IR-1).
-               DO 44  I = JFIRST, J-1
+               DO I = JFIRST, J-1
                   IF (XE(I) < X(IR-1))  GO TO 45
-   44          CONTINUE
+      END DO
 C              NOTE-- CANNOT DROP THROUGH HERE UNLESS THERE IS AN ERROR
 C                     IN CHFEV.
@@ -240,7 +240,7 @@
 C
 C              NOW FIND OUT HOW FAR TO BACK UP IN THE X-ARRAY.
-               DO 46  I = 1, IR-1
+               DO I = 1, IR-1
                   IF (XE(J) < X(I)) GO TO 47
-   46          CONTINUE
+      END DO
 C              NB-- CAN NEVER DROP THROUGH HERE, SINCE XE(J).LT.X(IR-1).
 C
Index: LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/pchsp.F	(revision 5086)
@@ -164,7 +164,7 @@
       IF ( N<2 )  GO TO 5001
       IF ( INCFD<1 )  GO TO 5002
-      DO 1  J = 2, N
+      DO J = 2, N
          IF ( X(J)<=X(J-1) )  GO TO 5003
-    1 CONTINUE
+      END DO
 C
       IBEG = IC(1)
@@ -181,8 +181,8 @@
 C  COMPUTE FIRST DIFFERENCES OF X SEQUENCE AND STORE IN WK(1,.). ALSO,
 C  COMPUTE FIRST DIVIDED DIFFERENCE OF DATA AND STORE IN WK(2,.).
-      DO 5  J=2,N
+      DO J=2,N
          WK(1,J) = X(J) - X(J-1)
          WK(2,J) = (F(1,J) - F(1,J-1))/WK(1,J)
-    5 CONTINUE
+      END DO
 C
 C  SET TO DEFAULT BOUNDARY CONDITIONS IF N IS TOO SMALL.
@@ -197,10 +197,10 @@
       ELSE IF (IBEG > 2)  THEN
 C        PICK UP FIRST IBEG POINTS, IN REVERSE ORDER.
-         DO 10  J = 1, IBEG
+         DO J = 1, IBEG
             INDEX = IBEG-J+1
 C           INDEX RUNS FROM IBEG DOWN TO 1.
             XTEMP(J) = X(INDEX)
             IF (J < IBEG)  STEMP(J) = WK(2,INDEX)
-   10    CONTINUE
+      END DO
 C                 --------------------------------
          D(1,1) = PCHDF (IBEG, XTEMP, STEMP, IERR)
@@ -214,10 +214,10 @@
       ELSE IF (IEND > 2)  THEN
 C        PICK UP LAST IEND POINTS.
-         DO 15  J = 1, IEND
+         DO J = 1, IEND
             INDEX = N-IEND+J
 C           INDEX RUNS FROM N+1-IEND UP TO N.
             XTEMP(J) = X(INDEX)
             IF (J < IEND)  STEMP(J) = WK(2,INDEX+1)
-   15    CONTINUE
+      END DO
 C                 --------------------------------
          D(1,N) = PCHDF (IEND, XTEMP, STEMP, IERR)
@@ -267,5 +267,5 @@
       NM1 = N-1
       IF (NM1 > 1)  THEN
-         DO 20 J=2,NM1
+         DO J=2,NM1
             IF (WK(2,J-1) == ZERO)  GO TO 5008
             G = -WK(1,J+1)/WK(2,J-1)
@@ -273,5 +273,5 @@
      *                  + THREE*(WK(1,J)*WK(2,J+1) + WK(1,J+1)*WK(2,J))
             WK(2,J) = G*WK(1,J-1) + TWO*(WK(1,J) + WK(1,J+1))
-   20    CONTINUE
+      END DO
       ENDIF
 C
@@ -324,8 +324,8 @@
 C
    30 CONTINUE
-      DO 40 J=NM1,1,-1
+      DO J=NM1,1,-1
          IF (WK(2,J) == ZERO)  GO TO 5008
          D(1,J) = (D(1,J) - WK(1,J)*D(1,J+1))/WK(2,J)
-   40 CONTINUE
+      END DO
 C --------------------(  END  CODING FROM CUBSPL )--------------------
 C
Index: LMDZ6/branches/Amaury_dev/libf/misc/ran1.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/ran1.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/ran1.F	(revision 5086)
@@ -21,9 +21,9 @@
         IX1=MOD(IA1*IX1+IC1,M1)
         IX3=MOD(IX1,M3)
-        DO 11 J=1,97
+        DO J=1,97
           IX1=MOD(IA1*IX1+IC1,M1)
           IX2=MOD(IA2*IX2+IC2,M2)
           R(J)=(REAL(IX1)+REAL(IX2)*RM2)*RM1
-11      CONTINUE
+      END DO
         IDUM=1
       ENDIF
Index: LMDZ6/branches/Amaury_dev/libf/misc/regr1_step_av_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/regr1_step_av_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/regr1_step_av_m.F90	(revision 5086)
@@ -73,5 +73,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it) = (vt(it) + (xt(it + 1) - left_edge) * vs(is)) &
@@ -79,5 +79,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr11_step_av
@@ -133,5 +133,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :) = (vt(it, :) + (xt(it + 1) - left_edge) * vs(is, :)) &
@@ -139,5 +139,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr12_step_av
@@ -194,5 +194,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :, :) = (vt(it, :, :) &
@@ -200,5 +200,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr13_step_av
@@ -256,5 +256,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :, :, :) = (vt(it, :, :, :) + (xt(it + 1) - left_edge) &
@@ -262,5 +262,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr14_step_av
Index: LMDZ6/branches/Amaury_dev/libf/misc/xermsg.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/xermsg.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/xermsg.F	(revision 5086)
@@ -328,7 +328,7 @@
       IF (LKNTRL > 0) THEN
          WRITE (TEMP, '(''ERROR NUMBER = '', I8)') NERR
-         DO 10 I=16,22
+         DO I=16,22
             IF (TEMP(I:I) /= ' ') GO TO 20
-   10    CONTINUE
+      END DO
 C
    20    CALL XERPRN (' *  ', -1, TEMP(1:15) // TEMP(I:23), 72)
Index: LMDZ6/branches/Amaury_dev/libf/misc/xerprn.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/xerprn.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/xerprn.F	(revision 5086)
@@ -92,7 +92,7 @@
 C
       N = I1MACH(4)
-      DO 10 I=1,NUNIT
+      DO I=1,NUNIT
          IF (IU(I) == 0) IU(I) = N
-   10 CONTINUE
+      END DO
 C
 C       LPREF IS THE LENGTH OF THE PREFIX.  THE PREFIX IS PLACED AT THE
@@ -117,8 +117,8 @@
       LENMSG = LEN(MESSG)
       N = LENMSG
-      DO 20 I=1,N
+      DO I=1,N
          IF (MESSG(LENMSG:LENMSG) /= ' ') GO TO 30
          LENMSG = LENMSG - 1
-   20 CONTINUE
+      END DO
    30 CONTINUE
 C
@@ -127,7 +127,7 @@
       IF (LENMSG == 0) THEN
          CBUFF(LPREF+1:LPREF+1) = ' '
-         DO 40 I=1,NUNIT
+         DO I=1,NUNIT
             WRITE(IU(I), '(A)') CBUFF(1:LPREF+1)
-   40    CONTINUE
+      END DO
          RETURN
       ENDIF
@@ -179,5 +179,5 @@
          LPIECE = MIN(LWRAP, LENMSG+1-NEXTC)
          IF (LPIECE < LENMSG+1-NEXTC) THEN
-            DO 52 I=LPIECE+1,2,-1
+            DO I=LPIECE+1,2,-1
                IF (MESSG(NEXTC+I-1:NEXTC+I-1) == ' ') THEN
                   LPIECE = I-1
@@ -185,5 +185,5 @@
                   GOTO 54
                ENDIF
-   52       CONTINUE
+      END DO
          ENDIF
    54    CBUFF(LPREF+1:LPREF+LPIECE) = MESSG(NEXTC:NEXTC+LPIECE-1)
@@ -202,5 +202,5 @@
          IDELTA = 0
          LPIECE = LWRAP
-         DO 56 I=LPIECE+1,2,-1
+         DO I=LPIECE+1,2,-1
             IF (MESSG(NEXTC+I-1:NEXTC+I-1) == ' ') THEN
                LPIECE = I-1
@@ -208,5 +208,5 @@
                GOTO 58
             ENDIF
-   56    CONTINUE
+      END DO
    58    CBUFF(LPREF+1:LPREF+LPIECE) = MESSG(NEXTC:NEXTC+LPIECE-1)
          NEXTC = NEXTC + LPIECE + IDELTA
@@ -223,7 +223,7 @@
 C       PRINT
 C
-      DO 60 I=1,NUNIT
+      DO I=1,NUNIT
          WRITE(IU(I), '(A)') CBUFF(1:LPREF+LPIECE)
-   60 CONTINUE
+      END DO
 C
       IF (NEXTC <= LENMSG) GO TO 50
Index: LMDZ6/branches/Amaury_dev/libf/misc/xersve.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/xersve.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/xersve.F	(revision 5086)
@@ -81,5 +81,5 @@
 C
          CALL XGETUA (LUN, NUNIT)
-         DO 20 KUNIT = 1,NUNIT
+         DO KUNIT = 1,NUNIT
             IUNIT = LUN(KUNIT)
             IF (IUNIT==0) IUNIT = I1MACH(4)
@@ -91,8 +91,8 @@
 C           Print body of table.
 C
-            DO 10 I = 1,NMSG
+            DO I = 1,NMSG
                WRITE (IUNIT,9010) LIBTAB(I), SUBTAB(I), MESTAB(I),
      *            NERTAB(I),LEVTAB(I),KOUNT(I)
-   10       CONTINUE
+      END DO
 C
 C           Print number of other errors.
@@ -100,5 +100,5 @@
             IF (KOUNTX/=0) WRITE (IUNIT,9020) KOUNTX
             WRITE (IUNIT,9030)
-   20    CONTINUE
+      END DO
 C
 C        Clear the error tables.
@@ -117,5 +117,5 @@
          SUB = SUBROU
          MES = MESSG
-         DO 30 I = 1,NMSG
+         DO I = 1,NMSG
             IF (LIB==LIBTAB(I) .AND. SUB==SUBTAB(I) .AND.
      *         MES==MESTAB(I) .AND. NERR==NERTAB(I) .AND.
@@ -125,5 +125,5 @@
                   RETURN
             ENDIF
-   30    CONTINUE
+      END DO
 C
          IF (NMSG<LENTAB) THEN
Index: LMDZ6/branches/Amaury_dev/libf/misc/xgetua.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/misc/xgetua.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/misc/xgetua.F	(revision 5086)
@@ -45,9 +45,9 @@
 C***FIRST EXECUTABLE STATEMENT  XGETUA
       N = J4SAVE(5,0,.FALSE.)
-      DO 30 I=1,N
+      DO I=1,N
          INDEX = I+4
          IF (I==1) INDEX = 3
          IUNITA(I) = J4SAVE(INDEX,0,.FALSE.)
-   30 CONTINUE
+      END DO
       RETURN
       END
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_conserv_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_conserv_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_conserv_m.F90	(revision 5086)
@@ -97,5 +97,5 @@
              vt(it) = vt(it) + (xs(is + 1) - xs(is)) * vs(is)
              is = is + 1
-          end do
+          END DO
           ! 1 <= is <= ns
           vt(it) = (vt(it) + mean_lin(xs(is), xt(it + 1)) * (xt(it + 1) &
@@ -105,5 +105,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   contains
@@ -170,5 +170,5 @@
              vt(it, :) = vt(it, :) + (xs(is + 1) - xs(is)) * vs(is, :)
              is = is + 1
-          end do
+          END DO
           ! 1 <= is <= ns
           vt(it, :) = (vt(it, :) + mean_lin(xs(is), xt(it + 1)) * (xt(it + 1) &
@@ -178,5 +178,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   contains
@@ -246,5 +246,5 @@
              vt(it, :, :) = vt(it, :, :) + (xs(is + 1) - xs(is)) * vs(is, :, :)
              is = is + 1
-          end do
+          END DO
           ! 1 <= is <= ns
           vt(it, :, :) = (vt(it, :, :) + mean_lin(xs(is), xt(it + 1)) &
@@ -254,5 +254,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   contains
@@ -324,5 +324,5 @@
                   * vs(is, :, :, :)
              is = is + 1
-          end do
+          END DO
           ! 1 <= is <= ns
           vt(it, :, :, :) = (vt(it, :, :, :) + mean_lin(xs(is), xt(it + 1)) &
@@ -332,5 +332,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   contains
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_lint_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_lint_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_lint_m.F90	(revision 5086)
@@ -53,5 +53,5 @@
        vt(it) = ((xs(is_b+1) - xt(it)) * vs(is_b) &
             + (xt(it) - xs(is_b)) * vs(is_b+1)) / (xs(is_b+1) - xs(is_b))
-    end do
+    END DO
 
   end function regr11_lint
@@ -92,5 +92,5 @@
        vt(it, :) = ((xs(is_b+1) - xt(it)) * vs(is_b, :) &
             + (xt(it) - xs(is_b)) * vs(is_b+1, :)) / (xs(is_b+1) - xs(is_b))
-    end do
+    END DO
 
   end function regr12_lint
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_step_av_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_step_av_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr1_step_av_m.F90	(revision 5086)
@@ -73,5 +73,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it) = (vt(it) + (xt(it + 1) - left_edge) * vs(is)) &
@@ -79,5 +79,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr11_step_av
@@ -133,5 +133,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :) = (vt(it, :) + (xt(it + 1) - left_edge) * vs(is, :)) &
@@ -139,5 +139,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr12_step_av
@@ -194,5 +194,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :, :) = (vt(it, :, :) &
@@ -200,5 +200,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr13_step_av
@@ -256,5 +256,5 @@
           is = is + 1
           left_edge = xs(is)
-       end do
+       END DO
        ! 1 <= is <= ns
        vt(it, :, :, :) = (vt(it, :, :, :) + (xt(it + 1) - left_edge) &
@@ -262,5 +262,5 @@
        if (xs(is + 1) == xt(it + 1)) is = is + 1
        ! 1 <= is <= ns .or. it == nt
-    end do
+    END DO
 
   end function regr14_step_av
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr3_lint_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr3_lint_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr3_lint_m.F90	(revision 5086)
@@ -53,5 +53,5 @@
        vt(:, :, it) = ((xs(is_b+1) - xt(it)) * vs(:, :, is_b) &
             + (xt(it) - xs(is_b)) * vs(:, :, is_b+1)) / (xs(is_b+1) - xs(is_b))
-    end do
+    END DO
 
   end function regr33_lint
@@ -94,5 +94,5 @@
             + (xt(it) - xs(is_b)) * vs(:, :, is_b+1, :)) &
             / (xs(is_b+1) - xs(is_b))
-    end do
+    END DO
 
   end function regr34_lint
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr_lat_time_climoz_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr_lat_time_climoz_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr_lat_time_climoz_m.F90	(revision 5086)
@@ -249,5 +249,5 @@
              do while (o3_in(j, 1, l, m) == missing_value)
                 j = j + 1
-             end do
+             END DO
              if (j > 1) o3_in(:j-1, :, l, m) = &
                   spread(o3_in(j, :, l, m), dim=1, ncopies=j-1)
@@ -257,5 +257,5 @@
              do while (o3_in(j, 1, l, m) == missing_value)
                 j = j - 1
-             end do
+             END DO
              if (j < n_lat) o3_in(j+1:, :, l, m) = &
                   spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j)
@@ -270,16 +270,16 @@
                 do while  (o3_in(j, k, l, m) /= missing_value .and. k < n_plev)
                    k = k + 1
-                end do
+                END DO
                 ! Replace missing values with the valid value at the
                 ! lowest level above missing values:
                 if (o3_in(j, k, l, m) == missing_value) &
                      o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m)
-             end do
-          end do
+             END DO
+          END DO
        else
           print *, "regr_lat_time_climoz: field ", m, &
                ", no missing value attribute"
        end if
-    end do
+    END DO
 
     call nf95_close(ncid_in)
@@ -322,5 +322,5 @@
        call nf95_put_var(ncid_out, varid_out(m), o3_out(nbp_lat:1:-1, :, :, m))
        ! (The order of "rlatu" is inverted in the output file)
-    end do
+    END DO
 
     call nf95_close(ncid_out)
Index: LMDZ6/branches/Amaury_dev/libf/obsolete/regr_pr_av_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/obsolete/regr_pr_av_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/obsolete/regr_pr_av_m.F90	(revision 5086)
@@ -99,5 +99,5 @@
           call handle_err("regr_pr_av nf90_get_var " // trim(name(i)), ncerr, &
                ncid)
-       end do
+       END DO
        
        ! Latitudes are in ascending order in the input file while
@@ -118,5 +118,5 @@
             slopes(v2(i, :, :), press_in_edg))
        ! (invert order of indices because "paprs" is in descending order)
-    end do
+    END DO
 
   end subroutine regr_pr_av
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/Ocean_skin/bulk_flux_m.F90	(revision 5086)
@@ -140,5 +140,5 @@
                   rain = null_array, qcol = rnl + hf + hlb - dels)
           end if
-       end do
+       END DO
     else
        tkt = 0.
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/StratAer/sulfate_aer_mod.F90	(revision 5086)
@@ -731,7 +731,7 @@
          JX=0
        ELSE
-         DO 10 I=1,N
+         DO I=1,N
            IF (X<XC(I)) GO TO 20
- 10      CONTINUE
+         END DO
          IER=1
  20      JX=I-1
@@ -756,7 +756,7 @@
          JX=0
        ELSE
-         DO 10 I=1,N
+         DO I=1,N
            IF (XT>X(I)) GO TO 20
- 10      CONTINUE
+         END DO
  20      JX=I
        ENDIF
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/concvl.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/concvl.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/concvl.F90	(revision 5086)
@@ -340,5 +340,5 @@
 ! em_wght(k)=wght_th(i,k)
 ! print*,'em_wght=',em_wght(k),wght_th(i,k)
-! end do
+! END DO
 ! END DO
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/icarus.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/icarus.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/icarus.F	(revision 5086)
@@ -364,5 +364,5 @@
       enddo 
 
-      do 12 ilev=1,nlev
+      do ilev=1,nlev
         do j=1,npoints 
          if (pfull(j,ilev) < 40000. .and.
@@ -375,12 +375,12 @@
            end if
         enddo
-12    continue
-
-      do 13 ilev=1,nlev
+      END DO
+
+      do ilev=1,nlev
         do j=1,npoints 
            if (at(j,ilev) > atmax(j) .and.
      &              ilev  >= itrop(j)) atmax(j)=at(j,ilev)
         enddo
-13    continue
+      END DO
 
       end if
@@ -391,10 +391,10 @@
               meantb(j) = 0.
 	      meantbclr(j) = 0. 
-          end do
+          END DO
       else
           do j=1,npoints
               meantb(j) = output_missing_value
        	      meantbclr(j) = output_missing_value
-          end do
+          END DO
       end if
       
@@ -466,5 +466,5 @@
  
       !initialize tau and albedocld to zero
-      do 15 ibox=1,ncol
+      do ibox=1,ncol
         do j=1,npoints 
             tau(j,ibox)=0.
@@ -474,5 +474,5 @@
           box_cloudy(j,ibox)=.false.
         enddo
-15    continue
+      END DO
 
       !compute total cloud optical depth for each column     
@@ -541,5 +541,5 @@
         if (ncolprint /= 0)
      &         write(6,*)  'ilev   pw (kg/m2)   tauwv(j)      dem_wv'
-        do 125 ilev=1,nlev
+        do ilev=1,nlev
           do j=1,npoints 
                !press and dpress are dyne/cm2 = Pascals *10
@@ -568,5 +568,5 @@
                enddo
              endif
-125     continue
+      END DO
 
         !initialize variables
@@ -740,5 +740,5 @@
             bb(j)=1/( exp(1307.27/skt(j)) - 1. )
             !bb(j)=5.67e-8*skt(j)**4
-          end do
+          END DO
 
         do ibox=1,ncol
@@ -751,6 +751,6 @@
      &         * trans_layers_above(j,ibox) 
             
-          end do
-        end do
+          END DO
+        END DO
 
         !calculate mean infrared brightness temperature
@@ -758,9 +758,9 @@
           do j=1,npoints 
             meantb(j) = meantb(j)+1307.27/(log(1.+(1./fluxtop(j,ibox))))
-	  end do
-        end do
+	  END DO
+        END DO
 	  do j=1, npoints
 	    meantb(j) = meantb(j) / real(ncol)
-	  end do        
+	  END DO
 
         if (ncolprint/=0) then
@@ -784,5 +784,5 @@
           write (6,'(8f7.2)') (meantb(j),ibox=1,ncolprint)
       
-          end do
+          END DO
       endif
     
@@ -925,5 +925,5 @@
 
       !compute cloud top pressure
-      do 30 ibox=1,ncol
+      do ibox=1,ncol
         !segregate according to optical thickness
         if (top_height == 1 .or. top_height == 3) then
@@ -933,5 +933,5 @@
             nmatch(j)=0
           enddo
-          do 29 k1=1,nlev-1
+          do k1=1,nlev-1
 	    if (top_height_direction == 2) then
 	      ilev = nlev - k1 
@@ -951,5 +951,5 @@
 	     end if                         
             enddo
-29        continue
+      END DO
 
           do j=1,npoints 
@@ -992,6 +992,6 @@
               levmatch(j,ibox)=ilev
               end if
-            end do
-          end do
+            END DO
+          END DO
         end if                            
           
@@ -1003,5 +1003,5 @@
         enddo
 
-30    continue
+      END DO
               
 !
@@ -1032,5 +1032,5 @@
 
       !reset frequencies
-      do 38 ilev=1,7
+      do ilev=1,7
       do 38 ilev2=1,7
         do j=1,npoints ! 
@@ -1042,4 +1042,5 @@
         enddo
 38    continue
+      END DO
 
       !reset variables need for averaging cloud properties
@@ -1060,5 +1061,5 @@
       boxarea = 1./real(ncol)
      
-      do 39 ibox=1,ncol
+      do ibox=1,ncol
         do j=1,npoints 
 
@@ -1166,5 +1167,5 @@
                        
         enddo ! j
-39    continue
+      END DO
       
       !compute mean cloud properties
@@ -1227,5 +1228,5 @@
      &           (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev) 
      &           ,(cchar(acc(ilev,ibox)+1),ilev=1,nlev) 
-             end do
+             END DO
              close(9)
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/math_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/math_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/math_lib.F90	(revision 5086)
@@ -273,5 +273,5 @@
        exit
     end if
-  end do
+  END DO
   
   if (lerror) then
@@ -316,5 +316,5 @@
     end if
     ilo = ilo + 1
-  end do
+  END DO
 
   ilo = max ( 2, ilo )
@@ -326,5 +326,5 @@
     end if
     ihi = ihi - 1
-  end do
+  END DO
   
   ihi = min ( ihi, ntab - 1 )
@@ -374,5 +374,5 @@
     syl = x2
  
-  end do
+  END DO
  
   result = sum1 &
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp_modis_simulator.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp_modis_simulator.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_cosp_modis_simulator.F90	(revision 5086)
@@ -176,7 +176,7 @@
               opticalThickness(i, j, k) = 0.   
             end if 
-          end do 
-        end do
-      end do
+          END DO
+        END DO
+      END DO
 
       !
@@ -197,7 +197,7 @@
           do i = 1, nSunlit
             if(subCols%frac_out(sunlit(i), j, k) == I_CVC) opticalThickness(i, j, k) = gridBox%dtau_c(sunlit(i), k)
-          end do 
-        end do
-      end do
+          END DO
+        END DO
+      END DO
 
       !
@@ -220,5 +220,5 @@
                                 retrievedPhase(i, :), retrievedCloudTopPressure(i, :),      & 
                                 retrievedTau(i, :), retrievedSize(i, :))
-     end do
+     END DO
      
       ! DJS2015: Call L3 modis simulator used by cospv2.0
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_modis_sim.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_modis_sim.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/mod_modis_sim.F90	(revision 5086)
@@ -331,5 +331,5 @@
         retrievedTau(i) = R_UNDEF 
       end if
-    end do
+    END DO
     where((retrievedSize(:) < 0.).and.(retrievedSize(:) /= R_UNDEF)) retrievedSize(:) = 1.0e-06*re_fill
 
@@ -802,5 +802,5 @@
         tauMask(:, :, i) = .false.
       end where
-    end do 
+    END DO
 
     do i = 1, numPressureHistogramBins 
@@ -811,5 +811,5 @@
         pressureMask(:, :, i) = .false.
       end where
-    end do 
+    END DO
     
     do i = 1, numPressureHistogramBins
@@ -817,6 +817,6 @@
         Optical_Thickness_vs_Cloud_Top_Pressure(:, j, i) = & 
           real(count(tauMask(:, :, j) .and. pressureMask(:, :, i), dim = 2)) / real(nSubcols)
-      end do 
-    end do 
+      END DO
+    END DO
     
   end subroutine modis_L3_simulator
@@ -851,5 +851,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
     cloud_top_pressure = totalProduct/totalTau
   end function cloud_top_pressure
@@ -877,5 +877,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
     weight_by_extinction = totalProduct/totalTau
   end function weight_by_extinction
@@ -1114,5 +1114,5 @@
     do i = 1, size(cloudIndicies)
       call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
-    end do 
+    END DO
                     
     call adding_doubling(Refl(:), Trans(:), Refl_tot, Trans_tot)  
@@ -1292,5 +1292,5 @@
           Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1 - Refl_cumulative(i-1) * Refl(i))
           Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1 - Refl_cumulative(i-1) * Refl(i))
-      end do
+      END DO
       
       Refl_tot = Refl_cumulative(size(Refl))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/optics_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/optics_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/optics_lib.F90	(revision 5086)
@@ -519,7 +519,4 @@
 
 !   // region from 0.045 microns to 167.0 microns - no temperature depend
-    do i=2,nwl
-      if(alam < wl(i)) continue
-    enddo
     x1=log(wl(i-1))
     x2=log(wl(i))
@@ -539,12 +536,12 @@
     if(tk > temref(1)) tk=temref(1)
     if(tk < temref(4)) tk=temref(4)
-    do 11 i=2,4
+    do i=2,4
       if(tk>=temref(i)) go to 12
-    11 continue
+    END DO
     12 lt1=i
     lt2=i-1
-    do 13 i=2,nwlt
+    do i=2,nwlt
       if(alam<=wlt(i)) go to 14
-    13 continue
+    END DO
     14 x1=log(wlt(i-1))
     x2=log(wlt(i))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/scops.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/scops.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/scops.F	(revision 5086)
@@ -161,5 +161,5 @@
       
       !loop over vertical levels
-      DO 200 ilev = 1,nlev
+      DO ilev = 1,nlev
                                   
 !     Initialise threshold
@@ -331,5 +331,5 @@
           endif
 
-200   CONTINUE    !loop over nlev
+      END DO    !loop over nlev
 
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/zeff.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/zeff.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/zeff.F90	(revision 5086)
@@ -113,5 +113,5 @@
       call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
         dg, xs1, xs2, dph, err)
-    end do
+    END DO
     
   else
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/cosp.F90	(revision 5086)
@@ -731,5 +731,5 @@
                                   modisRetrievedCloudTopPressure(i,:),                   &
                                   modisRetrievedTau(i,:),modisRetrievedSize(i,:))
-          end do
+          END DO
        endif
     endif
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/icarus.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/icarus.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/icarus.F90	(revision 5086)
@@ -155,5 +155,5 @@
                   (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev),&
                   (cchar(acc(ilev,ibox)+1),ilev=1,nlev) 
-          end do
+          END DO
           close(9)
 
@@ -322,5 +322,5 @@
        do ibox=1,ncol
           fluxtop(1:npoints,ibox) = fluxtop(1:npoints,ibox) + emsfc_lw*bb(1:npoints)*trans_layers_above(1:npoints,ibox) 
-       end do
+       END DO
 
        ! All Sky brightness temperature
@@ -445,5 +445,5 @@
                 levmatch(1:npoints,ibox)=ilev
              endwhere
-          end do
+          END DO
        end if
        where(tau(1:npoints,ibox) <= tauchk)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/math_lib.F90	(revision 5086)
@@ -209,5 +209,5 @@
           exit
        end if
-    end do
+    END DO
     
     if (lerror) then
@@ -244,5 +244,5 @@
        end if
        ilo = ilo + 1
-    end do
+    END DO
     
     ilo = max ( 2, ilo )
@@ -254,5 +254,5 @@
        end if
        ihi = ihi - 1
-    end do
+    END DO
     
     ihi = min ( ihi, ntab - 1 )
@@ -305,5 +305,5 @@
        syl = x2
        
-    end do
+    END DO
 
     result = sum1 + ca * ( b**3 - syl**3 ) / 3._wp &
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/modis_simulator.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/modis_simulator.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/modis_simulator.F90	(revision 5086)
@@ -222,5 +222,5 @@
           retrievedTau(i)              = R_UNDEF 
        end if
-    end do
+    END DO
     where((retrievedSize(1:nSubCols) < 0.).and.(retrievedSize(1:nSubCols) /= R_UNDEF)) &
          retrievedSize(1:nSubCols) = 1.0e-06_wp*re_fill
@@ -455,5 +455,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
 
     if (totalTau > 0._wp) then
@@ -489,5 +489,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
 
     if (totalTau > 0._wp) then
@@ -715,5 +715,5 @@
     do i = 1, size(cloudIndicies)
        call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
-    end do
+    END DO
     
     call adding_doubling(count(tau(1:nLevels) > 0),Refl(:), Trans(:), Refl_tot, Trans_tot)  
@@ -897,5 +897,5 @@
        Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i))
        Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1._wp - Refl_cumulative(i-1) * Refl(i))
-    end do
+    END DO
     
     Refl_tot = Refl_cumulative(size(Refl))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/optics_lib.F90	(revision 5086)
@@ -539,7 +539,4 @@
     if (alam < cutice) then
        ! Region from 0.045 microns to 167.0 microns - no temperature depend
-       do i=2,nwl
-          if(alam < wl(i)) continue
-       enddo
        x1  = log(wl(i-1))
        x2  = log(wl(i))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/quickbeam_optics.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/quickbeam_optics.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cosp2/quickbeam_optics.F90	(revision 5086)
@@ -1009,5 +1009,5 @@
           call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
                dg, xs1, xs2, dph, err)
-       end do
+       END DO
 
     else
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/cosp.F90	(revision 5086)
@@ -886,5 +886,5 @@
                                   modisRetrievedCloudTopPressure(i,:),                   &
                                   modisRetrievedTau(i,:),modisRetrievedSize(i,:))
-          end do
+          END DO
        endif
     endif
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/icarus.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/icarus.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/icarus.F90	(revision 5086)
@@ -155,5 +155,5 @@
                   (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev),&
                   (cchar(acc(ilev,ibox)+1),ilev=1,nlev) 
-          end do
+          END DO
           close(9)
 
@@ -322,5 +322,5 @@
        do ibox=1,ncol
           fluxtop(1:npoints,ibox) = fluxtop(1:npoints,ibox) + emsfc_lw*bb(1:npoints)*trans_layers_above(1:npoints,ibox) 
-       end do
+       END DO
 
        ! All Sky brightness temperature
@@ -445,5 +445,5 @@
                 levmatch(1:npoints,ibox)=ilev
              endwhere
-          end do
+          END DO
        end if
        where(tau(1:npoints,ibox) <= tauchk)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/lmdz_cosp_subsample_and_optics_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/lmdz_cosp_subsample_and_optics_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/lmdz_cosp_subsample_and_optics_mod.F90	(revision 5086)
@@ -365,6 +365,6 @@
              endif
              cospIN%g_vol_cloudsat(i,:,j)=g_vol(i,j)
-          end do
-       end do
+          END DO
+       END DO
        
        ! Loop over all subcolumns
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/math_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/math_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/math_lib.F90	(revision 5086)
@@ -209,5 +209,5 @@
           exit
        end if
-    end do
+    END DO
     
     if (lerror) then
@@ -244,5 +244,5 @@
        end if
        ilo = ilo + 1
-    end do
+    END DO
     
     ilo = max ( 2, ilo )
@@ -254,5 +254,5 @@
        end if
        ihi = ihi - 1
-    end do
+    END DO
     
     ihi = min ( ihi, ntab - 1 )
@@ -305,5 +305,5 @@
        syl = x2
        
-    end do
+    END DO
 
     result = sum1 + ca * ( b**3 - syl**3 ) / 3._wp &
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/modis_simulator.F90	(revision 5086)
@@ -222,5 +222,5 @@
           retrievedTau(i)              = R_UNDEF 
        end if
-    end do
+    END DO
     where((retrievedSize(1:nSubCols) < 0.).and.(retrievedSize(1:nSubCols) /= R_UNDEF)) &
          retrievedSize(1:nSubCols) = 1.0e-06_wp*re_fill
@@ -455,5 +455,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
 
     if (totalTau > 0._wp) then
@@ -489,5 +489,5 @@
       end if 
       if(totalTau >= tauLimit) exit
-    end do 
+    END DO
 
     if (totalTau > 0._wp) then
@@ -715,5 +715,5 @@
     do i = 1, size(cloudIndicies)
        call two_stream(tau(cloudIndicies(i)), g(cloudIndicies(i)), w0(cloudIndicies(i)), Refl(i), Trans(i))
-    end do
+    END DO
     
     call adding_doubling(count(tau(1:nLevels) > 0),Refl(:), Trans(:), Refl_tot, Trans_tot)  
@@ -897,5 +897,5 @@
        Refl_cumulative(i) = Refl_cumulative(i-1) + Refl(i)*(Tran_cumulative(i-1)**2)/(1._wp - Refl_cumulative(i-1) * Refl(i))
        Tran_cumulative(i) = (Tran_cumulative(i-1)*Tran(i)) / (1._wp - Refl_cumulative(i-1) * Refl(i))
-    end do
+    END DO
     
     Refl_tot = Refl_cumulative(size(Refl))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/optics_lib.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/optics_lib.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/optics_lib.F90	(revision 5086)
@@ -539,7 +539,4 @@
     if (alam < cutice) then
        ! Region from 0.045 microns to 167.0 microns - no temperature depend
-       do i=2,nwl
-          if(alam < wl(i)) continue
-       enddo
        x1  = log(wl(i-1))
        x2  = log(wl(i))
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/quickbeam_optics.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/quickbeam_optics.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cospv2/quickbeam_optics.F90	(revision 5086)
@@ -983,5 +983,5 @@
           call mieint(sizep(i), m0(i), one, dqv, qext(i), dqsc, qbsca(i), &
                dg, xs1, xs2, dph, err)
-       end do
+       END DO
 
     else
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cv30_routines.F90	(revision 5086)
@@ -1548,5 +1548,5 @@
             ! !!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
             ! !!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
-            ! !!!      end do
+            ! !!!      END DO
             elij(il, i, j) = altem
             elij(il, i, j) = amax1(0.0, elij(il,i,j))
@@ -2144,5 +2144,5 @@
             ! :            +tra(il,i,j)*(mp(il,i)-mp(il,i+1))
             ! trap(il,i,j)=trap(il,i,j)/mp(il,i)
-            ! end do
+            ! END DO
 
           ELSE
@@ -2161,5 +2161,5 @@
               ! do j=1,ntra
               ! trap(il,i,j)=trap(il,i+1,j)
-              ! end do
+              ! END DO
 
             END IF
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/cv3_routines.F90	(revision 5086)
@@ -2431,5 +2431,5 @@
 !!!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
 !!!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
-!!!!      end do
+!!!!      END DO
             elij(il, i, j) = altem
             elij(il, i, j) = max(0.0, elij(il,i,j))
@@ -3424,5 +3424,5 @@
 !AC!        endif ! (i.lt.inb(il) .and. lwork(il))
 !AC!       enddo
-!AC!      end do
+!AC!      END DO
 
 400 END DO
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1DUTILS_read_interp.h	(revision 5086)
@@ -2010,5 +2010,5 @@
           read (ilesfile,*) height(k),ugprof(k),vgprof(k),wfls(k),         &
      &                      dqtdxls(k),dqtdyls(k),dqtdtls(k),thlpcar(k)
-        end do
+        END DO
         do k=1,kmax
           if (height(k) .ne. height1(k)) then
@@ -2017,5 +2017,5 @@
             stop
           endif
-        end do
+        END DO
         close(ilesfile)
 
@@ -2035,5 +2035,5 @@
         do k=1,kmax
           read (ilesfile,*) height(k),(tracer(k,itrac),itrac=nt1,nt2)
-        end do
+        END DO
         close(ilesfile)
         endif
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/inlandsis/surf_inlandsis_mod.F90	(revision 5086)
@@ -470,5 +470,5 @@
                     depth = depth - dzsnSV(ikl,isl) / 2.
                     
-                end do
+                END DO
 
             END DO
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/isccp_cloud_types.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/isccp_cloud_types.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/isccp_cloud_types.F90	(revision 5086)
@@ -1027,5 +1027,5 @@
     ! write (6,'(a)') '100.*f:'
     ! write (6,'(8f7.2)') (100.*fluxtop(j,ibox),ibox=1,ncolprint)
-    ! end do
+    ! END DO
     ! endif
 
@@ -1569,5 +1569,5 @@
   ! &           (cchar_realtops(acc(ilev,ibox)+1),ilev=1,nlev)
   ! &           ,(cchar(acc(ilev,ibox)+1),ilev=1,nlev)
-  ! end do
+  ! END DO
   ! close(9)
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_wake.F90	(revision 5086)
@@ -1507,5 +1507,5 @@
     ! c      do i=1,klon
     ! c         print*,alpha(i)
-    ! c      end do
+    ! c      END DO
     ! cc
     DO k = 1, klev
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/o3_chem_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/o3_chem_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/o3_chem_m.F90	(revision 5086)
@@ -24,5 +24,5 @@
 
     integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
-    real, intent(in):: gmtime ! heure de la journée en fraction de jour
+    real, intent(in):: gmtime ! heure de la journ�e en fraction de jour
     real, intent(in):: t_seri(:, :) ! (klon, nbp_lev) temperature, in K
 
@@ -59,5 +59,5 @@
     real earth_long
     ! (longitude vraie de la Terre dans son orbite solaire, par
-    ! rapport au point vernal (21 mars), en degrés)
+    ! rapport au point vernal (21 mars), en degr�s)
 
     real pmu0(klon) ! mean of cosine of solar zenith angle during "pdtphys"
@@ -163,5 +163,5 @@
     do k =  nbp_lev - 1, 1, -1
        sigma_mass(:, k) = sigma_mass(:, k+1) + zmasse(:, k) * q(:, k)
-    end do
+    END DO
 
     o3_prod = c + b * q + a6_mass * sigma_mass
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/phyetat0_mod.F90	(revision 5086)
@@ -321,5 +321,5 @@
         IF (.NOT. found) THEN
            PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
-           PRINT*, "          Il prend donc la valeur de surface"
+           PRINT*, "          Il prEND DOnc la valeur de surface"
            tsoil(:, isoil, :)=ftsol(:, :)
         ENDIF
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/physiq_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/physiq_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/physiq_mod.F90	(revision 5086)
@@ -1705,5 +1705,4 @@
 998          CONTINUE
              CLOSE(98)
-             CONTINUE
              IF(nCFMIP>npCFMIP) THEN
                 print*,'nCFMIP > npCFMIP : augmenter npCFMIP et recompiler'
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/regr_lat_time_coefoz_m.F90	(revision 5086)
@@ -182,5 +182,5 @@
     do i_v = 1, n_o3_param
        call nf95_inq_varid(ncid_in, trim(name_in(i_v)), varid_in(i_v))
-    end do
+    END DO
 
     ! Create the output file and get the variable IDs:
@@ -225,5 +225,5 @@
             o3_par_out(nbp_lat:1:-1, :, :))
        ! (The order of "rlatu" is inverted in the output file)
-    end do
+    END DO
 
     call nf95_close(ncid_out)
@@ -309,5 +309,5 @@
             & varid_out(i))
        call handle_err_copy_att("standard_name")
-    end do
+    END DO
 
     ! Global attributes:
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_int_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_int_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_int_m.F90	(revision 5086)
@@ -98,5 +98,5 @@
                          v3(i, nbp_lev:1:-1))
        ! (invert order of indices because "pplay" is in descending order)
-    end do
+    END DO
 
   end subroutine regr_pr_int
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_o3_m.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_o3_m.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/regr_pr_o3_m.F90	(revision 5086)
@@ -77,5 +77,5 @@
             p3d(1, j, nbp_lev + 1:1:-1), o3_mob_regr(1, j, nbp_lev:1:-1))
        ! (invert order of indices because "p3d" is in descending order)
-    end do
+    END DO
 
     ! Other latitudes:
@@ -85,6 +85,6 @@
                p3d(i, j, nbp_lev + 1:1:-1), o3_mob_regr(i, j, nbp_lev:1:-1))
           ! (invert order of indices because "p3d" is in descending order)
-       end do
-    end do
+       END DO
+    END DO
 
     ! Duplicate pole values on all longitudes:
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F	(revision 5086)
@@ -27,7 +27,7 @@
       REAL zmin,zmax,zllu,zllm
       write(F1000,'(''(4x,'',I3,''(1H-))'')')im+3
-      DO 10001 i=1,200
+      DO i=1,200
       jline(1+(i-1)*5:5*i)='.    '
-10001 CONTINUE
+      END DO
 10002 zmin=z(1,1)
       imin=1
@@ -37,6 +37,6 @@
       jmax=1
       kzero=0
-      DO 10003 j=1,jm
-      DO 10005 i=1,im
+      DO j=1,jm
+      DO i=1,im
       IF(.NOT.( z(i,j)>zmax))GOTO 10007
       zmax=z(i,j)
@@ -50,7 +50,7 @@
       kzero=kzero+1
 10011 CONTINUE
-10005 CONTINUE
+      END DO
 10006 CONTINUE
-10003 CONTINUE
+      END DO
 10004 zsign=(sign(1.,zmin)*sign(1.,zmax)>0.)
       WRITE(*,*)'>>> dump2ds: ',trim(nom_z)
@@ -73,6 +73,6 @@
       zinf=.false.
       znan=.false.
-      DO 10017 j=1,jm
-      DO 10019 i=1,im
+      DO j=1,jm
+      DO i=1,im
       az=abs(z(i,j))
       IF(.NOT.( az==0.))GOTO 10021
@@ -97,8 +97,8 @@
       kchar(i)=32-kchar(i)
 10027 CONTINUE
-10019 CONTINUE
+      END DO
 10020 WRITE(*,'(1x,i3,''|'',1000a)')j+1,(iform(kchar(i)),i=1,im),'|','| 
      *'
-10017 CONTINUE
+      END DO
 10018 write(*,F1000)
       WRITE(*,'(5x,1000i1)')(mod(i,10),i=1,im)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90	(revision 5086)
@@ -353,10 +353,10 @@
   do k = 2, n
     w = w * (x - k)
-  end do
+  END DO
 else
   w = 1
   do k = 0, -n - 1
     y = y * (x + k)
-  end do
+  END DO
 end if
 gamma_res = w / y
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/set99.F
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/set99.F	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/set99.F	(revision 5086)
@@ -15,9 +15,9 @@
       NIL=0
       NHL=(N/2)-1
-      DO 10 K=NIL,NHL
+      DO K=NIL,NHL
       ANGLE=FLOAT(K)*DEL
       TRIGS(2*K+1)=COS(ANGLE)
       TRIGS(2*K+2)=SIN(ANGLE)
-   10 CONTINUE
+      END DO
 C
 C     FIND FACTORS OF N (8,6,5,4,3,2; ONLY ONE 8 ALLOWED)
Index: LMDZ6/branches/Amaury_dev/libf/phylmd/simu_airs.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmd/simu_airs.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmd/simu_airs.F90	(revision 5086)
@@ -30,5 +30,5 @@
         do i=1,N-1
         slope(i)=-(T(i+1)-T(i))/(alt(i+1)-alt(i))
-        end do
+        END DO
         slope(N)=slope(N-1)
 
@@ -60,5 +60,5 @@
         i=i+i_dir
         if (i<=1.or.i>=N) exit_flag=1
-        end do
+        END DO
 
         if (first_point<=0) P_tropo=65.4321
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/concvl.F90	(revision 5086)
@@ -415,5 +415,5 @@
 ! em_wght(k)=wght_th(i,k)
 ! print*,'em_wght=',em_wght(k),wght_th(i,k)
-! end do
+! END DO
 ! END DO
 
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv30_routines.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv30_routines.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv30_routines.F90	(revision 5086)
@@ -2032,5 +2032,5 @@
             ! !!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
             ! !!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
-            ! !!!      end do
+            ! !!!      END DO
             elij(il, i, j) = altem
             elij(il, i, j) = amax1(0.0, elij(il,i,j))
@@ -3155,5 +3155,5 @@
             ! :            +tra(il,i,j)*(mp(il,i)-mp(il,i+1))
             ! trap(il,i,j)=trap(il,i,j)/mp(il,i)
-            ! end do
+            ! END DO
 
           ELSE
@@ -3172,5 +3172,5 @@
               ! do j=1,ntra
               ! trap(il,i,j)=trap(il,i+1,j)
-              ! end do
+              ! END DO
 
             END IF
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/cv3_routines.F90	(revision 5086)
@@ -2965,5 +2965,5 @@
 !!!!      traent(il,i,j,k)=sij(il,i,j)*tra(il,i,k)
 !!!!     :      +(1.-sij(il,i,j))*tra(il,nk(il),k)
-!!!!      end do
+!!!!      END DO
             elij(il, i, j) = altem
             elij(il, i, j) = max(0.0, elij(il,i,j))
@@ -4485,5 +4485,5 @@
 !AC!        endif ! (i.lt.inb(il) .and. lwork(il))
 !AC!       enddo
-!AC!      end do
+!AC!      END DO
 
 #ifdef ISO
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/lmdz_wake.F90	(revision 5086)
@@ -127,10 +127,10 @@
   ! wdens_ref: initial number of wakes per unit area (3D) or per
   ! unit length (2D), at the beginning of each time step
-  ! Tgw    : 1 sur la période de onde de gravité
-  ! Cgw    : vitesse de propagation de onde de gravité
+  ! Tgw    : 1 sur la p�riode de onde de gravit�
+  ! Cgw    : vitesse de propagation de onde de gravit�
   ! LL     : distance entre 2 poches
 
   ! -------------------------------------------------------------------------
-  ! Déclaration de variables
+  ! D�claration de variables
   ! -------------------------------------------------------------------------
 
@@ -196,5 +196,5 @@
   ! -------------------
 
-  ! Variables à fixer
+  ! Variables � fixer
   INTEGER, SAVE                                         :: igout
   !$OMP THREADPRIVATE(igout)
@@ -383,7 +383,7 @@
   ! Configuration de coefgw,stark,wdens (22/02/06 by YU Jingmei)
 
-  ! coefgw : Coefficient pour les ondes de gravité
+  ! coefgw : Coefficient pour les ondes de gravit�
   ! stark : Coefficient k dans Cstar=k*sqrt(2*WAPE)
-  ! wdens : Densité surfacique de poche froide
+  ! wdens : Densit� surfacique de poche froide
   ! -------------------------------------------------------------------------
 
@@ -1088,6 +1088,6 @@
 
     ! cc nrlmd   Ajout d'un recalcul de wdens dans le cas d'un entrainement
-    ! négatif de ktop à kupper --------
-    ! cc           On calcule pour cela une densité wdens0 pour laquelle on
+    ! n�gatif de ktop � kupper --------
+    ! cc           On calcule pour cela une densit� wdens0 pour laquelle on
     ! aurait un entrainement nul ---
     !jyg<
@@ -1096,5 +1096,5 @@
     ! des descentes unsaturees. Nous faisons alors l'hypothese que la 
     ! convection profonde cree directement de nouvelles poches, sans passer 
-    ! par les thermiques. La nouvelle valeur de wdens est alors imposée.
+    ! par les thermiques. La nouvelle valeur de wdens est alors impos�e.
 
     DO i = 1, klon
@@ -1195,5 +1195,5 @@
       DO i = 1, klon
         IF (wk_adv(i)) THEN
-          ! cc nrlmd          Introduction du taux de mortalité des poches et
+          ! cc nrlmd          Introduction du taux de mortalit� des poches et
           ! test sur sigmaw_max=0.4
           ! cc         d_sigmaw(i) = gfl(i)*Cstar(i)*dtimesub
@@ -1308,5 +1308,5 @@
 
     ! c      DO i=1,klon
-    ! c        print*,'Pente entre 0 et kupper (référence)'
+    ! c        print*,'Pente entre 0 et kupper (r�f�rence)'
     ! c     $   	,omg(i,kupper(i)+1)/(pupper(i)-ph(i,1))
     ! c        print*,'Pente entre ktop et kupper'
@@ -1602,5 +1602,5 @@
 
 
-          ! Coefficient de répartition
+          ! Coefficient de r�partition
 
           crep(i, k) = crep_sol*(ph(i,kupper(i))-ph(i,k))/ &
@@ -1646,6 +1646,6 @@
 !
 
-          ! cc nrlmd          Prise en compte du taux de mortalité
-          ! cc               Définitions de entr, detr
+          ! cc nrlmd          Prise en compte du taux de mortalit�
+          ! cc               D�finitions de entr, detr
 !jyg<
 !!            detr(i, k) = 0.
@@ -1664,5 +1664,5 @@
 
 
-          ! ajout d'un effet onde de gravité -Tgw(k)*deltatw(k) 03/02/06 YU
+          ! ajout d'un effet onde de gravit� -Tgw(k)*deltatw(k) 03/02/06 YU
           ! Jingmei
 
@@ -1748,5 +1748,5 @@
     ! c      do i=1,klon
     ! c         print*,alpha(i)
-    ! c      end do
+    ! c      END DO
     ! cc
     DO k = 1, klev
@@ -2601,5 +2601,5 @@
         discrim = b*b - 4.*a*c
         ! print*, 'x, a, b, c, discrim', x, a, b, c, discrim
-        IF (a+b>=0.) THEN !! Condition suffisante pour la positivité de ovap
+        IF (a+b>=0.) THEN !! Condition suffisante pour la positivit� de ovap
           alpha1(i) = 1.
         ELSE
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyetat0_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyetat0_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/phyetat0_mod.F90	(revision 5086)
@@ -342,5 +342,5 @@
         IF (.NOT. found) THEN
            PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
-           PRINT*, "          Il prend donc la valeur de surface"
+           PRINT*, "          Il prEND DOnc la valeur de surface"
            tsoil(:, isoil, :)=ftsol(:, :)
         ENDIF
Index: LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90
===================================================================
--- LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/libf/phylmdiso/physiq_mod.F90	(revision 5086)
@@ -1869,5 +1869,4 @@
 998          CONTINUE
              CLOSE(98)
-             CONTINUE
              IF(nCFMIP>npCFMIP) THEN
                 print*,'nCFMIP > npCFMIP : augmenter npCFMIP et recompiler'
Index: LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_create_single.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_create_single.f90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_create_single.f90	(revision 5086)
@@ -47,6 +47,6 @@
           call nf95_put_att(ncid, varid_coord(i), coordinates(i)%attr_name(j), &
                coordinates(i)%attr_val(j))
-       end do
-    end do
+       END DO
+    END DO
 
     call nf95_def_var(ncid, name, NF90_FLOAT, dimids, varid)
Index: LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_find_coord.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_find_coord.f90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/tools/netcdf95/Datasets/nf95_find_coord.f90	(revision 5086)
@@ -96,5 +96,5 @@
           end if
        end if
-    end do
+    END DO
 
     if (found) then
Index: LMDZ6/branches/Amaury_dev/tools/netcdf95/Groups/nf95_inq_file_ncid.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/tools/netcdf95/Groups/nf95_inq_file_ncid.f90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/tools/netcdf95/Groups/nf95_inq_file_ncid.f90	(revision 5086)
@@ -31,5 +31,5 @@
        if (ncerr_local /= nf95_noerr) exit
        ncid_file = parent_ncid
-    end do
+    END DO
 
     if (ncerr_local == NF95_ENOGRP) then
Index: LMDZ6/branches/Amaury_dev/tools/netcdf95/Variables/nf95_gw_var.f90
===================================================================
--- LMDZ6/branches/Amaury_dev/tools/netcdf95/Variables/nf95_gw_var.f90	(revision 5082)
+++ LMDZ6/branches/Amaury_dev/tools/netcdf95/Variables/nf95_gw_var.f90	(revision 5086)
@@ -152,5 +152,5 @@
     do i = 1, 4
        call nf95_inquire_dimension(ncid, dimids(i), nclen=len_dim(i))
-    end do
+    END DO
 
     allocate(values(len_dim(1), len_dim(2), len_dim(3), len_dim(4)))
@@ -186,5 +186,5 @@
     do i = 1, 5
        call nf95_inquire_dimension(ncid, dimids(i), nclen=len_dim(i))
-    end do
+    END DO
 
     allocate(values(len_dim(1), len_dim(2), len_dim(3), len_dim(4), len_dim(5)))
@@ -323,5 +323,5 @@
     do i = 1, 4
        call nf95_inquire_dimension(ncid, dimids(i), nclen=len_dim(i))
-    end do
+    END DO
 
     allocate(values(len_dim(1), len_dim(2), len_dim(3), len_dim(4)))
