anchor.diff:        fdnx_const = -3.0d0/(8.0d0*pi)
   
         
  -!$OMP DO
  +ccc!$OMP DO
         nup     = dn1_in+ETA
         agrup   = agr1_in
    
  
positive.code:       rs_scale   = (0.75d0/pi)**onethird
        fdnx_const = -3.0d0/(8.0d0*pi)
        
  !$OMP DO
        n     = rho_in+ETA
        agr   = agr_in
        
  
negative.code: 
        fdn = t_agr*Hpbe_t
  
  CCC!$OMP END DO
        
        
        
  