anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine read_split(fault,nfault,numfn,numflt,ndof,numnp,kr,kw,
  -     & kp,idout,idsk,ierr,spfile,ofile,pfile)
  +      subroutine read_split(fault,nfault,numfn,numflt,numnp,kr,kw,
  +     & kp,idout,idsk,spfile,ofile,pfile,ierr,errstrng)
   c
   c...  reads and prints data on split nodes
   c
  
positive.code: c
  c...  open output files
  c
        if(idout.gt.0) open(kw,file=ofile,status="old",access="append")
        if(idsk.eq.0) open(kp,file=pfile,status="old",access="append")
        if(idsk.eq.1) open(kp,file=pfile,status="old",access="append",
       & form="unformatted")
  c
  c...  output results if requested
  c
        if(idsk.eq.0) write(kp,"(i7)") numfn
        if(idsk.eq.1) write(kp) numfn
        if(idout.gt.0.or.idsk.eq.0) then
          npage=50
          do i=1,numfn
            if((i.eq.1.or.mod(i,npage).eq.0).and.idout.gt.0)
       &     write(kw,3000)
            if(idout.gt.0) write(kw,4000) (nfault(j,i),j=1,3),
       &     (fault(j,i),j=1,ndof)
            if(idsk.eq.0) write(kp,5000) (nfault(j,i),j=1,2)
          end do
        end if
        if(idsk.eq.1.and.numfn.ne.0) write(kp) nfault
        if(idsk.eq.1.and.numfn.ne.0) write(kp) fault
  c
  c     compute total number of split nodes
  c
        if(numfn.ne.0) then
          do n=1,numnp
            do i=1,numfn
              if(n.eq.nfault(2,i)) then
                numflt=numflt+1
                goto 40
              end if
            end do
  
negative.code: c
   20   continue
          ierr=1
          errstrng="read_split"
          close(kr)
          return
  c
  