FILE FORMAT FOR GEOMETRY FILE

          READ (IFILE,ERR=63,END=63)

     1      N,M,((CORDS(J,K),K=1,2), ALFA(J),AO(J),J=1,N),
     2      ((NOPS(J,K),K=1,8),IMAT(J),TH(J),NFIXHS(J),J=1,M),
     +      (WIDTH(J),SS1(J),SS2(J),WIDS(J),J=1,N)

c
c     Part 1 of the geo file contains
c
c     1   N          Number of nodes on geo file                     INTEGER*4
c     2   M          Number of elements on "geo" file                INTEGER*4
c     3   CORD(J,1)  x-coordinate - node J                           REAL*8   ARRAY SIZE (NUMBER OF NODES, 2)
C     4   CORD(J,2)  y-coordinate - node J                           REAL*8   ARRAY SIZE (NUMBER OF NODES)
C     5   ALFA(J)    Boundary slope - node j                         REAL*4   ARRAY SIZE (NUMBER OF NODES)
C     6   AO(J)      Bed elevation - node j                          REAL*4   ARRAY SIZE (NUMBER OF NODES)
C     7   NOP(J,K)   Nodal connections, counterclockwise - element j INTEGER*8  ARRAY SIZE (NUMBER OF ELEMENTS,8)
C     8   IMAT(J)    Type number - element j                         INTEGER*2  ARRAY SIZE (NUMBER OF ELEMENTS)
C     9   TH(J)      Principal direction for eddy viscosity - element j  REAL*4  ARRAY SIZE (NUMBER OF ELEMENTS)
C    10   NFIXHS(J)  Ordering sequence for equation formation and solution - counter j  INTEGER*4  ARRAY SIZE (NUMBER OF ELEMENTS)
C    11   WIDTH(J)   Nodal bed width for 1-D elements                REAL*4   ARRAY SIZE (NUMBER OF NODES)
C    12   SS1(J)     Nodal left side slope for 1-D elements          REAL*4   ARRAY SIZE (NUMBER OF NODES)
C    13   SS2(J)     Nodal right side slope for 1-D elements         REAL*4   ARRAY SIZE (NUMBER OF NODES)
C    14   WIDS(J)    Nodal storage width for 1-D elements            REAL*4   ARRAY SIZE (NUMBER OF NODES)

The remainder of the file contains more information associated with 1-D elements 


C    11   WIDTH(J)   Nodal bed width ( for 1-D nodes) - node j
C    12   SS1(J)     Nodal left side slope ( for 1-D nodes) - node j
C    13   SS2(J)     Nodal right side slope ( for 1-D nodes) - node j
C    14   WIDS(J)    Off channel storage width - node j (see appendix to documentation)
