C   ---  Calculate real array pointers  ---
        J(1)  =  1
C   REAL    1, sys pt coords  
       RN(  1)  =  "X         "
        J(  2)  =  J(  1) + (M         )*(NSPACE    )*(1         )+0
C   REAL    2, jacobian       
       RN(  2)  =  "AJ        "
        J(  3)  =  J(  2) + (NSPACE    )*(NSPACE    )*(1         )+0
C   REAL    3, inverse jacobia
       RN(  3)  =  "AJINV     "
        J(  4)  =  J(  3) + (NSPACE    )*(NSPACE    )*(1         )+0
C   REAL    4, real pt average
       RN(  4)  =  "AVE       "
        J(  5)  =  J(  4) + (M     +1  )*(NSPACE+1  )*(1         )+1
C   REAL    5, b matrix       
       RN(  5)  =  "B         "
        J(  6)  =  J(  5) + (NRB       )*(NELFRE    )*(1         )+0
C   REAL    6, body force     
       RN(  6)  =  "BODY      "
        J(  7)  =  J(  6) + (NSPACE    )*(1         )*(1         )+0
C   REAL    7, el load vector 
       RN(  7)  =  "C         "
        J(  8)  =  J(  7) + (NELFRE    )*(1         )*(1         )+0
C   REAL    8, sys load vector
       RN(  8)  =  "CC        "
        J(  9)  =  J(  8) + (NDFREE    )*(1         )*(1         )+0
C   REAL    9, constrain coeff
       RN(  9)  =  "CEQ       "
        J( 10)  =  J(  9) + (MAXACT    )*(NUMCE     )*(1         )+0
C   REAL   10, el or pt coord 
       RN( 10)  =  "COORD     "
        J( 11)  =  J( 10) + (N         )*(NSPACE    )*(1         )+0
C   REAL   11, el or pt dof   
       RN( 11)  =  "D         "
        J( 12)  =  J( 11) + (NELFRE    )*(1         )*(1         )+0
C   REAL   12, old sys dof    
       RN( 12)  =  "DDOLD     "
        J( 13)  =  J( 12) + (NDFREE    )*(1         )*(1         )+0
C   REAL   13, global deriv h 
       RN( 13)  =  "DGH       "
        J( 14)  =  J( 13) + (NSPACE    )*(N         )*(NQP   +1  )+1
C   REAL   14, local deriv g  
       RN( 14)  =  "DLG       "
        J( 15)  =  J( 14) + (NPARM     )*(NGEOM     )*(NQP   +1  )+1
C   REAL   15, local deriv h  
       RN( 15)  =  "DLH       "
        J( 16)  =  J( 15) + (NSPACE    )*(N         )*(NQP   +1  )+1
C   REAL   16, constitutive ma
       RN( 16)  =  "E         "
        J( 17)  =  J( 16) + (NRB       )*(NRB       )*(1         )+0
C   REAL   17, matrix product 
       RN( 17)  =  "EB        "
        J( 18)  =  J( 17) + (NRB       )*(NELFRE    )*(1         )+0
C   REAL   18, real el propert
       RN( 18)  =  "ELPROP    "
        J( 19)  =  J( 18) + (NLPFLO+1  )*(1         )*(1         )+1
C   REAL   19, real prop all e
       RN( 19)  =  "FLTEL     "
        J( 20)  =  J( 19) + (NE        )*(NLPFLO+1  )*(1         )+1
C   REAL   20, real prop segme
       RN( 20)  =  "FLTBS     "
        J( 21)  =  J( 20) + (NSEG  +1  )*(NBSFLO+1  )*(1         )+1
C   REAL   21, real misc prop 
       RN( 21)  =  "FLTMIS    "
        J( 22)  =  J( 21) + (MISCFL+1  )*(1         )*(1         )+1
C   REAL   22, real prop all p
       RN( 22)  =  "FLTNP     "
        J( 23)  =  J( 22) + (M         )*(NNPFLO+1  )*(1         )+1
C   REAL   23, flux comps on e
       RN( 23)  =  "FLUX      "
        J( 24)  =  J( 23) + (NF    +1  )*(1         )*(1         )+1
C   REAL   24, all flux on nod
       RN( 24)  =  "FLUXBS    "
        J( 25)  =  J( 24) + (NSEG  +1  )*(NFLUX +1  )*(1         )+1
C   REAL   25, geom interpolat
       RN( 25)  =  "G         "
        J( 26)  =  J( 25) + (NGEOM     )*(NQP   +1  )*(1         )+1
C   REAL   26, gauss pts 1-d  
       RN( 26)  =  "GPT       "
        J( 27)  =  J( 26) + (NQP   +1  )*(1         )*(1         )+1
C   REAL   27, gauss wts 1-d  
       RN( 27)  =  "GWT       "
        J( 28)  =  J( 27) + (NQP   +1  )*(1         )*(1         )+1
C   REAL   28, solution interp
       RN( 28)  =  "H         "
        J( 29)  =  J( 28) + (N         )*(NQP   +1  )*(1         )+1
C   REAL   29, integral of h  
       RN( 29)  =  "HINTG     "
        J( 30)  =  J( 29) + (N         )*(1         )*(1         )+0
C   REAL   30, plotter data   
       RN( 30)  =  "PLTSET    "
        J( 31)  =  J( 30) + (NPLT  +1  )*(1         )*(1         )+1
C   REAL   31, real prop el no
       RN( 31)  =  "PRTLPT    "
        J( 32)  =  J( 31) + (N         )*(NNPFLO+1  )*(1         )+1
C   REAL   32, real mat numb p
       RN( 32)  =  "PRTMAT    "
        J( 33)  =  J( 32) + (NLPFLO+1  )*(NOMAT +1  )*(1         )+1
C   REAL   33, quadrature coor
       RN( 33)  =  "PT        "
        J( 34)  =  J( 33) + (NPARM     )*(NQP       )*(1         )+0
C   REAL   34, dof max min val
       RN( 34)  =  "RANGE     "
        J( 35)  =  J( 34) + (NG        )*(2         )*(1         )+0
C   REAL   35, el or edge sq m
       RN( 35)  =  "S         "
        J( 36)  =  J( 35) + (NELFRE    )*(NELFRE    )*(1         )+0
C   REAL   36, stress at el po
       RN( 36)  =  "SATPT     "
        J( 37)  =  J( 36) + (NRB+2     )*(N         )*(1         )+0
C   REAL   37, strain or grad 
       RN( 37)  =  "STRAIN    "
        J( 38)  =  J( 37) + (NRB+2     )*(1         )*(1         )+0
C   REAL   38, initial strain 
       RN( 38)  =  "STRAN0    "
        J( 39)  =  J( 38) + (NRB       )*(1         )*(1         )+0
C   REAL   39, stress + rms or
       RN( 39)  =  "STRESS    "
        J( 40)  =  J( 39) + (NRB+2     )*(1         )*(1         )+0
C   REAL   40, sys control dat
       RN( 40)  =  "SYSDAT    "
        J( 41)  =  J( 40) + (NSYS  +1  )*(1         )*(1         )+1
C   REAL   41, temporary work 
       RN( 41)  =  "TMP       "
        J( 42)  =  J( 41) + (NTMP  +1  )*(1         )*(1         )+1
C   REAL   42, values at corne
       RN( 42)  =  "VALC      "
        J( 43)  =  J( 42) + (NRB       )*(NC    +1  )*(1         )+1
C   REAL   43, values on edge 
       RN( 43)  =  "VALE      "
        J( 44)  =  J( 43) + (NRB       )*(NC    +1  )*(1         )+1
C   REAL   44, quadrature weig
       RN( 44)  =  "WT        "
        J( 45)  =  J( 44) + (NQP       )*(1         )*(1         )+0
C   REAL   45, xy ends of a li
       RN( 45)  =  "XPT       "
        J( 46)  =  J( 45) + (NSPACE    )*(2         )*(1         )+0
C   REAL   46, global of gauss
       RN( 46)  =  "XYZ       "
        J( 47)  =  J( 46) + (NSPACE    )*(1         )*(1         )+0
C   REAL   47, sys square matr
       RN( 47)  =  "SS        "
        J( 48)  =  J( 47) + (NCOEFF    )*(1         )*(1         )+0
