anchor.diff:        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
  
positive.code: c
  c...  output to UCD file
  c
        if(iucd.ne.0) then
          do ielf=1,nelfamily
            do l=1,ngauss
              npts=npts+ione
  
negative.code: c
        ielga=ielg
        ielgp=ielg
        floatlen=ifour
        call fill(sout,zero,3*nstatesmax)
  c
  c...  output to human-readable ascii file
  