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:       else
          cstep="prest"
        end if
  c
  c...  write mesh info
  c
        filenm=ucdroot(i1:i2)//".mesh.time."//cstep//".inp"
        open(kucd,file=filenm,status="new")
        write(kucd,"(7i7)") nnvals,(ival(i),i=1,nnvals)
        do i=1,nnvals
          write(kucd,"(a9)") dout(i)
        end do
  c
  c...  write nodal displacements
  c
        do i=1,numnp
          write(kucd,"(i7,6(2x,1pe15.8))") i,(d(j,i),j=1,ndof),
       &   (rmult*deld(j,i),j=1,ndof)
        end do
        close(kucd)
  c
  c...  write split node displacements, if there are any
  c
        if(numfn.ne.0) then
          filenm=ucdroot(i1:i2)//".mesh.split.time."//cstep//".inp"
          open(kucd,file=filenm,status="new")
          write(kucd,"(7i7)") nnvals,(ival(i),i=1,nnvals)
          do i=1,nnvals
            write(kucd,"(a9)") dout(i)
          end do
          do i=1,numfn
            write(kucd,"(2i7,6(2x,1pe15.8))") nfault(1,i),nfault(2,i),
       &     (tfault(j,i),j=1,ndof),(rmult*dfault(j,i),j=1,ndof)
          end do
          close(kucd)
        end if
  c
  c...  write slippery node displacements, if there are any
  c
        if(numsn.ne.0) then
          filenm=ucdroot(i1:i2)//".mesh.slip.time."//cstep//".inp"
          open(kucd,file=filenm,status="new")
          write(kucd,"(7i7)") nnvals,(ival(i),i=1,nnvals)
          do i=1,nnvals
            write(kucd,"(a9)") dout(i)
          end do
          do i=1,numsn
            inode=idslp(i)
            write(kucd,"(i7,6(2x,1pe15.8))") idslp(inode),
       &     (dx(j,inode),j=1,ndof),(rmult*deldx(j,inode),j=1,ndof)
          end do
          close(kucd)
        end if
        return
        end
  c
  c version
  c $Id: write_ucd_node_vals.f,v 1.4 2005/04/14 01:01:22 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.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     Either binary or ascii UCD files are written, depending on the
  c     value of the iucd parameter.
  c
        include "implicit.inc"
  c
  