anchor.diff:  c
   c...  open input file and define scaling factors
   c
  -      ierr=0
  +      ierr=izero
  +      nconcforce=izero
         open(kr,file=bcfile,status="old",err=20)
         scale(1)=zero
         scale(2)=dscale
  
positive.code:           ihist(j)=imode/10
            itype(j)=imode-10*ihist(j)
            bond(j,n)=scale(itype(j)+1)*bond(j,n)
          end do
        end do
        close(kr)
  
negative.code: c
  c...  subroutine arguments
  c
        integer numnp,numbc,neq,nconcforce,kr,kw,idout,ierr
        integer ibond(ndof,numnp),id(ndof,numnp)
        double precision bond(ndof,numnp)
        double precision dscale,vscale,fscale
  