anchor.diff:       & gauss,sh,shj,ngauss,ngaussdim,nen,                               ! eltype
        & s,nee,                                                           ! stiff
        & getshape,bmatrix,                                                ! bbar
  -     & idout,kto,kw,                                                    ! ioinf
        & ierr)                                                            ! errcode
   c
         if(ierr.ne.izero) return
  
positive.code:      & gauss,sh,shj,ngauss,ngaussdim,nen,nee,                           ! eltype
       & skew,nskdim,numrot,                                              ! skew
       & getshape,bmatrix,                                                ! bbar
       & idout,kto,kw,                                                    ! ioinfo
       & ierr)                                                            ! errcode
  c
  c...  subroutine to form the elemental stiffness matrix
  
negative.code: c...  subroutine arguments
  c
        integer nsd,ndof,numnp,nstr,nddmat,iel,ngauss,ngaussdim,nen,nee
        integer nskdim,numrot,ierr
        integer ien(nen),lm(ndof,nen),iddmat(nstr,nstr)
        double precision x(nsd,numnp),s(neemax*neemax)
        double precision stemp(neemax*neemax),dmat(nddmat,ngaussdim)
  