Fix hyper_cleaning psi initialization only covering k=0 plane

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

Author: sbryngelson

src/pre_process/m_start_up.fpp

@@ -767,7 +767,7 @@ contains
 
         real(wp), intent(inout) :: start, finish
 
-        integer :: j, k
+        integer :: j, k, l
 
         ! Setting up the grid and the initial condition. If the grid is read in from
         ! preexisting grid data files, it is checked for consistency. If the grid is
@@ -787,10 +787,12 @@ contains
 
         ! hard-coded psi
         if (hyper_cleaning) then
-            do j = 0, m
+            do l = 0, p
                 do k = 0, n
-                    q_cons_vf(psi_idx)%sf(j, k, 0) = 1d-2*exp(-(x_cc(j)**2 + y_cc(k)**2)/(2.0*0.05**2))
-                    q_prim_vf(psi_idx)%sf(j, k, 0) = q_cons_vf(psi_idx)%sf(j, k, 0)
+                    do j = 0, m
+                        q_cons_vf(psi_idx)%sf(j, k, l) = 1d-2*exp(-(x_cc(j)**2 + y_cc(k)**2 + z_cc(l)**2)/(2.0*0.05**2))
+                        q_prim_vf(psi_idx)%sf(j, k, l) = q_cons_vf(psi_idx)%sf(j, k, l)
+                    end do
                 end do
             end do
         end if