anchor.diff:  c        ************
   c        **** up ****
   c        ************
  -         n     = 2.0d0*nup
  -         agr   = 2.0d0*agrup
  -
  -         n_onethird = (3.0d0*n/pi)**onethird
  +         n     = nup
  +         agr   = agrup
  +
  +         n_onethird = (6.0d0*n/pi)**onethird
            ex_lda     = -0.75d0*n_onethird
   
  -         kf = (3.0d0*pi*pi*n)**onethird
  +         kf = (6.0d0*pi*pi*n)**onethird
            s  = agr/(2.0d0*kf*n)
            P0 = 1.0d0 + (MU/KAPPA)*s*s
   
            F   = (1.0d0 + KAPPA - KAPPA/P0)
  -         Fs  = 2.0d0*MU/(P0*P0)*s
  +         Fs  = 2.0d0*MU/(P0*P0)*s*
   
            exup = ex_lda*F
            fnxup = fourthird*(exup - ex_lda*Fs*s)
  -         fdnxup = fdnx_const*Fs
  -
  -c         exup = ex_lda
  -c         fnxup = fourthird*(exup)
  -c         fdnxup = 0.0d0
  -
  -			
  +         fdnxup = ex_lda*Fs*0.5d0/kf
  +
   c        **************
   c        **** down ****
   c        **************
  -         n     = 2.0d0*ndn
  -         agr   = 2.0d0*agrdn
  -
  -         n_onethird = (3.0d0*n/pi)**onethird
  +         n     = ndn
  +         agr   = agrdn
  +
  +         n_onethird = (6.0d0*n/pi)**onethird
            ex_lda     = -0.75d0*n_onethird
   
  -         kf = (3.0d0*pi*pi*n)**onethird
  +         kf = (6.0d0*pi*pi*n)**onethird
            s  = agr/(2.0d0*kf*n)
            P0 = 1.0d0 + (MU/KAPPA)*s*s
   
  
positive.code:          s  = agr/(2.0d0*kf*n)
           P0 = 1.0d0 + (MU/KAPPA)*s*s
  
  c        if (n.gt.DNS_CUT) then
  c           F   = (1.0d0 + KAPPA - KAPPA/P0)
  c           Fs  = 2.0d0*MU/(P0*P0)*s
  c        else
  c           F   = 1.0d0
  c           Fs  = 0.0d0
  c        end if
           F   = (1.0d0 + KAPPA - KAPPA/P0)
           Fs  = 2.0d0*MU/(P0*P0)*s
  
  
negative.code: *           x_parameter: scale parameter for exchange
  *           c_parameter: scale parameter for correlation
  *
  *   Exit - xce(*)  : PBE96 energy density/ divided by electronic density ?
  *        - fn(*,2) : d(n*xce)/dnup, d(n*xce)/dndn
  *        - fdn(*,3): d(n*xce)/d|grad nup|, d(n*xce)/d|grad ndn|
  *                    d(n*xce)/d|grad n|
  