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"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/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"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/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"
 
# 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"
 
# 126 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 197 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 211 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 236 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 247 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 249 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 260 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 310 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 320 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 356 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 366 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 373 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 379 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 385 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 391 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 397 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 403 "/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"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/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.
# 52 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 72 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 80 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 104 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 120 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 145 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 158 "/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) then
# 83 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or."n == 0  p > 0  patch_ib(patch_id)%radius <= 0._wp", 'in circle IB patch ' // trim(istr))
# 83 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
 
    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) then
# 95 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or."n == 0  p > 0  patch_ib(patch_id)%length_x <= 0._wp  patch_ib(patch_id)%length_y <= 0._wp", 'in ellipse IB patch ' // trim(istr))
# 95 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 97 "/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) then
# 108 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.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", 'in airfoil IB patch ' // trim(istr))
# 108 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 113 "/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)%length_z)) then
# 124 "/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  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)%length_z)", 'in 3d airfoil IB patch ' // trim(istr))
# 124 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 130 "/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. patch_ib(patch_id)%length_x <= 0._wp .or. patch_ib(patch_id)%length_y <= 0._wp) then
# 141 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or..or."n == 0  p > 0  patch_ib(patch_id)%length_x <= 0._wp  patch_ib(patch_id)%length_y <= 0._wp", 'in rectangle IB patch ' // trim(istr))
# 141 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 143 "/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. patch_ib(patch_id)%radius <= 0._wp) then
# 154 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..or."n == 0  p == 0  patch_ib(patch_id)%radius <= 0._wp", 'in sphere IB patch ' // trim(istr))
# 154 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
 
    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. 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
# 166 "/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)%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))
# 166 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 168 "/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. (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
# 179 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.or..and..and..or."p == 0  (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))
# 179 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 182 "/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
# 183 "/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))
# 183 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 189 "/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
# 205 "/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]')
# 205 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 207 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
 
    if (stl_models(mid)%model_filepath == dflt_char) then
# 208 "/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)//')')
# 208 "/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
# 210 "/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) // ')')
# 210 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 212 "/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)%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
# 222 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
        call s_prohibit_abort(.not..or..not..or..not..or..not."( 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))
# 222 "/home/runner/work/MFC/MFC/src/pre_process/m_check_ib_patches.fpp"
    end if
# 225 "/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