anchor.diff:       & x,numnp,                                                         ! global
        & s,stemp,                                                         ! stiff
        & dmat,ien,lm,lmx,infiel,iddmat,ndmatsz,numelt,nconsz,             ! elemnt
  -     & infmat,numat,                                                    ! materl
  +     & infmat,infmatmod,numat,                                          ! materl
        & gauss,sh,shj,infetype,                                           ! eltype
        & skew,numrot,                                                     ! skew
        & nfault,dfault,tfault,numfn,                                      ! split
  
positive.code: c
        integer neq,numnp,ndmatsz,numelt,nconsz,numat,numrot,numfn,ierr
        integer ien(nconsz),lm(ndof,nconsz),lmx(ndof,nconsz)
        integer infiel(6,numelt),iddmat(nstr,nstr),infmat(6,numat)
        integer infetype(4,netypes),nfault(3,numfn)
        character errstrng*(*)
        double precision bdisp(neq),x(nsd,numnp),s(neemax*neemax)
        double precision stemp(neemax*neemax),dfault(ndof,numfn)
  
negative.code: c
        include "ndimens.inc"
        include "nshape.inc"
        include "materials.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  