m_compute_cbc.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!>
!! @file
!! @brief CBC computation module
# 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"
# 5 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp" 2
# 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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 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"
 
# 95 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 100 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 115 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 165 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 219 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 233 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 271 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 282 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 332 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 342 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 361 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 378 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 388 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 395 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 401 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 407 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 413 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 419 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 425 "/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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 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
# 6 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp" 2
 
!> @brief Characteristic boundary condition (CBC) computations for subsonic inflow, outflow, and slip walls
module m_compute_cbc
 
    use m_global_parameters
 
    implicit none
 
    private; public :: s_compute_slip_wall_l, s_compute_nonreflecting_subsonic_buffer_l, &
        & s_compute_nonreflecting_subsonic_inflow_l, s_compute_nonreflecting_subsonic_outflow_l, &
        & s_compute_force_free_subsonic_outflow_l, s_compute_constant_pressure_subsonic_outflow_l, s_compute_supersonic_inflow_l, &
        & s_compute_supersonic_outflow_l
 
contains
    !> Base L1 calculation
    function f_base_l1(lambda, rho, c, dpres_ds, dvel_ds) result(L1)
 
 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 23 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 24 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), dimension(3), intent(in) :: lambda
        real(wp), intent(in)               :: rho, c, dpres_ds
# 29 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_dims), intent(in) :: dvel_ds
# 31 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp) :: l1
        l1 = lambda(1)*(dpres_ds - rho*c*dvel_ds(dir_idx(1)))
 
    end function f_base_l1
 
    !> Fill density L variables
    subroutine s_fill_density_l(L, lambda_factor, lambda2, c, mf, dalpha_rho_ds, dpres_ds)
 
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 39 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: L
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: lambda_factor, lambda2, c
        real(wp), intent(in) :: dpres_ds
        integer              :: i
 
        do i = 2, eqn_idx%mom%beg
            l(i) = lambda_factor*lambda2*(c*c*dalpha_rho_ds(i - 1) - mf(i - 1)*dpres_ds)
        end do
 
    end subroutine s_fill_density_l
 
    !> Fill velocity L variables
    subroutine s_fill_velocity_l(L, lambda_factor, lambda2, dvel_ds)
 
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 63 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 67 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: L
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 72 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_dims), intent(in) :: dvel_ds
# 74 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: lambda_factor, lambda2
        integer              :: i
 
        do i = eqn_idx%mom%beg + 1, eqn_idx%mom%end
            l(i) = lambda_factor*lambda2*dvel_ds(dir_idx(i - eqn_idx%cont%end))
        end do
 
    end subroutine s_fill_velocity_l
 
    !> Fill advection L variables
    subroutine s_fill_advection_l(L, lambda_factor, lambda2, dadv_ds)
 
 
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 86 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 87 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: L
# 92 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: lambda_factor, lambda2
        integer              :: i
 
        do i = eqn_idx%E, eqn_idx%adv%end - 1
            l(i) = lambda_factor*lambda2*dadv_ds(i - eqn_idx%mom%end)
        end do
 
    end subroutine s_fill_advection_l
 
    !> Fill chemistry L variables
    subroutine s_fill_chemistry_l(L, lambda_factor, lambda2, dYs_ds)
 
 
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 109 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 110 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: L
# 115 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), dimension(1), intent(in) :: dYs_ds
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in)                             :: lambda_factor, lambda2
        integer                                          :: i
 
        if (.not. chemistry) return
 
        do i = eqn_idx%species%beg, eqn_idx%species%end
            l(i) = lambda_factor*lambda2*dys_ds(i - eqn_idx%species%beg + 1)
        end do
 
    end subroutine s_fill_chemistry_l
 
    !> Slip wall CBC (Thompson 1990, pg. 451)
    subroutine s_compute_slip_wall_l(lambda, L, rho, c, dpres_ds, dvel_ds)
 
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_slip_wall_L
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_dims), intent(in) :: dvel_ds
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: rho, c, dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        l(2:eqn_idx%adv%end - 1) = 0._wp
        l(eqn_idx%adv%end) = l(1)
 
    end subroutine s_compute_slip_wall_l
 
    !> Nonreflecting subsonic buffer CBC (Thompson 1987, pg. 13)
    subroutine s_compute_nonreflecting_subsonic_buffer_l(lambda, L, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dYs_ds)
 
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_nonreflecting_subsonic_buffer_L
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 155 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 160 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 162 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 167 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
            real(wp), dimension(num_dims), intent(in)   :: dvel_ds
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 171 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), dimension(1), intent(in) :: dys_ds
# 172 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in)                             :: rho, c
        real(wp), intent(in)                             :: dpres_ds
        real(wp)                                         :: lambda_factor
 
        lambda_factor = (5.e-1_wp - 5.e-1_wp*sign(1._wp, lambda(1)))
        l(1) = lambda_factor*lambda(1)*(dpres_ds - rho*c*dvel_ds(dir_idx(1)))
 
        lambda_factor = (5.e-1_wp - 5.e-1_wp*sign(1._wp, lambda(2)))
        call s_fill_density_l(l, lambda_factor, lambda(2), c, mf, dalpha_rho_ds, dpres_ds)
        call s_fill_velocity_l(l, lambda_factor, lambda(2), dvel_ds)
        call s_fill_advection_l(l, lambda_factor, lambda(2), dadv_ds)
        call s_fill_chemistry_l(l, lambda_factor, lambda(2), dys_ds)
 
        lambda_factor = (5.e-1_wp - 5.e-1_wp*sign(1._wp, lambda(3)))
        l(eqn_idx%adv%end) = lambda_factor*lambda(3)*(dpres_ds + rho*c*dvel_ds(dir_idx(1)))
 
    end subroutine s_compute_nonreflecting_subsonic_buffer_l
 
    !> Nonreflecting subsonic inflow CBC (Thompson 1990, pg. 455)
    subroutine s_compute_nonreflecting_subsonic_inflow_l(lambda, L, rho, c, dpres_ds, dvel_ds)
 
 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_nonreflecting_subsonic_inflow_L
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 193 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 194 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 199 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 201 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 204 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_dims), intent(in) :: dvel_ds
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: rho, c, dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        l(2:eqn_idx%adv%end) = 0._wp
        if (chemistry) l(eqn_idx%species%beg:eqn_idx%species%end) = 0._wp
 
    end subroutine s_compute_nonreflecting_subsonic_inflow_l
 
    !> Nonreflecting subsonic outflow CBC (Thompson 1990, pg. 454)
    subroutine s_compute_nonreflecting_subsonic_outflow_l(lambda, L, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dYs_ds)
 
 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_nonreflecting_subsonic_outflow_L
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 217 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 218 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 225 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 230 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
            real(wp), dimension(num_dims), intent(in)   :: dvel_ds
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 234 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), dimension(1), intent(in) :: dys_ds
# 235 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in)                             :: rho, c
        real(wp), intent(in)                             :: dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        call s_fill_density_l(l, 1._wp, lambda(2), c, mf, dalpha_rho_ds, dpres_ds)
        call s_fill_velocity_l(l, 1._wp, lambda(2), dvel_ds)
        call s_fill_advection_l(l, 1._wp, lambda(2), dadv_ds)
        call s_fill_chemistry_l(l, 1._wp, lambda(2), dys_ds)
        l(eqn_idx%adv%end) = 0._wp
 
    end subroutine s_compute_nonreflecting_subsonic_outflow_l
 
    !> Force-free subsonic outflow CBC (Thompson 1990, pg. 454)
    subroutine s_compute_force_free_subsonic_outflow_l(lambda, L, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
 
 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_force_free_subsonic_outflow_L
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 250 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 251 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 256 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 258 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 263 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
            real(wp), dimension(num_dims), intent(in)   :: dvel_ds
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 267 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: rho, c
        real(wp), intent(in) :: dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        call s_fill_density_l(l, 1._wp, lambda(2), c, mf, dalpha_rho_ds, dpres_ds)
        call s_fill_velocity_l(l, 1._wp, lambda(2), dvel_ds)
        call s_fill_advection_l(l, 1._wp, lambda(2), dadv_ds)
        l(eqn_idx%adv%end) = l(1) + 2._wp*rho*c*lambda(2)*dvel_ds(dir_idx(1))
 
    end subroutine s_compute_force_free_subsonic_outflow_l
 
    !> Constant pressure subsonic outflow CBC (Thompson 1990, pg. 455)
    subroutine s_compute_constant_pressure_subsonic_outflow_l(lambda, L, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds)
 
 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_constant_pressure_subsonic_outflow_L
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 281 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 282 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 287 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 289 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 294 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
            real(wp), dimension(num_dims), intent(in)   :: dvel_ds
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in) :: rho, c
        real(wp), intent(in) :: dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        call s_fill_density_l(l, 1._wp, lambda(2), c, mf, dalpha_rho_ds, dpres_ds)
        call s_fill_velocity_l(l, 1._wp, lambda(2), dvel_ds)
        call s_fill_advection_l(l, 1._wp, lambda(2), dadv_ds)
        l(eqn_idx%adv%end) = -l(1)
 
    end subroutine s_compute_constant_pressure_subsonic_outflow_l
 
    !> Supersonic inflow CBC (Thompson 1990, pg. 453)
    subroutine s_compute_supersonic_inflow_l(L)
 
 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_supersonic_inflow_L
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 312 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 313 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 316 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 318 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        l(1:eqn_idx%adv%end) = 0._wp
        if (chemistry) l(eqn_idx%species%beg:eqn_idx%species%end) = 0._wp
 
    end subroutine s_compute_supersonic_inflow_l
 
    !> Supersonic outflow CBC (Thompson 1990, pg. 453)
    subroutine s_compute_supersonic_outflow_l(lambda, L, rho, c, mf, dalpha_rho_ds, dpres_ds, dvel_ds, dadv_ds, dYs_ds)
 
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#ifdef _CRAYFTN
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#if MFC_OpenACC
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#else
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!DIR$ INLINEALWAYS s_compute_supersonic_outflow_L
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenACC
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$acc routine seq 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#elif MFC_OpenMP
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
!$omp declare target device_type(any) 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
#endif
# 327 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
 
        real(wp), dimension(3), intent(in) :: lambda
# 332 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(sys_size), intent(inout) :: l
# 334 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
            real(wp), dimension(num_fluids), intent(in) :: mf, dalpha_rho_ds
            real(wp), dimension(num_dims), intent(in)   :: dvel_ds
            real(wp), dimension(num_fluids), intent(in) :: dadv_ds
# 343 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), dimension(1), intent(in) :: dys_ds
# 344 "/home/runner/work/MFC/MFC/src/simulation/m_compute_cbc.fpp"
        real(wp), intent(in)                             :: rho, c
        real(wp), intent(in)                             :: dpres_ds
 
        l(1) = f_base_l1(lambda, rho, c, dpres_ds, dvel_ds)
        call s_fill_density_l(l, 1._wp, lambda(2), c, mf, dalpha_rho_ds, dpres_ds)
        call s_fill_velocity_l(l, 1._wp, lambda(2), dvel_ds)
        call s_fill_advection_l(l, 1._wp, lambda(2), dadv_ds)
        call s_fill_chemistry_l(l, 1._wp, lambda(2), dys_ds)
        l(eqn_idx%adv%end) = lambda(3)*(dpres_ds + rho*c*dvel_ds(dir_idx(1)))
 
    end subroutine s_compute_supersonic_outflow_l
 
end module m_compute_cbc