m_conduction.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!>
!! @file
!! @brief Contains module m_conduction
 
# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 197 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 211 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 236 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 247 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 249 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 260 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 310 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 320 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 356 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 366 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 373 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 379 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 385 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 391 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 397 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 403 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
! New line at end of file is required for FYPP
# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 15 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 16 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
 
# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
! GPU parallel region (scalar reductions, maxval/minval)
# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! GPU parallel loop over threads (most common GPU macro)
# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Required closing for GPU_PARALLEL_LOOP
# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Mark routine for device compilation
# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Declare device-resident data
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Inner loop within a GPU parallel region
# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Scoped GPU data region
# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Host code with device pointers (for MPI with GPU buffers)
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Allocate device memory (unscoped)
# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Free device memory
# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Atomic operation on device
# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! End atomic capture block
# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Copy data between host and device
# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Synchronization barrier
# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Import GPU library module (openacc or omp_lib)
# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for AMD compiler
# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-Cray compilers
# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for Cray compiler
# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-NVIDIA compilers
# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
! example see misc/nvidia_uvm/bind.sh.
# 52 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 72 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 80 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 104 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 120 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 145 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 158 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 6 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.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_conduction.fpp" 2
 
!> @brief Fourier heat conduction, div(k grad T), as a face-centered source flux on the energy equation. The term has no
!! cross-derivatives, so the direction-split face difference below is exact for it.
module m_conduction
 
    use m_global_parameters
 
    implicit none
 
    private; public :: s_compute_conduction_source_flux, s_compute_conduction_axis_source
 
    type(int_bounds_info) :: isc1, isc2, isc3
 
# 19 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 19 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc declare create(isc1, isc2, isc3) 
# 19 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 19 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp declare target (isc1, isc2, isc3) 
# 19 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
    integer, dimension(3) :: offsets_c
 
# 21 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 21 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc declare create(offsets_c) 
# 21 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 21 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp declare target (offsets_c) 
# 21 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
 
contains
 
    !> Accumulate -k*dT/dx_idir into the energy source flux at each idir-normal face.
    subroutine s_compute_conduction_source_flux(idir, q_prim_qp, q_T_sf, flux_src_vf, irx, iry, irz)
 
        integer, intent(in)                                    :: idir
        type(scalar_field), dimension(sys_size), intent(in)    :: q_prim_qp
        type(scalar_field), intent(in)                         :: q_t_sf
        type(scalar_field), dimension(sys_size), intent(inout) :: flux_src_vf
        type(int_bounds_info), intent(in)                      :: irx, iry, irz
        real(wp)                                               :: k_face, dt_dxi, grid_spacing, alpha_face
        integer                                                :: x, y, z, i
 
        isc1 = irx; isc2 = iry; isc3 = irz
        offsets_c = 0
        offsets_c(idir) = 1
 
 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc update device(isc1, isc2, isc3, offsets_c) 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp target update to(isc1, isc2, isc3, offsets_c) 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
 
 
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(k_face, dT_dxi, grid_spacing, alpha_face, i)
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) &
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp& private(k_face, dT_dxi, grid_spacing, alpha_face, i)
# 42 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
        do z = isc3%beg, isc3%end
            do y = isc2%beg, isc2%end
                do x = isc1%beg, isc1%end
                    select case (idir)
                    case (1)
                        grid_spacing = x_cc(x + 1) - x_cc(x)
                    case (2)
                        grid_spacing = y_cc(y + 1) - y_cc(y)
                    case (3)
                        grid_spacing = z_cc(z + 1) - z_cc(z)
                        ! 3D cylindrical z is the azimuth: dz is in radians. The r**2 carries both metric
                        ! factors of (1/r**2) d2T/dtheta2, since m_rhs divides this flux by dz alone.
                        if (grid_geometry == 3) grid_spacing = y_cc(y)**2*grid_spacing
                    end select
 
                    ! Volume-fraction-weighted face conductivity. Raw cell-centered alphas over- and
                    ! undershoot near interfaces, so clamp the face average as the viscous path does.
                    k_face = 0._wp
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc loop seq 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
                    do i = 1, num_fluids
                        alpha_face = 0.5_wp*(q_prim_qp(eqn_idx%adv%beg + i - 1)%sf(x, y, &
                                             & z) + q_prim_qp(eqn_idx%adv%beg + i - 1)%sf(x + offsets_c(1), y + offsets_c(2), &
                                             & z + offsets_c(3)))
                        alpha_face = min(max(alpha_face, 0._wp), 1._wp)
                        k_face = k_face + alpha_face*fluid_k_therm(i)
                    end do
 
                    dt_dxi = (q_t_sf%sf(x + offsets_c(1), y + offsets_c(2), z + offsets_c(3)) - q_t_sf%sf(x, y, z))/grid_spacing
 
                    flux_src_vf(eqn_idx%E)%sf(x, y, z) = flux_src_vf(eqn_idx%E)%sf(x, y, z) - k_face*dt_dxi
                end do
            end do
        end do
 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc end parallel loop
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp end target teams loop
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
 
    end subroutine s_compute_conduction_source_flux
 
    !> Cell-centered -k*dT/dr for the axis cell of a cylindrical grid. The generic geometric source in m_rhs uses face values of
    !! flux_src(E) in this role; the axis cell has no face pair, so it reads tau_Re_vf(E) instead. Accumulates: call after
    !! s_compute_viscous_stress_cylindrical_boundary, which zeroes tau_Re_vf(mom%beg:E); without viscosity this routine does that
    !! zeroing itself.
    subroutine s_compute_conduction_axis_source(q_prim_vf, q_T_sf, tau_Re_vf, ix, iy, iz)
 
        type(scalar_field), dimension(sys_size), intent(in)      :: q_prim_vf
        type(scalar_field), intent(in)                           :: q_t_sf
        type(scalar_field), dimension(1:sys_size), intent(inout) :: tau_re_vf
        type(int_bounds_info), intent(in)                        :: ix, iy, iz
        real(wp)                                                 :: k_cell, dt_dr, dr_l, dr_r, alpha_cell
        integer                                                  :: j, k, l, i
 
        isc1 = ix; isc2 = iy; isc3 = iz
 
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc update device(isc1, isc2, isc3) 
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp target update to(isc1, isc2, isc3) 
# 94 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
 
        if (.not. viscous) then
 
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i)
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i)
# 97 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
            do l = isc3%beg, isc3%end
                do k = isc2%beg, isc2%end
                    do j = isc1%beg, isc1%end
 
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc loop seq 
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
                        do i = eqn_idx%mom%beg, eqn_idx%E
                            tau_re_vf(i)%sf(j, k, l) = 0._wp
                        end do
                    end do
                end do
            end do
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc end parallel loop
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp end target teams loop
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
        end if
 
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(k_cell, dT_dr, dr_l, dr_r, alpha_cell, i)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) &
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp& private(k_cell, dT_dr, dr_l, dr_r, alpha_cell, i)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
        do l = isc3%beg, isc3%end
            do k = isc2%beg + 1, isc2%end - 1
                do j = isc1%beg, isc1%end
                    k_cell = 0._wp
 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc loop seq 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
                    do i = 1, num_fluids
                        alpha_cell = min(max(q_prim_vf(eqn_idx%adv%beg + i - 1)%sf(j, k, l), 0._wp), 1._wp)
                        k_cell = k_cell + alpha_cell*fluid_k_therm(i)
                    end do
 
                    ! Distance-weighted central difference. The axisymmetric grid halves the cell at
                    ! r = 0, so the plain form would be only first order over exactly this stencil.
                    dr_l = y_cc(k) - y_cc(k - 1)
                    dr_r = y_cc(k + 1) - y_cc(k)
                    dt_dr = ((q_t_sf%sf(j, k, l) - q_t_sf%sf(j, k - 1, l))*dr_r/dr_l + (q_t_sf%sf(j, k + 1, l) - q_t_sf%sf(j, k, &
                             & l))*dr_l/dr_r)/(dr_l + dr_r)
 
                    tau_re_vf(eqn_idx%E)%sf(j, k, l) = tau_re_vf(eqn_idx%E)%sf(j, k, l) - k_cell*dt_dr
                end do
            end do
        end do
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#if defined(MFC_OpenACC)
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$acc end parallel loop
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#elif defined(MFC_OpenMP)
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
!$omp end target teams loop
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_conduction.fpp"
#endif
 
    end subroutine s_compute_conduction_axis_source
 
end module m_conduction