Testing if this causes the seg fault

Author: danieljvickers

src/simulation/m_compute_levelset.fpp

@@ -33,6 +33,8 @@ contains
 
         integer :: i, patch_id, patch_geometry
 
+        print *, "num_gps: ", num_gps
+
         !  3D Patch Geometries
         if (p > 0) then
 

src/simulation/m_ib_patches.fpp

@@ -1035,8 +1035,10 @@ contains
                 stl_bounding_boxes(patch_id, 3, 1:3) = [bbox%min(3), (bbox%min(3) + bbox%max(3))/2._wp, bbox%max(3)]
 
                 models(patch_id)%model = model
-                models(patch_id)%boundary_v = boundary_v
-                models(patch_id)%boundary_edge_count = boundary_edge_count
+                if (p == 0) then                                                                                                                       
+                    models(patch_id)%boundary_v = boundary_v                                                                                           
+                    models(patch_id)%boundary_edge_count = boundary_edge_count                                                                         
+                 end if  
             end if
         end do
 
@@ -1101,6 +1103,8 @@ contains
             end if
         end block
 
+        print *, "Finished reading models"
+
     end subroutine s_instantiate_STL_models
 
     !> Subroutine that computes a rotation matrix for converting to the rotating frame of the boundary