fixed formatting issue

Author: Daniel Vickers

src/common/m_constants.fpp

@@ -24,7 +24,7 @@ module m_constants
     integer, parameter  :: num_fluids_max = 10                 !< Maximum number of fluids in the simulation
     integer, parameter  :: num_probes_max = 10                 !< Maximum number of flow probes in the simulation
     integer, parameter  :: num_patches_max = 10                !< Maximum number of IC patches
-    integer, parameter  :: num_ib_patches_max = 300000         !< Maximum number of immersed boundary patches (patch_ib)
+    integer, parameter  :: num_ib_patches_max = 2050000        !< Maximum number of immersed boundary patches (patch_ib)
     integer, parameter  :: num_local_ibs_max = 2000            !< Maximum number of immersed boundary patches (patch_ib)
     integer, parameter  :: num_bc_patches_max = 10             !< Maximum number of boundary condition patches
     integer, parameter  :: max_2d_fourier_modes = 10           !< Max Fourier mode index for 2D modal patch (geometry 13)

src/simulation/m_ibm.fpp

@@ -998,8 +998,7 @@ contains
                                     call s_compute_viscous_stress_tensor(viscous_stress_div_2, q_prim_vf, dynamic_viscosity, i, &
                                                                          & j, k + 1)
                                     viscous_stress_div(3,1:3) = (viscous_stress_div_2(3,1:3) - viscous_stress_div_1(3, &
-                                                       & 1:3))/(2._wp*dz) &
-                                                       &  ! get z derivative of the third-row of viscous stress tensor
+                                                       & 1:3))/(2._wp*dz)
                                     local_force_contribution(1:3) = local_force_contribution(1:3) + viscous_stress_div(3, &
                                                              & 1:3)  ! add the z components of the divergence to the force
                                 end if