anchor.diff:  c
    30   continue
           ierr=3
  +        errstrng="read_diff"
           close(kr)
  +        return
  +c
  +c...  error opening output file
  +c
  + 40   continue
  +        ierr=2
  +        errstrng="read_diff"
  +        close(kw)
  +        return
  +c
  +c...  error writing to output file
  +c
  + 50   continue
  +        ierr=4
  +        errstrng="read_diff"
  +        close(kw)
           return
   c
    6000 format(//' differential forces on slippery nodes'//
  
positive.code:       end
  c
  c version
  c $Id: read_diff.f,v 1.1 2004/04/14 21:18:30 willic3 Exp $
  c
  c Generated automatically by Fortran77Mill on Wed May 21 14:15:03 2003
  c
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine read_diff(diforc,nslip,idhist,numslp,numdif,numnp,
       & kr,kw,idout,difile,ofile,ierr,errstrng)
  c
  c...  reads and prints differential forces applied to slippery nodes
  c
  c     Error codes:
  c         0:  No error
  c         1:  Error opening input file (if numdif.ne.zero)
  c         2:  Error opening output file (if numdif.ne.zero)
  c         3:  Read error
  c         4:  Read error
  c       109:  Differential force applied to non-slippery node
  c
        include "implicit.inc"
  c
  c...  parameter definitions
  c
        include "ndimens.inc"
        include "nconsts.inc"
        include "rconsts.inc"
  c
  c...  subroutine arguments
  c
        integer numslp,numdif,numnp,kr,kw,idout,ierr
        integer nslip(nsdim,numslp),idhist(numnp)
        double precision diforc(ndof,numnp)
        character difile*(*),ofile*(*),errstrng*(*)
  c
  c...  included dimension and type statements
  c
        include "labeld_dim.inc"
  c
  c...  intrinsic functions
  c
  