anchor.diff:  c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
   c
   c
  -      subroutine read_stateout(istatout,kr,kw,kp,idout,idsk,stfile,
  -     & ofile,pfile,ierr,errstrng)
  +      subroutine read_stateout(istatout,nstatout,kr,kw,kp,idout,idsk,
  +     & stfile,ofile,pfile,ierr,errstrng)
   c
   c...     reads integer array indicating which state variables
   c        to output.  It is assumed that the same set of variables are
  
positive.code: c      7-12:   Total strain
  c      13-18:  Viscous strain
  c      18-24:  Plastic strain
  c
  c     Error codes:
  c         0:  No error
  
negative.code: c     The istatout array specifies output options for each individual
  c     state variable.  At present there are a maximum of 24 possible
  c     state variables, and this number may increase with the addition
  c     of new material models.  There are three types of state variable
  c     output:
  c
  c           1  Total accumulated values for the current time step
  c           2  Incremental values from the previous step to the current
  c           3  Rates computed from the previous step to the current
  c
  c      Present state variables occur in groups of 6, corresponding to
  c      the number of stress/strain components, although this may change
  