Seely-Smith prob 188 page 306               
    3    1    1    3    1    0    0    0    0    0    3    2    4    1    1    1
    0    0    0    4    0    0    0    1    0    0    8    0    0    1    0    0
-----
The exact reaction is 2.8522e5, single elem reaction = 2.9325e5
Error =  2.8 %
         1         1    3.86  
         2         0    4.93  
         3         0    6.0   
    1    1    2    3
         1         1     0.0040
         1  30.e6    0.25      2.2e-6   000.

    M   NE   NG    N Spac  Seg  LBN Iter Curv  Rhs ISAY  NRB  NQP Shap LTYP MODE
 Nfix Nflo Lfix Lflo Mfix Mflo Nhom Lhom Nwrt Lwrt Tap1 Tap2  NGF NoCL Sfix Sflo
         1         2         3         4         5         6         7         8
12345678901234567890123456789012345678901234567890123456789012345678901234567890
      (20A4) TITLE, A general description of the problem.
 
      Control parameters. (16I5)
      M      = Number of nodal points in the system.
      NE     = Number of elements in the system.
      NG     = Number of parameters per node.
      N      = Number of nodes per element.
      NSPACE = Dimension of solution space.
      NSEG   = No. of elem boundary segments with given flux.
      LBN    = Number of nodes on an elem. boundary segment.
      NITER  = Number of iterations to be run. (If any.)
      NCURVE = No. of contours per parameter.
      INRHS  = Initial system column matrix flag.
             = 0, Not to be input. (Default)
             > 0, Will be input by CARD TYPE 16.
      ISAY   = Number of user remark cards to be read and
               printed. See CARD TYPE 4.
      NRB    = Number of rows in element B matrix. (Default=1)
      NQP    = Number of quadrature points. (Default=0)
      LSHAPE = Shape 1-line 2-tri 3-quad 4-hex -5 tet
      NLTYPE = Number of element types. (Default=1)
      MODE   = Storage 1-banded 2-skyline    (default=1)
 
      Control parameters. (16I5)
      NNPFIX = Number of integer properties per node.
      NNPFLO = Number of real properties per node.
      NLPFIX = Number of integer proterties per elem.
      NLPFLO = Number of real properties per elem.
      MISCFX = Number of miscellaneous integer properties.
      MISCFL = Number of miscellaneous real properties
      NHOMO  = Nodal properties flag.
            > 0, Properties are homogeneous. Only
                 values for the first node will be input.
            = 0, Non-homogeneous properties. (Default)
      LHOMO  = Element properties flag.
            > 0, Properties are homogeneous. Only
                 values for the first element will be input.
            = 0, Non-homogeneous properties. (Default)
      NPTWRT = Nodal parameter print flag.
            = 0, List calculated values by nodes. (Default)
            Otherwise omit nodal list.
      LEMWRT = Nodal parameter print flag.
            = 0, List node and calculated nodal values for
            each element. Otherwise omit list.
      NTAPE1 = I/O Unit for post solution calculation data.
            = 0, No post solution calculations. (Default)
            > 0, Unit to be utilized.
      NTAPE2 = Extra scratch unit when > 0. (Note: the above
            units are for the analyst's problem dependent data.)
      NGF    = Number of generalized flux components (NSEG>0)
      NULCOL = Element column matrix flag.
            = 0, Problem requires the matrix. (Default)
            > 0, The matrix is always a null vector.
      NBSFIX = Number of segment integer properties.
      NBSFLO = Number of segment real properties.
 
      User comments. Utilized only when ISAY > 0 . (20A4)
 
      Nodal point data. (2I10,(6F10.0))
      J      = System node number.
      IBC(J) = Parameters constraint indicators at
               node J. (Right justified) Digit K equals
               constraint type of parameter K.
      X(J,K) = Spatial coordinates (K=1,NSPACE) of node J
 
      Element Connectivity Data  (I5,(15I5))
      If NLTYPE = 1:
        J          = Element number.
        NODES(J,K) = Nodal incidences (K=1,N) of elem. J
      If NLTYPE > 1:
        J          = Element number.
        LTYPE(J)   = Element type number
        NODES(J,K) = Nodal incidences (K=1,N) of elem. J
 
      Essential boundary conditions. (2I10,F10.0)
      NODE1 = system nodal point number.
      IPAR1 = Constrainted parameter number, 1.le.IPAR1.le.NG
      A1    = Value assigned to the parameter.
 
      Nodal properties (I10,(7I10))
      J          = Nodal point number.
      NPFIX(J,K) = Fixed point properties. K=1,NNPFIX
 
      Nodal properties (I10,(7F10.0))
      J           = Nodal point number.
      FLTNP(J,K) = Floating point properties. K=1,NNPFLO
 
      Element properties (I5,(7I10))
      J          = Element number.
      LPFIX(J,K) = Fixed point properties. K=1,NLPFIX
 
      Element properties (I5,(7F10.0))
      J          = Element number.
      FLTEL(J,K) = Floating point properties. K=1,NLPFLO
 
      Segment properties (I5,(7I10))
      J           = Segment number.
      NBSPFX(J,K) = Fixed point properties. K=1,NBSFIX
 
      Segment properties (I5,(7F10.0))
      J          = Segment number.
      FLTBS(J,K) = Floating point properties. K=1,NBSFLO
 
      Integer properties. (8I10.0)
      MISFIX(K) = Misc. fixed point properties.K=1,MISCFX
 
      Real properties. (8F10.0)
      FLTMIS(K) = Misc. floating point properties.K=1,MISCFL
 
      External point sources (2I10,F10.0)
      NODE  = System node number of the point.
      IPARM = Parameter number associated with the value.
      VALUE = Specified value to be added to column matrix.
 
      Segment nodes (16I5)
      LNODE(K) = List of nodes on the segment. K=1,LBN.
      Flux component values (8F10.4)
      (FLUX(K,L),L=1,NGF) = The NGF components of flux at 
      each of the above nodes. K=1,LBN.
 
      ***** END OF DATA INSTRUCTIONS *****
