m_body_forces.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!>
!! @file
!! @brief Contains module m_body_forces
 
# 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"
! 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"
! 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"
 
# 207 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 232 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 243 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 245 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 283 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 293 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 303 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 329 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 346 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 358 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 364 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 370 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 376 "/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"
! 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"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 213 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 241 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 266 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 275 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 291 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 301 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 308 "/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
 
# 21 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 104 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 119 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 171 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 182 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 204 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 214 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 225 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 236 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 246 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 252 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 264 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 270 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 272 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 273 "/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.
# 63 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 153 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 167 "/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_body_forces.fpp" 2
 
!> @brief Computes gravitational and user-defined body force source terms for the momentum equations
module m_body_forces
 
    use m_derived_types        !< definitions of the derived types
 
    use m_global_parameters    !< definitions of the global parameters
 
    use m_variables_conversion
 
    use m_nvtx
 
! $:USE_GPU_MODULE()
 
    implicit none
 
    private; 
    public :: s_compute_body_forces_rhs, &
              s_initialize_body_forces_module, &
              s_finalize_body_forces_module
 
    real(wp), allocatable, dimension(:, :, :) :: rhom
 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc declare create(rhoM) 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp declare target (rhoM) 
# 28 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
 
contains
 
    !> This subroutine initializes the module global array of mixture
    !! densities in each grid cell
    impure subroutine s_initialize_body_forces_module
 
        ! Simulation is at least 2D
        if (n > 0) then
            ! Simulation is 3D
            if (p > 0) then
#ifdef MFC_DEBUG
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    block
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        use iso_fortran_env, only: output_unit
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        print *, 'm_body_forces.fpp:40: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, -buff_size:buff_size + n, -buff_size:buff_size + p))'
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        call flush (output_unit)
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    end block
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    allocate (rhom(-buff_size:buff_size + m,                     -buff_size:buff_size + n,                     -buff_size:buff_size + p))
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc enter data create(rhoM) 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target enter data map(always,alloc:rhoM) 
# 40 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
                ! Simulation is 2D
            else
#ifdef MFC_DEBUG
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    block
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        use iso_fortran_env, only: output_unit
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        print *, 'm_body_forces.fpp:45: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, -buff_size:buff_size + n, 0:0))'
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        call flush (output_unit)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    end block
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    allocate (rhom(-buff_size:buff_size + m,                     -buff_size:buff_size + n,                     0:0))
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc enter data create(rhoM) 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target enter data map(always,alloc:rhoM) 
# 45 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 48 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
            end if
            ! Simulation is 1D
        else
#ifdef MFC_DEBUG
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    block
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        use iso_fortran_env, only: output_unit
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        print *, 'm_body_forces.fpp:51: ', '@:ALLOCATE(rhoM(-buff_size:buff_size + m, 0:0, 0:0))'
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        call flush (output_unit)
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    end block
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    allocate (rhom(-buff_size:buff_size + m,                 0:0,                 0:0))
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc enter data create(rhoM) 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target enter data map(always,alloc:rhoM) 
# 51 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 54 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        end if
 
    end subroutine s_initialize_body_forces_module
 
    !> This subroutine computes the acceleration at time t
    subroutine s_compute_acceleration(t)
 
        real(wp), intent(in) :: t
 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
            if (bf_x) then
                accel_bf(1) = g_x+k_x*sin(w_x*t - p_x)
            end if
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
            if (bf_y) then
                accel_bf(2) = g_y+k_y*sin(w_y*t - p_y)
            end if
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
            if (bf_z) then
                accel_bf(3) = g_z+k_z*sin(w_z*t - p_z)
            end if
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc update device(accel_bf) 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target update to(accel_bf) 
# 69 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
 
    end subroutine s_compute_acceleration
 
    !> This subroutine calculates the mixture density at each cell
    !! center
    !! param q_cons_vf Conservative variable
    subroutine s_compute_mixture_density(q_cons_vf)
 
        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
        integer :: i, j, k, l !< standard iterators
 
 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l) 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l) 
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 81 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
        do l = 0, p
            do k = 0, n
                do j = 0, m
                    rhom(j, k, l) = 0._wp
                    do i = 1, num_fluids
                        rhom(j, k, l) = rhom(j, k, l) + &
                                        q_cons_vf(contxb + i - 1)%sf(j, k, l)
                    end do
                end do
            end do
        end do
 
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc end parallel loop
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp end target teams loop
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 93 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
 
    end subroutine s_compute_mixture_density
 
    !> This subroutine calculates the source term due to body forces
    !! so the system can be advanced in time
    !! @param q_cons_vf Conservative variables
    !! @param q_prim_vf Primitive variables
    !! @param rhs_vf Right-hand side accumulator
    subroutine s_compute_body_forces_rhs(q_prim_vf, q_cons_vf, rhs_vf)
 
        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
        type(scalar_field), dimension(sys_size), intent(inout) :: rhs_vf
 
        integer :: i, j, k, l !< Loop variables
 
        call s_compute_acceleration(mytime)
        call s_compute_mixture_density(q_cons_vf)
 
 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc parallel loop collapse(4) gang vector default(present) private(i, j, k, l) 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(4) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k, l) 
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 113 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
        do i = momxb, e_idx
            do l = 0, p
                do k = 0, n
                    do j = 0, m
                        rhs_vf(i)%sf(j, k, l) = 0._wp
                    end do
                end do
            end do
        end do
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc end parallel loop
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp end target teams loop
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
 
        if (bf_x) then ! x-direction body forces
 
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l) 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l) 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
            do l = 0, p
                do k = 0, n
                    do j = 0, m
                        rhs_vf(momxb)%sf(j, k, l) = rhs_vf(momxb)%sf(j, k, l) + &
                                                    rhom(j, k, l)*accel_bf(1)
                        rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
                                                    q_cons_vf(momxb)%sf(j, k, l)*accel_bf(1)
                    end do
                end do
            end do
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc end parallel loop
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp end target teams loop
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
        end if
 
        if (bf_y) then ! y-direction body forces
 
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l) 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l) 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
            do l = 0, p
                do k = 0, n
                    do j = 0, m
                        rhs_vf(momxb + 1)%sf(j, k, l) = rhs_vf(momxb + 1)%sf(j, k, l) + &
                                                        rhom(j, k, l)*accel_bf(2)
                        rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
                                                    q_cons_vf(momxb + 1)%sf(j, k, l)*accel_bf(2)
                    end do
                end do
            end do
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc end parallel loop
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp end target teams loop
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 154 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
        end if
 
        if (bf_z) then ! z-direction body forces
 
 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l) 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(j, k, l) 
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 159 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
            do l = 0, p
                do k = 0, n
                    do j = 0, m
                        rhs_vf(momxe)%sf(j, k, l) = rhs_vf(momxe)%sf(j, k, l) + &
                                                    rhom(j, k, l)*accel_bf(3)
                        rhs_vf(e_idx)%sf(j, k, l) = rhs_vf(e_idx)%sf(j, k, l) + &
                                                    q_cons_vf(momxe)%sf(j, k, l)*accel_bf(3)
                    end do
                end do
            end do
 
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc end parallel loop
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp end target teams loop
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 170 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
 
        end if
 
    end subroutine s_compute_body_forces_rhs
 
    !> @brief Deallocates module variables used for body force computations.
    impure subroutine s_finalize_body_forces_module
 
#ifdef MFC_DEBUG
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    block
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        use iso_fortran_env, only: output_unit
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        print *, 'm_body_forces.fpp:179: ', '@:DEALLOCATE(rhoM)'
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
        call flush (output_unit)
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    end block
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#if defined(MFC_OpenACC)
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$acc exit data delete(rhoM) 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#elif defined(MFC_OpenMP)
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
!$omp target exit data map(release:rhoM) 
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
#endif
# 179 "/home/runner/work/MFC/MFC/src/simulation/m_body_forces.fpp"
    deallocate (rhom)
 
    end subroutine s_finalize_body_forces_module
 
end module m_body_forces