Apply ffmt: remove blank lines after #else/#elif continuations

sbryngelson · sbryngelson · commit 023a5b473cf3 · 2026-03-21T15:16:03.000-04:00
diff --git a/src/common/include/3dHardcodedIC.fpp b/src/common/include/3dHardcodedIC.fpp
@@ -140,7 +140,6 @@
 
         q_prim_vf(E_idx)%sf(i, j, k) = p_th*rcut*xcut + p_am
     case (303) ! 3D Multijet
-
         eps_smooth = 3.0_wp
         ux_th = 10*sqrt(1.4*0.4)
         ux_am = 2.5*sqrt(1.4*0.4)
diff --git a/src/common/m_mpi_common.fpp b/src/common/m_mpi_common.fpp
@@ -334,7 +334,6 @@ contains
             call MPI_REDUCE(Rc_min_loc, Rc_min_glb, 1, mpi_p, MPI_MIN, 0, MPI_COMM_WORLD, ierr)
         end if
 #else
-
         icfl_max_glb = icfl_max_loc
 
         if (viscous) then
@@ -1046,7 +1045,6 @@ contains
                         end if
                     end do
                 else
-
                     if (cyl_coord .and. p > 0) then
                         ! Implement pencil processor blocking if using cylindrical coordinates so that all cells in azimuthal
                         ! direction are stored on a single processor. This is necessary for efficient application of Fourier filter
@@ -1080,7 +1078,6 @@ contains
                             end if
                         end do
                     else
-
                         ! Initial estimate of optimal processor topology
                         num_procs_x = 1
                         num_procs_y = 1
@@ -1206,7 +1203,6 @@ contains
 
                 ! 2D Cartesian Processor Topology
             else
-
                 ! Initial estimate of optimal processor topology
                 num_procs_x = 1
                 num_procs_y = num_procs
@@ -1320,7 +1316,6 @@ contains
 
             ! 1D Cartesian Processor Topology
         else
-
             ! Optimal processor topology
             num_procs_x = num_procs
 
@@ -1428,20 +1423,17 @@ contains
                     call MPI_SENDRECV(dx(m - buff_size + 1), buff_size, mpi_p, bc_x%end, 0, dx(-buff_size), buff_size, mpi_p, &
                                       & bc_x%beg, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the beginning only
-
                     ! Send/receive buffer to/from bc_x%beg/bc_x%beg
                     call MPI_SENDRECV(dx(0), buff_size, mpi_p, bc_x%beg, 1, dx(-buff_size), buff_size, mpi_p, bc_x%beg, 0, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 end if
             else ! PBC at the end
-
                 if (bc_x%beg >= 0) then ! PBC at the end and beginning
 
                     ! Send/receive buffer to/from bc_x%beg/bc_x%end
                     call MPI_SENDRECV(dx(0), buff_size, mpi_p, bc_x%beg, 1, dx(m + 1), buff_size, mpi_p, bc_x%end, 1, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the end only
-
                     ! Send/receive buffer to/from bc_x%end/bc_x%end
                     call MPI_SENDRECV(dx(m - buff_size + 1), buff_size, mpi_p, bc_x%end, 0, dx(m + 1), buff_size, mpi_p, &
                                       & bc_x%end, 1, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
@@ -1459,20 +1451,17 @@ contains
                     call MPI_SENDRECV(dy(n - buff_size + 1), buff_size, mpi_p, bc_y%end, 0, dy(-buff_size), buff_size, mpi_p, &
                                       & bc_y%beg, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the beginning only
-
                     ! Send/receive buffer to/from bc_y%beg/bc_y%beg
                     call MPI_SENDRECV(dy(0), buff_size, mpi_p, bc_y%beg, 1, dy(-buff_size), buff_size, mpi_p, bc_y%beg, 0, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 end if
             else ! PBC at the end
-
                 if (bc_y%beg >= 0) then ! PBC at the end and beginning
 
                     ! Send/receive buffer to/from bc_y%beg/bc_y%end
                     call MPI_SENDRECV(dy(0), buff_size, mpi_p, bc_y%beg, 1, dy(n + 1), buff_size, mpi_p, bc_y%end, 1, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the end only
-
                     ! Send/receive buffer to/from bc_y%end/bc_y%end
                     call MPI_SENDRECV(dy(n - buff_size + 1), buff_size, mpi_p, bc_y%end, 0, dy(n + 1), buff_size, mpi_p, &
                                       & bc_y%end, 1, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
@@ -1482,7 +1471,6 @@ contains
 
             ! MPI Communication in z-direction
         else
-
             if (pbc_loc == -1) then ! PBC at the beginning
 
                 if (bc_z%end >= 0) then ! PBC at the beginning and end
@@ -1491,20 +1479,17 @@ contains
                     call MPI_SENDRECV(dz(p - buff_size + 1), buff_size, mpi_p, bc_z%end, 0, dz(-buff_size), buff_size, mpi_p, &
                                       & bc_z%beg, 0, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the beginning only
-
                     ! Send/receive buffer to/from bc_z%beg/bc_z%beg
                     call MPI_SENDRECV(dz(0), buff_size, mpi_p, bc_z%beg, 1, dz(-buff_size), buff_size, mpi_p, bc_z%beg, 0, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 end if
             else ! PBC at the end
-
                 if (bc_z%beg >= 0) then ! PBC at the end and beginning
 
                     ! Send/receive buffer to/from bc_z%beg/bc_z%end
                     call MPI_SENDRECV(dz(0), buff_size, mpi_p, bc_z%beg, 1, dz(p + 1), buff_size, mpi_p, bc_z%end, 1, &
                                       & MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
                 else ! PBC at the end only
-
                     ! Send/receive buffer to/from bc_z%end/bc_z%end
                     call MPI_SENDRECV(dz(p - buff_size + 1), buff_size, mpi_p, bc_z%end, 0, dz(p + 1), buff_size, mpi_p, &
                                       & bc_z%end, 1, MPI_COMM_WORLD, MPI_STATUS_IGNORE, ierr)
diff --git a/src/common/m_phase_change.fpp b/src/common/m_phase_change.fpp
@@ -203,7 +203,6 @@ contains
                             ! correcting the vapor partial density
                             q_cons_vf(vp + contxb - 1)%sf(j, k, l) = mixM*rM
                         else
-
                             ! returning partial pressures to what they were from the homogeneous solver liquid
                             q_cons_vf(lp + contxb - 1)%sf(j, k, l) = m1
 
@@ -604,7 +603,6 @@ contains
             ! assigning Saturation temperature
             TSat = 0.0_wp
         else
-
             ! calculating initial estimate for temperature in the TSat procedure. I will also use this variable to iterate over the
             ! Newton's solver
             TSat = TSIn
diff --git a/src/common/m_variables_conversion.fpp b/src/common/m_variables_conversion.fpp
@@ -152,7 +152,6 @@ contains
                 pres = (energy - 0.5_wp*(mom**2._wp)/rho - pi_inf - qv - E_e)/gamma
             end if
         #:else
-
             Y_rs(:) = rhoYks(:)/rho
             e_Per_Kg = energy/rho
             Pdyn_Per_Kg = dyn_p/rho
@@ -418,7 +417,6 @@ contains
 
                 ! Simulation is 2D
             else
-
                 allocate (rho_sf(-buff_size:m + buff_size, -buff_size:n + buff_size, 0:0))
                 allocate (gamma_sf(-buff_size:m + buff_size, -buff_size:n + buff_size, 0:0))
                 allocate (pi_inf_sf(-buff_size:m + buff_size, -buff_size:n + buff_size, 0:0))
@@ -427,7 +425,6 @@ contains
 
             ! Simulation is 1D
         else
-
             allocate (rho_sf(-buff_size:m + buff_size, 0:0, 0:0))
             allocate (gamma_sf(-buff_size:m + buff_size, 0:0, 0:0))
             allocate (pi_inf_sf(-buff_size:m + buff_size, 0:0, 0:0))
diff --git a/src/post_process/m_data_output.fpp b/src/post_process/m_data_output.fpp
@@ -204,7 +204,6 @@ contains
 
             ! Generating Binary Directory Tree
         else
-
             ! Creating the directory associated with the local process
             dbdir = trim(case_dir) // '/binary'
 
@@ -451,7 +450,6 @@ contains
 
             ! Binary Database Format
         else
-
             ! Generating the relative path to the formatted database slave file, that is to be opened for the current time-step,
             ! t_step
             write (file_loc, '(A,I0,A)') '/', t_step, '.dat'
@@ -645,7 +643,6 @@ contains
                 ! One-dimensional local grid data is written to the formatted database slave file. In addition, the local grid data
                 ! is put together by the root process and written to the master file.
             else
-
                 if (precision == 1) then
                     write (dbfile) real(x_cb, sp)
                 else
@@ -805,7 +802,6 @@ contains
 
             ! Binary Database Format
         else
-
             ! Writing the name of the flow variable and its data, associated with the local processor, to the formatted database
             ! slave file
             if (precision == 1) then
diff --git a/src/post_process/m_global_parameters.fpp b/src/post_process/m_global_parameters.fpp
@@ -892,7 +892,6 @@ contains
             ! Allocating the grid variables, only used for the 1D simulations, and containing the defragmented computational domain
             ! grid data
         else
-
             allocate (x_root_cb(-1:m_root))
             allocate (x_root_cc(0:m_root))
 
diff --git a/src/post_process/m_mpi_proxy.fpp b/src/post_process/m_mpi_proxy.fpp
@@ -217,7 +217,6 @@ contains
             call MPI_GATHERV(maxval(y_cb), 1, mpi_p, spatial_extents(4, 0), recvcounts, 4*displs, mpi_p, 0, MPI_COMM_WORLD, ierr)
             ! Simulation is 1D
         else
-
             ! For 1D, recvcounts/displs are sized for grid defragmentation (m+1 per rank), not for scalar gathers. Use MPI_GATHER
             ! instead.
 
@@ -246,7 +245,6 @@ contains
 
             ! Binary database format
         else
-
             call MPI_GATHERV(x_cb(0), m + 1, mpi_p, x_root_cb(0), recvcounts, displs, mpi_p, 0, MPI_COMM_WORLD, ierr)
 
             if (proc_rank == 0) x_root_cb(-1) = x_cb(-1)
diff --git a/src/pre_process/m_grid.f90 b/src/pre_process/m_grid.f90
@@ -85,7 +85,6 @@ impure subroutine s_generate_serial_grid
                 y_cb(i - 1) = y_domain%beg + dy*real(2*i - 1, wp)
             end do
         else
-
             dy = (y_domain%end - y_domain%beg)/real(n + 1, wp)
 
             do i = 0, n
diff --git a/src/pre_process/m_icpp_patches.fpp b/src/pre_process/m_icpp_patches.fpp
@@ -145,7 +145,6 @@ contains
 
             ! 1D Patch Geometries
         else
-
             do i = 1, num_patches
                 if (proc_rank == 0) then
                     print *, 'Processing patch', i
diff --git a/src/simulation/m_bubbles_EL.fpp b/src/simulation/m_bubbles_EL.fpp
@@ -586,7 +586,6 @@ contains
                 gas_p(k, 1) = myPb
                 gas_mv(k, 1) = myMass_v
             else
-
                 ! Radial acceleration from bubble models
                 intfc_dveldt(k, stage) = f_rddot(myRho, myPinf, myR, myV, myR0, myPb, myPbdot, dmalf, dmntait, dmBtait, &
                              & dm_bub_adv_src, dm_divu, myCson)
diff --git a/src/simulation/m_bubbles_EL_kernels.fpp b/src/simulation/m_bubbles_EL_kernels.fpp
@@ -230,7 +230,6 @@ contains
                                                       & **(3._wp*(strength_idx + 1._wp))
                 end do
             else
-
                 !> 2D cartesian function:
                 ! We smear particles considering a virtual depth (lag_params%charwidth)
                 theta = 0._wp
diff --git a/src/simulation/m_cbc.fpp b/src/simulation/m_cbc.fpp
@@ -398,7 +398,6 @@ contains
 
                     ! Computing CBC4 Coefficients
                 else
-
                     fd_coef_${XYZ}$ (:, cbc_loc_in) = 0._wp
                     fd_coef_${XYZ}$ (0, &
                                      & cbc_loc_in) = -50._wp/(25._wp*ds(0) + 2._wp*ds(1) - 1.e1_wp*ds(2) + 1.e1_wp*ds(3) &
@@ -479,7 +478,6 @@ contains
 
             ! Associating CBC Coefficients in z-direction
         else
-
             ! fd_coef => fd_coef_z; if (weno_order > 1) pi_coef => pi_coef_z
 
             if (cbc_loc_in == -1) then
@@ -948,7 +946,6 @@ contains
                                                         & i) + ds(0)*dadv_dt(i - E_idx))
                             end do
                         else
-
                             $:GPU_LOOP(parallelism='[seq]')
                             do i = advxb, advxe
                                 flux_rs${XYZ}$_vf_l(-1, k, r, i) = flux_rs${XYZ}$_vf_l(0, k, r, i) + ds(0)*dadv_dt(i - E_idx)
@@ -1158,7 +1155,6 @@ contains
 
             ! Reshaping Inputted Data in z-direction
         else
-
             $:GPU_PARALLEL_LOOP(private='[i, j, k, r]', collapse=4)
             do i = 1, sys_size
                 do r = is3%beg, is3%end
@@ -1349,7 +1345,6 @@ contains
 
             ! Reshaping Outputted Data in z-direction
         else
-
             $:GPU_PARALLEL_LOOP(private='[i, j, k, r]', collapse=4)
             do i = 1, flux_cbc_index
                 do r = is3%beg, is3%end
diff --git a/src/simulation/m_fftw.fpp b/src/simulation/m_fftw.fpp
@@ -74,7 +74,6 @@ contains
 
 #if defined(MFC_GPU)
         rank = 1; istride = 1; ostride = 1
-
         allocate (gpu_fft_size(1:rank), iembed(1:rank), oembed(1:rank))
 
         gpu_fft_size(1) = real_size;
diff --git a/src/simulation/m_global_parameters.fpp b/src/simulation/m_global_parameters.fpp
@@ -905,7 +905,6 @@ contains
 
             ! Volume Fraction Model
         else
-
             ! Annotating structure of the state and flux vectors belonging to the system of equations defined by the selected number
             ! of spatial dimensions and the volume fraction model
             if (model_eqns == 2) then
diff --git a/src/simulation/m_ibm.fpp b/src/simulation/m_ibm.fpp
@@ -698,7 +698,6 @@ contains
                         interp_coeffs(:,:, 1) = eta(:,:, 1)/buf
                     end if
                 else
-
                     if (ib_markers%sf(i, j, k + 1) /= 0) alpha(1, 1, 2) = 0._wp
                     if (ib_markers%sf(i + 1, j, k + 1) /= 0) alpha(2, 1, 2) = 0._wp
                     if (ib_markers%sf(i, j + 1, k + 1) /= 0) alpha(1, 2, 2) = 0._wp
diff --git a/src/simulation/m_rhs.fpp b/src/simulation/m_rhs.fpp
@@ -1449,7 +1449,6 @@ contains
                 $:END_GPU_PARALLEL_LOOP()
             end if
         else if (idir == 2) then ! y-direction
-
             if (surface_tension) then
                 $:GPU_PARALLEL_LOOP(private='[j, k, l]', collapse=3)
                 do l = 0, p
@@ -1504,7 +1503,6 @@ contains
                 end do
                 $:END_GPU_PARALLEL_LOOP()
             else
-
                 if ((surface_tension .or. viscous) .or. chem_params%diffusion) then
                     $:GPU_PARALLEL_LOOP(private='[i, j, k, l]', collapse=3)
                     do l = 0, p
@@ -1567,7 +1565,6 @@ contains
                         $:END_GPU_PARALLEL_LOOP()
                     end if
                 else
-
                     $:GPU_PARALLEL_LOOP(private='[i, j, k, l]', collapse=3)
                     do l = 0, p
                         do k = 0, n
@@ -1584,7 +1581,6 @@ contains
                 end if
             end if
         else if (idir == 3) then ! z-direction
-
             if (surface_tension) then
                 $:GPU_PARALLEL_LOOP(private='[j, k, l]', collapse=3)
                 do l = 0, p
@@ -1695,7 +1691,6 @@ contains
                                            & recon_dir, is1, is2, is3)
                     end if
                 else
-
                     call s_${SCHEME}$ (v_vf(iv%beg:iv%end), vL_x(:,:,:, iv%beg:iv%end), vL_y(:,:,:,:), vL_z(:,:,:,:), vR_x(:,:,:, &
                                        & iv%beg:iv%end), vR_y(:,:,:,:), vR_z(:,:,:,:), recon_dir, is1, is2, is3)
                 end if
diff --git a/src/simulation/m_riemann_solvers.fpp b/src/simulation/m_riemann_solvers.fpp
@@ -2624,7 +2624,6 @@ contains
                                         R3V2Lbar = mom_sp_rs${XYZ}$_vf(j, k, l, 3)
                                         R3V2Rbar = mom_sp_rs${XYZ}$_vf(j + 1, k, l, 3)
                                     else
-
                                         PbwR3Lbar = 0._wp
                                         PbwR3Rbar = 0._wp
 
@@ -3940,7 +3939,6 @@ contains
 
             ! Population of Buffers in z-direction
         else
-
             if (bc_z%beg == BC_RIEMANN_EXTRAP) then ! Riemann state extrap. BC at beginning
                 $:GPU_PARALLEL_LOOP(collapse=3)
                 do i = 1, sys_size
@@ -4137,7 +4135,6 @@ contains
 
             ! Reshaping Inputted Data in z-direction
         else
-
             if (viscous .or. (surface_tension) .or. dummy) then
                 $:GPU_PARALLEL_LOOP(collapse=4)
                 do i = momxb, E_idx
diff --git a/src/simulation/m_start_up.fpp b/src/simulation/m_start_up.fpp
@@ -531,7 +531,6 @@ contains
                 if (ib) then
                     call s_initialize_mpi_data(q_cons_vf, ib_markers)
                 else
-
                     call s_initialize_mpi_data(q_cons_vf)
                 end if
 
diff --git a/src/simulation/m_viscous.fpp b/src/simulation/m_viscous.fpp
diff --git a/src/simulation/m_weno.fpp b/src/simulation/m_weno.fpp
