Updated pressure communication for multi-node case

Author: danieljvickers

src/simulation/m_ibm.fpp

@@ -998,7 +998,9 @@ contains
     ! by Archana Sridhar and Jesse Capecelatro
     subroutine s_compute_ib_forces(pressure)
 
-        real(wp), dimension(0:m, 0:n, 0:p), intent(in) :: pressure
+        real(wp), dimension(idwbuff(1)%beg:idwbuff(1)%end, &
+                    idwbuff(2)%beg:idwbuff(2)%end, &
+                    idwbuff(3)%beg:idwbuff(3)%end), intent(in) :: pressure
 
         integer :: i, j, k, ib_idx
         real(wp), dimension(num_ibs, 3) :: forces, torques
@@ -1061,8 +1063,6 @@ contains
             patch_ib(i)%torque(:) = matmul(patch_ib(i)%rotation_matrix_inverse, torques(i, :)) ! torques must be computed in the local coordinates of the IB
         end do
 
-        print *, forces(1, 1:2)
-
     end subroutine s_compute_ib_forces
 
     !> Subroutine to deallocate memory reserved for the IBM module

src/simulation/m_time_steppers.fpp

@@ -611,7 +611,7 @@ contains
             if (ib) then
                 ! check if any IBMS are moving, and if so, update the markers, ghost points, levelsets, and levelset norms
                 if (moving_immersed_boundary_flag) then
-                    call s_compute_ib_forces(q_prim_vf(E_idx)%sf(0:m, 0:n, 0:p)) ! compute the force and torque on the IB from the fluid
+                    call s_compute_ib_forces(q_prim_vf(E_idx)%sf) ! compute the force and torque on the IB from the fluid
 
                     do i = 1, num_ibs
                         if (s == 1) then