anchor.diff:  c
   c...  output to UCD file
   c
  -      if(iucd.ne.0) then
  +      if(iucd.eq.1) then
           do ielf=1,nelfamily
             do l=1,ngauss
               npts=npts+ione
  
positive.code: c
        ielga=ielg
        ielgp=ielg
        call fill(sout,zero,3*nstatesmax)
  c
  c...  output to human-readable ascii file
  
negative.code:       integer nelfamily,nstate,ielg,ngauss,nstep,nstatestot,nout,npts
        integer idout,idsk,iucd,kw,kp,kucd
        integer istatout(nstatesmax,3),istatoutc(nstatestot),nstatout(3)
        integer ibyteg(3*nstatesmax)
        double precision state(nstate,ngauss,nelfamily)
        double precision dstate(nstate,ngauss,nelfamily)
        double precision statemin(3*nstatesmax),statemax(3*nstatesmax)
        double precision delti
  c
  c...  external routines
  