anchor.diff:  c
         ielga=ielg
         ielgp=ielg
  -      ielgu=ielg
         call fill(sout,zero,3*nstatesmax)
   c
   c...  output to human-readable ascii file
  
positive.code:       end
  c
  c version
  c $Id: write_state_cmp.f,v 1.4 2005/05/03 18:42:05 willic3 Exp $
  c
  c End of file 
  
negative.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) npts,(sout(istatoutc(j)),j=1,nstatestot)
            end do
          end do
        end if
  c
  