m_checker.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
!>
!!@file
!!@brief Contains module m_checker
 
# 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"
 
# 167 "/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"
 
# 167 "/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"
 
# 167 "/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.
# 55 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 75 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 83 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 107 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 123 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 148 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 161 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 6 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/case.fpp" 1
! This file exists so that Fypp can be run without generating case.fpp files for
! each target. This is useful when generating documentation, for example. This
! should also let MFC be built with CMake directly, without invoking mfc.sh.
 
! For pre-process.
# 8 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
 
! For moving immersed boundaries in simulation
# 12 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
# 7 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp" 2
 
!> @brief Validates simulation input parameters for consistency and supported configurations
module m_checker
 
    use m_global_parameters
    use m_mpi_proxy
    use m_helper
    use m_helper_basic
    use m_constants, only: model_eqns_5eq, riemann_solver_hll, riemann_solver_hllc, riemann_solver_hlld, recon_type_weno, &
        & recon_type_muscl, muscl_order_first_order, riemann_solver_lax_friedrichs, wave_speeds_pressure
 
    implicit none
 
    private; public :: s_check_inputs
 
contains
 
    !> Checks compatibility of parameters in the input file. Used by the simulation stage
    impure subroutine s_check_inputs
 
        call s_check_inputs_compilers
 
        if (igr) then
            call s_check_inputs_nvidia_uvm
        else
            if (recon_type == recon_type_weno) then
                call s_check_inputs_weno
            else if (recon_type == recon_type_muscl) then
                call s_check_inputs_muscl
            end if
        end if
 
        call s_check_inputs_time_stepping
 
        call s_check_inputs_hll_non_conservative
 
        call s_check_inputs_hypo_branch
 
    if (chemistry .and. chem_params%reaction_substeps < 0) then
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."chemistry  chem_params%reaction_substeps < 0", "chem_params%reaction_substeps must be >= 0 (0 = reaction source in the flow RHS; > 0 = operator-split sub-stepping)")
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 47 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    if (chemistry .and. igr .and. chem_params%reaction_substeps > 0) then
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."chemistry  igr  chem_params%reaction_substeps > 0", "operator-split reaction sub-stepping (reaction_substeps > 0) is not supported with igr: the reactor reads the post-flow (rho, e, T) state, which the IGR update path does not guarantee")
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    if (chemistry .and. chem_params%adap_substeps .and. chem_params%reaction_substeps < 1) then
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."chemistry  chem_params%adap_substeps  chem_params%reaction_substeps < 1", "chem_params%adap_substeps requires reaction_substeps >= 1 (the operator-split floor)")
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 53 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    if (chemistry .and. chem_params%adap_substeps  .and. chem_params%reaction_substeps_max < chem_params%reaction_substeps) then
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."chemistry  chem_params%adap_substeps   chem_params%reaction_substeps_max < chem_params%reaction_substeps", "chem_params%reaction_substeps_max must be >= reaction_substeps when adap_substeps = T")
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 57 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
        ! Chemistry with Euler bubbles is not currently supported: the IBM image-point
        ! interpolation selects the bubbles/QBMM branch before the chemistry branch, so the
        ! species state (Ys_IP) is not carried when both are enabled. Disallow until implemented.
    if (chemistry .and. (bubbles_euler .or. qbmm)) then
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..or."chemistry  (bubbles_euler  qbmm)", "chemistry is not currently supported with Euler bubbles (bubbles_euler/qbmm)")
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    if (ib_state_wrt .and. .not. ib) then
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."ib_state_wrt   ib", "ib_state_wrt requires ib to be enabled")
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (many_ib_patch_parallelism .and. .not. ib) then
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."many_ib_patch_parallelism   ib", "many_ib_patch_parallelism requires ib to be enabled")
# 65 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
 
    if (bf_spatial_support .and. (n == 0 .or. p /= 0)) then
# 67 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..or."bf_spatial_support  (n == 0  p /= 0)", "bf_spatial_support is implemented for 2D only (it forces mom%beg and mom%beg+1)")
# 67 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
        ! Condensed-phase reactive burn assumes exactly two fluids (reactant=1, product=2) that share the
        ! stiffened-gas EOS and differ only in qv; violating these silently corrupts the mass/energy balance.
    if (reactive_burn .and. num_fluids /= 2) then
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  num_fluids /= 2", "reactive_burn requires num_fluids = 2 (reactant then product)")
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (reactive_burn .and. .not. f_approx_equal(fluid_pp(1)%gamma, fluid_pp(2)%gamma)) then
# 73 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."reactive_burn   f_approx_equal(fluid_pp(1)%gamma, fluid_pp(2)%gamma)", "reactive_burn requires fluid_pp(1)%gamma == fluid_pp(2)%gamma (reactant and product share the EOS)")
# 73 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 75 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. .not. f_approx_equal(fluid_pp(1)%pi_inf, fluid_pp(2)%pi_inf)) then
# 75 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."reactive_burn   f_approx_equal(fluid_pp(1)%pi_inf, fluid_pp(2)%pi_inf)", "reactive_burn requires fluid_pp(1)%pi_inf == fluid_pp(2)%pi_inf (reactant and product share the EOS)")
# 75 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 77 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. fluid_pp(1)%qv <= fluid_pp(2)%qv) then
# 77 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  fluid_pp(1)%qv <= fluid_pp(2)%qv", "reactive_burn requires fluid_pp(1)%qv > fluid_pp(2)%qv (reactant releases energy on conversion to product)")
# 77 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 79 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. rburn%pref <= 0._wp) then
# 79 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  rburn%pref <= 0._wp", "reactive_burn requires rburn%pref > 0 (it normalizes the pressure drive (p - rburn%pign)/rburn%pref and is used as a divisor)")
# 79 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        ! The rate uses rburn%k, rburn%pign, rburn%n directly; each defaults to the sentinel dflt_real,
        ! so an unset value silently produces spurious ignition (pign), NaN via drive**n (n), or a
        ! backward reaction (k). Require each to be set to a physical value.
    if (reactive_burn .and. rburn%k <= 0._wp) then
# 84 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  rburn%k <= 0._wp", "reactive_burn requires rburn%k > 0 (rate coefficient [1/s]; unset defaults to a negative sentinel that runs the reaction backward)")
# 84 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. f_is_default(rburn%pign)) then
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  f_is_default(rburn%pign)", "reactive_burn requires rburn%pign to be set (ignition pressure threshold [Pa]; unset defaults to a negative sentinel, so the reactant ignites everywhere from t = 0)")
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. rburn%n < 0._wp) then
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  rburn%n < 0._wp", "reactive_burn requires rburn%n >= 0 (pressure-drive exponent; unset defaults to a negative sentinel, so drive**n overflows to Inf and the field goes NaN)")
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. model_eqns /= 2 .and. model_eqns /= 3) then
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."reactive_burn  model_eqns /= 2  model_eqns /= 3", "reactive_burn requires model_eqns = 2 or 3 (the 5-equation pressure-equilibrium or 6-equation multi-fluid model)")
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 92 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. rburn%ta < 0._wp) then
# 92 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."reactive_burn  rburn%ta < 0._wp", "reactive_burn requires rburn%ta >= 0 (activation temperature [K]; 0 disables the Arrhenius factor)")
# 92 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (reactive_burn .and. rburn%ta > 0._wp .and. fluid_pp(1)%cv <= 0._wp) then
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."reactive_burn  rburn%ta > 0._wp  fluid_pp(1)%cv <= 0._wp", "reactive_burn with rburn%ta > 0 requires fluid_pp(1)%cv > 0 (the reactant temperature T = (p + pi_inf)/((gamma - 1) cv rho) needs a physical heat capacity; cv = 0 silently disables the Arrhenius factor)")
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 96 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
        if (ib .and. chemistry) then
            call s_check_inputs_ib_injection
        end if
 
        if (num_particle_clouds > 0) then
            call s_check_inputs_particle_clouds
        end if
 
        if (synthetic_turbulence) then
            call s_check_inputs_synthetic_turbulence
        end if
 
    end subroutine s_check_inputs
 
    !> Checks constraints on compiler options
    impure subroutine s_check_inputs_compilers
 
#if !defined(MFC_OpenACC) && !(defined(__PGI) || defined(_CRAYFTN))
    if (rdma_mpi) then
# 115 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("rdma_mpi", "Unsupported value of rdma_mpi for the current compiler")
# 115 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
#endif
 
    end subroutine s_check_inputs_compilers
 
    !> Checks constraints on WENO scheme parameters
    impure subroutine s_check_inputs_weno
 
        character(len=5) :: numstr  !< for int to string conversion
 
        call s_int_to_str(num_stcls_min*weno_order, numstr)
    if (m + 1 < num_stcls_min*weno_order) then
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("m + 1 < num_stcls_min*weno_order", "m must be greater than or equal to (num_stcls_min*weno_order - 1), whose value is " // trim(numstr))
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (n + 1 < min(1, n)*num_stcls_min*weno_order) then
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("n + 1 < min(1, n)*num_stcls_min*weno_order", "For 2D simulation, n must be greater than or equal to (num_stcls_min*weno_order - 1), whose value is "  // trim(numstr))
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (p + 1 < min(1, p)*num_stcls_min*weno_order) then
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("p + 1 < min(1, p)*num_stcls_min*weno_order", "For 3D simulation, p must be greater than or equal to (num_stcls_min*weno_order - 1), whose value is "  // trim(numstr))
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    end subroutine s_check_inputs_weno
 
    !> Validate that the grid resolution is sufficient for the MUSCL reconstruction order
    impure subroutine s_check_inputs_muscl
 
        character(len=5) :: numstr  !< for int to string conversion
 
        call s_int_to_str(num_stcls_min*muscl_order, numstr)
    if (m + 1 < num_stcls_min*muscl_order) then
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("m + 1 < num_stcls_min*muscl_order", "m must be greater than or equal to (num_stcls_min*muscl_order - 1), whose value is " // trim(numstr))
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 145 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (n + 1 < min(1, n)*num_stcls_min*muscl_order) then
# 145 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("n + 1 < min(1, n)*num_stcls_min*muscl_order", "For 2D simulation, n must be greater than or equal to (num_stcls_min*muscl_order - 1), whose value is "  // trim(numstr))
# 145 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (p + 1 < min(1, p)*num_stcls_min*muscl_order) then
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("p + 1 < min(1, p)*num_stcls_min*muscl_order", "For 3D simulation, p must be greater than or equal to (num_stcls_min*muscl_order - 1), whose value is "  // trim(numstr))
# 148 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 151 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (muscl_order == muscl_order_first_order .and. int_comp > 0) then
# 151 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."muscl_order == muscl_order_first_order  int_comp > 0", "int_comp requires muscl_order >= 2 (muscl_order=1 leaves the reconstruction workspace uninitialised)")
# 151 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 153 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    end subroutine s_check_inputs_muscl
 
    !> Checks constraints on time stepping parameters
    impure subroutine s_check_inputs_time_stepping
 
        if (.not. cfl_dt) then
    if (dt <= 0) then
# 160 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("dt <= 0", "")
# 160 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
        end if
 
    end subroutine s_check_inputs_time_stepping
 
    !> Validate NVIDIA unified virtual memory configuration parameters
    impure subroutine s_check_inputs_nvidia_uvm
 
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
    if (nv_uvm_igr_temps_on_gpu > 3 .or. nv_uvm_igr_temps_on_gpu < 0) then
# 169 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.or."nv_uvm_igr_temps_on_gpu > 3  nv_uvm_igr_temps_on_gpu < 0", "nv_uvm_igr_temps_on_gpu must be in the range [0, 3]")
# 169 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 171 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (nv_uvm_igr_temps_on_gpu == 3 .and. igr_iter_solver == 2) then
# 171 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."nv_uvm_igr_temps_on_gpu == 3  igr_iter_solver == 2", "nv_uvm_igr_temps_on_gpu must be in the range [0, 2] for igr_iter_solver == 2")
# 171 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 173 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
#endif
 
    end subroutine s_check_inputs_nvidia_uvm
 
    impure subroutine s_check_inputs_hll_non_conservative
 
    if ((riemann_solver == riemann_solver_hll) .and. hll_u_interface .and. cyl_coord .and. p > 0) then
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..and."(riemann_solver == riemann_solver_hll)  hll_u_interface  cyl_coord  p > 0", "HLL Method 2 is not supported for 3D cylindrical geometry")
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (alt_soundspeed .and. riemann_solver == riemann_solver_hll .and. (.not. hll_u_interface) .and. cyl_coord  .and. p == 0) then
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..not..and..and."alt_soundspeed  riemann_solver == riemann_solver_hll  ( hll_u_interface)  cyl_coord   p == 0", "alt_soundspeed with HLL Method 1 is not supported for 2D axisymmetric geometry")
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 183 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (alt_soundspeed .and. riemann_solver == riemann_solver_hll .and. cyl_coord .and. p > 0) then
# 183 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..and."alt_soundspeed  riemann_solver == riemann_solver_hll  cyl_coord  p > 0", "alt_soundspeed with HLL is not currently supported for 3D cylindrical geometry")
# 183 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 185 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hll_u_interface .and. riemann_solver /= riemann_solver_hll) then
# 185 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hll_u_interface  riemann_solver /= riemann_solver_hll", "hll_u_interface requires the HLL Riemann solver (riemann_solver = 1)")
# 185 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 187 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hll_u_interface .and. mhd) then
# 187 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hll_u_interface  mhd", "HLL Method 2 does not support MHD (the MHD path zeroes the shared interface-velocity trace)")
# 187 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (surface_tension .and. riemann_solver == riemann_solver_hll .and. (.not. hll_u_interface)) then
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..not."surface_tension  riemann_solver == riemann_solver_hll  ( hll_u_interface)", "surface_tension requires a shared interface-velocity representation (HLL Method 2 or HLLC)")
# 189 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 191 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (surface_tension .and. riemann_solver == riemann_solver_lax_friedrichs) then
# 191 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."surface_tension  riemann_solver == riemann_solver_lax_friedrichs", "surface_tension requires a shared interface-velocity representation (HLL Method 2 or HLLC)")
# 191 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    end subroutine s_check_inputs_hll_non_conservative
 
    impure subroutine s_check_inputs_hypo_branch
 
    if (hypoelasticity .and. cyl_coord .and. p > 0) then
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  cyl_coord  p > 0", "3D cylindrical hypoelasticity is not supported")
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
 
        ! Hypoelasticity solver restrictions
    if (hypoelasticity .and. model_eqns /= model_eqns_5eq) then
# 201 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypoelasticity  model_eqns /= model_eqns_5eq", "hypoelasticity requires model_eqns = 2")
# 201 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. riemann_solver /= riemann_solver_hll .and. riemann_solver /= riemann_solver_hllc  .and. riemann_solver /= riemann_solver_hlld) then
# 202 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..and."hypoelasticity  riemann_solver /= riemann_solver_hll  riemann_solver /= riemann_solver_hllc   riemann_solver /= riemann_solver_hlld", "hypoelasticity requires HLL (1), HLLC (2), or HLLD (4) Riemann solver")
# 202 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 205 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. n == 0) then
# 205 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  n == 0", "HLLD hypoelasticity requires at least 2D (n must be > 0)")
# 205 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. num_fluids /= 2) then
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  num_fluids /= 2", "HLLD hypoelasticity currently requires exactly 2 fluid components")
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. mhd) then
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypoelasticity  mhd", "MHD and hypoelasticity cannot be enabled together")
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. bubbles_euler) then
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypoelasticity  bubbles_euler", "Hypoelasticity does not support Euler-Euler bubbles")
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. viscous) then
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  viscous", "HLLD hypoelasticity does not support viscous effects (the dual-pass omits the viscous source term)")
# 211 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. surface_tension) then
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  surface_tension", "HLLD hypoelasticity does not support surface tension (the dual-pass omits the surface-tension source term)")
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. cont_damage) then
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  cont_damage", "HLLD hypoelasticity does not support continuum damage (the dual-pass does not damage-scale the shear modulus)")
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. riemann_solver == riemann_solver_hlld .and. chemistry) then
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."hypoelasticity  riemann_solver == riemann_solver_hlld  chemistry", "HLLD hypoelasticity does not support chemistry")
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 219 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (riemann_solver == riemann_solver_hlld .and. (.not. mhd) .and. (.not. hypoelasticity)) then
# 219 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not..and..not."riemann_solver == riemann_solver_hlld  ( mhd)  ( hypoelasticity)", "HLLD is only available for MHD or hypoelasticity")
# 219 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 221 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
        ! Feature flag prerequisites
    if (riemann_hypo_adc .and. .not. hypoelasticity) then
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."riemann_hypo_ADC   hypoelasticity", "riemann_hypo_ADC requires hypoelasticity = T")
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (riemann_hypo_adc .and. riemann_solver /= riemann_solver_hllc .and. riemann_solver /= riemann_solver_hlld) then
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and."riemann_hypo_ADC  riemann_solver /= riemann_solver_hllc  riemann_solver /= riemann_solver_hlld", "riemann_hypo_ADC only applies to hypo HLLC/HLLD")
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 226 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (riemann_hypo_adc .and. (bubbles_euler .or. surface_tension .or. chemistry .or. cont_damage)) then
# 226 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..or..or..or."riemann_hypo_ADC  (bubbles_euler  surface_tension  chemistry  cont_damage)", "riemann_hypo_ADC does not support bubbles, surface tension, chemistry, or continuum damage (the ADC HLL blend omits their flux components)")
# 226 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypo_hll_interface_rhs .and. .not. hypoelasticity) then
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..not."hypo_hll_interface_rhs   hypoelasticity", "hypo_hll_interface_rhs requires hypoelasticity = T")
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypo_hll_interface_rhs .and. riemann_solver /= riemann_solver_hll) then
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypo_hll_interface_rhs  riemann_solver /= riemann_solver_hll", "hypo_hll_interface_rhs requires HLL Riemann solver (riemann_solver = 1)")
# 229 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 231 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (alt_soundspeed .and. riemann_solver == riemann_solver_hlld .and. .not. hypoelasticity) then
# 231 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..not."alt_soundspeed  riemann_solver == riemann_solver_hlld   hypoelasticity", "alt_soundspeed with HLLD requires hypoelasticity = T")
# 231 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (alt_soundspeed .and. num_fluids /= 2) then
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."alt_soundspeed  num_fluids /= 2", "alt_soundspeed requires exactly 2 fluid components (the Kapila K coefficient is a two-fluid closure)")
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 235 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (cont_damage .and. alt_soundspeed) then
# 235 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."cont_damage  alt_soundspeed", "Continuum damage does not support alt_soundspeed")
# 235 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. igr) then
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypoelasticity  igr", "Hypoelasticity is not compatible with IGR")
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. fd_order /= 1 .and. fd_order /= 2 .and. fd_order /= 4) then
# 237 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..and..and."hypoelasticity  fd_order /= 1  fd_order /= 2  fd_order /= 4", "hypoelasticity requires fd_order to be set to 1, 2, or 4 (the finite-difference coefficients are initialized unconditionally)")
# 237 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 239 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (hypoelasticity .and. wave_speeds == wave_speeds_pressure) then
# 239 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."hypoelasticity  wave_speeds == wave_speeds_pressure", "Pressure-based wave speeds (wave_speeds = 2) omit the elastic longitudinal speed and are not supported with hypoelasticity")
# 239 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 241 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (riemann_solver == riemann_solver_hlld .and. wave_speeds == wave_speeds_pressure) then
# 241 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."riemann_solver == riemann_solver_hlld  wave_speeds == wave_speeds_pressure", "HLLD uses its own direct wave-speed estimates; wave_speeds = 2 is not supported")
# 241 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 243 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (riemann_solver == riemann_solver_hlld .and. low_mach > 0) then
# 243 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."riemann_solver == riemann_solver_hlld  low_Mach > 0", "low_Mach corrections are not implemented for HLLD")
# 243 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (riemann_hypo_adc .and. low_mach > 0) then
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."riemann_hypo_ADC  low_Mach > 0", "riemann_hypo_ADC does not support low_Mach (the ADC HLL reference flux uses pre-correction velocities)")
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (riemann_hypo_adc .and. adc_kappa <= 0._wp) then
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."riemann_hypo_ADC  ADC_kappa <= 0._wp", "ADC_kappa must be positive")
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
    if (hypoelasticity .and. (cfl_const_dt .or. cfl_adap_dt)) then
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and..or."hypoelasticity  (cfl_const_dt  cfl_adap_dt)", "Automatic CFL time stepping uses the acoustic sound speed only and does not bound the elastic characteristic speed; set dt explicitly for hypoelastic runs")
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 249 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
 
    end subroutine s_check_inputs_hypo_branch
 
    !> Validates that each burning immersed-boundary patch injects a species index within the mechanism. inj_species indexes the
    !! image-point mass-fraction array Ys_IP(1:num_species) in m_ibm; an out-of-range value is an out-of-bounds write (silent
    !! corruption). Only reachable with chemistry.
    impure subroutine s_check_inputs_ib_injection
 
        integer :: i
 
        do i = 1, num_ibs
    if (patch_ib(i)%inj_species > num_species) then
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("patch_ib(i)%inj_species > num_species", "patch_ib inj_species must be <= num_species (it indexes the image-point species mass fractions; an out-of-range value writes out of bounds)")
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 262 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        end do
 
    end subroutine s_check_inputs_ib_injection
 
    !> Checks that each active particle cloud has a valid packing_method specified
    impure subroutine s_check_inputs_particle_clouds
 
        integer          :: i
        character(len=5) :: idxstr
 
        do i = 1, num_particle_clouds
            call s_int_to_str(i, idxstr)
    if (particle_cloud(i)%packing_method == dflt_int) then
# 274 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("particle_cloud(i)%packing_method == dflt_int", "particle_cloud("//trim(idxstr)  //")%packing_method must be specified (1 = rejection sampling, 2 = lattice)")
# 274 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 277 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (particle_cloud(i)%packing_method /= 1 .and. particle_cloud(i)%packing_method /= 2) then
# 277 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.and."particle_cloud(i)%packing_method /= 1  particle_cloud(i)%packing_method /= 2", "particle_cloud("//trim(idxstr) //")%packing_method must be 1 (rejection sampling) or 2 (lattice)")
# 277 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        end do
 
    end subroutine s_check_inputs_particle_clouds
 
    !> Checks that each active synthetic-turbulence forcing zone has a fully specified position and a positive size in every active
    !! dimension
    impure subroutine s_check_inputs_synthetic_turbulence
 
        integer          :: i, d
        character(len=5) :: idxstr
 
    if (num_turbulent_sources <= 0) then
# 290 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("num_turbulent_sources <= 0", "num_turbulent_sources must be > 0 when synthetic_turbulence is enabled")
# 290 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
 
        do i = 1, num_turbulent_sources
            call s_int_to_str(i, idxstr)
            do d = 1, num_dims
    if (f_is_default(turb_pos(i, d))) then
# 295 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort("f_is_default(turb_pos(i, d))", "turb_pos("//trim(idxstr)  //",:) must be specified for all num_dims when synthetic_turbulence is enabled")
# 295 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    if (f_is_default(synth_l(i, d)) .or. synth_l(i, d) <= 0._wp) then
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
        call s_prohibit_abort(.or."f_is_default(synth_L(i, d))  synth_L(i, d) <= 0._wp", "synth_L("//trim(idxstr)//",:) must be positive for all num_dims when synthetic_turbulence is enabled")
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
    end if
# 300 "/home/runner/work/MFC/MFC/src/simulation/m_checker.fpp"
            end do
        end do
 
    end subroutine s_check_inputs_synthetic_turbulence
 
end module m_checker