anchor.diff:    NVTX_BEG('dt_evolve_omp_KB',1)
     call timelog_begin(LOG_DT_EVOLVE)
   
  -!$acc data pcopy(zu, ekr_omp)
  +!$acc data pcopy(zu, ekr_omp, vloc)
   
     NVTX_BEG('nonlocal part?',2)
   !Constructing nonlocal part ! sato
  
positive.code: ! yabana
    select case(functional)
    case('VS98','TPSS','TBmBJ')
  !$acc update self(zu, ekr_omp)
  
  !$omp parallel do private(ik,ib)
    do ikb=1,NKB
  
negative.code:   do ixyz=1,3
      kAc(:,ixyz)=kAc0(:,ixyz)+0.5*(Ac_tot(iter,ixyz)+Ac_tot(iter+1,ixyz))
    enddo
  !$acc update device(kAc)
  
  !$acc update device(zu, vloc)
    end select
  ! yabana
  
  