@@ -86,10 +86,8 @@ contains
8686 call s_ib_cylinder(i, ib_markers_sf)
8787 call s_cylinder_levelset(i, levelset, levelset_norm)
8888 elseif (patch_ib(i)%geometry == 11 ) then
89- #ifdef MFC_PRE_PROCESS
9089 call s_ib_3D_airfoil(i, ib_markers_sf)
9190 call s_3D_airfoil_levelset(i, levelset, levelset_norm)
92- #endif
9391 ! STL+ IBM patch
9492 elseif (patch_ib(i)%geometry == 12 ) then
9593 call s_ib_model(i, ib_markers_sf, levelset, levelset_norm)
@@ -110,10 +108,8 @@ contains
110108 call s_ib_rectangle(i, ib_markers_sf)
111109 call s_rectangle_levelset(i, levelset, levelset_norm)
112110 elseif (patch_ib(i)%geometry == 4 ) then
113- #ifdef MFC_PRE_PROCESS
114111 call s_ib_airfoil(i, ib_markers_sf)
115112 call s_airfoil_levelset(i, levelset, levelset_norm)
116- #endif
117113 ! STL+ IBM patch
118114 elseif (patch_ib(i)%geometry == 5 ) then
119115 call s_ib_model(i, ib_markers_sf, levelset, levelset_norm)
@@ -277,7 +273,7 @@ contains
277273 end if
278274 if (y_act >= y0) then
279275 k = 1
280- do while (airfoil_grid_u(k)%x < x_act)
276+ do while (airfoil_grid_u(k)%x < x_act .and. k < = Np )
281277 k = k + 1
282278 end do
283279 if (f_approx_equal(airfoil_grid_u(k)%x, x_act)) then
@@ -771,7 +767,6 @@ contains
771767 !> The STL patch is a 2/3D geometry that is imported from an STL file.
772768 !! @param patch_id is the patch identifier
773769 !! @param ib_markers_sf Array to track patch ids
774- !! @param ib True if this patch is an immersed boundary
775770 !! @param STL_levelset STL levelset
776771 !! @param STL_levelset_norm STL levelset normals
777772 subroutine s_ib_model(patch_id, ib_markers_sf, STL_levelset, STL_levelset_norm)
0 commit comments