anchor.diff:  c
   c...  subroutine arguments
   c
  -      integer nstatestot,kucd
  +      integer nstatestot,kucd,iucd
         integer istatoutc(nstatestot)
   c
   c... included dimension and type statements
   c
         include "labels_dim.inc"
   c
  +c...  local constants
  +c
  +      character stress*6,strain*4,strate*9
  +      data stress,strain,strate/"Pascal","None","1/seconds"/
  +c
   c...  local variables
   c
  -      integer i
  +      integer i,ibyte,indl,indu
         integer iout(3*nstatesmax)
  +      character nlabels*1024,nunits*1024
   c
   c...  included variable definitions
   c
  
positive.code:         iout(i)=ione
        end do
  c
        write(kucd,"(100i5)") nstatestot,(iout(i),i=1,nstatestot)
  c
        do i=1,nstatestot
          write(kucd,"(a11)") labels(istatoutc(i))
        end do
  c
        return
        end
  
negative.code: c~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  c
  c
        subroutine write_ucd_header(istatoutc,nstatestot,kucd,iucd)
  c
  c...subroutine to write the header for a UCD file containing state
  c   variables
  