anchor.diff:  
   !$OMP DO
         do i=1,n2ft3d
  -       n        = rho_in(i) + ETA
  -       agr      = agr_in(i)
  -       n_thrd 	= n**thrd
  +       n      = rho_in(i) + ETA
  +       agr    = agr_in(i)
  +       n_thrd = n**thrd
   
          if (rho_in(i).lt.DNS_CUT) then
            xe     = 0.0d0
  
positive.code:  
  
  
         xce(i) 	= x_parameter*xe   + c_parameter*ce
         fn(i)	= x_parameter*fdnx  + c_parameter*fdnc
         fdn(i) 	= x_parameter*fdagrx + c_parameter*fdagrc
  
        end do
  !$OMP END DO
  
negative.code:       real*8 DNS_CUT, ETA
        parameter (DNS_CUT        =      1.0d-18)
        parameter (ETA            =      1.0d-20)
  
        real*8 tolrho,minagr
        parameter (tolrho=2.0e-8,minagr=1.0e-10)
  ****** Becke constants *****************************
        real*8 lda_c,beta
        parameter (beta = 0.0042d0)
  