NOARCLDA for no arc length damping
    NREST for number of time steps between restart file generation
    NSMD for nonlinear solution method which can be
      BFGS for BFGS (default)
      BROY for Broyden
      DFP for Davidon-Fletcher-Powell
      MDAV for modified Davidon
      MNEWT for modified Newton
      MARC for modified constant arc length
      MCLS for modified constant arc length with line search
      MABFGS for modified constant arc length with BFGS
      MABROY for modified constant arc length with Broyden
      MADFP for modified constant arc length with DFP
      MAMBFGS for modified constant arc length with modified BFGS
      MADAV for modified constant arc length with Davidon
      NEWT for full Newton
      NEWTLS for full Newton with line search
    NSTEPS for number of time steps
    NUNLOAD for number of unloading steps in modified arc length method
    OPNIT for optimal number of iterations per step
    PRLIS for linear iterative solver print-out flag with values
      0 for no inner loop information
      1 for time step convergence information
      2 for iteration norm information
      3 for residual input, solution output
      4 for residual/solution each iteration
    RETOL to set the convergence tolerance on residual norm
    SEGS to set the shell element geometric stiffness flag
    SFOR for shell element formulation which can be
      HL for Hughes-Liu
      YASE for YASE
      BT for Belytschko-Tsay (DYNA3D-like)
    SHIFT for frequency shift
    SSDM for shell surface strain data dumps
    STIFCORE for stiffness matrix storage option
    SW3 for toggle the default sense switch number 3
    SW6 for toggle the default sense switch number 6
    SW7 for toggle the default sense switch number 7
    TEO for thermal effects option
    XVEL for load due to x-angular velocity
    YVEL for load due to y-angular velocity
    ZVEL for load due to z-angular velocity
  LSNKOPTS: LSNIKE3D analysis options.
    ANAL with the options
      STAT for static
      DYN for dynamic analysis
      DYNS for dynamic analysis with stresses initialized statically
    ARCL for desired arc length
    ARCLCM for the arc length constraint method which can be
      CRISFIELD
      RAMM
    ARCLDAMP for arc length damping
    AUTO for automatic time step control invoked
    BRSTIF for brick element geometric stiffness included
    BWMO for bandwidth minimization on/off switch
    DCTOL for displacement convergence tolerance
    DELT for time step
    DISPNODE for node number for displacement controlled arc length method
    DISPDIR for direction of displacement at node for arc length control
      1 is for the global x-direction
      2 is for the global y-direction
      3 is for the global z-direction
    DSFCT for force vs. displacement table displacement scale factor
    DTROC for type of force recovery resulting from enforced boundary conditions
    ECTOL for energy convergence tolerance
    FRQSHT for frequency shift, cycles per unit time
    FSFCT for force vs. displacement table force scale factor
    ICNT1 for initial stiffness matrix formulation frequency
    ICNT2 for equilibrium iteration frequency
    ICNVRG for displacement criteria to check convergence
    IDYNRST for number of materials to initialize stresses in LSDYNA
    IFB for output force balance table
    IFDOF for degree-of-freedom for node in force vs. displacement table
    IFFD for output force vs. displacement table
    IFPRI for force dump print interval
    INFNM1 for interface number to end first phase
    INFNM2 for interface number to end second phase
    INITSTR for read initial element stresses
    IOBUF for buffer size (words) for element data I/O
    IORO for Quasi-Newton update vectors storage
    IPLT for TUARUS dump interval
    IPRT for print dump interval
      positive number causes printing each interval excluding step 0 0 causes
      printing each step including step 0 negative number causes printing each
      interval including step 0
    ISLVDT for linear equation solver option
    LAUTO for parameter to control of increase in the time step
    LCLOD1 for load curve for first phase
    LCLOD2 for load curve for second phase
    LSTOL for line search convergence tolerance
    MNSS for minimum allowable step size
    MSRF for maximum number of reform/time steps
    MUNLOAD for arc length unloading method which can be
      BFGS for BFGS
      BROY for Broyden
      DFP for Davidon-Fletcher-Powell
      DAV for Davidon
      MNEWT for modified Newton
    MEMSOL for percentage of memory for the solver
    MXNRE for maximum number of retries allowable per step
    MXSS for maximum allowable step size
    NBEI for number of steps between equilibrium iterations
    NBSR for number of steps between matrix reformations
    NCYCSP for number of cycles for monotonic pressure after reversal
    NDPT for node number for force vs. displacement table
    NECHO for nodal input data echo flag
    NEIG for number of eigenvectors
    NGRAV grav_x load_curve grav_y load_curve grav_z load_curve
    NIBSR for maximum number of equilibrium iter./matrix reform.
    NIP1 for first Newmark integration parameter
    NIP2 for second Newmark integration parameter
    NPSD for discrete element input flag
    NREST for number of time steps between restart file generation
    NSMD for nonlinear solution method which can be
      BFGS for BFGS (default)
      BROY for Broyden
      DFP for Davidon-Fletcher-Powell
      MDAV for modified Davidon
      MNEWT for modified Newton
      MARC for modified constant arc length
      MCLS for modified constant arc length with line search
      MABFGS for modified constant arc length with BFGS
      MABROY for modified constant arc length with Broyden
      MADFP for modified constant arc length with DFP
      MAMBFGS for modified constant arc length with modified BFGS
      MADAV for modified constant arc length with Davidon
      NEWT for full Newton
      NEWTLS for full Newton with line search
    NSTEPS for number of time steps
    NUNLOAD for number of unloading steps in modified arc length method
    NUSRLDP for number of user loading subroutine parameters
    N1DECHO for beam element input data echo flag
    N2DECHO for shell element input data echo flag
    N3DECHO for solid element input data echo flag
    OPNIT for optimal number of iterations per step
    PENLDC for load curve defining multiple vs. time
    PENSTAB for penalty stiffness multiplier for artificial stabilization
    SCLMIN for minimum pressure load scale factor
    SCLMAX for maximum pressure load scale factor
    SPRCNT1 for percentage of nodes in contact to end first phase
    SPRCNT2 for percentage of nodes in contact to end second phase
    SSDM for shell surface strain data dumps
    STIFCORE for stiffness matrix storage option
    TARGET for target strain rate
    TEO for thermal effects option
    XVEL for load due to x-angular velocity
    YVEL for load due to y-angular velocity
    ZVEL for load due to z-angular velocity
  TZ3DOPTS: TOPAZ3D analysis options.                                           T
    IUNIT for temperature units with the two exclusive options (enclosure 
     radiation data only)
      Kelvin
      Rankine
    BWMO for bandwidth and profile minimization selection with
     the three exclusive options
      no minimization (default)
      minimization
      minimization, nodal destination vector is read from input file
    PHASE to perform phase change calculations
    SOLVER to select the analysis solver
      FISSLE out of core profile solver
      ACTCOL in core profile solver
      DACTCOL nonsymmetric profile solver
      DSCG diagonally scaled conjugate gradient
      ICCG incomplete cholesky conjugate gradient
      Nonsymmetric diagonally scaled conjugate gradient
      Diagonally scaled bi-conjugate gradient
      Incomplete LU bi-conjugate gradient 
    RADIATION to select the radiation calculation type
     with the five exclusive options
      VIEW for view factors
      FBEXCHNG for familied binary exchange factors
      FAEXCHNG for familied ascii exchange factors
      BBEXCHNG for black body exchange factors
      NFXCHNG for non-familied ascii exchange factors 
    select type of analysis from the two exclusive options
      STEADY for steady state or
      TRANSIENT with the options
        FIXED for fixed time step with the options
          TIMIN for the initial problem time
          TIMEND for the final problem time
          DELT for the time step size
        VARIABLE for variable time step with the options
          TIMIN for the initial problem time
          TIMEND for the final problem time
          DELT for the initial time step size
          DTMIN for the minimum time step
          DTMAX for the maximum time step
          DCMX for the desired maximum temperature change in each time step
           above which the time step will be decreased
          MFTS for the modification factor for increasing/decreasing time step
        IPRT for the number of time steps between printed data output
        IPLT for the number of steps between plotted data output
        SBRF for the number of time steps between dump file generation
        TIPAR for the time integration parameter
    select type of problem with the two exclusive options
      LINEAR for a linear problem or
      NONLINEAR for a nonlinear problem with the options
        NBSR for the number of time steps between conductance matrix
         reformations
        NBEI for the number of time steps between equilibrium iterations
        MSRF for the maximum number of conductance matrix reformations per time
         step
        NIBSR for the maximum number of equilibrium iterations permitted per
         conductance matrix reformation
        DCTOL for the convergence tolerance for equilibrium iterations
        TZRELAX for the relaxation parameter
    SBC for the Stefan-Boltzman constant
    RCTOL for the radiosity convergence tolerance
    MRDI for the maximum number of radiosity iterations
    RBAND to specify a set of wavelength breakpoints
    EMISSIVITY to specify multiple numbered 2D emissivity curves with the same
     number of points as the number of wavelength breakpoints
ELEMENT commands:
  BIND: Hughes-Liu beam user defined integration points.
  BSD: global beam element cross section definition.
   A cross section is referenced by every beam.
   Each beam is defined by two end nodes and either a third node or
   a coordinate triple to define the orientation of the cross-sectional
   local coordinate system of the beam.
   Throughout TG, the local coordinate system for the beam element
   is referred to by x, y, and z.
   They correspond to,respectively, the longitudinal axis
   of the beam stretched between the first two nodes, the orthogonal
   component from the first node to the third node/coordinate, and
   the third direction orthogonal to the first two directions.
   The thickness parameters STHI, STHI1, STHI2, TTHI, TTHI1, and TTHI2
   are the distances from the center to the boundary of the element.
   They are the half thicknesses used by DYNA3D and NIKE3D.
   Other parameters are used to define specific dimensions for the
   different types of cross sections.
   When there is a possibility of specifying a quantity differently
   at each end, there will be one option which assigns the same value
   to both ends, and then two additional options, one for each end.
   These two options are distinguished by the suffix of 1 and 2, respectively.
   When multiple beams are generated from information at two ends, some
   quantities such as thickness are linearly interpolated.
   Other quantities, such as cross-sectional area and moments are not
   interpolated.
   In this latter case, the values found at the first cross section
   are broadcast throughout the string of beams.
   Some of the parameters in the cross section definition are also definable
   within the definition of the individual beams.
   The value of these parameters within the beam definition take
   precedence over this beam cross section definition.
   The parameters in this cross section definition are code dependent.
   The parameters available for each code are listed.
   The STHI, STHI1, STHI2, TTHI, TTHI1, TTHI2, AREA, SCAREA, IXX, IYY,
   IZZ and LDP options can be
   used any time.
    ABAQUS: the *BEAM SECTION card is generated using CSTYPE with
      7 which defines a Pipe cross section,
      8 which defines a Box cross section,
      9 which defines a Circle cross section,
      10 which defines a I-Beam cross section,
      11 which defines a Rectangle cross section,
      12 which defines a Hexagon cross section,
      13 which defines a Elbow cross section,
      14 which defines a Trapezoid cross section,
      15 which defines a L-Section cross section,
      16 which defines a Arbitrary cross section.
    ANSYS: 
      BAN4 defines the Elastic Beam element 4,
      BAN8 defines the Spar element 8,
      BAN10 defines the Tension/Compression Spar element 10,
      BAN24 defines the Thin Walled Plastic element 24,
      BAN33 defines the Thermal Bar element 33,
      BAN44 defines the Tapered Unsymmetrical element 44.
    DYNA3D: default cross section properties are defined within the
     material definition. To override the default material definition
     use STHI, STHI1, STHI2, TTHI, TTHI1, TTHI2, AREA, SCAREA, IXX,
     IYY, IZZ;
    LSDYNA: the following cross section options are available
     for LSDYNA beam elements:
      Thicknesses
      Area and Inertia
      Discrete Beams
      W-section
      C-section
      angle section
      T-section
      rectangular tubing
      Z-section
      trapezoidal section
    MARC: the following cross section options are available for MARC elements: