PROGRAM INPDAT
C
C  THIS PROGRAM READS N POINTS FROM A DATA FILE AND STORES THEM IN 
C  3 ARRAYS X, Y, Z.
C
      INTEGER NMAX, U
      PARAMETER (NMAX=1000, U=20)
      REAL X(NMAX), Y(NMAX), Z(NMAX)

C  OPEN THE DATA FILE
      OPEN (U, FILE='POINTS.DAT', STATUS='OLD')

C  READ THE NUMBER OF POINTS
      READ(U,*) N
      IF (N.GT.NMAX) THEN
         WRITE(*,*) 'ERROR: N = ', N, 'IS LARGER THAN NMAX =', NMAX
         GOTO 9999
      ENDIF

C  LOOP OVER THE DATA POINTS
      DO 10 I= 1, N
         READ(U,*) X(I), Y(I), Z(I)
   10 CONTINUE
  100 FORMAT (3(F10.4))

C  CLOSE THE FILE
      CLOSE (U)

C  NOW WE SHOULD PROCESS THE DATA SOMEHOW...
C  (MISSING PART)

 9999 STOP
      END