      SUBROUTINE  POSTEL (N, NSPACE, NELFRE, NRB, NQP, NGEOM,
     1    NPARM, NNPFIX, NNPFLO, MISCFX, MISCFL, NLPFIX,
     2    NLPFLO, COORD, S, C, H, DGH, B, E, EB, STRAIN,
     3    STRAN0, STRESS, BODY, PT, WT, XYZ, DLH, G, DLG,
     4    AJ, AJINV, HINTG, D, PRTLPT, FLTMIS, ELPROP, 
     5    PRTMAT, MISFIX, LTSHAP, LPROP, LPPROP, NTAPE1, 
     6    NTAPE2, NTAPE3, NTAPE4, NTAPE5, IT, NITER, IE, NE, 
     7    LNODE, NG, USEREL, USERPT )
C     * * * * * * * * * * * * * * * * * * * * * * * * * * * *
C           ELEMENT LEVEL POST-SOLUTION CALCULATIONS
C     * * * * * * * * * * * * * * * * * * * * * * * * * * * *
CDP   IMPLICIT   REAL*8  (A-H,O-Z)
C        ALWAYS USED
      DIMENSION  COORD(N,NSPACE), S(NELFRE,NELFRE), C(NELFRE)
C
C        USUALLY USED
C     DIMENSION  H(N,0:NQP), DGH(NSPACE,N,0:NQP), B(NRB,NELFRE), 
      DIMENSION  H(N), DGH(NSPACE,N), B(NRB,NELFRE), 
     1           E(NRB,NRB), EB(NRB,NELFRE), STRAIN(NRB+2), 
     2           STRAN0(NRB), STRESS(NRB+2), BODY(NSPACE)
C
C        OPTIONAL FOR NUMERICAL INTEGRATION
      DIMENSION  PT(NPARM,0:NQP), WT(0:NQP), DLH(NSPACE,N), 
     1           G(NGEOM), DLG(NPARM,NGEOM), AJ(NSPACE,NSPACE), 
     2           AJINV(NSPACE,NSPACE), HINTG(N), LNODE(N),
     3           XYZ(NSPACE)
C
C        OPTIONAL PROPERTY AND SOLUTION VALUES
      DIMENSION  D(NELFRE), PRTLPT(N,0:NNPFLO), FLTMIS(0:MISCFL), 
     1           ELPROP(0:NLPFLO), PRTMAT(0:NLPFLO), 
     2           MISFIX(0:MISCFX), LPROP(0:NLPFIX), 
     3           LPPROP(0:NNPFIX)
C
C        OPTIONAL USER APPLICATION AT NODE OR ELEMENT
      DIMENSION  USERPT(NG), USEREL(NG,N)
C
C                     VARIABLES:
C     AJ      = JACOBIAN
C     AJINV   = JACOBIAN INVERSE
C     B       = STRAIN-DISPLACEMENT (GRADIENT) MATRIX
C     BODY    = BODY FORCE VECTOR
C     C       = ELEMENT COLUMN MATRIX
C     COORD   = SPATIAL COORDINATES OF ELEMENT'S NODES
C     D       = NODAL PARAMETERS ASSOCIATED WITH AN ELEMENT
C     DGH     = GLOBAL DERIVATIVES INTERPOLATION FUNCTIONS
C     DLG     = LOCAL DERIVATIVES GEOMETRIC INTERPOLATION
C     DLH     = LOCAL DERIVATIVES INTERPOLATION FUNCTIONS
C     E       = CONSTITUTIVE MATRIX
C     EB      = PRODUCT OF E*B
C     ELPROP  = ELEMENT ARRAY OF FLOATING PT PROPERTIES
C     FLTMIS  = SYSTEM STORAGE OF FLOATING PT MISC PROP
C     G       = GEOMETRIC INTERPOLATION FUNCTIONS
C     H       = SOLUTION INTERPOLATION FUNCTIONS
C     HINTG   = INTEGRAL OF INTERPOLATION FUNCTIONS
C     IE      = CURRENT ELEMENT NUMBER
C     IT      = CURRENT ITERATION NUMBER
C     LNODE   = ELEMENT TOPOLOGY LIST
C     LPPROP  = INTEGER PROPERTIES AT EACH ELEMENT NODE
C     LPROP   = ARRAY INTEGER POINT ELEMENT PROPERTIES
C     MISFIX  = MISCELLANEOUS INTEGER SYSTEM PROPERTIES
C     N       = NUMBER OF NODES PER ELEMENT
C     NE      = TOTAL NUMBER OF ELEMENTS
C     NELFRE  = NUMBER OF DEGREES OF FREEDOM PER ELEMENT
C     NG      = NUMBER OF GENERALIZED UNKNOWNS PER NODE
C     NGEOM   = NUMBER OF GEOMETRY NODES
C     NITER   = MAX NUMBER OF ITERATIONS
C     NMAT    = NUMBER OF MATERIAL TYPES
C     NPARM   = DIMENSION OF PARAMWETRIC SPACE
C     NQP     = NUMBER OF QUADRATURE POINTS
C     NRB     = NUMBER OF ROWS IN B AND E MATRICES
C     NSPACE  = DIMENSION OF SPACE
C     NTAPE1  = UNIT FOR POST SOLUTION MATRICES STORAGE
C     NTAPE2,3,4 = OPTIONAL UNITS FOR USER (USED WHEN > 0)
C     PRTLPT  = REAL PROPERTIES AT ELEMENT NODES
C     PRTMAT  = REAL ELEM PROPERTIES BASED ON MATERIAL NUMBER
C     PT      = QUADRATURE COORDINATES
C     S       = ELEMENT SQUARE MATRIX
C     STRAIN  = STRAIN OR GRADIENT VECTOR
C     STRAN0  = INITIAL STRAIN OR GRADIENT VECTOR
C     STRESS  = STRESS VECTOR
C     USEREL  = (USER CHOICE) ELEMENT APPLICATION RESULT
C     USERPT  = (USER CHOICE) NODAL APPLICATION RESULT
C     WT      = QUADRATURE WEIGHTS
C     XYZ     = SPACE COORDINATES AT A POINT
C     ....................................................
C      *** POSTEL PROBLEM DEPENDENT STATEMENTS FOLLOW ***
C     ....................................................
      RETURN 
      END
