anchor.diff:       &     ngauss,ngaussdim,nen,nee,                                    ! eltype
        &     skew,nskdim,numrot,                                          ! skew
        &     getshape,bmatrix,                                            ! bbar
  -     &     idout,kto,kw,                                                ! ioinfo
        &     ierr)                                                        ! errcode
   c
             if(ierr.ne.izero) return
  
positive.code:      & gauss,sh,shj,infetype,netypes,                                   ! eltype
       & skew,nskdim,numrot,                                              ! skew
       & getshape,bmatrix,                                                ! bbar
       & idout,kto,kw,                                                    ! ioinfo
       & ierr)                                                            ! errcode
  c
  c...program to compute forces due to kinematic boundary conditions
  
negative.code: c...  subroutine arguments
  c
        integer neq,nsd,ndof,numnp,nstr,nddmat,ndmatsz,numelt,nconsz,numat
        integer netypes,nskdim,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)
  