anchor.diff:  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
  
positive.code: c     velocities.
  c     Note also that the split node output is not really a UCD file,
  c     since it contains both node and element information.
  c
        include "implicit.inc"
  c
  
negative.code: c
  c
        subroutine write_ucd_node_vals(d,deld,tfault,dfault,nfault,numfn,
       & dx,deldx,idslp,numsn,deltp,nstep,numnp,kucd,iucd,ucdroot,
       & iprestress)
  c
  c...  Specialized routine to output displacement info for SCEC
  c     benchmarks.
  