* -------------------------------------------------------------------- *
      PROGRAM LATTICE
*
C     Developes a lattice from input of:
C     Xmin, Ymin, Xinc, Yinc
C     The output is in ascending order of X,Y-pairs
*
      REAL XMIN, YMIN, XINC, YINC, XVAL, YVAL
      INTEGER ROWS, COLS, X, Y, ID, DIMENS
*
      OPEN (UNIT=99, FILE='coord.dat', STATUS='NEW')
*
      PRINT 1
    1 FORMAT (1X,' Enter value (m) for X-min (lower left): ')
      READ*, XMIN
      PRINT 2
    2 FORMAT ('+',' Enter value (m) for Y-min (lower left): ')
      READ*, YMIN
      PRINT 3
    3 FORMAT ('+',' Enter value (m) for X-Increment:        ')
      READ*, XINC
      PRINT 4
    4 FORMAT ('+',' Enter value (m) for Y-Increment:        ')
      READ*, YINC
      PRINT 5
    5 FORMAT ('+',' Enter number of (vertical) columns:     ')
      READ*, COLS
      PRINT 6
    6 FORMAT ('+',' Enter number of (horizontal) rows:      ')
      READ*, ROWS
      XVAL = XMIN
      ID   = 1
      DO 50 X = 1,COLS
         YVAL = YMIN
         DO 40 Y = 1,ROWS
            WRITE (99,20) ID, XVAL, YVAL
   20       FORMAT (1X,I7,2(',',2X,F11.1))
            ID = ID + 1
            YVAL = YVAL + YINC
   40    CONTINUE
         XVAL = XVAL + XINC
   50 CONTINUE
      WRITE(99,60)
   60 FORMAT (1X)
      DIMENS = REAL(ROWS) * REAL(COLS)
      PRINT 70, DIMENS
   70 FORMAT ('0','ÚÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ¿',/,
     +         1X,'³ Output has ',    I11,' pairs of coordinates ³',/,
     +         1X,'³ ÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄ ³',/,
     +         1X,'³        Output is in file COORD.DAT          ³',/,
     +         1X,'ÀÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÄÙ')
      END
* -------------------------------------------------------------------- *