Simplify hyper_cleaning init and use size() for IB broadcast

- Remove redundant if(p>0) branch in psi initialization: z_cc(0)=0
  in 2D, so the 3D formula produces identical results in both cases.
- Replace hardcoded 3 with size(patch_ib(i)%VAR) in MPI_BCAST for
  vel/angular_vel/angles, consistent with other IB broadcasts.

Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>

Author: sbryngelson

src/pre_process/m_mpi_proxy.fpp

@@ -118,7 +118,7 @@ contains
             end if
             ! Broadcast IB variables
             #:for VAR in ['vel', 'angular_vel', 'angles']
-                call MPI_BCAST(patch_ib(i)%${VAR}$, 3, mpi_p, 0, MPI_COMM_WORLD, ierr)
+                call MPI_BCAST(patch_ib(i)%${VAR}$, size(patch_ib(i)%${VAR}$), mpi_p, 0, MPI_COMM_WORLD, ierr)
             #:endfor
             call MPI_BCAST(patch_ib(i)%geometry, 1, MPI_INTEGER, 0, MPI_COMM_WORLD, ierr)
             call MPI_BCAST(patch_ib(i)%model_filepath, len(patch_ib(i)%model_filepath), MPI_CHARACTER, 0, MPI_COMM_WORLD, ierr)

src/pre_process/m_start_up.fpp

@@ -790,11 +790,7 @@ contains
             do l = 0, p
                 do k = 0, n
                     do j = 0, m
-                        if (p > 0) then
-                            q_cons_vf(psi_idx)%sf(j, k, l) = 1.0e-2_wp*exp(-(x_cc(j)**2 + y_cc(k)**2 + z_cc(l)**2)/(2.0_wp*0.05_wp**2))
-                        else
-                            q_cons_vf(psi_idx)%sf(j, k, l) = 1.0e-2_wp*exp(-(x_cc(j)**2 + y_cc(k)**2)/(2.0_wp*0.05_wp**2))
-                        end if
+                        q_cons_vf(psi_idx)%sf(j, k, l) = 1.0e-2_wp*exp(-(x_cc(j)**2 + y_cc(k)**2 + z_cc(l)**2)/(2.0_wp*0.05_wp**2))
                         q_prim_vf(psi_idx)%sf(j, k, l) = q_cons_vf(psi_idx)%sf(j, k, l)
                     end do
                 end do