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