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:     call timelog_end(LOG_OTHER)
  
    enddo RTiteratopm !end of RT iteraction========================
    etime2=MPI_WTIME()
    call timelog_disable_verbose
  
  
negative.code:     open(943,file=trim(directory)//'Ac_M'//trim(cMacro_x)//'.out',position = position_option)
    end if
  
  !$acc enter data copyin(ik_table,ib_table)
  !$acc enter data copyin(lapx,lapy,lapz)
  !$acc enter data copyin(nabx,naby,nabz)
  !$acc enter data copyin(modx,mody,modz)
  !$acc enter data copyin(zJxyz,zKxyz)
  !$acc enter data copyin(uV,iuV)
  
  !$acc enter data create(kAc)
  
    call timelog_reset
    call timelog_enable_verbose
    etime1=MPI_WTIME()
  !$acc enter data copyin(zu)
    RTiteratopm : do iter=entrance_iter+1,Nt ! sato
  
      call dt_evolve_Ac ! sato
  