m_check_ib_patches.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
!>
!! @file
!! @brief Contains module m_check_ib_patches
 
!> @brief Validates geometry parameters and constraints for immersed boundary patches
 
# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 206 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 231 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 242 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 284 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 340 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 347 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 353 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 359 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 365 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 371 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 377 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
! New line at end of file is required for FYPP
# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 145 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
 
# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
! GPU parallel region (scalar reductions, maxval/minval)
# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! GPU parallel loop over threads (most common GPU macro)
# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Required closing for GPU_PARALLEL_LOOP
# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Mark routine for device compilation
# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Declare device-resident data
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Inner loop within a GPU parallel region
# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Scoped GPU data region
# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Host code with device pointers (for MPI with GPU buffers)
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Allocate device memory (unscoped)
# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Free device memory
# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Atomic operation on device
# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! End atomic capture block
# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Copy data between host and device
# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Synchronization barrier
# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Import GPU library module (openacc or omp_lib)
# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for AMD compiler
# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-Cray compilers
# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for Cray compiler
# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-NVIDIA compilers
# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
! example see misc/nvidia_uvm/bind.sh. NVIDIA unified memory page placement hint
# 57 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 77 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 85 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 109 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 125 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 150 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 163 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 8 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp" 2
 
module m_check_ib_patches
 
    use m_derived_types
    use m_global_parameters
    use m_mpi_proxy
    use m_data_output
#ifdef MFC_MPI
    use mpi  !< message passing interface (mpi) module
#endif
 
    use m_compile_specific
    use m_helper_basic
    use m_helper
 
    implicit none
 
    private
    public :: s_check_ib_patches
 
    character(len=10) :: istr
 
contains
 
    !> Validate the geometry parameters of all active and inactive immersed boundary patches.
    impure subroutine s_check_ib_patches
 
        integer :: i
 
        do i = 1, num_patches_max
            if (i <= num_ibs) then
                call s_int_to_str(i, istr)
    if (patch_ib(i)%geometry == dflt_int) then
# 40 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort("patch_ib(i)%geometry == dflt_int", "IB patch undefined. patch_ib("//trim(istr)//")%geometry must be set.")
# 40 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
 
                ! Constraints on the geometric initial condition patch parameters
                if (patch_ib(i)%geometry == 2) then
                    call s_check_circle_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 3) then
                    call s_check_rectangle_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 8) then
                    call s_check_sphere_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 9) then
                    call s_check_cuboid_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 4) then
                    call s_check_airfoil_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 11) then
                    call s_check_3d_airfoil_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 10) then
                    call s_check_cylinder_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 5 .or. patch_ib(i)%geometry == 12) then
                    call s_check_model_ib_patch_geometry(i)
                else if (patch_ib(i)%geometry == 6) then
                    call s_check_ellipse_ib_patch_geometry(i)
                else
                    call s_prohibit_abort("Invalid IB patch", &
                                          & "patch_ib(" // trim(istr) // ")%geometry must be " // "2-4, 8-10, 11 or 12.")
                end if
            else
    if (patch_ib(i)%geometry /= dflt_int) then
# 66 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort("patch_ib(i)%geometry /= dflt_int", "Inactive IB patch defined. " // "patch_ib(" // trim(istr)  // ")%geometry must not be set for inactive patches.")
# 66 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 69 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
                call s_check_inactive_ib_patch_geometry(i)
            end if
        end do
 
    end subroutine s_check_ib_patches
 
    !> Verify that the geometric parameters of the circle patch have been consistently inputted.
 
    impure subroutine s_check_circle_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p > 0 .or. patch_ib(patch_id)%radius <= 0._wp .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid)) then
# 83 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or."n == 0  p > 0  patch_ib(patch_id)%radius <= 0._wp  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)", 'in circle IB patch ' // trim(istr))
# 83 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 85 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_circle_ib_patch_geometry
 
    !> Verify that the geometric parameters of the ellipse patch have been consistently inputted.
 
    impure subroutine s_check_ellipse_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p > 0 .or. patch_ib(patch_id)%length_x <= 0._wp .or. patch_ib(patch_id)%length_y <= 0._wp  .or. f_is_default(patch_ib(patch_id)%x_centroid) .or. f_is_default(patch_ib(patch_id)%y_centroid)) then
# 96 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or."n == 0  p > 0  patch_ib(patch_id)%length_x <= 0._wp  patch_ib(patch_id)%length_y <= 0._wp   f_is_default(patch_ib(patch_id)%x_centroid)  f_is_default(patch_ib(patch_id)%y_centroid)", 'in ellipse IB patch ' // trim(istr))
# 96 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 99 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_ellipse_ib_patch_geometry
 
    !> Verify that the geometric parameters of the airfoil patch have been consistently inputted.
 
    impure subroutine s_check_airfoil_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p > 0 .or. patch_ib(patch_id)%airfoil_id <= 0  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%c <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%p <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%t <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%m <= 0._wp .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid)) then
# 110 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or..or..or..or."n == 0  p > 0  patch_ib(patch_id)%airfoil_id <= 0   ib_airfoil(patch_ib(patch_id)%airfoil_id)%c <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%p <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%t <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%m <= 0._wp  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)", 'in airfoil IB patch ' // trim(istr))
# 110 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 116 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_airfoil_ib_patch_geometry
 
    !> Verify that the geometric parameters of the 3D airfoil patch have been consistently inputted.
 
    impure subroutine s_check_3d_airfoil_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p == 0 .or. patch_ib(patch_id)%airfoil_id <= 0  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%c <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%p <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%t <= 0._wp  .or. ib_airfoil(patch_ib(patch_id)%airfoil_id)%m <= 0._wp .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid) .or. f_is_default(patch_ib(patch_id)%z_centroid)  .or. f_is_default(patch_ib(patch_id)%length_z)) then
# 127 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or..or..or..or..or..or."n == 0  p == 0  patch_ib(patch_id)%airfoil_id <= 0   ib_airfoil(patch_ib(patch_id)%airfoil_id)%c <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%p <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%t <= 0._wp   ib_airfoil(patch_ib(patch_id)%airfoil_id)%m <= 0._wp  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)  f_is_default(patch_ib(patch_id)%z_centroid)   f_is_default(patch_ib(patch_id)%length_z)", 'in 3d airfoil IB patch ' // trim(istr))
# 127 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 134 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_3d_airfoil_ib_patch_geometry
 
    !> Verify that the geometric parameters of the rectangle patch have been consistently inputted.
 
    impure subroutine s_check_rectangle_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p > 0 .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid) .or. patch_ib(patch_id)%length_x <= 0._wp  .or. patch_ib(patch_id)%length_y <= 0._wp) then
# 145 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or."n == 0  p > 0  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)  patch_ib(patch_id)%length_x <= 0._wp   patch_ib(patch_id)%length_y <= 0._wp", 'in rectangle IB patch ' // trim(istr))
# 145 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 148 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_rectangle_ib_patch_geometry
 
    !> Verify that the geometric parameters of the sphere patch have been consistently inputted.
 
    impure subroutine s_check_sphere_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p == 0 .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid) .or. f_is_default(patch_ib(patch_id)%z_centroid)  .or. patch_ib(patch_id)%radius <= 0._wp) then
# 159 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or."n == 0  p == 0  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)  f_is_default(patch_ib(patch_id)%z_centroid)   patch_ib(patch_id)%radius <= 0._wp", 'in sphere IB patch ' // trim(istr))
# 159 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 162 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_sphere_ib_patch_geometry
 
    !> Verify that the geometric parameters of the cuboid patch have been consistently inputted.
 
    impure subroutine s_check_cuboid_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (n == 0 .or. p == 0 .or. f_is_default(patch_ib(patch_id)%x_centroid)  .or. f_is_default(patch_ib(patch_id)%y_centroid) .or. f_is_default(patch_ib(patch_id)%z_centroid)  .or. patch_ib(patch_id)%length_x <= 0._wp .or. patch_ib(patch_id)%length_y <= 0._wp  .or. patch_ib(patch_id)%length_z <= 0._wp) then
# 173 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..or..or..or."n == 0  p == 0  f_is_default(patch_ib(patch_id)%x_centroid)   f_is_default(patch_ib(patch_id)%y_centroid)  f_is_default(patch_ib(patch_id)%z_centroid)   patch_ib(patch_id)%length_x <= 0._wp  patch_ib(patch_id)%length_y <= 0._wp   patch_ib(patch_id)%length_z <= 0._wp", 'in cuboid IB patch ' // trim(istr))
# 173 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 177 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_cuboid_ib_patch_geometry
 
    !> Verify that the geometric parameters of the cylinder patch have been consistently inputted.
 
    impure subroutine s_check_cylinder_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if (p == 0 .or. f_is_default(patch_ib(patch_id)%x_centroid) .or. f_is_default(patch_ib(patch_id)%y_centroid)  .or. f_is_default(patch_ib(patch_id)%z_centroid) .or. (patch_ib(patch_id)%length_x <= 0._wp  .and. patch_ib(patch_id)%length_y <= 0._wp .and. patch_ib(patch_id)%length_z <= 0._wp)  .or. patch_ib(patch_id)%radius <= 0._wp) then
# 188 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or..or..and..and..or."p == 0  f_is_default(patch_ib(patch_id)%x_centroid)  f_is_default(patch_ib(patch_id)%y_centroid)   f_is_default(patch_ib(patch_id)%z_centroid)  (patch_ib(patch_id)%length_x <= 0._wp   patch_ib(patch_id)%length_y <= 0._wp  patch_ib(patch_id)%length_z <= 0._wp)   patch_ib(patch_id)%radius <= 0._wp", 'in cylinder IB patch ' // trim(istr))
# 188 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 192 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    if ((patch_ib(patch_id)%length_x > 0._wp .and. ((.not. f_is_default(patch_ib(patch_id)%length_y))  .or. (.not. f_is_default(patch_ib(patch_id)%length_z)))) .or. (patch_ib(patch_id)%length_y > 0._wp  .and. ((.not. f_is_default(patch_ib(patch_id)%length_x))  .or. (.not. f_is_default(patch_ib(patch_id)%length_z)))) .or. (patch_ib(patch_id)%length_z > 0._wp  .and. ((.not. f_is_default(patch_ib(patch_id)%length_x))  .or. (.not. f_is_default(patch_ib(patch_id)%length_y))))) then
# 193 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.and..not..or..not..or..and..not..or..not..or..and..not..or..not."(patch_ib(patch_id)%length_x > 0._wp  (( f_is_default(patch_ib(patch_id)%length_y))   ( f_is_default(patch_ib(patch_id)%length_z))))  (patch_ib(patch_id)%length_y > 0._wp   (( f_is_default(patch_ib(patch_id)%length_x))   ( f_is_default(patch_ib(patch_id)%length_z))))  (patch_ib(patch_id)%length_z > 0._wp   (( f_is_default(patch_ib(patch_id)%length_x))   ( f_is_default(patch_ib(patch_id)%length_y))))", 'in cylinder IB patch ' // trim(istr))
# 193 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 199 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_cylinder_ib_patch_geometry
 
    !> Verify that the geometric parameters of the model patch have been consistently inputted.
 
    impure subroutine s_check_model_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
        integer             :: mid
        character(len=10)   :: midstr
 
        call s_int_to_str(patch_id, istr)
 
        mid = patch_ib(patch_id)%model_id
        call s_int_to_str(mid, midstr)
 
    if (mid <= 0 .or. mid > num_stl_models) then
# 215 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or."mid <= 0  mid > num_stl_models", 'patch_ib('//trim(istr)//')%model_id='//trim(midstr)//' must be in [1, num_stl_models]')
# 215 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 217 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    if (stl_models(mid)%model_filepath == dflt_char) then
# 218 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort("stl_models(mid)%model_filepath == dflt_char", 'Empty model file path for stl_models('//trim(midstr)//')')
# 218 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
 
    if (stl_models(mid)%model_scale(1) <= 0._wp .or. stl_models(mid)%model_scale(2) <= 0._wp  .or. stl_models(mid)%model_scale(3) <= 0._wp) then
# 220 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or."stl_models(mid)%model_scale(1) <= 0._wp  stl_models(mid)%model_scale(2) <= 0._wp   stl_models(mid)%model_scale(3) <= 0._wp", 'Negative scale in stl_models(' // trim(midstr) // ')')
# 220 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 222 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_model_ib_patch_geometry
 
    !> Verify that inactive IB patch geometry parameters remain at defaults
    impure subroutine s_check_inactive_ib_patch_geometry(patch_id)
 
        integer, intent(in) :: patch_id
 
        call s_int_to_str(patch_id, istr)
 
    if ((.not. f_is_default(patch_ib(patch_id)%x_centroid)) .or. (.not. f_is_default(patch_ib(patch_id)%y_centroid))  .or. (.not. f_is_default(patch_ib(patch_id)%z_centroid))  .or. (.not. f_is_default(patch_ib(patch_id)%length_x)) .or. (.not. f_is_default(patch_ib(patch_id)%length_y))  .or. (.not. f_is_default(patch_ib(patch_id)%length_z)) .or. (.not. f_is_default(patch_ib(patch_id)%radius))) then
# 232 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.not..or..not..or..not..or..not..or..not..or..not..or..not."( f_is_default(patch_ib(patch_id)%x_centroid))  ( f_is_default(patch_ib(patch_id)%y_centroid))   ( f_is_default(patch_ib(patch_id)%z_centroid))   ( f_is_default(patch_ib(patch_id)%length_x))  ( f_is_default(patch_ib(patch_id)%length_y))   ( f_is_default(patch_ib(patch_id)%length_z))  ( f_is_default(patch_ib(patch_id)%radius))", 'in inactive IB patch ' // trim(istr))
# 232 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 237 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    end subroutine s_check_inactive_ib_patch_geometry
 
end module m_check_ib_patches