m_bubbles_EL.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!>
!! @file
!! @brief Contains module @ref m_bubbles_el "m_bubbles_EL"
 
# 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_bubbles_EL.fpp" 2
 
!> @brief Tracks Lagrangian bubbles and couples their dynamics to the Eulerian flow via volume averaging
module m_bubbles_el
 
    use m_global_parameters             !< definitions of the global parameters
 
    use m_mpi_proxy                     !< message passing interface (mpi) module proxy
 
    use m_bubbles_el_kernels            !< definitions of the kernel functions
 
    use m_bubbles                       !< general bubble dynamics procedures
 
    use m_variables_conversion          !< state variables type conversion procedures
 
    use m_compile_specific
 
    use m_boundary_common
 
    use m_helper_basic         !< functions to compare floating point numbers
 
    use m_sim_helpers
 
    use m_helper
 
    implicit none
 
    !(nBub)
    integer, allocatable, dimension(:, :) :: lag_id                 !< Global and local IDs
    real(wp), allocatable, dimension(:) :: bub_r0            !< Initial bubble radius
    real(wp), allocatable, dimension(:) :: rmax_stats        !< Maximum radius
    real(wp), allocatable, dimension(:) :: rmin_stats        !< Minimum radius
 
# 37 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 37 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(lag_id, bub_R0, Rmax_stats, Rmin_stats) 
# 37 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 37 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (lag_id, bub_R0, Rmax_stats, Rmin_stats) 
# 37 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
    real(wp), allocatable, dimension(:) :: gas_mg            !< Bubble's gas mass
    real(wp), allocatable, dimension(:) :: gas_betat         !< heatflux model (Preston et al., 2007)
    real(wp), allocatable, dimension(:) :: gas_betac         !< massflux model (Preston et al., 2007)
    real(wp), allocatable, dimension(:) :: bub_dphidt        !< subgrid velocity potential (Maeda & Colonius, 2018)
 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(gas_mg, gas_betaT, gas_betaC, bub_dphidt) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (gas_mg, gas_betaT, gas_betaC, bub_dphidt) 
# 43 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
    !(nBub, 1 -> actual val or 2 -> temp val)
    real(wp), allocatable, dimension(:, :) :: gas_p          !< Pressure in the bubble
    real(wp), allocatable, dimension(:, :) :: gas_mv         !< Vapor mass in the bubble
    real(wp), allocatable, dimension(:, :) :: intfc_rad      !< Bubble radius
    real(wp), allocatable, dimension(:, :) :: intfc_vel      !< Velocity of the bubble interface
 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(gas_p, gas_mv, intfc_rad, intfc_vel) 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (gas_p, gas_mv, intfc_rad, intfc_vel) 
# 50 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
    !(nBub, 1-> x or 2->y or 3 ->z, 1 -> actual or 2 -> temporal val)
    real(wp), allocatable, dimension(:, :, :) :: mtn_pos     !< Bubble's position
    real(wp), allocatable, dimension(:, :, :) :: mtn_posprev !< Bubble's previous position
    real(wp), allocatable, dimension(:, :, :) :: mtn_vel     !< Bubble's velocity
    real(wp), allocatable, dimension(:, :, :) :: mtn_s       !< Bubble's computational cell position in real format
 
# 56 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 56 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(mtn_pos, mtn_posPrev, mtn_vel, mtn_s) 
# 56 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 56 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (mtn_pos, mtn_posPrev, mtn_vel, mtn_s) 
# 56 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
    !(nBub, 1-> x or 2->y or 3 ->z, time-stage)
    real(wp), allocatable, dimension(:, :) :: intfc_draddt   !< Time derivative of bubble's radius
    real(wp), allocatable, dimension(:, :) :: intfc_dveldt   !< Time derivative of bubble's interface velocity
    real(wp), allocatable, dimension(:, :) :: gas_dpdt       !< Time derivative of gas pressure
    real(wp), allocatable, dimension(:, :) :: gas_dmvdt      !< Time derivative of the vapor mass in the bubble
    real(wp), allocatable, dimension(:, :, :) :: mtn_dposdt  !< Time derivative of the bubble's position
    real(wp), allocatable, dimension(:, :, :) :: mtn_dveldt  !< Time derivative of the bubble's velocity
 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(intfc_draddt, intfc_dveldt, gas_dpdt, gas_dmvdt, mtn_dposdt, mtn_dveldt) 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (intfc_draddt, intfc_dveldt, gas_dpdt, gas_dmvdt, mtn_dposdt, mtn_dveldt) 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
    integer, private :: lag_num_ts                                  !<  Number of time stages in the time-stepping scheme
 
 
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(lag_num_ts) 
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (lag_num_ts) 
# 68 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
    integer :: nbubs                            !< Number of bubbles in the local domain
    real(wp) :: rmax_glb, rmin_glb       !< Maximum and minimum bubbe size in the local domain
    !< Projection of the lagrangian particles in the Eulerian framework
    type(scalar_field), dimension(:), allocatable :: q_beta
    integer :: q_beta_idx                       !< Size of the q_beta vector field
 
 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc declare create(nBubs, Rmax_glb, Rmin_glb, q_beta, q_beta_idx) 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target (nBubs, Rmax_glb, Rmin_glb, q_beta, q_beta_idx) 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
contains
 
    !> Initializes the lagrangian subgrid bubble solver
        !! @param q_cons_vf Initial conservative variables
    impure subroutine s_initialize_bubbles_el_module(q_cons_vf)
 
        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
 
        integer :: nbubs_glb, i
 
        ! Setting number of time-stages for selected time-stepping scheme
        lag_num_ts = time_stepper
 
        ! Allocate space for the Eulerian fields needed to map the effect of the bubbles
        if (lag_params%solver_approach == 1) then
            ! One-way coupling
            q_beta_idx = 3
        elseif (lag_params%solver_approach == 2) then
            ! Two-way coupling
            q_beta_idx = 4
            if (p == 0) then
                !Subgrid noise model for 2D approximation
                q_beta_idx = 6
            end if
        else
            call s_mpi_abort('Please check the lag_params%solver_approach input')
        end if
 
 
# 106 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 106 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(lag_num_ts, q_beta_idx) 
# 106 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 106 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(lag_num_ts, q_beta_idx) 
# 106 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
#ifdef MFC_DEBUG
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:108: ', '@:ALLOCATE(q_beta(1:q_beta_idx))'
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (q_beta(1:q_beta_idx))
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(q_beta) 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:q_beta) 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
        do i = 1, q_beta_idx
#ifdef MFC_DEBUG
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:111: ', '@:ALLOCATE(q_beta(i)%sf(idwbuff(1)%beg:idwbuff(1)%end, idwbuff(2)%beg:idwbuff(2)%end, idwbuff(3)%beg:idwbuff(3)%end))'
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (q_beta(i)%sf(idwbuff(1)%beg:idwbuff(1)%end,                 idwbuff(2)%beg:idwbuff(2)%end,                 idwbuff(3)%beg:idwbuff(3)%end))
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(q_beta(i)%sf) 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:q_beta(i)%sf) 
# 111 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 114 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        end do
 
        do i = 1, q_beta_idx
#ifdef _CRAYFTN
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#ifdef MFC_DEBUG
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:117: ', '@:ACC_SETUP_SFs(q_beta(i))'
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data copyin(q_beta(i)) 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(to:q_beta(i)) 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
            if (associated(q_beta(i)%sf)) then
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data copyin(q_beta(i)%sf) 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(to:q_beta(i)%sf) 
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
            end if
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 117 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
        end do
 
        ! Allocating space for lagrangian variables
        nbubs_glb = lag_params%nBubs_glb
 
#ifdef MFC_DEBUG
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:123: ', '@:ALLOCATE(lag_id(1:nBubs_glb, 1:2))'
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (lag_id(1:nbubs_glb, 1:2))
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(lag_id) 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:lag_id) 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:124: ', '@:ALLOCATE(bub_R0(1:nBubs_glb))'
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (bub_r0(1:nbubs_glb))
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(bub_R0) 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:bub_R0) 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:125: ', '@:ALLOCATE(Rmax_stats(1:nBubs_glb))'
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (rmax_stats(1:nbubs_glb))
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(Rmax_stats) 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:Rmax_stats) 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:126: ', '@:ALLOCATE(Rmin_stats(1:nBubs_glb))'
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (rmin_stats(1:nbubs_glb))
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(Rmin_stats) 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:Rmin_stats) 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:127: ', '@:ALLOCATE(gas_mg(1:nBubs_glb))'
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_mg(1:nbubs_glb))
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_mg) 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_mg) 
# 127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:128: ', '@:ALLOCATE(gas_betaT(1:nBubs_glb))'
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_betat(1:nbubs_glb))
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_betaT) 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_betaT) 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:129: ', '@:ALLOCATE(gas_betaC(1:nBubs_glb))'
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_betac(1:nbubs_glb))
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_betaC) 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_betaC) 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:130: ', '@:ALLOCATE(bub_dphidt(1:nBubs_glb))'
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (bub_dphidt(1:nbubs_glb))
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(bub_dphidt) 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:bub_dphidt) 
# 130 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:131: ', '@:ALLOCATE(gas_p(1:nBubs_glb, 1:2))'
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_p(1:nbubs_glb, 1:2))
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_p) 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_p) 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:132: ', '@:ALLOCATE(gas_mv(1:nBubs_glb, 1:2))'
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_mv(1:nbubs_glb, 1:2))
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_mv) 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_mv) 
# 132 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:133: ', '@:ALLOCATE(intfc_rad(1:nBubs_glb, 1:2))'
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (intfc_rad(1:nbubs_glb, 1:2))
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(intfc_rad) 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:intfc_rad) 
# 133 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:134: ', '@:ALLOCATE(intfc_vel(1:nBubs_glb, 1:2))'
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (intfc_vel(1:nbubs_glb, 1:2))
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(intfc_vel) 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:intfc_vel) 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:135: ', '@:ALLOCATE(mtn_pos(1:nBubs_glb, 1:3, 1:2))'
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_pos(1:nbubs_glb, 1:3, 1:2))
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_pos) 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_pos) 
# 135 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:136: ', '@:ALLOCATE(mtn_posPrev(1:nBubs_glb, 1:3, 1:2))'
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_posprev(1:nbubs_glb, 1:3, 1:2))
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_posPrev) 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_posPrev) 
# 136 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:137: ', '@:ALLOCATE(mtn_vel(1:nBubs_glb, 1:3, 1:2))'
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_vel(1:nbubs_glb, 1:3, 1:2))
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_vel) 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_vel) 
# 137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:138: ', '@:ALLOCATE(mtn_s(1:nBubs_glb, 1:3, 1:2))'
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_s(1:nbubs_glb, 1:3, 1:2))
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_s) 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_s) 
# 138 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:139: ', '@:ALLOCATE(intfc_draddt(1:nBubs_glb, 1:lag_num_ts))'
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (intfc_draddt(1:nbubs_glb, 1:lag_num_ts))
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(intfc_draddt) 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:intfc_draddt) 
# 139 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:140: ', '@:ALLOCATE(intfc_dveldt(1:nBubs_glb, 1:lag_num_ts))'
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (intfc_dveldt(1:nbubs_glb, 1:lag_num_ts))
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(intfc_dveldt) 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:intfc_dveldt) 
# 140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:141: ', '@:ALLOCATE(gas_dpdt(1:nBubs_glb, 1:lag_num_ts))'
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_dpdt(1:nbubs_glb, 1:lag_num_ts))
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_dpdt) 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_dpdt) 
# 141 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:142: ', '@:ALLOCATE(gas_dmvdt(1:nBubs_glb, 1:lag_num_ts))'
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (gas_dmvdt(1:nbubs_glb, 1:lag_num_ts))
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(gas_dmvdt) 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:gas_dmvdt) 
# 142 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:143: ', '@:ALLOCATE(mtn_dposdt(1:nBubs_glb, 1:3, 1:lag_num_ts))'
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_dposdt(1:nbubs_glb, 1:3, 1:lag_num_ts))
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_dposdt) 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_dposdt) 
# 143 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
#ifdef MFC_DEBUG
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:144: ', '@:ALLOCATE(mtn_dveldt(1:nBubs_glb, 1:3, 1:lag_num_ts))'
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    allocate (mtn_dveldt(1:nbubs_glb, 1:3, 1:lag_num_ts))
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc enter data create(mtn_dveldt) 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target enter data map(always,alloc:mtn_dveldt) 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
        if (adap_dt .and. f_is_default(adap_dt_tol)) adap_dt_tol = dflt_adap_dt_tol
 
        ! Starting bubbles
        call s_read_input_bubbles(q_cons_vf)
 
    end subroutine s_initialize_bubbles_el_module
 
    !> The purpose of this procedure is to obtain the initial bubbles' information
        !! @param q_cons_vf Conservative variables
    impure subroutine s_read_input_bubbles(q_cons_vf)
 
        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
 
        real(wp), dimension(8) :: inputbubble
        real(wp) :: qtime
        integer :: id, bub_id, save_count
        integer :: i, ios
        logical :: file_exist, indomain
 
        character(LEN=path_len + 2*name_len) :: path_d_dir !<
 
        ! Initialize number of particles
        bub_id = 0
        id = 0
 
        ! Read the input lag_bubble file or restart point
        if (cfl_dt) then
            save_count = n_start
            qtime = n_start*t_save
        else
            save_count = t_step_start
            qtime = t_step_start*dt
        end if
 
        if (save_count == 0) then
            if (proc_rank == 0) print *, 'Reading lagrange bubbles input file.'
            inquire (file='input/lag_bubbles.dat', exist=file_exist)
            if (file_exist) then
                open (94, file='input/lag_bubbles.dat', form='formatted', iostat=ios)
                do while (ios == 0)
                    read (94, *, iostat=ios) (inputbubble(i), i=1, 8)
                    if (ios /= 0) cycle
                    indomain = particle_in_domain_physical(inputbubble(1:3))
                    id = id + 1
                    if (id > lag_params%nBubs_glb .and. proc_rank == 0) then
                        call s_mpi_abort("Current number of bubbles is larger than nBubs_glb")
                    end if
                    if (indomain) then
                        bub_id = bub_id + 1
                        call s_add_bubbles(inputbubble, q_cons_vf, bub_id)
                        lag_id(bub_id, 1) = id      !global ID
                        lag_id(bub_id, 2) = bub_id  !local ID
                        nbubs = bub_id              ! local number of bubbles
                    end if
                end do
                close (94)
            else
                call s_mpi_abort("Initialize the lagrange bubbles in input/lag_bubbles.dat")
            end if
        else
            if (proc_rank == 0) print *, 'Restarting lagrange bubbles at save_count: ', save_count
            call s_restart_bubbles(bub_id, save_count)
        end if
 
        print *, " Lagrange bubbles running, in proc", proc_rank, "number:", bub_id, "/", id
 
 
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(bubbles_lagrange, lag_params) 
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(bubbles_lagrange, lag_params) 
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(lag_id, bub_R0, Rmax_stats, Rmin_stats, gas_mg, gas_betaT, gas_betaC, bub_dphidt, gas_p, gas_mv, intfc_rad, intfc_vel, mtn_pos, mtn_posPrev, mtn_vel, mtn_s, intfc_draddt, intfc_dveldt, gas_dpdt, gas_dmvdt, mtn_dposdt, mtn_dveldt, nBubs) 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(lag_id, bub_R0, Rmax_stats, Rmin_stats, gas_mg, gas_betaT, gas_betaC, bub_dphidt, gas_p, gas_mv, intfc_rad, intfc_vel, mtn_pos, mtn_posPrev, mtn_vel, mtn_s, intfc_draddt, intfc_dveldt, gas_dpdt, gas_dmvdt, mtn_dposdt, mtn_dveldt, nBubs) 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 219 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        rmax_glb = min(dflt_real, -dflt_real)
        rmin_glb = max(dflt_real, -dflt_real)
 
# 222 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 222 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(Rmax_glb, Rmin_glb) 
# 222 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 222 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(Rmax_glb, Rmin_glb) 
# 222 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
 
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(dx, dy, dz, x_cb, x_cc, y_cb, y_cc, z_cb, z_cc) 
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(dx, dy, dz, x_cb, x_cc, y_cb, y_cc, z_cb, z_cc) 
# 224 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
        !Populate temporal variables
        call s_transfer_data_to_tmp()
        call s_smear_voidfraction()
 
        if (lag_params%write_bubbles) call s_write_lag_particles(qtime)
 
        if (save_count == 0) then
            ! Create ./D directory
            write (path_d_dir, '(A,I0,A,I0)') trim(case_dir)//'/D'
            call my_inquire(path_d_dir, file_exist)
            if (.not. file_exist) call s_create_directory(trim(path_d_dir))
            call s_write_restart_lag_bubbles(save_count) ! Needed for post_processing
            call s_write_void_evol(qtime)
        end if
 
    end subroutine s_read_input_bubbles
 
    !> The purpose of this procedure is to obtain the information of the bubbles when starting fresh
        !! @param inputBubble Bubble information
        !! @param q_cons_vf Conservative variables
        !! @param bub_id Local id of the bubble
    impure subroutine s_add_bubbles(inputBubble, q_cons_vf, bub_id)
 
        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
        real(wp), dimension(8), intent(in) :: inputbubble
        integer, intent(in) :: bub_id
        integer :: i
 
        real(wp) :: pliq, volparticle, concvap, totalmass, kparticle, cpparticle
        real(wp) :: omegan_local, peg, pet, rhol, pcrit, qv, gamma, pi_inf, dynp
        integer, dimension(3) :: cell
        real(wp), dimension(2) :: re
        real(wp) :: massflag, heatflag, re_trans, im_trans
 
        massflag = 0._wp
        heatflag = 0._wp
        if (lag_params%massTransfer_model) massflag = 1._wp
        if (lag_params%heatTransfer_model) heatflag = 1._wp
 
        bub_r0(bub_id) = inputbubble(7)
        rmax_stats(bub_id) = min(dflt_real, -dflt_real)
        rmin_stats(bub_id) = max(dflt_real, -dflt_real)
        bub_dphidt(bub_id) = 0._wp
        intfc_rad(bub_id, 1) = inputbubble(7)
        intfc_vel(bub_id, 1) = inputbubble(8)
        mtn_pos(bub_id, 1:3, 1) = inputbubble(1:3)
        mtn_posprev(bub_id, 1:3, 1) = mtn_pos(bub_id, 1:3, 1)
        mtn_vel(bub_id, 1:3, 1) = inputbubble(4:6)
 
        if (cyl_coord .and. p == 0) then
            mtn_pos(bub_id, 2, 1) = sqrt(mtn_pos(bub_id, 2, 1)**2._wp + &
                                         mtn_pos(bub_id, 3, 1)**2._wp)
            !Storing azimuthal angle (-Pi to Pi)) into the third coordinate variable
            mtn_pos(bub_id, 3, 1) = atan2(inputbubble(3), inputbubble(2))
            mtn_posprev(bub_id, 1:3, 1) = mtn_pos(bub_id, 1:3, 1)
        end if
 
        cell = -buff_size
        call s_locate_cell(mtn_pos(bub_id, 1:3, 1), cell, mtn_s(bub_id, 1:3, 1))
 
        ! Check if the bubble is located in the ghost cell of a symmetric, or wall boundary
        if ((any(bc_x%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(1) < 0) .or. &
            (any(bc_x%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(1) > m) .or. &
            (any(bc_y%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(2) < 0) .or. &
            (any(bc_y%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(2) > n)) then
            call s_mpi_abort("Lagrange bubble is in the ghost cells of a symmetric or wall boundary.")
        end if
 
        if (p > 0) then
            if ((any(bc_z%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(3) < 0) .or. &
                (any(bc_z%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. cell(3) > p)) then
                call s_mpi_abort("Lagrange bubble is in the ghost cells of a symmetric or wall boundary.")
            end if
        end if
 
        call s_convert_to_mixture_variables(q_cons_vf, cell(1), cell(2), cell(3), &
                                            rhol, gamma, pi_inf, qv, re)
        dynp = 0._wp
        do i = 1, num_dims
            dynp = dynp + 0.5_wp*q_cons_vf(contxe + i)%sf(cell(1), cell(2), cell(3))**2/rhol
        end do
        pliq = (q_cons_vf(e_idx)%sf(cell(1), cell(2), cell(3)) - dynp - pi_inf)/gamma
        if (pliq < 0) print *, "Negative pressure", proc_rank, &
            q_cons_vf(e_idx)%sf(cell(1), cell(2), cell(3)), pi_inf, gamma, pliq, cell, dynp
 
        ! Initial particle pressure
        gas_p(bub_id, 1) = pliq + 2._wp*(1._wp/web)/bub_r0(bub_id)
        if (.not. f_approx_equal((1._wp/web), 0._wp)) then
            pcrit = pv - 4._wp*(1._wp/web)/(3._wp*sqrt(3._wp*gas_p(bub_id, 1)*bub_r0(bub_id)**3._wp/(2._wp*(1._wp/web))))
            pref = gas_p(bub_id, 1)
        else
            pcrit = 0._wp
        end if
 
        ! Initial particle mass
        volparticle = 4._wp/3._wp*pi*bub_r0(bub_id)**3._wp ! volume
        gas_mv(bub_id, 1) = pv*volparticle*(1._wp/(r_v*tw))*(massflag) ! vapermass
        gas_mg(bub_id) = (gas_p(bub_id, 1) - pv*(massflag))*volparticle*(1._wp/(r_g*tw)) ! gasmass
        if (gas_mg(bub_id) <= 0._wp) then
            call s_mpi_abort("The initial mass of gas inside the bubble is negative. Check the initial conditions.")
        end if
        totalmass = gas_mg(bub_id) + gas_mv(bub_id, 1) ! totalmass
 
        ! Bubble natural frequency
        concvap = gas_mv(bub_id, 1)/(gas_mv(bub_id, 1) + gas_mg(bub_id))
        omegan_local = (3._wp*(gas_p(bub_id, 1) - pv*(massflag)) + 4._wp*(1._wp/web)/bub_r0(bub_id))/rhol
        if (pv*(massflag) > gas_p(bub_id, 1)) then
            call s_mpi_abort("Lagrange bubble initially located in a region with pressure below the vapor pressure.")
        end if
        omegan_local = sqrt(omegan_local/bub_r0(bub_id)**2._wp)
 
        cpparticle = concvap*cp_v + (1._wp - concvap)*cp_g
        kparticle = concvap*k_vl + (1._wp - concvap)*k_gl
 
        ! Mass and heat transfer coefficients (based on Preston 2007)
        pet = totalmass/volparticle*cpparticle*bub_r0(bub_id)**2._wp*omegan_local/kparticle
        call s_transcoeff(1._wp, pet, re_trans, im_trans)
        gas_betat(bub_id) = re_trans*(heatflag)*kparticle
 
        peg = bub_r0(bub_id)**2._wp*omegan_local/vd
        call s_transcoeff(1._wp, peg, re_trans, im_trans)
        gas_betac(bub_id) = re_trans*(massflag)*vd
 
        if (gas_mg(bub_id) <= 0._wp) then
            call s_mpi_abort("Negative gas mass in the bubble, check if the bubble is in the domain.")
        end if
 
    end subroutine s_add_bubbles
 
    !> The purpose of this procedure is to obtain the information of the bubbles from a restart point.
        !! @param bub_id Local ID of the particle
        !! @param save_count File identifier
    impure subroutine s_restart_bubbles(bub_id, save_count)
 
        integer, intent(inout) :: bub_id, save_count
 
        character(LEN=path_len + 2*name_len) :: file_loc
        real(wp) :: file_time, file_dt
        integer :: file_num_procs, file_tot_part, tot_part
 
#ifdef MFC_MPI
        real(wp), dimension(20) :: inputvals
        integer, dimension(MPI_STATUS_SIZE) :: status
        integer(kind=MPI_OFFSET_KIND) :: disp
        integer :: view
 
        integer, dimension(3) :: cell
        logical :: indomain, particle_file, file_exist
 
        integer, dimension(2) :: gsizes, lsizes, start_idx_part
        integer :: ifile, ierr, tot_data, id
        integer :: i
 
        integer, dimension(:), allocatable :: proc_bubble_counts
        real(wp), dimension(1:1, 1:lag_io_vars) :: dummy
        dummy = 0._wp
 
        ! Construct file path
        write (file_loc, '(A,I0,A)') 'lag_bubbles_', save_count, '.dat'
        file_loc = trim(case_dir)//'/restart_data'//trim(mpiiofs)//trim(file_loc)
 
        ! Check if file exists
        inquire (file=trim(file_loc), exist=file_exist)
        if (.not. file_exist) then
            call s_mpi_abort('Restart file '//trim(file_loc)//' does not exist!')
        end if
 
        if (.not. parallel_io) return
 
        if (proc_rank == 0) then
            call mpi_file_open(mpi_comm_self, file_loc, mpi_mode_rdonly, &
                               mpi_info_int, ifile, ierr)
 
            call mpi_file_read(ifile, file_tot_part, 1, mpi_integer, status, ierr)
            call mpi_file_read(ifile, file_time, 1, mpi_p, status, ierr)
            call mpi_file_read(ifile, file_dt, 1, mpi_p, status, ierr)
            call mpi_file_read(ifile, file_num_procs, 1, mpi_integer, status, ierr)
 
            call mpi_file_close(ifile, ierr)
        end if
 
        call mpi_bcast(file_tot_part, 1, mpi_integer, 0, mpi_comm_world, ierr)
        call mpi_bcast(file_time, 1, mpi_p, 0, mpi_comm_world, ierr)
        call mpi_bcast(file_dt, 1, mpi_p, 0, mpi_comm_world, ierr)
        call mpi_bcast(file_num_procs, 1, mpi_integer, 0, mpi_comm_world, ierr)
 
        allocate (proc_bubble_counts(file_num_procs))
 
        if (proc_rank == 0) then
            call mpi_file_open(mpi_comm_self, file_loc, mpi_mode_rdonly, &
                               mpi_info_int, ifile, ierr)
 
            ! Skip to processor counts position
            disp = int(sizeof(file_tot_part) + 2*sizeof(file_time) + sizeof(file_num_procs), &
                       mpi_offset_kind)
            call mpi_file_seek(ifile, disp, mpi_seek_set, ierr)
            call mpi_file_read(ifile, proc_bubble_counts, file_num_procs, mpi_integer, status, ierr)
 
            call mpi_file_close(ifile, ierr)
        end if
 
        call mpi_bcast(proc_bubble_counts, file_num_procs, mpi_integer, 0, mpi_comm_world, ierr)
 
        ! Set time variables from file
        mytime = file_time
        dt = file_dt
 
        bub_id = proc_bubble_counts(proc_rank + 1)
 
        start_idx_part(1) = 0
        do i = 1, proc_rank
            start_idx_part(1) = start_idx_part(1) + proc_bubble_counts(i)
        end do
 
        start_idx_part(2) = 0
        lsizes(1) = bub_id
        lsizes(2) = lag_io_vars
 
        gsizes(1) = file_tot_part
        gsizes(2) = lag_io_vars
 
        if (bub_id > 0) then
 
            allocate (mpi_io_data_lag_bubbles(bub_id, 1:lag_io_vars))
 
            call mpi_type_create_subarray(2, gsizes, lsizes, start_idx_part, &
                                          mpi_order_fortran, mpi_p, view, ierr)
            call mpi_type_commit(view, ierr)
 
            call mpi_file_open(mpi_comm_world, file_loc, mpi_mode_rdonly, &
                               mpi_info_int, ifile, ierr)
 
            ! Skip extended header
            disp = int(sizeof(file_tot_part) + 2*sizeof(file_time) + sizeof(file_num_procs) + &
                       file_num_procs*sizeof(proc_bubble_counts(1)), mpi_offset_kind)
            call mpi_file_set_view(ifile, disp, mpi_p, view, 'native', mpi_info_int, ierr)
 
            call mpi_file_read_all(ifile, mpi_io_data_lag_bubbles, &
                                   lag_io_vars*bub_id, mpi_p, status, ierr)
 
            call mpi_file_close(ifile, ierr)
            call mpi_type_free(view, ierr)
 
            nbubs = bub_id
 
            do i = 1, bub_id
                lag_id(i, 1) = int(mpi_io_data_lag_bubbles(i, 1))
                mtn_pos(i, 1:3, 1) = mpi_io_data_lag_bubbles(i, 2:4)
                mtn_posprev(i, 1:3, 1) = mpi_io_data_lag_bubbles(i, 5:7)
                mtn_vel(i, 1:3, 1) = mpi_io_data_lag_bubbles(i, 8:10)
                intfc_rad(i, 1) = mpi_io_data_lag_bubbles(i, 11)
                intfc_vel(i, 1) = mpi_io_data_lag_bubbles(i, 12)
                bub_r0(i) = mpi_io_data_lag_bubbles(i, 13)
                rmax_stats(i) = mpi_io_data_lag_bubbles(i, 14)
                rmin_stats(i) = mpi_io_data_lag_bubbles(i, 15)
                bub_dphidt(i) = mpi_io_data_lag_bubbles(i, 16)
                gas_p(i, 1) = mpi_io_data_lag_bubbles(i, 17)
                gas_mv(i, 1) = mpi_io_data_lag_bubbles(i, 18)
                gas_mg(i) = mpi_io_data_lag_bubbles(i, 19)
                gas_betat(i) = mpi_io_data_lag_bubbles(i, 20)
                gas_betac(i) = mpi_io_data_lag_bubbles(i, 21)
                cell = -buff_size
                call s_locate_cell(mtn_pos(i, 1:3, 1), cell, mtn_s(i, 1:3, 1))
            end do
 
            deallocate (mpi_io_data_lag_bubbles)
 
        else
            nbubs = 0
 
            call mpi_type_contiguous(0, mpi_p, view, ierr)
            call mpi_type_commit(view, ierr)
 
            call mpi_file_open(mpi_comm_world, file_loc, mpi_mode_rdonly, &
                               mpi_info_int, ifile, ierr)
 
            ! Skip extended header
            disp = int(sizeof(file_tot_part) + 2*sizeof(file_time) + sizeof(file_num_procs) + &
                       file_num_procs*sizeof(proc_bubble_counts(1)), mpi_offset_kind)
            call mpi_file_set_view(ifile, disp, mpi_p, view, 'native', mpi_info_int, ierr)
 
            call mpi_file_read_all(ifile, dummy, 0, mpi_p, status, ierr)
 
            call mpi_file_close(ifile, ierr)
            call mpi_type_free(view, ierr)
        end if
 
        if (proc_rank == 0) then
            write (*, '(A,I0,A,I0)') 'Read ', file_tot_part, ' particles from restart file at t_step = ', save_count
            write (*, '(A,E15.7,A,E15.7)') 'Restart time = ', mytime, ', dt = ', dt
        end if
 
        deallocate (proc_bubble_counts)
#endif
 
    end subroutine s_restart_bubbles
 
    !>  Contains the bubble dynamics subroutines.
        !! @param q_prim_vf Primitive variables
        !! @param stage Current stage in the time-stepper algorithm
    subroutine s_compute_bubble_el_dynamics(q_prim_vf, stage)
 
        type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
        integer, intent(in) :: stage
 
        real(wp) :: myVapFlux
        real(wp) :: preterm1, term2, paux, pint, Romega, term1_fac
        real(wp) :: myR_m, mygamma_m, myPb, myMass_n, myMass_v
        real(wp) :: myR, myV, myBeta_c, myBeta_t, myR0, myPbdot, myMvdot
        real(wp) :: myPinf, aux1, aux2, myCson, myRho
        real(wp) :: gamma, pi_inf, qv
# 540 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
            real(wp), dimension(num_fluids) :: myalpha_rho, myalpha
# 542 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        real(wp), dimension(2) :: Re
        integer, dimension(3) :: cell
 
        integer :: adap_dt_stop_max, adap_dt_stop !< Fail-safe exit if max iteration count reached
        real(wp) :: dmalf, dmntait, dmBtait, dm_bub_adv_src, dm_divu !< Dummy variables for unified subgrid bubble subroutines
 
        integer :: i, k, l
 
        call nvtxstartrange("LAGRANGE-BUBBLE-DYNAMICS")
 
        ! Subgrid p_inf model based on Maeda and Colonius (2018).
        if (lag_params%pressure_corrector) then
            ! Calculate velocity potentials (valid for one bubble per cell)
 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k, cell, paux, preterm1, term2, Romega, myR0, myR, myV, myPb, pint, term1_fac) 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k, cell, paux, preterm1, term2, Romega, myR0, myR, myV, myPb, pint, term1_fac) 
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 555 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            do k = 1, nbubs
                call s_get_pinf(k, q_prim_vf, 2, paux, cell, preterm1, term2, romega)
                myr0 = bub_r0(k)
                myr = intfc_rad(k, 2)
                myv = intfc_vel(k, 2)
                mypb = gas_p(k, 2)
                pint = f_cpbw_km(myr0, myr, myv, mypb)
                pint = pint + 0.5_wp*myv**2._wp
                if (lag_params%cluster_type == 2) then
                    bub_dphidt(k) = (paux - pint) + term2
                    ! Accounting for the potential induced by the bubble averaged over the control volume
                    ! Note that this is based on the incompressible flow assumption near the bubble.
                    term1_fac = 3._wp/2._wp*(myr*(romega**2._wp - myr**2._wp))/(romega**3._wp - myr**3._wp)
                    bub_dphidt(k) = bub_dphidt(k)/(1._wp - term1_fac)
                end if
            end do
 
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 572 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        end if
 
        ! Radial motion model
        adap_dt_stop_max = 0
 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k, i, myalpha_rho, myalpha, Re, cell, myVapFlux, preterm1, term2, paux, pint, Romega, term1_fac, myR_m, mygamma_m, myPb, myMass_n, myMass_v, myR, myV, myBeta_c, myBeta_t, myR0, myPbdot, myMvdot, myPinf, aux1, aux2, myCson, myRho, gamma, pi_inf, qv, dmalf, dmntait, dmBtait, dm_bub_adv_src, dm_divu, adap_dt_stop) reduction(MAX:adap_dt_stop_max) copy(adap_dt_stop_max) copyin(stage) 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k, i, myalpha_rho, myalpha, Re, cell, myVapFlux, preterm1, term2, paux, pint, Romega, term1_fac, myR_m, mygamma_m, myPb, myMass_n, myMass_v, myR, myV, myBeta_c, myBeta_t, myR0, myPbdot, myMvdot, myPinf, aux1, aux2, myCson, myRho, gamma, pi_inf, qv, dmalf, dmntait, dmBtait, dm_bub_adv_src, dm_divu, adap_dt_stop) reduction(MAX:adap_dt_stop_max) map(tofrom:adap_dt_stop_max) map(to:stage) 
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 577 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 580 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        do k = 1, nbubs
            ! Keller-Miksis model
 
            ! Current bubble state
            mypb = gas_p(k, 2)
            mymass_n = gas_mg(k)
            mymass_v = gas_mv(k, 2)
            myr = intfc_rad(k, 2)
            myv = intfc_vel(k, 2)
            mybeta_c = gas_betac(k)
            mybeta_t = gas_betat(k)
            myr0 = bub_r0(k)
 
            ! Vapor and heat fluxes
            call s_vflux(myr, myv, mypb, mymass_v, k, myvapflux, mymass_n, mybeta_c, myr_m, mygamma_m)
            mypbdot = f_bpres_dot(myvapflux, myr, myv, mypb, mymass_v, k, mybeta_t, myr_m, mygamma_m)
            mymvdot = 4._wp*pi*myr**2._wp*myvapflux
 
            ! Obtaining driving pressure
            call s_get_pinf(k, q_prim_vf, 1, mypinf, cell, aux1, aux2)
 
            ! Obtain liquid density and computing speed of sound from pinf
            call s_compute_species_fraction(q_prim_vf, cell(1), cell(2), cell(3), myalpha_rho, myalpha)
            call s_convert_species_to_mixture_variables_acc(myrho, gamma, pi_inf, qv, myalpha, &
                                                            myalpha_rho, re)
            call s_compute_cson_from_pinf(q_prim_vf, mypinf, cell, myrho, gamma, pi_inf, mycson)
 
            ! Adaptive time stepping
            adap_dt_stop = 0
 
            if (adap_dt) then
 
                call s_advance_step(myrho, mypinf, myr, myv, myr0, mypb, mypbdot, dmalf, &
                                    dmntait, dmbtait, dm_bub_adv_src, dm_divu, &
                                    k, mymass_v, mymass_n, mybeta_c, &
                                    mybeta_t, mycson, adap_dt_stop)
 
                ! Update bubble state
                intfc_rad(k, 1) = myr
                intfc_vel(k, 1) = myv
                gas_p(k, 1) = mypb
                gas_mv(k, 1) = mymass_v
 
            else
 
                ! Radial acceleration from bubble models
                intfc_dveldt(k, stage) = f_rddot(myrho, mypinf, myr, myv, myr0, &
                                                 mypb, mypbdot, dmalf, dmntait, dmbtait, &
                                                 dm_bub_adv_src, dm_divu, &
                                                 mycson)
                intfc_draddt(k, stage) = myv
                gas_dmvdt(k, stage) = mymvdot
                gas_dpdt(k, stage) = mypbdot
 
            end if
 
            adap_dt_stop_max = max(adap_dt_stop_max, adap_dt_stop)
 
        end do
 
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 639 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
        if (adap_dt .and. adap_dt_stop_max > 0) call s_mpi_abort("Adaptive time stepping failed to converge.")
 
        ! Bubbles remain in a fixed position
 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(2) gang vector default(present) private(k, l) copyin(stage) 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(2) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k, l) map(to:stage) 
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 644 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        do k = 1, nbubs
            do l = 1, 3
                mtn_dposdt(k, l, stage) = 0._wp
                mtn_dveldt(k, l, stage) = 0._wp
            end do
        end do
 
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 651 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
        call nvtxendrange
 
    end subroutine s_compute_bubble_el_dynamics
 
    !>  The purpose of this subroutine is to obtain the bubble source terms based on Maeda and Colonius (2018)
        !!      and add them to the RHS scalar field.
        !! @param q_cons_vf Conservative variables
        !! @param q_prim_vf Conservative variables
        !! @param rhs_vf Time derivative of the conservative variables
    subroutine s_compute_bubbles_el_source(q_cons_vf, q_prim_vf, rhs_vf)
 
        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
        type(scalar_field), dimension(sys_size), intent(inout) :: q_prim_vf
        type(scalar_field), dimension(sys_size), intent(inout) :: rhs_vf
 
        integer :: i, j, k, l
 
        if (.not. adap_dt) call s_smear_voidfraction()
 
        if (lag_params%solver_approach == 2) then
 
            ! (q / (1 - beta)) * d(beta)/dt source
            if (p == 0) then
 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(4) gang vector default(present) private(i, j, k, l) 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.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) 
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 0, p
                    do j = 0, n
                        do i = 0, m
                            do l = 1, e_idx
                                if (q_beta(1)%sf(i, j, k) > (1._wp - lag_params%valmaxvoid)) then
                                    rhs_vf(l)%sf(i, j, k) = rhs_vf(l)%sf(i, j, k) + &
                                                            q_cons_vf(l)%sf(i, j, k)*(q_beta(2)%sf(i, j, k) + &
                                                                                      q_beta(5)%sf(i, j, k))
 
                                end if
                            end do
                        end do
                    end do
                end do
 
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 691 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            else
 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(4) gang vector default(present) private(i, j, k, l) 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.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) 
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 693 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 0, p
                    do j = 0, n
                        do i = 0, m
                            do l = 1, e_idx
                                if (q_beta(1)%sf(i, j, k) > (1._wp - lag_params%valmaxvoid)) then
                                    rhs_vf(l)%sf(i, j, k) = rhs_vf(l)%sf(i, j, k) + &
                                                            q_cons_vf(l)%sf(i, j, k)/q_beta(1)%sf(i, j, k)* &
                                                            q_beta(2)%sf(i, j, k)
                                end if
                            end do
                        end do
                    end do
                end do
 
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 707 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            end if
 
            do l = 1, num_dims
 
                call s_gradient_dir(q_prim_vf(e_idx)%sf, q_beta(3)%sf, l)
 
                ! (q / (1 - beta)) * d(beta)/dt source
 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 0, p
                    do j = 0, n
                        do i = 0, m
                            if (q_beta(1)%sf(i, j, k) > (1._wp - lag_params%valmaxvoid)) then
                                rhs_vf(contxe + l)%sf(i, j, k) = rhs_vf(contxe + l)%sf(i, j, k) - &
                                                                 (1._wp - q_beta(1)%sf(i, j, k))/ &
                                                                 q_beta(1)%sf(i, j, k)* &
                                                                 q_beta(3)%sf(i, j, k)
                            end if
                        end do
                    end do
                end do
 
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 728 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
                !source in energy
 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 731 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = idwbuff(3)%beg, idwbuff(3)%end
                    do j = idwbuff(2)%beg, idwbuff(2)%end
                        do i = idwbuff(1)%beg, idwbuff(1)%end
                            q_beta(3)%sf(i, j, k) = q_prim_vf(e_idx)%sf(i, j, k)*q_prim_vf(contxe + l)%sf(i, j, k)
                        end do
                    end do
                end do
 
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 739 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
                call s_gradient_dir(q_beta(3)%sf, q_beta(4)%sf, l)
 
                ! (beta / (1 - beta)) * d(Pu)/dl source
 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 744 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 0, p
                    do j = 0, n
                        do i = 0, m
                            if (q_beta(1)%sf(i, j, k) > (1._wp - lag_params%valmaxvoid)) then
                                rhs_vf(e_idx)%sf(i, j, k) = rhs_vf(e_idx)%sf(i, j, k) - &
                                                            q_beta(4)%sf(i, j, k)*(1._wp - q_beta(1)%sf(i, j, k))/ &
                                                            q_beta(1)%sf(i, j, k)
                            end if
                        end do
                    end do
                end do
 
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 756 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            end do
 
        end if
 
    end subroutine s_compute_bubbles_el_source
 
    !>  This procedure computes the speed of sound from a given driving pressure
        !! @param q_prim_vf Primitive variables
        !! @param pinf Driving pressure
        !! @param cell Bubble cell
        !! @param rhol Liquid density
        !! @param gamma Liquid specific heat ratio
        !! @param pi_inf Liquid stiffness
        !! @param cson Calculated speed of sound
    subroutine s_compute_cson_from_pinf(q_prim_vf, pinf, cell, rhol, gamma, pi_inf, cson)
 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#ifdef _CRAYFTN
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if MFC_OpenACC
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc routine seq 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenMP
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target device_type(any) 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#else
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!DIR$ INLINEALWAYS s_compute_cson_from_pinf
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenACC
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc routine seq 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenMP
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target device_type(any) 
# 772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
        real(wp), intent(in) :: pinf, rhol, gamma, pi_inf
        integer, dimension(3), intent(in) :: cell
        real(wp), intent(out) :: cson
 
        real(wp) :: E, H
# 784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
            real(wp), dimension(num_dims) :: vel
# 786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        integer :: i
 
        vel(:) = 0._wp
 
# 789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
        do i = 1, num_dims
            vel(i) = q_prim_vf(i + contxe)%sf(cell(1), cell(2), cell(3))
        end do
        e = gamma*pinf + pi_inf + 0.5_wp*rhol*dot_product(vel, vel)
        h = (e + pinf)/rhol
        cson = sqrt((h - 0.5_wp*dot_product(vel, vel))/gamma)
 
    end subroutine s_compute_cson_from_pinf
 
    !>  The purpose of this subroutine is to smear the effect of the bubbles in the Eulerian framework
    subroutine s_smear_voidfraction()
 
        integer :: i, j, k, l
 
        call nvtxstartrange("BUBBLES-LAGRANGE-KERNELS")
 
 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(4) gang vector default(present) private(i, j, k, l) 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.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) 
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 806 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        do i = 1, q_beta_idx
            do l = idwbuff(3)%beg, idwbuff(3)%end
                do k = idwbuff(2)%beg, idwbuff(2)%end
                    do j = idwbuff(1)%beg, idwbuff(1)%end
                        q_beta(i)%sf(j, k, l) = 0._wp
                    end do
                end do
            end do
        end do
 
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 816 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
        call s_smoothfunction(nbubs, intfc_rad, intfc_vel, &
                              mtn_s, mtn_pos, q_beta)
 
        !Store 1-beta
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(j, k, l) 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.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) 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        do l = idwbuff(3)%beg, idwbuff(3)%end
            do k = idwbuff(2)%beg, idwbuff(2)%end
                do j = idwbuff(1)%beg, idwbuff(1)%end
                    q_beta(1)%sf(j, k, l) = 1._wp - q_beta(1)%sf(j, k, l)
                    ! Limiting void fraction given max value
                    q_beta(1)%sf(j, k, l) = max(q_beta(1)%sf(j, k, l), &
                                                1._wp - lag_params%valmaxvoid)
                end do
            end do
        end do
 
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 833 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
        call nvtxendrange
 
    end subroutine s_smear_voidfraction
 
    !> The purpose of this procedure is obtain the bubble driving pressure p_inf
        !! @param bub_id Particle identifier
        !! @param q_prim_vf  Primitive variables
        !! @param ptype 1: p at infinity, 2: averaged P at the bubble location
        !! @param f_pinfl Driving pressure
        !! @param cell Bubble cell
        !! @param preterm1 Pre-computed term 1
        !! @param term2 Computed term 2
        !! @param Romega Control volume radius
    subroutine s_get_pinf(bub_id, q_prim_vf, ptype, f_pinfl, cell, preterm1, term2, Romega)
 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#ifdef _CRAYFTN
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if MFC_OpenACC
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc routine seq 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenMP
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target device_type(any) 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#else
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!DIR$ INLINEALWAYS s_get_pinf
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenACC
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc routine seq 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif MFC_OpenMP
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp declare target device_type(any) 
# 849 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 851 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        integer, intent(in) :: bub_id, ptype
        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
        real(wp), intent(out) :: f_pinfl
        integer, dimension(3), intent(out) :: cell
        real(wp), intent(out), optional :: preterm1, term2, Romega
 
        real(wp), dimension(3) :: scoord, psi
        real(wp) :: dc, vol, aux
        real(wp) :: volgas, term1, Rbeq, denom
        real(wp) :: charvol, charpres, charvol2, charpres2
        integer, dimension(3) :: cellaux
        integer :: i, j, k
        integer :: smearGrid, smearGridz
        logical :: celloutside
 
        scoord = mtn_s(bub_id, 1:3, 2)
        f_pinfl = 0._wp
 
        !< Find current bubble cell
        cell(:) = int(scoord(:))
 
# 872 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 872 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 872 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 872 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 872 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
        do i = 1, num_dims
            if (scoord(i) < 0._wp) cell(i) = cell(i) - 1
        end do
 
        if ((lag_params%cluster_type == 1)) then
            !< Getting p_cell in terms of only the current cell by interpolation
 
            !< Getting the cell volulme as Omega
            if (p > 0) then
                vol = dx(cell(1))*dy(cell(2))*dz(cell(3))
            else
                if (cyl_coord) then
                    vol = dx(cell(1))*dy(cell(2))*y_cc(cell(2))*2._wp*pi
                else
                    vol = dx(cell(1))*dy(cell(2))*lag_params%charwidth
                end if
            end if
 
            !< Obtain bilinear interpolation coefficients, based on the current location of the bubble.
            psi(1) = (scoord(1) - real(cell(1)))*dx(cell(1)) + x_cb(cell(1) - 1)
            if (cell(1) == (m + buff_size)) then
                cell(1) = cell(1) - 1
                psi(1) = 1._wp
            else if (cell(1) == (-buff_size)) then
                psi(1) = 0._wp
            else
                if (psi(1) < x_cc(cell(1))) cell(1) = cell(1) - 1
                psi(1) = abs((psi(1) - x_cc(cell(1)))/(x_cc(cell(1) + 1) - x_cc(cell(1))))
            end if
 
            psi(2) = (scoord(2) - real(cell(2)))*dy(cell(2)) + y_cb(cell(2) - 1)
            if (cell(2) == (n + buff_size)) then
                cell(2) = cell(2) - 1
                psi(2) = 1._wp
            else if (cell(2) == (-buff_size)) then
                psi(2) = 0._wp
            else
                if (psi(2) < y_cc(cell(2))) cell(2) = cell(2) - 1
                psi(2) = abs((psi(2) - y_cc(cell(2)))/(y_cc(cell(2) + 1) - y_cc(cell(2))))
            end if
 
            if (p > 0) then
                psi(3) = (scoord(3) - real(cell(3)))*dz(cell(3)) + z_cb(cell(3) - 1)
                if (cell(3) == (p + buff_size)) then
                    cell(3) = cell(3) - 1
                    psi(3) = 1._wp
                else if (cell(3) == (-buff_size)) then
                    psi(3) = 0._wp
                else
                    if (psi(3) < z_cc(cell(3))) cell(3) = cell(3) - 1
                    psi(3) = abs((psi(3) - z_cc(cell(3)))/(z_cc(cell(3) + 1) - z_cc(cell(3))))
                end if
            else
                psi(3) = 0._wp
            end if
 
            !< Perform bilinear interpolation
            if (p == 0) then  !2D
                f_pinfl = q_prim_vf(e_idx)%sf(cell(1), cell(2), cell(3))*(1._wp - psi(1))*(1._wp - psi(2))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2), cell(3))*psi(1)*(1._wp - psi(2))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2) + 1, cell(3))*psi(1)*psi(2)
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1), cell(2) + 1, cell(3))*(1._wp - psi(1))*psi(2)
            else              !3D
                f_pinfl = q_prim_vf(e_idx)%sf(cell(1), cell(2), cell(3))*(1._wp - psi(1))*(1._wp - psi(2))*(1._wp - psi(3))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2), cell(3))*psi(1)*(1._wp - psi(2))*(1._wp - psi(3))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2) + 1, cell(3))*psi(1)*psi(2)*(1._wp - psi(3))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1), cell(2) + 1, cell(3))*(1._wp - psi(1))*psi(2)*(1._wp - psi(3))
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1), cell(2), cell(3) + 1)*(1._wp - psi(1))*(1._wp - psi(2))*psi(3)
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2), cell(3) + 1)*psi(1)*(1._wp - psi(2))*psi(3)
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1) + 1, cell(2) + 1, cell(3) + 1)*psi(1)*psi(2)*psi(3)
                f_pinfl = f_pinfl + q_prim_vf(e_idx)%sf(cell(1), cell(2) + 1, cell(3) + 1)*(1._wp - psi(1))*psi(2)*psi(3)
            end if
 
            !R_Omega
            dc = (3._wp*vol/(4._wp*pi))**(1._wp/3._wp)
 
        else if (lag_params%cluster_type >= 2) then
            ! Bubble dynamic closure from Maeda and Colonius (2018)
 
            ! Include the cell that contains the bubble (mapCells+1+mapCells)
            smeargrid = mapcells - (-mapcells) + 1
            smeargridz = smeargrid
            if (p == 0) smeargridz = 1
 
            charvol = 0._wp
            charpres = 0._wp
            charvol2 = 0._wp
            charpres2 = 0._wp
            vol = 0._wp
 
 
# 963 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 963 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 963 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 963 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 963 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
            do i = 1, smeargrid
 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
                do j = 1, smeargrid
 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
                    do k = 1, smeargridz
                        cellaux(1) = cell(1) + i - (mapcells + 1)
                        cellaux(2) = cell(2) + j - (mapcells + 1)
                        cellaux(3) = cell(3) + k - (mapcells + 1)
                        if (p == 0) cellaux(3) = 0
 
                        !< check if the current cell is outside the computational domain or not (including ghost cells)
                        celloutside = .false.
                        if (num_dims == 2) then
                            if ((cellaux(1) < -buff_size) .or. (cellaux(2) < -buff_size)) then
                                celloutside = .true.
                            end if
                            if (cyl_coord .and. y_cc(cellaux(2)) < 0._wp) then
                                celloutside = .true.
                            end if
                            if ((cellaux(2) > n + buff_size) .or. (cellaux(1) > m + buff_size)) then
                                celloutside = .true.
                            end if
                        else
                            if ((cellaux(3) < -buff_size) .or. (cellaux(1) < -buff_size) .or. (cellaux(2) < -buff_size)) then
                                celloutside = .true.
                            end if
 
                            if ((cellaux(3) > p + buff_size) .or. (cellaux(2) > n + buff_size) .or. (cellaux(1) > m + buff_size)) then
                                celloutside = .true.
                            end if
                        end if
                        if (.not. celloutside) then
                            if (cyl_coord .and. (p == 0) .and. (y_cc(cellaux(2)) < 0._wp)) then
                                celloutside = .true.
                            end if
                        end if
 
                        if (.not. celloutside) then
                            !< Obtaining the cell volulme
                            if (p > 0) then
                                vol = dx(cellaux(1))*dy(cellaux(2))*dz(cellaux(3))
                            else
                                if (cyl_coord) then
                                    vol = dx(cellaux(1))*dy(cellaux(2))*y_cc(cellaux(2))*2._wp*pi
                                else
                                    vol = dx(cellaux(1))*dy(cellaux(2))*lag_params%charwidth
                                end if
                            end if
                            !< Update values
                            charvol = charvol + vol
                            charpres = charpres + q_prim_vf(e_idx)%sf(cellaux(1), cellaux(2), cellaux(3))*vol
                            charvol2 = charvol2 + vol*q_beta(1)%sf(cellaux(1), cellaux(2), cellaux(3))
                            charpres2 = charpres2 + q_prim_vf(e_idx)%sf(cellaux(1), cellaux(2), cellaux(3)) &
                                        *vol*q_beta(1)%sf(cellaux(1), cellaux(2), cellaux(3))
                        end if
 
                    end do
                end do
            end do
 
            f_pinfl = charpres2/charvol2
            vol = charvol
            dc = (3._wp*abs(vol)/(4._wp*pi))**(1._wp/3._wp)
 
        end if
 
        if (lag_params%pressure_corrector) then
 
            !Valid if only one bubble exists per cell
            volgas = intfc_rad(bub_id, 2)**3._wp
            denom = intfc_rad(bub_id, 2)**2._wp
            term1 = bub_dphidt(bub_id)*intfc_rad(bub_id, 2)**2._wp
            term2 = intfc_vel(bub_id, 2)*intfc_rad(bub_id, 2)**2._wp
 
            rbeq = volgas**(1._wp/3._wp) !surrogate bubble radius
            aux = dc**3._wp - rbeq**3._wp
            term2 = term2/denom
            term2 = 3._wp/2._wp*term2**2._wp*rbeq**3._wp*(1._wp - rbeq/dc)/aux
            preterm1 = 3._wp/2._wp*rbeq*(dc**2._wp - rbeq**2._wp)/(aux*denom)
 
            !Control volume radius
            if (ptype == 2) romega = dc
 
            ! Getting p_inf
            if (ptype == 1) then
                f_pinfl = f_pinfl + preterm1*term1 + term2
            end if
 
        end if
 
    end subroutine s_get_pinf
 
    !>  This subroutine updates the Lagrange variables using the tvd RK time steppers.
        !!      The time derivative of the bubble variables must be stored at every stage to avoid precision errors.
        !! @param stage Current tvd RK stage
    impure subroutine s_update_lagrange_tdv_rk(stage)
 
        integer, intent(in) :: stage
 
        integer :: k
 
        if (time_stepper == 1) then ! 1st order TVD RK
 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1066 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            do k = 1, nbubs
                !u{1} = u{n} +  dt * RHS{n}
                intfc_rad(k, 1) = intfc_rad(k, 1) + dt*intfc_draddt(k, 1)
                intfc_vel(k, 1) = intfc_vel(k, 1) + dt*intfc_dveldt(k, 1)
                mtn_pos(k, 1:3, 1) = mtn_pos(k, 1:3, 1) + dt*mtn_dposdt(k, 1:3, 1)
                mtn_vel(k, 1:3, 1) = mtn_vel(k, 1:3, 1) + dt*mtn_dveldt(k, 1:3, 1)
                gas_p(k, 1) = gas_p(k, 1) + dt*gas_dpdt(k, 1)
                gas_mv(k, 1) = gas_mv(k, 1) + dt*gas_dmvdt(k, 1)
            end do
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
            call s_transfer_data_to_tmp()
            call s_write_void_evol(mytime)
            if (lag_params%write_bubbles_stats) call s_calculate_lag_bubble_stats()
 
            if (lag_params%write_bubbles) then
 
# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update host(gas_p, gas_mv, intfc_rad, intfc_vel) 
# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update from(gas_p, gas_mv, intfc_rad, intfc_vel) 
# 1083 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
                call s_write_lag_particles(mytime)
            end if
 
        elseif (time_stepper == 2) then ! 2nd order TVD RK
            if (stage == 1) then
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 1, nbubs
                    !u{1} = u{n} +  dt * RHS{n}
                    intfc_rad(k, 2) = intfc_rad(k, 1) + dt*intfc_draddt(k, 1)
                    intfc_vel(k, 2) = intfc_vel(k, 1) + dt*intfc_dveldt(k, 1)
                    mtn_pos(k, 1:3, 2) = mtn_pos(k, 1:3, 1) + dt*mtn_dposdt(k, 1:3, 1)
                    mtn_vel(k, 1:3, 2) = mtn_vel(k, 1:3, 1) + dt*mtn_dveldt(k, 1:3, 1)
                    gas_p(k, 2) = gas_p(k, 1) + dt*gas_dpdt(k, 1)
                    gas_mv(k, 2) = gas_mv(k, 1) + dt*gas_dmvdt(k, 1)
                end do
 
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1099 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
            elseif (stage == 2) then
 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1102 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 1, nbubs
                    !u{1} = u{n} + (1/2) * dt * (RHS{n} + RHS{1})
                    intfc_rad(k, 1) = intfc_rad(k, 1) + dt*(intfc_draddt(k, 1) + intfc_draddt(k, 2))/2._wp
                    intfc_vel(k, 1) = intfc_vel(k, 1) + dt*(intfc_dveldt(k, 1) + intfc_dveldt(k, 2))/2._wp
                    mtn_pos(k, 1:3, 1) = mtn_pos(k, 1:3, 1) + dt*(mtn_dposdt(k, 1:3, 1) + mtn_dposdt(k, 1:3, 2))/2._wp
                    mtn_vel(k, 1:3, 1) = mtn_vel(k, 1:3, 1) + dt*(mtn_dveldt(k, 1:3, 1) + mtn_dveldt(k, 1:3, 2))/2._wp
                    gas_p(k, 1) = gas_p(k, 1) + dt*(gas_dpdt(k, 1) + gas_dpdt(k, 2))/2._wp
                    gas_mv(k, 1) = gas_mv(k, 1) + dt*(gas_dmvdt(k, 1) + gas_dmvdt(k, 2))/2._wp
                end do
 
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1112 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
                call s_transfer_data_to_tmp()
                call s_write_void_evol(mytime)
                if (lag_params%write_bubbles_stats) call s_calculate_lag_bubble_stats()
 
                if (lag_params%write_bubbles) then
 
# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update host(gas_p, gas_mv, intfc_rad, intfc_vel) 
# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update from(gas_p, gas_mv, intfc_rad, intfc_vel) 
# 1119 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
                    call s_write_lag_particles(mytime)
                end if
 
            end if
 
        elseif (time_stepper == 3) then ! 3rd order TVD RK
            if (stage == 1) then
 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1127 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 1, nbubs
                    !u{1} = u{n} +  dt * RHS{n}
                    intfc_rad(k, 2) = intfc_rad(k, 1) + dt*intfc_draddt(k, 1)
                    intfc_vel(k, 2) = intfc_vel(k, 1) + dt*intfc_dveldt(k, 1)
                    mtn_pos(k, 1:3, 2) = mtn_pos(k, 1:3, 1) + dt*mtn_dposdt(k, 1:3, 1)
                    mtn_vel(k, 1:3, 2) = mtn_vel(k, 1:3, 1) + dt*mtn_dveldt(k, 1:3, 1)
                    gas_p(k, 2) = gas_p(k, 1) + dt*gas_dpdt(k, 1)
                    gas_mv(k, 2) = gas_mv(k, 1) + dt*gas_dmvdt(k, 1)
                end do
 
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1137 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
            elseif (stage == 2) then
 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1140 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 1, nbubs
                    !u{2} = u{n} + (1/4) * dt * [RHS{n} + RHS{1}]
                    intfc_rad(k, 2) = intfc_rad(k, 1) + dt*(intfc_draddt(k, 1) + intfc_draddt(k, 2))/4._wp
                    intfc_vel(k, 2) = intfc_vel(k, 1) + dt*(intfc_dveldt(k, 1) + intfc_dveldt(k, 2))/4._wp
                    mtn_pos(k, 1:3, 2) = mtn_pos(k, 1:3, 1) + dt*(mtn_dposdt(k, 1:3, 1) + mtn_dposdt(k, 1:3, 2))/4._wp
                    mtn_vel(k, 1:3, 2) = mtn_vel(k, 1:3, 1) + dt*(mtn_dveldt(k, 1:3, 1) + mtn_dveldt(k, 1:3, 2))/4._wp
                    gas_p(k, 2) = gas_p(k, 1) + dt*(gas_dpdt(k, 1) + gas_dpdt(k, 2))/4._wp
                    gas_mv(k, 2) = gas_mv(k, 1) + dt*(gas_dmvdt(k, 1) + gas_dmvdt(k, 2))/4._wp
                end do
 
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1150 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            elseif (stage == 3) then
 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1152 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
                do k = 1, nbubs
                    !u{n+1} = u{n} + (2/3) * dt * [(1/4)* RHS{n} + (1/4)* RHS{1} + RHS{2}]
                    intfc_rad(k, 1) = intfc_rad(k, 1) + (2._wp/3._wp)*dt*(intfc_draddt(k, 1)/4._wp + intfc_draddt(k, 2)/4._wp + intfc_draddt(k, 3))
                    intfc_vel(k, 1) = intfc_vel(k, 1) + (2._wp/3._wp)*dt*(intfc_dveldt(k, 1)/4._wp + intfc_dveldt(k, 2)/4._wp + intfc_dveldt(k, 3))
                    mtn_pos(k, 1:3, 1) = mtn_pos(k, 1:3, 1) + (2._wp/3._wp)*dt*(mtn_dposdt(k, 1:3, 1)/4._wp + mtn_dposdt(k, 1:3, 2)/4._wp + mtn_dposdt(k, 1:3, 3))
                    mtn_vel(k, 1:3, 1) = mtn_vel(k, 1:3, 1) + (2._wp/3._wp)*dt*(mtn_dveldt(k, 1:3, 1)/4._wp + mtn_dveldt(k, 1:3, 2)/4._wp + mtn_dveldt(k, 1:3, 3))
                    gas_p(k, 1) = gas_p(k, 1) + (2._wp/3._wp)*dt*(gas_dpdt(k, 1)/4._wp + gas_dpdt(k, 2)/4._wp + gas_dpdt(k, 3))
                    gas_mv(k, 1) = gas_mv(k, 1) + (2._wp/3._wp)*dt*(gas_dmvdt(k, 1)/4._wp + gas_dmvdt(k, 2)/4._wp + gas_dmvdt(k, 3))
                end do
 
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1162 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
                call s_transfer_data_to_tmp()
                call s_write_void_evol(mytime)
                if (lag_params%write_bubbles_stats) call s_calculate_lag_bubble_stats()
 
                if (lag_params%write_bubbles) then
 
# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update host(gas_p, gas_mv, gas_mg, intfc_rad, intfc_vel) 
# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update from(gas_p, gas_mv, gas_mg, intfc_rad, intfc_vel) 
# 1169 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
                    call s_write_lag_particles(mytime)
                end if
 
            end if
 
        end if
 
    end subroutine s_update_lagrange_tdv_rk
 
    !> This subroutine returns the computational coordinate of the cell for the given position.
          !! @param pos Input coordinates
          !! @param cell Computational coordinate of the cell
          !! @param scoord Calculated particle coordinates
    subroutine s_locate_cell(pos, cell, scoord)
 
        real(wp), dimension(3), intent(in) :: pos
        real(wp), dimension(3), intent(out) :: scoord
        integer, dimension(3), intent(inout) :: cell
 
        integer :: i
 
        do while (pos(1) < x_cb(cell(1) - 1))
            cell(1) = cell(1) - 1
        end do
 
        do while (pos(1) > x_cb(cell(1)))
            cell(1) = cell(1) + 1
        end do
 
        do while (pos(2) < y_cb(cell(2) - 1))
            cell(2) = cell(2) - 1
        end do
 
        do while (pos(2) > y_cb(cell(2)))
            cell(2) = cell(2) + 1
        end do
 
        if (p > 0) then
            do while (pos(3) < z_cb(cell(3) - 1))
                cell(3) = cell(3) - 1
            end do
            do while (pos(3) > z_cb(cell(3)))
                cell(3) = cell(3) + 1
            end do
        end if
 
        ! The numbering of the cell of which left boundary is the domain boundary is 0.
        ! if comp.coord of the pos is s, the real coordinate of s is
        ! (the coordinate of the left boundary of the Floor(s)-th cell)
        ! + (s-(int(s))*(cell-width).
        ! In other words,  the coordinate of the center of the cell is x_cc(cell).
 
        !coordinates in computational space
        scoord(1) = cell(1) + (pos(1) - x_cb(cell(1) - 1))/dx(cell(1))
        scoord(2) = cell(2) + (pos(2) - y_cb(cell(2) - 1))/dy(cell(2))
        scoord(3) = 0._wp
        if (p > 0) scoord(3) = cell(3) + (pos(3) - z_cb(cell(3) - 1))/dz(cell(3))
        cell(:) = int(scoord(:))
        do i = 1, num_dims
            if (scoord(i) < 0._wp) cell(i) = cell(i) - 1
        end do
 
    end subroutine s_locate_cell
 
    !> This subroutine transfer data into the temporal variables.
    impure subroutine s_transfer_data_to_tmp()
 
        integer :: k
 
 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) 
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1239 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        do k = 1, nbubs
            gas_p(k, 2) = gas_p(k, 1)
            gas_mv(k, 2) = gas_mv(k, 1)
            intfc_rad(k, 2) = intfc_rad(k, 1)
            intfc_vel(k, 2) = intfc_vel(k, 1)
            mtn_pos(k, 1:3, 2) = mtn_pos(k, 1:3, 1)
            mtn_posprev(k, 1:3, 2) = mtn_posprev(k, 1:3, 1)
            mtn_vel(k, 1:3, 2) = mtn_vel(k, 1:3, 1)
            mtn_s(k, 1:3, 2) = mtn_s(k, 1:3, 1)
        end do
 
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1250 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
    end subroutine s_transfer_data_to_tmp
 
    !> The purpose of this procedure is to determine if the global coordinates of the bubbles
        !!      are present in the current MPI processor (including ghost cells).
        !! @param pos_part Spatial coordinates of the bubble
    function particle_in_domain(pos_part)
 
        logical :: particle_in_domain
        real(wp), dimension(3), intent(in) :: pos_part
 
        ! 2D
        if (p == 0 .and. cyl_coord .neqv. .true.) then
            ! Defining a virtual z-axis that has the same dimensions as y-axis
            ! defined in the input file
            particle_in_domain = ((pos_part(1) < x_cb(m + buff_size)) .and. (pos_part(1) >= x_cb(-buff_size - 1)) .and. &
                                  (pos_part(2) < y_cb(n + buff_size)) .and. (pos_part(2) >= y_cb(-buff_size - 1)) .and. &
                                  (pos_part(3) < lag_params%charwidth/2._wp) .and. (pos_part(3) >= -lag_params%charwidth/2._wp))
        else
            ! cyl_coord
            particle_in_domain = ((pos_part(1) < x_cb(m + buff_size)) .and. (pos_part(1) >= x_cb(-buff_size - 1)) .and. &
                                  (abs(pos_part(2)) < y_cb(n + buff_size)) .and. (abs(pos_part(2)) >= max(y_cb(-buff_size - 1), 0._wp)))
        end if
 
        ! 3D
        if (p > 0) then
            particle_in_domain = ((pos_part(1) < x_cb(m + buff_size)) .and. (pos_part(1) >= x_cb(-buff_size - 1)) .and. &
                                  (pos_part(2) < y_cb(n + buff_size)) .and. (pos_part(2) >= y_cb(-buff_size - 1)) .and. &
                                  (pos_part(3) < z_cb(p + buff_size)) .and. (pos_part(3) >= z_cb(-buff_size - 1)))
        end if
 
        ! For symmetric and wall boundary condition
        if (any(bc_x%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/))) then
            particle_in_domain = (particle_in_domain .and. (pos_part(1) >= x_cb(-1)))
        end if
        if (any(bc_x%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/))) then
            particle_in_domain = (particle_in_domain .and. (pos_part(1) < x_cb(m)))
        end if
        if (any(bc_y%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. (.not. cyl_coord)) then
            particle_in_domain = (particle_in_domain .and. (pos_part(2) >= y_cb(-1)))
        end if
        if (any(bc_y%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/)) .and. (.not. cyl_coord)) then
            particle_in_domain = (particle_in_domain .and. (pos_part(2) < y_cb(n)))
        end if
        if (p > 0) then
            if (any(bc_z%beg == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/))) then
                particle_in_domain = (particle_in_domain .and. (pos_part(3) >= z_cb(-1)))
            end if
            if (any(bc_z%end == (/bc_reflective, bc_char_slip_wall, bc_slip_wall, bc_no_slip_wall/))) then
                particle_in_domain = (particle_in_domain .and. (pos_part(3) < z_cb(p)))
            end if
        end if
 
    end function particle_in_domain
 
    !> The purpose of this procedure is to determine if the lagrangian bubble is located in the
        !!       physical domain. The ghost cells are not part of the physical domain.
        !! @param pos_part Spatial coordinates of the bubble
    function particle_in_domain_physical(pos_part)
 
        logical :: particle_in_domain_physical
        real(wp), dimension(3), intent(in) :: pos_part
 
        particle_in_domain_physical = ((pos_part(1) < x_cb(m)) .and. (pos_part(1) >= x_cb(-1)) .and. &
                                       (pos_part(2) < y_cb(n)) .and. (pos_part(2) >= y_cb(-1)))
 
        if (p > 0) then
            particle_in_domain_physical = (particle_in_domain_physical .and. (pos_part(3) < z_cb(p)) .and. (pos_part(3) >= z_cb(-1)))
        end if
 
    end function particle_in_domain_physical
 
    !> The purpose of this procedure is to calculate the gradient of a scalar field along the x, y and z directions
        !!      following a second-order central difference considering uneven widths
        !! @param q Input scalar field
        !! @param dq Output gradient of q
        !! @param dir Gradient spatial direction
    subroutine s_gradient_dir(q, dq, dir)
 
        real(stp), dimension(idwbuff(1)%beg:, idwbuff(2)%beg:, idwbuff(3)%beg:), intent(inout) :: q, dq
        integer, intent(in) :: dir
 
        integer :: i, j, k
 
        if (dir == 1) then
            ! Gradient in x dir.
 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1337 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            do k = 0, p
                do j = 0, n
                    do i = 0, m
                        dq(i, j, k) = q(i, j, k)*(dx(i + 1) - dx(i - 1)) &
                                      + q(i + 1, j, k)*(dx(i) + dx(i - 1)) &
                                      - q(i - 1, j, k)*(dx(i) + dx(i + 1))
                        dq(i, j, k) = dq(i, j, k)/ &
                                      ((dx(i) + dx(i - 1))*(dx(i) + dx(i + 1)))
                    end do
                end do
            end do
 
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1349 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        elseif (dir == 2) then
            ! Gradient in y dir.
 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1352 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            do k = 0, p
                do j = 0, n
                    do i = 0, m
                        dq(i, j, k) = q(i, j, k)*(dy(j + 1) - dy(j - 1)) &
                                      + q(i, j + 1, k)*(dy(j) + dy(j - 1)) &
                                      - q(i, j - 1, k)*(dy(j) + dy(j + 1))
                        dq(i, j, k) = dq(i, j, k)/ &
                                      ((dy(j) + dy(j - 1))*(dy(j) + dy(j + 1)))
                    end do
                end do
            end do
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        elseif (dir == 3) then
            ! Gradient in z dir.
 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(i, j, k) 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(i, j, k) 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
            do k = 0, p
                do j = 0, n
                    do i = 0, m
                        dq(i, j, k) = q(i, j, k)*(dz(k + 1) - dz(k - 1)) &
                                      + q(i, j, k + 1)*(dz(k) + dz(k - 1)) &
                                      - q(i, j, k - 1)*(dz(k) + dz(k + 1))
                        dq(i, j, k) = dq(i, j, k)/ &
                                      ((dz(k) + dz(k - 1))*(dz(k) + dz(k + 1)))
                    end do
                end do
            end do
 
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1379 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        end if
 
    end subroutine s_gradient_dir
 
    !> Subroutine that writes on each time step the changes of the lagrangian bubbles.
        !!  @param qtime Current time
    impure subroutine s_write_lag_particles(qtime)
 
        real(wp), intent(in) :: qtime
        integer :: k
 
        character(LEN=path_len + 2*name_len) :: file_loc
        logical :: file_exist
 
        character(LEN=25) :: fmt
 
        write (file_loc, '(A,I0,A)') 'lag_bubble_evol_', proc_rank, '.dat'
        file_loc = trim(case_dir)//'/D/'//trim(file_loc)
        inquire (file=trim(file_loc), exist=file_exist)
 
        if (precision == 1) then
            fmt = "(A16,A14,8A16)"
        else
            fmt = "(A24,A14,8A24)"
        end if
 
        if (.not. file_exist) then
            open (11, file=trim(file_loc), form='formatted', position='rewind')
            write (11, fmt) 'currentTime', 'particleID', 'x', 'y', 'z', &
                'coreVaporMass', 'coreVaporConcentration', 'radius', 'interfaceVelocity', &
                'corePressure'
        else
            open (11, file=trim(file_loc), form='formatted', position='append')
        end if
 
        if (precision == 1) then
            fmt = "(F16.8,I14,8F16.8)"
        else
            fmt = "(F24.16,I14,8F24.16)"
        end if
 
        ! Cycle through list
        do k = 1, nbubs
            write (11, fmt) &
                qtime, &
                lag_id(k, 1), &
                mtn_pos(k, 1, 1), &
                mtn_pos(k, 2, 1), &
                mtn_pos(k, 3, 1), &
                gas_mv(k, 1), &
                gas_mv(k, 1)/(gas_mv(k, 1) + gas_mg(k)), &
                intfc_rad(k, 1), &
                intfc_vel(k, 1), &
                gas_p(k, 1)
        end do
 
        close (11)
 
    end subroutine s_write_lag_particles
 
    !>  Subroutine that writes some useful statistics related to the volume fraction
            !!       of the particles (void fraction) in the computatioational domain
            !!       on each time step.
            !!  @param qtime Current time
    impure subroutine s_write_void_evol(qtime)
 
        real(wp), intent(in) :: qtime
        real(wp) :: volcell, voltot
        real(wp) :: lag_void_max, lag_void_avg, lag_vol
        real(wp) :: void_max_glb, void_avg_glb, vol_glb
 
        integer :: i, j, k
 
        character(LEN=path_len + 2*name_len) :: file_loc
        logical :: file_exist
 
        if (proc_rank == 0) then
            write (file_loc, '(A)') 'voidfraction.dat'
            file_loc = trim(case_dir)//'/D/'//trim(file_loc)
            inquire (file=trim(file_loc), exist=file_exist)
            if (.not. file_exist) then
                open (12, file=trim(file_loc), form='formatted', position='rewind')
                !write (12, *) 'currentTime, averageVoidFraction, ', &
                !    'maximumVoidFraction, totalParticlesVolume'
                !write (12, *) 'The averageVoidFraction value does ', &
                !    'not reflect the real void fraction in the cloud since the ', &
                !    'cells which do not have bubbles are not accounted'
            else
                open (12, file=trim(file_loc), form='formatted', position='append')
            end if
        end if
 
        lag_void_max = 0._wp
        lag_void_avg = 0._wp
        lag_vol = 0._wp
 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop collapse(3) gang vector default(present) private(volcell) reduction(+:lag_vol, lag_void_avg) reduction(MAX:lag_void_max) copy(lag_vol, lag_void_avg, lag_void_max) 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) collapse(3) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(volcell) reduction(+:lag_vol, lag_void_avg) reduction(MAX:lag_void_max) map(tofrom:lag_vol, lag_void_avg, lag_void_max) 
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1475 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
        do k = 0, p
            do j = 0, n
                do i = 0, m
                    lag_void_max = max(lag_void_max, 1._wp - q_beta(1)%sf(i, j, k))
                    call s_get_char_vol(i, j, k, volcell)
                    if ((1._wp - q_beta(1)%sf(i, j, k)) > 5.0d-11) then
                        lag_void_avg = lag_void_avg + (1._wp - q_beta(1)%sf(i, j, k))*volcell
                        lag_vol = lag_vol + volcell
                    end if
                end do
            end do
        end do
 
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1488 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
#ifdef MFC_MPI
        if (num_procs > 1) then
            call s_mpi_allreduce_max(lag_void_max, void_max_glb)
            lag_void_max = void_max_glb
            call s_mpi_allreduce_sum(lag_vol, vol_glb)
            lag_vol = vol_glb
            call s_mpi_allreduce_sum(lag_void_avg, void_avg_glb)
            lag_void_avg = void_avg_glb
        end if
#endif
        voltot = lag_void_avg
        ! This voidavg value does not reflect the real void fraction in the cloud
        ! since the cell which does not have bubbles are not accounted
        if (lag_vol > 0._wp) lag_void_avg = lag_void_avg/lag_vol
 
        if (proc_rank == 0) then
            write (12, '(6X,4e24.8)') &
                qtime, &
                lag_void_avg, &
                lag_void_max, &
                voltot
            close (12)
        end if
 
    end subroutine s_write_void_evol
 
    !>  Subroutine that writes the restarting files for the particles in the lagrangian solver.
        !!  @param t_step Current time step
    impure subroutine s_write_restart_lag_bubbles(t_step)
 
        ! Generic string used to store the address of a particular file
        integer, intent(in) :: t_step
 
        character(LEN=path_len + 2*name_len) :: file_loc
        logical :: file_exist
        integer :: bub_id, tot_part
        integer :: i, k
 
#ifdef MFC_MPI
        ! For Parallel I/O
        integer :: ifile, ierr
        integer, dimension(MPI_STATUS_SIZE) :: status
        integer(KIND=MPI_OFFSET_KIND) :: disp
        integer :: view
        integer, dimension(2) :: gsizes, lsizes, start_idx_part
        integer, allocatable :: proc_bubble_counts(:)
        real(wp), dimension(1:1, 1:lag_io_vars) :: dummy
        dummy = 0._wp
 
        bub_id = 0._wp
        if (nbubs /= 0) then
            do k = 1, nbubs
                if (particle_in_domain_physical(mtn_pos(k, 1:3, 1))) then
                    bub_id = bub_id + 1
                end if
            end do
        end if
 
        if (.not. parallel_io) return
 
        allocate (proc_bubble_counts(num_procs))
 
        lsizes(1) = bub_id
        lsizes(2) = lag_io_vars
 
        ! Total number of particles
        call mpi_allreduce(bub_id, tot_part, 1, mpi_integer, &
                           mpi_sum, mpi_comm_world, ierr)
 
        call mpi_allgather(bub_id, 1, mpi_integer, proc_bubble_counts, 1, mpi_integer, &
                           mpi_comm_world, ierr)
 
        ! Calculate starting index for this processor's particles
        call mpi_exscan(lsizes(1), start_idx_part(1), 1, mpi_integer, mpi_sum, mpi_comm_world, ierr)
        if (proc_rank == 0) start_idx_part(1) = 0
        start_idx_part(2) = 0
 
        gsizes(1) = tot_part
        gsizes(2) = lag_io_vars
 
        write (file_loc, '(A,I0,A)') 'lag_bubbles_', t_step, '.dat'
        file_loc = trim(case_dir)//'/restart_data'//trim(mpiiofs)//trim(file_loc)
 
        ! Clean up existing file
        if (proc_rank == 0) then
            inquire (file=trim(file_loc), exist=file_exist)
            if (file_exist) then
                call mpi_file_delete(file_loc, mpi_info_int, ierr)
            end if
        end if
 
        call mpi_barrier(mpi_comm_world, ierr)
 
        if (proc_rank == 0) then
            call mpi_file_open(mpi_comm_self, file_loc, &
                               ior(mpi_mode_wronly, mpi_mode_create), &
                               mpi_info_int, ifile, ierr)
 
            ! Write header using MPI I/O for consistency
            call mpi_file_write(ifile, tot_part, 1, mpi_integer, status, ierr)
            call mpi_file_write(ifile, mytime, 1, mpi_p, status, ierr)
            call mpi_file_write(ifile, dt, 1, mpi_p, status, ierr)
            call mpi_file_write(ifile, num_procs, 1, mpi_integer, status, ierr)
            call mpi_file_write(ifile, proc_bubble_counts, num_procs, mpi_integer, status, ierr)
 
            call mpi_file_close(ifile, ierr)
        end if
 
        call mpi_barrier(mpi_comm_world, ierr)
 
        if (bub_id > 0) then
            allocate (mpi_io_data_lag_bubbles(max(1, bub_id), 1:lag_io_vars))
 
            i = 1
            do k = 1, nbubs
                if (particle_in_domain_physical(mtn_pos(k, 1:3, 1))) then
                    mpi_io_data_lag_bubbles(i, 1) = real(lag_id(k, 1))
                    mpi_io_data_lag_bubbles(i, 2:4) = mtn_pos(k, 1:3, 1)
                    mpi_io_data_lag_bubbles(i, 5:7) = mtn_posprev(k, 1:3, 1)
                    mpi_io_data_lag_bubbles(i, 8:10) = mtn_vel(k, 1:3, 1)
                    mpi_io_data_lag_bubbles(i, 11) = intfc_rad(k, 1)
                    mpi_io_data_lag_bubbles(i, 12) = intfc_vel(k, 1)
                    mpi_io_data_lag_bubbles(i, 13) = bub_r0(k)
                    mpi_io_data_lag_bubbles(i, 14) = rmax_stats(k)
                    mpi_io_data_lag_bubbles(i, 15) = rmin_stats(k)
                    mpi_io_data_lag_bubbles(i, 16) = bub_dphidt(k)
                    mpi_io_data_lag_bubbles(i, 17) = gas_p(k, 1)
                    mpi_io_data_lag_bubbles(i, 18) = gas_mv(k, 1)
                    mpi_io_data_lag_bubbles(i, 19) = gas_mg(k)
                    mpi_io_data_lag_bubbles(i, 20) = gas_betat(k)
                    mpi_io_data_lag_bubbles(i, 21) = gas_betac(k)
                    i = i + 1
                end if
            end do
 
            call mpi_type_create_subarray(2, gsizes, lsizes, start_idx_part, &
                                          mpi_order_fortran, mpi_p, view, ierr)
            call mpi_type_commit(view, ierr)
 
            call mpi_file_open(mpi_comm_world, file_loc, &
                               ior(mpi_mode_wronly, mpi_mode_create), &
                               mpi_info_int, ifile, ierr)
 
            ! Skip header (written by rank 0)
            disp = int(sizeof(tot_part) + 2*sizeof(mytime) + sizeof(num_procs) + &
                       num_procs*sizeof(proc_bubble_counts(1)), mpi_offset_kind)
            call mpi_file_set_view(ifile, disp, mpi_p, view, 'native', mpi_info_int, ierr)
 
            call mpi_file_write_all(ifile, mpi_io_data_lag_bubbles, &
                                    lag_io_vars*bub_id, mpi_p, status, ierr)
 
            call mpi_file_close(ifile, ierr)
 
            deallocate (mpi_io_data_lag_bubbles)
 
        else
            call mpi_type_contiguous(0, mpi_p, view, ierr)
            call mpi_type_commit(view, ierr)
 
            call mpi_file_open(mpi_comm_world, file_loc, &
                               ior(mpi_mode_wronly, mpi_mode_create), &
                               mpi_info_int, ifile, ierr)
 
            ! Skip header (written by rank 0)
            disp = int(sizeof(tot_part) + 2*sizeof(mytime) + sizeof(num_procs) + &
                       num_procs*sizeof(proc_bubble_counts(1)), mpi_offset_kind)
            call mpi_file_set_view(ifile, disp, mpi_p, view, 'native', mpi_info_int, ierr)
 
            call mpi_file_write_all(ifile, dummy, 0, mpi_p, status, ierr)
 
            call mpi_file_close(ifile, ierr)
        end if
 
        deallocate (proc_bubble_counts)
 
#endif
 
    end subroutine s_write_restart_lag_bubbles
 
    !>  This procedure calculates the maximum and minimum radius of each bubble.
    subroutine s_calculate_lag_bubble_stats()
 
        integer :: k
 
 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc parallel loop gang vector default(present) private(k) reduction(MAX:Rmax_glb) reduction(MIN:Rmin_glb) copy(Rmax_glb, Rmin_glb) 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target teams loop defaultmap(firstprivate:scalar) bind(teams,parallel) defaultmap(tofrom:aggregate) defaultmap(tofrom:allocatable) defaultmap(tofrom:pointer) private(k) reduction(MAX:Rmax_glb) reduction(MIN:Rmin_glb) map(tofrom:Rmax_glb, Rmin_glb) 
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1674 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1676 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        do k = 1, nbubs
            rmax_glb = max(rmax_glb, intfc_rad(k, 1)/bub_r0(k))
            rmin_glb = min(rmin_glb, intfc_rad(k, 1)/bub_r0(k))
            rmax_stats(k) = max(rmax_stats(k), intfc_rad(k, 1)/bub_r0(k))
            rmin_stats(k) = min(rmin_stats(k), intfc_rad(k, 1)/bub_r0(k))
        end do
 
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc end parallel loop
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp end target teams loop
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1682 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
 
    end subroutine s_calculate_lag_bubble_stats
 
    !>  Subroutine that writes the maximum and minimum radius of each bubble.
    impure subroutine s_write_lag_bubble_stats()
 
        integer :: k
        character(LEN=path_len + 2*name_len) :: file_loc
 
        character(len=20) :: fmt
 
        write (file_loc, '(A,I0,A)') 'stats_lag_bubbles_', proc_rank, '.dat'
        file_loc = trim(case_dir)//'/D/'//trim(file_loc)
 
 
# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update host(Rmax_glb, Rmin_glb) 
# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update from(Rmax_glb, Rmin_glb) 
# 1697 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
        if (precision == 1) then
            fmt = "(A10,A14,5A16)"
        else
            fmt = "(A10,A14,5A24)"
        end if
 
        open (13, file=trim(file_loc), form='formatted', position='rewind')
        write (13, fmt) 'proc_rank', 'particleID', 'x', 'y', 'z', 'Rmax_glb', 'Rmin_glb'
 
        if (precision == 1) then
            fmt = "(I10,I14,5F16.8)"
        else
            fmt = "(I10,I14,5F24.16)"
        end if
 
        do k = 1, nbubs
            write (13, fmt) &
                proc_rank, &
                lag_id(k, 1), &
                mtn_pos(k, 1, 1), &
                mtn_pos(k, 2, 1), &
                mtn_pos(k, 3, 1), &
                rmax_stats(k), &
                rmin_stats(k)
        end do
 
        close (13)
 
    end subroutine s_write_lag_bubble_stats
 
    !> The purpose of this subroutine is to remove one specific particle if dt is too small.
          !! @param bub_id Particle id
    impure subroutine s_remove_lag_bubble(bub_id)
 
        integer, intent(in) :: bub_id
 
        integer :: i
 
 
# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc loop seq 
# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1737 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
        do i = bub_id, nbubs - 1
            lag_id(i, 1) = lag_id(i + 1, 1)
            bub_r0(i) = bub_r0(i + 1)
            rmax_stats(i) = rmax_stats(i + 1)
            rmin_stats(i) = rmin_stats(i + 1)
            gas_mg(i) = gas_mg(i + 1)
            gas_betat(i) = gas_betat(i + 1)
            gas_betac(i) = gas_betac(i + 1)
            bub_dphidt(i) = bub_dphidt(i + 1)
            gas_p(i, 1:2) = gas_p(i + 1, 1:2)
            gas_mv(i, 1:2) = gas_mv(i + 1, 1:2)
            intfc_rad(i, 1:2) = intfc_rad(i + 1, 1:2)
            intfc_vel(i, 1:2) = intfc_vel(i + 1, 1:2)
            mtn_pos(i, 1:3, 1:2) = mtn_pos(i + 1, 1:3, 1:2)
            mtn_posprev(i, 1:3, 1:2) = mtn_posprev(i + 1, 1:3, 1:2)
            mtn_vel(i, 1:3, 1:2) = mtn_vel(i + 1, 1:3, 1:2)
            mtn_s(i, 1:3, 1:2) = mtn_s(i + 1, 1:3, 1:2)
            intfc_draddt(i, 1:lag_num_ts) = intfc_draddt(i + 1, 1:lag_num_ts)
            intfc_dveldt(i, 1:lag_num_ts) = intfc_dveldt(i + 1, 1:lag_num_ts)
            gas_dpdt(i, 1:lag_num_ts) = gas_dpdt(i + 1, 1:lag_num_ts)
            gas_dmvdt(i, 1:lag_num_ts) = gas_dmvdt(i + 1, 1:lag_num_ts)
        end do
 
        nbubs = nbubs - 1
 
# 1762 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1762 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc update device(nBubs) 
# 1762 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1762 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target update to(nBubs) 
# 1762 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
 
    end subroutine s_remove_lag_bubble
 
    !> The purpose of this subroutine is to deallocate variables
    impure subroutine s_finalize_lagrangian_solver()
 
        integer :: i
 
        do i = 1, q_beta_idx
#ifdef MFC_DEBUG
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1772: ', '@:DEALLOCATE(q_beta(i)%sf)'
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(q_beta(i)%sf) 
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:q_beta(i)%sf) 
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1772 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (q_beta(i)%sf)
        end do
#ifdef MFC_DEBUG
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1774: ', '@:DEALLOCATE(q_beta)'
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(q_beta) 
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:q_beta) 
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1774 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (q_beta)
 
        !Deallocating space
#ifdef MFC_DEBUG
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1777: ', '@:DEALLOCATE(lag_id)'
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(lag_id) 
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:lag_id) 
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1777 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (lag_id)
#ifdef MFC_DEBUG
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1778: ', '@:DEALLOCATE(bub_R0)'
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(bub_R0) 
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:bub_R0) 
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1778 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (bub_r0)
#ifdef MFC_DEBUG
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1779: ', '@:DEALLOCATE(Rmax_stats)'
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(Rmax_stats) 
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:Rmax_stats) 
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1779 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (rmax_stats)
#ifdef MFC_DEBUG
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1780: ', '@:DEALLOCATE(Rmin_stats)'
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(Rmin_stats) 
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:Rmin_stats) 
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1780 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (rmin_stats)
#ifdef MFC_DEBUG
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1781: ', '@:DEALLOCATE(gas_mg)'
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_mg) 
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_mg) 
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1781 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_mg)
#ifdef MFC_DEBUG
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1782: ', '@:DEALLOCATE(gas_betaT)'
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_betaT) 
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_betaT) 
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1782 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_betat)
#ifdef MFC_DEBUG
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1783: ', '@:DEALLOCATE(gas_betaC)'
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_betaC) 
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_betaC) 
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1783 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_betac)
#ifdef MFC_DEBUG
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1784: ', '@:DEALLOCATE(bub_dphidt)'
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(bub_dphidt) 
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:bub_dphidt) 
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1784 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (bub_dphidt)
#ifdef MFC_DEBUG
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1785: ', '@:DEALLOCATE(gas_p)'
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_p) 
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_p) 
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1785 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_p)
#ifdef MFC_DEBUG
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1786: ', '@:DEALLOCATE(gas_mv)'
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_mv) 
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_mv) 
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1786 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_mv)
#ifdef MFC_DEBUG
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1787: ', '@:DEALLOCATE(intfc_rad)'
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(intfc_rad) 
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:intfc_rad) 
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1787 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (intfc_rad)
#ifdef MFC_DEBUG
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1788: ', '@:DEALLOCATE(intfc_vel)'
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(intfc_vel) 
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:intfc_vel) 
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1788 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (intfc_vel)
#ifdef MFC_DEBUG
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1789: ', '@:DEALLOCATE(mtn_pos)'
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_pos) 
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_pos) 
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1789 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_pos)
#ifdef MFC_DEBUG
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1790: ', '@:DEALLOCATE(mtn_posPrev)'
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_posPrev) 
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_posPrev) 
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1790 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_posprev)
#ifdef MFC_DEBUG
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1791: ', '@:DEALLOCATE(mtn_vel)'
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_vel) 
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_vel) 
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1791 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_vel)
#ifdef MFC_DEBUG
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1792: ', '@:DEALLOCATE(mtn_s)'
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_s) 
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_s) 
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1792 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_s)
#ifdef MFC_DEBUG
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1793: ', '@:DEALLOCATE(intfc_draddt)'
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(intfc_draddt) 
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:intfc_draddt) 
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1793 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (intfc_draddt)
#ifdef MFC_DEBUG
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1794: ', '@:DEALLOCATE(intfc_dveldt)'
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(intfc_dveldt) 
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:intfc_dveldt) 
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1794 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (intfc_dveldt)
#ifdef MFC_DEBUG
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1795: ', '@:DEALLOCATE(gas_dpdt)'
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_dpdt) 
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_dpdt) 
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1795 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_dpdt)
#ifdef MFC_DEBUG
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1796: ', '@:DEALLOCATE(gas_dmvdt)'
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(gas_dmvdt) 
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:gas_dmvdt) 
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1796 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (gas_dmvdt)
#ifdef MFC_DEBUG
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1797: ', '@:DEALLOCATE(mtn_dposdt)'
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_dposdt) 
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_dposdt) 
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1797 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_dposdt)
#ifdef MFC_DEBUG
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    block
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        use iso_fortran_env, only: output_unit
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        print *, 'm_bubbles_EL.fpp:1798: ', '@:DEALLOCATE(mtn_dveldt)'
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
        call flush (output_unit)
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    end block
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
 
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#if defined(MFC_OpenACC)
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$acc exit data delete(mtn_dveldt) 
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#elif defined(MFC_OpenMP)
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
!$omp target exit data map(release:mtn_dveldt) 
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
#endif
# 1798 "/home/runner/work/MFC/MFC/src/simulation/m_bubbles_EL.fpp"
    deallocate (mtn_dveldt)
 
    end subroutine s_finalize_lagrangian_solver
 
end module m_bubbles_el