>       PARAMETER  ( NREACT = 20, CUTOFF = 1.E-6 )
>       LOGICAL    FACT, BACK
>       DIMENSION  TITLE(15), RARRAY(MAXR), RNAME(NUMR),         
>      1  X(M,NSPACE), AJ(NSPACE,NSPACE), AJINV(NSPACE,NSPACE),
>      2  AVE(0:M,NRB+2), B(NRB,NELFRE), BODY(NSPACE), C(NELFRE),
>      3  CC(NDFREE), CEQ(MAXACT,NUMCE), COORD(N,NSPACE), D(NELFRE),
>      4  DDOLD(NDFREE), DGH(NSPACE,N,0:NQP), DLG(NPARM,NGEOM,0:NQP),
>      5  DLH(NSPACE,N,0:NQP), E(NRB,NRB), EB(NRB,NELFRE), 
>      6  ELPROP(0:NLPFLO), FLTEL(NE,0:NLPFLO), FLTBS(0:NSEG,0:NBSFLO), 
>      7  FLTMIS(0:MISCFL), FLTNP(M,0:NNPFLO), FLUX(0:NGF), 
>      8  FLUXBS(0:NSEG,0:NFLUX), G(NGEOM,0:NQP), GPT(0:NQP), GWT(0:NQP)
>       DIMENSION  H(N,0:NQP), HINTG(N), PLTSET(0:NPLT), 
>      1  PRTLPT(N,0:NNPFLO), PRTMAT(0:NLPFLO,0:NOMAT), PT(NPARM,NQP), 
>      2  RANGE(NG,2), S(NELFRE,NELFRE), SATPT(NRB+2,N), 
>      3  STRAIN(NRB+2),STRAN0(NRB), STRESS(NRB+2), SYSDAT(0:NSYS), 
>      4  TMP(0:NTMP), VALC(NRB,0:NC), VALE(NRB,0:NC), WT(NQP), 
>      5  XPT(NSPACE,2), XYZ(NSPACE), DD(NDFREE), USEREL(NG,N),
>      6  USERPT(NG), SS(1)
>       DIMENSION  IARRAY(MAXI), INAME(NUMI), IBC(M), KODES(NG), 
>      1  NODES(NE,N), NRES(MAXTYP), LTYPE(NE), IADD(0:M), 
>      2  IDIAG(NDFREE), INDEX(NELFRE), LFIRST(M), LLAST(M), LNODE(N), 
>      3  LPFIX(NE,0:NLPFIX), NBSPFX(0:NSEG,0:NBSFIX), LPPROP(0:NNPFIX), 
>      4  LPROP(0:NLPFIX), LTDATA(6,NLTYPE), MISFIX(0:MISCFX), 
>      5  NDXC(MAXACT,NUMCE), NODEF(0:NSEG,0:LBN), NRANGE(NG,2), 
>      6  NREQ(MAXTYP), NPFIX(M,0:NNPFIX), LHIGH(NELFRE)
