anchor.diff:           fnxup = fourthird*(exup - ex_lda*Fs*s)
            fdnxup = fdnx_const*Fs
   
  +c         exup = ex_lda
  +c         fnxup = fourthird*(exup)
  +c         fdnxup = 0.0d0
  +
  +			
   c        **************
   c        **** down ****
   c        **************
            n     = 2.0d0*ndn
            agr   = 2.0d0*agrdn
   
  -	 n13   = n**onethird         
  -	 n13   = max(n13,1.d-14)
  -         n_onethird = n13*((3.0d0/pi)**onethird)
  +         n_onethird = (3.0d0*n/pi)**onethird
            ex_lda     = -0.75d0*n_onethird
   
  -         kf = n13*((3.0d0*pi*pi)**onethird)
  -         s  = agr/(2.0d0*kf*(n13**3))
  +         kf = (3.0d0*pi*pi*n)**onethird
  +         s  = agr/(2.0d0*kf*n)
            P0 = 1.0d0 + (MU/KAPPA)*s*s
   
            F   = (1.0d0 + KAPPA - KAPPA/P0)
  
negative.code: 
           ex   = ex_lda*F
           fnx  = fourthird*(ex - ex_lda*Fs*s) 
           fdnx = fdnx_const*Fs
  
  			
  
  