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:       end
  c
  c version
  c $Id: formes_ss.f,v 1.2 2004/06/16 20:24:46 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  