anchor.diff:  c
   c...compute current split node displacements
   c
  -      if(numfn.ne.0) then
  +      if(numfn.ne.izero) then
           call loadf(fault,dfault,histry,deltp,nfault,nstep,numfn,nhist,
        &   lastep,ierr,errstrng)
           if(ierr.ne.izero) return
  
positive.code: c
          if(ierr.ne.izero) return
  c
        else if(lgdefp.eq.0.and.ibbarp.eq.1) then
          call matinit_drv(
       &   alnz,ja,nnz,neq,                                               ! sparse
       &   x,d,iwink,wink,numnp,nwink,                                    ! global
  
negative.code: c
  c...  local variables
  c
        integer igroup,naxstp,nfirst
        double precision time,tminmax
        logical*4 fulout,skc
  c
  