Index: LMDZ6/trunk/libf/filtrez/inifgn.F
===================================================================
--- LMDZ6/trunk/libf/filtrez/inifgn.F	(revision 5079)
+++ LMDZ6/trunk/libf/filtrez/inifgn.F	(revision 5084)
@@ -28,42 +28,39 @@
       pi = 2.* ASIN(1.)
 C
-      DO i=1,iim
+      DO 5 i=1,iim
        dlonu(i)=  xprimu( i )
        dlonv(i)=  xprimv( i )
-      END DO
+   5  CONTINUE
 
-      DO i=1,iim
+      DO 12 i=1,iim
       sddv(i)   = SQRT(dlonv(i))
       sddu(i)   = SQRT(dlonu(i))
       unsddu(i) = 1./sddu(i)
       unsddv(i) = 1./sddv(i)
-      END DO
+  12  CONTINUE
 C
-      DO j=1,iim
-      DO i=1,iim
+      DO 17 j=1,iim
+      DO 17 i=1,iim
       vec(i,j)     = 0.
       vec1(i,j)    = 0.
       eignfnv(i,j) = 0.
       eignfnu(i,j) = 0.
-      END DO
-      END DO
+  17  CONTINUE
 c
 c
       eignfnv(1,1)    = -1.
       eignfnv(iim,1)  =  1.
-      DO i=1,imm1
+      DO 20 i=1,imm1
       eignfnv(i+1,i+1)= -1.
       eignfnv(i,i+1)  =  1.
-      END DO
-      DO j=1,iim
-      DO i=1,iim
+  20  CONTINUE
+      DO 25 j=1,iim
+      DO 25 i=1,iim
       eignfnv(i,j) = eignfnv(i,j)/(sddu(i)*sddv(j))
-      END DO
-      END DO
-      DO j=1,iim
-      DO i=1,iim
+  25  CONTINUE
+      DO 30 j=1,iim
+      DO 30 i=1,iim
       eignfnu(i,j) = -eignfnv(j,i)
-      END DO
-      END DO
+  30  CONTINUE
 c
 #ifdef CRAY