--- trunk/bibio/formcoord.f	2013/11/15 18:45:49	76
+++ trunk/bibio/formcoord.f90	2014/03/05 14:38:41	81
@@ -1,54 +1,55 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/bibio/formcoord.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
-!
-      subroutine formcoord(unit,n,x,a,rev,text)
-      implicit none
-      integer n,unit,ndec
-      logical rev
-      real x(n),a
-      character*4 text
-
-      integer i,id,i1,i2,in
-      real dx,dxmin
-
-      if(rev) then
-         id=-1
-         i1=n
-         i2=n-1
-         in=1
-         write(unit,3000) text(1:1)
-      else
-         id=1
-         i1=1
-         i2=2
-         in=n
-      endif
-
-      if (n.lt.2) then
-         ndec=1
-         write(unit,1000) text,n,x(1)*a
-      else
-         dxmin=abs(x(2)-x(1))
-         do i=2,n-1
-            dx=abs(x(i+1)-x(i))
-            if (dx.lt.dxmin) dxmin=dx
-         enddo
-
-         ndec=-log10(dxmin)+2
-         if(mod(n,6).eq.1) then
-            write(unit,1000) text,n,x(i1)*a
-            write(unit,2000) (x(i)*a,i=i2,in,id)
-         else
-            write(unit,1000) text,n
-            write(unit,2000) (x(i)*a,i=i1,in,id)
-         endif
-      endif
-
-1000  format(a4,2x,i4,' LEVELS',43x,f12.2)
-2000  format(6f12.2)
-c1000  format(a4,2x,i4,' LEVELS',43x,f12.<ndec>)
-c2000  format(6f12.<ndec>)
-3000  format('FORMAT ',a1,'REV')
-      return
 
-      end
+! $Header: /home/cvsroot/LMDZ4/libf/bibio/formcoord.F,v 1.1.1.1 2004/05/19
+! 12:53:05 lmdzadmin Exp $
+
+SUBROUTINE formcoord(unit, n, x, a, rev, text)
+  IMPLICIT NONE
+  INTEGER n, unit, ndec
+  LOGICAL rev
+  REAL x(n), a
+  CHARACTER *4 text
+
+  INTEGER i, id, i1, i2, in
+  REAL dx, dxmin
+
+  IF (rev) THEN
+    id = -1
+    i1 = n
+    i2 = n - 1
+    in = 1
+    WRITE (unit, 3000) text(1:1)
+  ELSE
+    id = 1
+    i1 = 1
+    i2 = 2
+    in = n
+  END IF
+
+  IF (n<2) THEN
+    ndec = 1
+    WRITE (unit, 1000) text, n, x(1)*a
+  ELSE
+    dxmin = abs(x(2)-x(1))
+    DO i = 2, n - 1
+      dx = abs(x(i+1)-x(i))
+      IF (dx<dxmin) dxmin = dx
+    END DO
+
+    ndec = -log10(dxmin) + 2
+    IF (mod(n,6)==1) THEN
+      WRITE (unit, 1000) text, n, x(i1)*a
+      WRITE (unit, 2000)(x(i)*a, i=i2, in, id)
+    ELSE
+      WRITE (unit, 1000) text, n
+      WRITE (unit, 2000)(x(i)*a, i=i1, in, id)
+    END IF
+  END IF
+
+1000 FORMAT (A4, 2X, I4, ' LEVELS', 43X, F12.2)
+2000 FORMAT (6F12.2)
+  ! 1000  format(a4,2x,i4,' LEVELS',43x,f12.<ndec>)
+  ! 2000  format(6f12.<ndec>)
+3000 FORMAT ('FORMAT ', A1, 'REV')
+  RETURN
+
+END SUBROUTINE formcoord