anchor.diff:  c
         integer neq,numnp,ndmatsz,numelt,nconsz,numat,numrot,ierr
         integer ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz)
  -      integer infiel(6,numelt),iddmat(nstr,nstr),infmat(6,numat)
  -      integer infetype(4,netypes)
  +      integer infiel(6,numelt),iddmat(nstr,nstr),infmat(3,numat)
  +      integer infmatmod(5,nmatmodmax),infetype(4,netypes)
         character errstrng*(*)
         double precision bdisp(neq),x(nsd,numnp),d(ndof,numnp)
         double precision deld(ndof,numnp),s(neemax*neemax)
  
positive.code: c...  loop over material groups
  c
        do imat=1,numat
          nmatel=infmat(2,imat)
          imatvar=infmat(4,imat)
          ngtest=0
          if(imatvar.eq.izero) ngtest=ngaussmax
          do ielg=matgpt,matgpt+nmatel-1
  
negative.code:      & x,d,deld,numnp,                                                  ! global
       & s,stemp,                                                         ! stiff
       & dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,             ! elemnt
       & infmat,infmatmod,numat,                                          ! materl
       & gauss,sh,shj,infetype,                                           ! eltype
       & skew,numrot,                                                     ! skew
       & getshape,bmatrix,                                                ! bbar
  