1,10c1,10
<       SUBROUTINE  ASYMBL (NG, NCOEFF, MODE, IDIAG, NODES, SS, CC, 
<      1    M, NE, NDFREE, NITER, LPTEST, LHOMO, NHOMO, NULCOL, N, 
<      2    NSPACE, NELFRE, NRB, NQP, NGEOM, NPARM, NNPFIX, NNPFLO, 
<      3    MISCFX, MISCFL, NLPFIX, NLPFLO, LNODE, INDEX, X, DDOLD, 
<      4    COORD, S, C, H, DGH, B, E, EB, STRAIN, STRAN0, STRESS, 
<      5    BODY, PT, WT, XYZ, DLH, G, DLG, AJ, AJINV, HINTG, D, 
<      6    PRTLPT, FLTNP, FLTEL, FLTMIS, ELPROP, PRTLPT, PRTMAT, 
<      7    MISFIX, NPFIX, LPFIX, LPROP, LPPROP, NTAPE1, NTAPE2, 
<      8    NTAPE3, NTAPE4, NTAPE5, LTYPE, NLTYPE, LTDATA, LSHAPE, 
<      9    GPT, GWT ) 
---
>       SUBROUTINE  ASYMBL ( NG, NCOEFF, MODE, IDIAG, NODES, SS, 
>      1    CC, M, NE, NDFREE, NITER, LPTEST, LHOMO, 
>      2    NHOMO, NULCOL, N, NSPACE, NELFRE, NRB, NQP, NGEOM,
>      3    NPARM, NNPFIX, NNPFLO, MISCFX, MISCFL, NLPFIX, NLPFLO, 
>      4    LNODE, INDEX, X, DDOLD, COORD, S, C, H, DGH, B, E, EB, 
>      5    STRAIN, STRAN0, STRESS, BODY, PT, WT, XYZ, DLH, G, DLG,
>      6    AJ, AJINV, HINTG, D, PRTLPT, FLTNP, FLTEL, FLTMIS,
>      7    ELPROP, PRTLPT, PRTMAT, MISFIX,         NPFIX, LPFIX, 
>      8    LPROP, LPPROP, NTAPE1, NTAPE2, NTAPE3, NTAPE4, NTAPE5,
>      9    LTYPE, NLTYPE, LTDATA, LSHAPE, gpt, gwt ) 
19c19,20
< C       JUST PASSING THROUGH: SYSTEM DATA
---
> C     JUST PASSING THROUGH
> C         SYSTEM DATA
25,37c26,39
< C      FOR USE IN ELSQ, ELCOL, OR ELPOST:
<       DIMENSION  COORD(N,NSPACE), S(NELFRE,NELFRE), C(NELFRE), 
<      1           H(N,0:NQP), DGH(NSPACE,N), B(NRB,NELFRE), 
<      2           E(NRB,NRB), EB(NRB,NELFRE), STRAIN(NRB+2), 
<      3           STRAN0(NRB), STRESS(NRB+2), BODY(NSPACE),
<      4           PT(NPARM,0:NQP), WT(0:NQP), DLH(NSPACE,N,0:NQP), 
<      5           G(NGEOM,0:NQP), DLG(NPARM,NGEOM,0:NQP), 
<      6           AJ(NSPACE,NSPACE), AJINV(NSPACE,NSPACE), 
<      7           HINTG(N,0:NQP+1), GPT(0:NQP), GWT(0:NQP), 
<      8           XYZ(NSPACE), D(NELFRE), FLTMIS(0:MISCFL), 
<      9           ELPROP(0:NLPFLO), PRTMAT(0:NLPFLO), 
<      1           MISFIX(0:MISCFX), LPROP(0:NLPFIX), 
<      2           LPPROP(0:NNPFIX)
---
> C        FOR USE IN ELSQ OR ELCOL:
>       DIMENSION  COORD(N,NSPACE), S(NELFRE,NELFRE)
>       DIMENSION  C(NELFRE), H(N,0:NQP), DGH(NSPACE,N), B(NRB,NELFRE), 
>      1           E(NRB,NRB), EB(NRB,NELFRE), STRAIN(NRB+2), 
>      2           STRAN0(NRB), STRESS(NRB+2), BODY(NSPACE)
>       DIMENSION  PT(NPARM,0:NQP), WT(0:NQP), DLH(NSPACE,N,0:NQP), 
>      1           G(NGEOM,0:NQP), DLG(NPARM,NGEOM,0:NQP), 
>      2           AJ(NSPACE,NSPACE), AJINV(NSPACE,NSPACE), 
>      3           HINTG(N,0:NQP+1), gpt(0:nqp), gwt(0:nqp), 
>      4           xyz(nspace)
>       DIMENSION  D(NELFRE), FLTMIS(0:MISCFL), ELPROP(0:NLPFLO), 
>      1           PRTMAT(0:NLPFLO), MISFIX(0:MISCFX), 
>      2                             LPROP(0:NLPFIX), 
>      3           LPPROP(0:NNPFIX)
68c70
< C     LTQP    = NUMBER OF QUADRATURE PTS FOR ELEMENT TYPE
---
> C     LTQP    = NUMBER OF QUADRATURE POINTS FOR ELEM TYPE
117,119c119,120
<           IF ( LTQP .GT. 0 )  THEN
<             IF ( LTQP .GT. NQP ) STOP 'LTQP > NQP IN ASYMBL'
<             CALL  GETQD (LTSHAP, LTQP, NSPACE, GPT, GWT, PT, WT)
---
>           IF ( LTQP .GT. 0 ) CALL  GETQD (LTSHAP, LTQP, NSPACE, 
>      1                                    GPT, GWT, PT, WT)
121d121
<         ENDIF
128,136c128,135
<      1        NHOMO, NULCOL, LTN, NSPACE, LTFREE, NRB, LTQP, 
<      2        LTGEOM, LTPARM, NNPFIX, NNPFLO, MISCFX, MISCFL, 
<      3        NLPFIX, NLPFLO, LNODE, INDEX, X, DDOLD, COORD, S, 
<      4        C, H, DGH, B, E, EB, STRAIN, STRAN0, STRESS, BODY, 
<      5        PT, WT, XYZ, DLH, G, DLG, AJ, AJINV, HINTG, D, 
<      6        PRTLPT, FLTNP, FLTEL, FLTMIS, ELPROP, PRTLPT, 
<      7        PRTMAT, MISFIX, NPFIX, LPFIX, LPROP, LPPROP, 
<      8        NTAPE1, NTAPE2, NTAPE3, NTAPE4, NTAPE5, LT, 
<      9        LTSHAP, LTUSER, NG )
---
>      1    NHOMO, NULCOL, LTN, NSPACE, LTFREE, NRB, LTQP, LTGEOM,
>      2    LTPARM, NNPFIX, NNPFLO, MISCFX, MISCFL, NLPFIX, NLPFLO, 
>      3    LNODE, INDEX, X, DDOLD, COORD, S, C, H, DGH, B, E, EB, 
>      4    STRAIN, STRAN0, STRESS, BODY, PT, WT, XYZ, DLH, G, DLG,
>      5    AJ, AJINV, HINTG, D, PRTLPT, FLTNP, FLTEL, FLTMIS,
>      6    ELPROP, PRTLPT, PRTMAT, MISFIX, NPFIX, LPFIX, 
>      7    LPROP, LPPROP, NTAPE1, NTAPE2, NTAPE3, NTAPE4, NTAPE5,
>      8    LT, LTSHAP, LTUSER, NG )
140,141c139
< c         CALL  SKYSTR (NCOEFF, NDFREE, LTFREE, INDEX, IDIAG, 
< c    1                  S, SS)
---
> c         CALL  SKYSTR (NCOEFF,NDFREE,LTFREE,INDEX,IDIAG,S,SS)
147,148c145
<         IF ( NULCOL .EQ. 0 )  
<      1       CALL STORCL (NDFREE, LTFREE, INDEX, C, CC)
---
>         IF ( NULCOL .EQ. 0 )  CALL STORCL (NDFREE,LTFREE,INDEX,C,CC)
