anchor.diff:        kAc(:,1)=kAc0(:,1)
         kAc(:,2)=kAc0(:,2)+Ac_new_m(2,ix_m,iy_m)
         kAc(:,3)=kAc0(:,3)+Ac_new_m(3,ix_m,iy_m)
  +!$acc update device(kAc)
   ! sato ---------------------------------------
         call timelog_end(LOG_OTHER)
   
  
positive.code:       kAc(:,1)=kAc0(:,1)
        kAc(:,2)=kAc0(:,2)+(Ac_new_m(2,ix_m,iy_m)+Ac_m(2,ix_m,iy_m))/2d0
        kAc(:,3)=kAc0(:,3)+(Ac_new_m(3,ix_m,iy_m)+Ac_m(3,ix_m,iy_m))/2d0
  ! sato ---------------------------------------
        call timelog_end(LOG_OTHER)
        
  
negative.code: 
        javt(iter,:)=jav(:)
        if (MD_option == 'Y') then
  !$acc update self(zu)
          call Ion_Force_omp(Rion_update,'RT')
          if (iter/Nstep_write*Nstep_write == iter) then ! sato
            call Total_Energy_omp(Rion_update,'RT')
          end if
        else
          if (iter/Nstep_write*Nstep_write == iter) then ! sato
  !$acc update self(zu)
            call Total_Energy_omp(Rion_update,'RT')
            call Ion_Force_omp(Rion_update,'RT')
          end if
  