anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer numfn,numsn,nstep,numnp,kucd,iprestress
  +      integer numfn,numsn,nstep,numnp,kucd,iucd,iprestress
         integer nfault(3,numfn),idslp(numsn)
         double precision deltp
         double precision d(ndof,numnp),deld(ndof,numnp)
  
positive.code: c
  c
        subroutine write_ucd_node_vals(d,deld,tfault,dfault,nfault,numfn,
       & dx,deldx,idslp,numsn,deltp,nstep,numnp,kucd,ucdroot,iprestress)
  c
  c...  Specialized routine to output displacement info for SCEC
  c     benchmarks.
  
negative.code: c...  local variables
  c
        integer i,j,i1,i2,inode
        integer ibyte,intlen,floatlen,len
        double precision rmult
        double precision vmin(6),vmax(6)
        character filenm*200,cstep*5
        character nlabels*1024,nunits*1024
  c
  cdebug      write(6,*) "Hello from write_ucd_node_vals!"
  c
        if(iucd.eq.izero) return
  c
  c...  set up labels and byte-counting info for binary output
  c
        nlabels="X-Displacement.Y-Displacement.Z-Displacement."//
       &        "X-Velocity.Y-Velocity.Z-Velocity."
        nlabels(1024:1024)="0"
        nunits="meters.meters.meters."//
       &       "meters/second.meters/second.meters/second."
        nunits(1024:1024)="0"
        ibyte=ione
        len=ione
        intlen=ifour
        floatlen=ifour
  c
  c...  determine filename
  c
        if(deltp.ne.zero) then
          rmult=one/deltp
  