Resolved some prohibit statements with airfoils from an untested example

Author: danieljvickers

examples/2D_ibm_airfoil/case.py

@@ -90,10 +90,11 @@
             "patch_ib(1)%geometry": 4,
             "patch_ib(1)%x_centroid": 1.0e-03,
             "patch_ib(1)%y_centroid": 3.0e-03,
-            "patch_ib(1)%c": 1.0e-03,
-            "patch_ib(1)%t": 0.15,
-            "patch_ib(1)%p": 0.4,
-            "patch_ib(1)%m": 0.02,
+            "patch_ib(1)%airfoil_id": 1,
+            "ib_airfoil(1)%c": 1.0e-03,
+            "ib_airfoil(1)%t": 0.15,
+            "ib_airfoil(1)%p": 0.4,
+            "ib_airfoil(1)%m": 0.02,
             "patch_ib(1)%angles(3)": -0.5235987756,  # 30 degrees clockwise rotation, in radians
             # Fluids Physical Parameters
             # Use the same stiffness as the air bubble

src/pre_process/m_global_parameters.fpp

@@ -93,10 +93,9 @@ module m_global_parameters
 
     !> @name Immersed Boundaries
     !> @{
-    type(ib_patch_parameters), dimension(num_ib_patches_max_namelist) :: patch_ib        !< Immersed boundary patch parameters
-    type(ib_airfoil_parameters), dimension(num_ib_airfoils_max)       :: ib_airfoil      !< Per-airfoil NACA parameters
-    type(ib_stl_parameters), dimension(num_stl_models_max)            :: stl_models      !< Per-STL model parameters (namelist)
-    integer                                                           :: num_stl_models  !< Number of STL models
+    type(ib_patch_parameters), dimension(num_ib_patches_max_namelist) :: patch_ib    !< Immersed boundary patch parameters
+    type(ib_airfoil_parameters), dimension(num_ib_airfoils_max)       :: ib_airfoil  !< Per-airfoil NACA parameters
+    type(ib_stl_parameters), dimension(num_stl_models_max)            :: stl_models  !< Per-STL model parameters (namelist)
     !> @}
 
     !> @name Non-polytropic bubble gas compression

src/simulation/m_global_parameters.fpp

@@ -268,7 +268,6 @@ module m_global_parameters
     type(ib_airfoil_parameters), dimension(num_ib_airfoils_max) :: ib_airfoil  !< Per-airfoil NACA user inputs (namelist)
     type(ib_airfoil_grid), dimension(num_ib_airfoils_max) :: ib_airfoil_grids  !< Per-airfoil computed surface grids
     type(ib_stl_parameters), dimension(num_stl_models_max) :: stl_models  !< Per-STL model parameters (namelist)
-    integer :: num_stl_models  !< Number of STL models
 
     $:GPU_DECLARE(create='[ib, num_ibs, patch_ib, ib_airfoil, ib_airfoil_grids]')
     $:GPU_DECLARE(create='[ib_coefficient_of_friction]')

src/simulation/m_ibm.fpp

@@ -1529,6 +1529,8 @@ contains
     !> Finalize the IBM module
     impure subroutine s_finalize_ibm_module()
 
+        integer :: i
+
         @:DEALLOCATE(ib_markers%sf)
         @:DEALLOCATE(ib_gbl_idx_lookup)
         do i = 1, num_ib_airfoils_max

src/simulation/m_start_up.fpp

@@ -819,7 +819,6 @@ contains
         if (bubbles_euler) call s_initialize_bubbles_EE_module()
         if (ib) then
             call s_initialize_ibm_module()
-            call s_initialize_ib_airfoils()
         end if
         if (qbmm) call s_initialize_qbmm_module()
 
@@ -885,6 +884,7 @@ contains
                     call s_generate_particle_beds(particle_bed_ibs)
                 end if
                 call s_instantiate_STL_models()
+                call s_initialize_ib_airfoils()
                 call s_reduce_ib_patch_array(particle_bed_ibs)
                 deallocate (particle_bed_ibs)
             end block