Merge

Author: Anand

src/simulation/m_ibm.fpp

@@ -18,6 +18,8 @@ module m_ibm
 
     use m_helper
 
+    use m_helper_basic         !< Functions to compare floating point numbers
+
     use m_constants
 
     implicit none
@@ -164,7 +166,7 @@ contains
         type(ghost_point) :: gp
         type(ghost_point) :: innerp
 
-        !$acc parallel loop gang vector private(physical_loc, dyn_pres, alpha_rho_IP, alpha_IP, pres_IP, vel_IP, vel_g, vel_norm_IP, r_IP, v_IP, pb_IP, mv_IP, nmom_IP, presb_IP, massv_IP, rho, gamma, pi_inf, Re_K, G_K, Gs, gp, innerp, norm, buf, j, k, l, q, coeff)
+        !$acc parallel loop gang vector private(physical_loc, dyn_pres, alpha_rho_IP, alpha_IP, pres_IP, vel_IP, vel_g, vel_norm_IP, r_IP, v_IP, pb_IP, mv_IP, nmom_IP, presb_IP, massv_IP, rho, gamma, pi_inf, Re_K, G_K, Gs, gp, innerp, norm, buf, j, k, l, q)
         do i = 1, num_gps
 
             gp = ghost_points(i)
@@ -224,7 +226,7 @@ contains
                                                                     alpha_rho_IP, Re_K)
                 end if
             end if
-            
+
             ! Calculate velocity of ghost cell
             if (gp%slip) then
                 norm(1:3) = levelset_norm%sf(gp%loc(1), gp%loc(2), gp%loc(3), gp%ib_patch_id, 1:3)
@@ -382,7 +384,7 @@ contains
                     bound = p
                 end if
 
-                if (norm(dim) == 0) then
+                if (f_approx_equal(norm(dim), 0._wp)) then
                     ghost_points(q)%ip_grid(dim) = ghost_points(q)%loc(dim)
                 else
                     if (norm(dim) > 0) then