anchor.diff:    do ixyz=1,3
       kAc(:,ixyz)=kAc0(:,ixyz)+0.5*(Ac_tot(iter,ixyz)+Ac_tot(iter+1,ixyz))
     enddo
  -
  -!$acc update device(zu)
  +!$acc update device(kAc)
  +
  +!$acc update device(zu, vloc)
     end select
   ! yabana
   
  
positive.code:   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 end data
  
  
negative.code: ! yabana
    select case(functional)
    case('VS98','TPSS','TBmBJ')
  !$acc update self(zu, ekr_omp, vloc)
  
  !$omp parallel do private(ik,ib)
    do ikb=1,NKB
  