anchor.diff:  c       write(*,*) "pbe96:",i,ec,fnc,fdnc
   
   
  -!$OMP END DO
  +ccc!$OMP END DO
   
         return
         end
  
positive.code:       fdnx_const = -3.0d0/(8.0d0*pi)
  
        
  !$OMP DO
        nup     = dn1_in+ETA
        agrup   = agr1_in
   
  
negative.code:       rs_scale   = (0.75d0/pi)**onethird
        fdnx_const = -3.0d0/(8.0d0*pi)
        
  ccc!$OMP DO
        n     = rho_in+ETA
        agr   = agr_in
        
  