Merged with master

Author: danieljvickers

src/simulation/m_data_output.fpp

@@ -1062,6 +1062,8 @@ contains
 
         integer, intent(in) :: time_step
 
+        $:GPU_UPDATE(host='[patch_ib(1:num_ibs)]')
+
         if (parallel_io) then
             call s_write_parallel_ib_state(time_step)
         else

src/simulation/m_ibm.fpp

@@ -936,9 +936,9 @@ contains
             end do
         end if
 
-        $:GPU_PARALLEL_LOOP(private='[ib_idx, ib_idx_temp, encoded_ib_idx, fluid_idx, radial_vector, local_force_contribution, &
-                            & cell_volume, local_torque_contribution, dynamic_viscosity, viscous_stress]', copy='[forces, &
-                            & torques]', copyin='[dynamic_viscosities]', collapse=3)
+        $:GPU_PARALLEL_LOOP(private='[i, j, k, l, ib_idx, ib_idx_temp, encoded_ib_idx, fluid_idx, radial_vector, &
+                            & local_force_contribution, cell_volume, local_torque_contribution, dynamic_viscosity, &
+                            & viscous_stress]', copy='[forces, torques]', copyin='[dynamic_viscosities]', collapse=3)
         do i = 0, m
             do j = 0, n
                 do k = 0, p

src/simulation/m_start_up.fpp

@@ -1494,6 +1494,13 @@ contains
                     call MPI_SENDRECV(end_val, 1, mpi_p, send_neighbor, ${TAG + 1}$, recv_val, 1, mpi_p, recv_neighbor, &
                                       & ${TAG + 1}$, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                     end_val = recv_val
+
+                    ! protect from looping back around on yourself multiple times
+                    if (f_approx_equal(beg_val, ${X}$_cb(${DIM}$)) .or. f_approx_equal(end_val, ${X}$_cb(-1))) then
+                        beg_val = -huge(0._wp)
+                        end_val = huge(0._wp)
+                        exit
+                    end if
                 end do
                 neighbor_domain_${X}$%beg = beg_val
                 neighbor_domain_${X}$%end = end_val