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
