* ------------------------------------------------------------------
* -------------- ex.f ----------------------------------------------
* -   writes some formatted dummy data, with x and y axis scaling --
* ------------------------------------------------------------------

       REAL c (10,10), xax(10), yax(10)
       INTEGER i, j, k


* ----- define the axes --------------------------------------------

       DO k = 1, 10
          yax(k) = FLOAT(k) * 1000D0
          xax(k) = FLOAT(k) * 500D0
       ENDDO


* ----- create some dummy data -------------------------------------

       DO i = 1, 10
          DO j = 1, 10
             c(i, j) = 1D0 / (1D0 + 
     &       ABS (FLOAT(i - 5) * FLOAT(j - 5)) )
          ENDDO
       ENDDO

* ----- write the data as fomratted output -------------------------

       OPEN(23, FILE = 'test.dat', STATUS = 'UNKNOWN')

* ----- first, write the x-axis ------------------------------------

       WRITE(23, 100) (xax(i), i = 1, 10)

100    FORMAT( 12X, 10(F12.4) )

* ----- now the data; include the y-axis ---------------------------

       DO i = 1, 10
          WRITE(23, 101) yax(i),(c(i,j), j = 1, 10)
       ENDDO

101    FORMAT( 11(F12.4) )

       CLOSE (23)

       END