anchor.diff:        call MPI_BARRIER(MPI_COMM_WORLD,ierr)
         write(*,*) Myrank,'iter =',iter
         iter_now=iter
  +!$acc update self(zu)
         call prep_Reentrance_write
         go to 1
       end if
  
positive.code: 
        javt(iter,:)=jav(:)
        if (MD_option == 'Y') then
          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
            call Total_Energy_omp(Rion_update,'RT')
            call Ion_Force_omp(Rion_update,'RT')
          end if
  
negative.code:     call timelog_end(LOG_OTHER)
  
    enddo RTiteratopm !end of RT iteraction========================
  !$acc exit data copyout(zu)
    etime2=MPI_WTIME()
    call timelog_disable_verbose
  
  