Minor fixes and formatting

Author: Mark Zhang

examples/2D_IGR_forward_facing_step/case.py

@@ -68,23 +68,23 @@
             "patch_icpp(1)%pres": p0,
             "patch_icpp(1)%alpha_rho(1)": rho0,
             "patch_icpp(1)%alpha(1)": 1.0,
-            # Patch: Slip rectangle with rouded corner
+            # Patch: Slip rectangle with rounded corner
             "patch_ib(1)%geometry": 3,
             "patch_ib(1)%x_centroid": (9 + rc / 2) * h,
             "patch_ib(1)%y_centroid": 0.5 * h,
             "patch_ib(1)%length_x": (12 - rc) * h,
             "patch_ib(1)%length_y": h,
             "patch_ib(1)%slip": "T",
             "patch_ib(2)%geometry": 3,
-            "patch_ib(2)%x_centroid": (3 + rc/2) * h,
-            "patch_ib(2)%y_centroid": (0.5 - rc/2) * h,
+            "patch_ib(2)%x_centroid": (3 + rc / 2) * h,
+            "patch_ib(2)%y_centroid": (0.5 - rc / 2) * h,
             "patch_ib(2)%length_x": rc * h,
             "patch_ib(2)%length_y": (1 - rc) * h,
             "patch_ib(2)%slip": "T",
             "patch_ib(3)%geometry": 2,
-            "patch_ib(3)%x_centroid": (3 + rc)*h,
-            "patch_ib(3)%y_centroid": (1 - rc)*h,
-            "patch_ib(3)%radius" : rc * h,
+            "patch_ib(3)%x_centroid": (3 + rc) * h,
+            "patch_ib(3)%y_centroid": (1 - rc) * h,
+            "patch_ib(3)%radius": rc * h,
             "patch_ib(3)%slip": "T",
             # Fluids Physical Parameters
             "fluid_pp(1)%gamma": 1.0 / (gam_a - 1.0),

src/simulation/m_ibm.fpp

@@ -970,6 +970,9 @@ contains
                     if (igr) then
                         ! For IGR, we will need to perform operations on
                         ! the conservative variables instead
+                        if (.not. present(q_cons_vf)) then
+                            stop 'q_cons_vf is required inside s_interpolate_image_point when running with IGR'
+                        end if
                         alpha_sum = 0._wp
                         dyn_pres = 0._wp
                         if (num_fluids == 1) then
@@ -1030,7 +1033,7 @@ contains
                                 alpha_sum = alpha_sum + coeff*q_cons_vf(E_idx + l)%sf(i, j, k)
                             end do
                             alpha_rho_IP(num_fluids) = alpha_rho_IP(num_fluids) + coeff*q_cons_vf(num_fluids)%sf(i, j, k)
-                            alpha_IP(num_fluids) = alpha_IP(num_fluids) + coeff*(1._wp - alpha_sum)
+                            alpha_IP(num_fluids) = alpha_IP(num_fluids) + (1._wp - alpha_sum)
                         end if
                     else
                         pres_IP = pres_IP + coeff* &

src/simulation/m_igr.fpp

@@ -15,7 +15,7 @@ module m_igr
 
     use m_boundary_common
 
-    use m_ibm, only: ib, ib_markers, s_interpolate_sigma_igr
+    use m_ibm, only: ib_markers, s_interpolate_sigma_igr
 
     implicit none