Small updates to code for strictness and updating case file to be well posed

Author: Daniel Vickers

examples/3D_mibm_sphere_head_on_collision/case.py

@@ -92,9 +92,9 @@
             # Patch: Constant Tube filled with air
             # Specify the cylindrical air tube grid geometry
             "patch_icpp(1)%geometry": 9,
-            "patch_icpp(1)%x_centroid": 0.5 * domain_size,
+            "patch_icpp(1)%x_centroid": 0.0,
             "patch_icpp(1)%y_centroid": 0.5 * domain_size,
-            "patch_icpp(1)%z_centroid": 0.5 * domain_size,
+            "patch_icpp(1)%z_centroid": 0.0,
             "patch_icpp(1)%length_x": domain_size,
             "patch_icpp(1)%length_y": domain_size,
             "patch_icpp(1)%length_z": domain_size,

src/simulation/m_collisions.fpp

@@ -222,8 +222,6 @@ contains
         end do
         $:END_GPU_PARALLEL_LOOP()
 
-        print *, forces(1,1:3)
-
     end subroutine s_apply_wall_collision_forces_soft_sphere
 
     !> uses ghost-point/image-point information to determine if it is possible if two IBs are colliding, effectively an optimized

src/simulation/m_ibm.fpp

@@ -915,7 +915,7 @@ contains
 
         $:GPU_PARALLEL_LOOP(private='[ib_idx, fluid_idx, radial_vector, local_force_contribution, cell_volume, &
                             & local_torque_contribution, dynamic_viscosity, viscous_stress_div, viscous_stress_div_1, &
-                            & viscous_stress_div_2, dx, dy, dz]', copy='[forces, torques]', copyin='[ib_markers, patch_ib, &
+                            & viscous_stress_div_2, dx, dy, dz]', copy='[forces, torques]', copyin='[patch_ib, &
                             & dynamic_viscosities]', collapse=3)
         do i = 0, m
             do j = 0, n

src/simulation/m_time_steppers.fpp

@@ -771,6 +771,7 @@ contains
             end if
         end do
 
+        $:GPU_UPDATE(device='[patch_ib]')
         call s_update_mib(num_ibs)
 
         call nvtxEndRange