anchor.diff:    do ixyz=1,3
       kAc(:,ixyz)=kAc0(:,ixyz)+Ac_tot(iter,ixyz)
     enddo
  -
  -!$omp parallel do
  -  do i=1,NL
  -    Vloc(i)=Vh(i)+Vpsl(i)+Vexc(i)
  -  end do
  +!$acc update device(kAc)
  +
  +!$acc kernels pcopy(Vloc) pcopyin(Vh,Vpsl,Vexc)
  +#ifndef _OPENACC
  +!$omp parallel do
  +#endif
  +  do i=1,NL
  +    Vloc(i)=Vh(i)+Vpsl(i)+Vexc(i)
  +  end do
  +!$acc end kernels
   
   !$acc end data
   
  
positive.code:   NVTX_BEG('dt_evolve_omp_KB',1)
    call timelog_begin(LOG_DT_EVOLVE)
  
  !$acc data pcopy(zu, ekr_omp)
  
    NVTX_BEG('nonlocal part?',2)
  !Constructing nonlocal part ! sato
  
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
  
  