anchor.diff:    end select
   ! yabana
   
  +  NVTX_BEG('dt_evolve_omp_KB_MS(): dt_evolve_hpsi',3)
     call dt_evolve_hpsi
  -
  +  NVTX_END()
  +
  +  NVTX_BEG('dt_evolve_omp_KB_MS(): psi_rho_RT',4)
     call psi_rho_RT
  +  NVTX_END()
  +
  +  NVTX_BEG('dt_evolve_omp_KB_MS(): Hartree',5)
     call Hartree
  -! yabana
  +  NVTX_END()
  +
  +! yabana
  +  NVTX_BEG('dt_evolve_omp_KB_MS(): Hartree',5)
     call Exc_Cor('RT')
  -! yabana
  -
  +  NVTX_END()
  +! yabana
  +
  +#ifdef _OPENACC
  +!$acc kernels pcopy(Vloc) pcopyin(Vh,Vpsl,Vexc)
  +#else
   !$omp parallel do
  +#endif
     do i=1,NL
       Vloc(i)=Vh(i)+Vpsl(i)+Vexc(i)
     end do
  +!$acc end kernels
  +
  +!$acc end data
   
     call timelog_end(LOG_DT_EVOLVE)
  +  NVTX_END()
   
     return
   End Subroutine dt_evolve_omp_KB_MS
  
positive.code:   end do
    end do
    end do
  #ifndef _OPENACC
  !$omp end parallel
  #endif
  !$acc end kernels
    NVTX_END()
  
  ! yabana
  
negative.code:   end do
    end do
    end do
  #ifdef _OPENACC
  !$acc end kernels
  #else
  !$omp end parallel
  #endif
  
  ! yabana
    select case(functional)
    case('VS98','TPSS','TBmBJ')
  !$acc update self(zu, ekr_omp, vloc)
  
  
  !$omp parallel do private(ik,ib)
  