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:             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
  c...  local variables
  c
        integer ielga,ielgp,ielf,j,l
        double precision sout(3*nstatesmax)
  c
  cdebug      integer idb
  