diff --git a/examples/2D_kelvin_helmholtz/case.py b/examples/2D_kelvin_helmholtz/case.py
index 3a12034696..e6873f8145 100644
--- a/examples/2D_kelvin_helmholtz/case.py
+++ b/examples/2D_kelvin_helmholtz/case.py
@@ -2,7 +2,6 @@
 import json
 import math
 
-
 eps = 1e-6
 time_end = 1.0
 time_save = time_end / 100.0
diff --git a/examples/2D_richtmyer_meshkov/case.py b/examples/2D_richtmyer_meshkov/case.py
index efb30f17f4..fed6529c55 100644
--- a/examples/2D_richtmyer_meshkov/case.py
+++ b/examples/2D_richtmyer_meshkov/case.py
@@ -2,7 +2,6 @@
 import json
 import math
 
-
 mu = 1.0e-4
 lambd = 1.0
 time_end = 15.0
diff --git a/examples/2D_viscous_shock_tube/case.py b/examples/2D_viscous_shock_tube/case.py
index e61f34cb97..7c760e696b 100644
--- a/examples/2D_viscous_shock_tube/case.py
+++ b/examples/2D_viscous_shock_tube/case.py
@@ -2,7 +2,6 @@
 import json
 import math
 
-
 mu = 1.0e-3
 time_end = 1.0
 time_save = time_end / 100.0
diff --git a/examples/nD_perfect_reactor/analyze.py b/examples/nD_perfect_reactor/analyze.py
index ab0c6c0a94..b2bd4e1fa8 100644
--- a/examples/nD_perfect_reactor/analyze.py
+++ b/examples/nD_perfect_reactor/analyze.py
@@ -6,7 +6,6 @@
 import mfc.viz
 from case import dt, Tend, SAVE_COUNT, sol
 
-
 case = mfc.viz.Case(".", dt)
 
 sns.set_theme(style=mfc.viz.generate_cpg_style())
diff --git a/examples/nD_perfect_reactor/export.py b/examples/nD_perfect_reactor/export.py
index f60f55e3fe..112622a32f 100644
--- a/examples/nD_perfect_reactor/export.py
+++ b/examples/nD_perfect_reactor/export.py
@@ -5,7 +5,6 @@
 import mfc.viz
 from case import dt, NS, Tend, SAVE_COUNT, sol
 
-
 case = mfc.viz.Case(".", dt)
 
 for name in tqdm(sol.species_names, desc="Loading Variables"):
diff --git a/src/simulation/m_ibm.fpp b/src/simulation/m_ibm.fpp
index 0b8ec1ab81..2c5d133060 100644
--- a/src/simulation/m_ibm.fpp
+++ b/src/simulation/m_ibm.fpp
@@ -1179,7 +1179,7 @@ contains
         elseif (patch_ib(ib_marker)%geometry == 8) then ! sphere
             patch_ib(ib_marker)%moment = 0.4*patch_ib(ib_marker)%mass*(patch_ib(ib_marker)%radius)**2
 
-        else ! we do not have an analytic moment of inertia calculation and need to approximate it directly
+        else ! we do not have an analytic moment of inertia calculation and need to approximate it directly via a sum
             count = 0
             moment = 0._wp
             cell_volume = (x_cc(1) - x_cc(0))*(y_cc(1) - y_cc(0)) ! computed without grid stretching. Update in the loop to perform with stretching