VortexX Nozzle Thermal Stress  
    7    3    1    3    1    0    0    0    0    0    5    2    4    1    1    1
    0    0    0    4    0    0    0    0    0    0    8    0    0    0    0    0
                         ---
Nozzle, Brass, Tungsten, Brass thickness = 0.001                                
Radial max min: 57,000 0 : 53,000 53,000 : 15,000 -4,000                        
Hoop   max min: 108,000 66,000 : 228,000 228,000 : -46,000 -7,000               
(Reactions should be zero, possible model bug)
         1         0    0.2   
         2         0    0.4125 
         3         0    0.625 
         4         0    0.6255
         5         0    0.626
         6         0    1.627
         7         0    2.627
    1    1    2    3
    2    3    4    5
    3    5    6    7
         1  90.e6    0.24      2.2e-6   -1400.
         2  15.e6    0.10      11.0-6   -1400. 
         3 100.e6    0.24      4.4e-6   -1400.

    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 *****
