fix: grow patch_ib on non-rank-0 before MPI broadcast of IB patches

When particle beds are used, rank 0 grows patch_ib to num_ib_patches_max
inside s_generate_particle_beds. Non-rank-0 ranks only have the
namelist-sized allocation (num_ib_patches_max_namelist). The num_ibs
scalar is broadcast first, so we can check and grow before the per-patch
MPI_BCAST loop accesses patch_ib(i) for i > num_ib_patches_max_namelist.

Author: sbryngelson

src/simulation/m_mpi_proxy.fpp

@@ -192,6 +192,19 @@ contains
             #:endfor
         end do
 
+        ! Non-rank-0 processes may have patch_ib sized to num_ib_patches_max_namelist while rank 0 grew it
+        ! for particle beds. Grow here before receiving the broadcast entries.
+        if (proc_rank /= 0 .and. num_ibs > size(patch_ib)) then
+            block
+                type(ib_patch_parameters), allocatable :: tmp(:)
+                integer                                :: n
+                n = size(patch_ib)
+                call move_alloc(patch_ib, tmp)
+                allocate (patch_ib(num_ib_patches_max))
+                patch_ib(1:n) = tmp
+            end block
+        end if
+
         do i = 1, num_ibs
             #:for VAR in [ 'radius', 'length_x', 'length_y', 'length_z', &
                 & 'x_centroid', 'y_centroid', 'z_centroid', 'c', 'm', 'p', 't', 'theta', 'slip', 'mass', &