anchor.diff:  
         
         do i=1,n2ft3d
  -         nup     = dn_in(i,1)
  +         nup     = dn_in(i,1)+ETA
            agrup   = agr_in(i,1)
  -         ndn     = dn_in(i,2)
  +
  + 
  +         ndn     = dn_in(i,2)+ETA
            agrdn   = agr_in(i,2)
  +
  +         
  + 
   c        ****************************************************************
   c        ***** calculate polarized Exchange energies and potentials *****
   c        ****************************************************************
  
positive.code: 
  
  *        **** calculate Wigner radius ****
           rs    = rs_scale/n13
           rss   = dsqrt(rs)
  
  *        **** calculate n*drs/dn ****
           rs_n = onethirdm*rs
  
  
  
negative.code: *
  * $Id: pbe96.f,v 1.12 2007-03-13 03:51:31 d3p708 Exp $
  *
  
  *    ************************************
  