anchor.diff:        end
   c
   c version
  -c $Id: write_state_cmp.f,v 1.4 2005/05/03 18:42:05 willic3 Exp $
  +c $Id: write_state_cmp.f,v 1.5 2005/05/13 15:16:51 willic3 Exp $
   c
   c End of file 
  
positive.code:             call daxpy(nstatesmax,delti,sout(nstatesmax+ione),ione,
       &       sout(2*nstatesmax+ione),ione)
  cdebug            write(6,*) "sout:",(sout(idb),idb=1,3*nstatesmax)
              write(kucd,720) ielgu,(sout(istatoutc(j)),j=1,nstatestot)
            end do
            ielgu=ielgu+ione
          end do
        end if
  c
  
negative.code: c
        ielga=ielg
        ielgp=ielg
        call fill(sout,zero,3*nstatesmax)
  c
  c...  output to human-readable ascii file
  