m_global_parameters.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!>
!! @file
!! @brief Contains module m_global_parameters
 
# 1 "/home/runner/work/MFC/MFC/src/common/include/case.fpp" 1
! This file exists so that Fypp can be run without generating case.fpp files for
! each target. This is useful when generating documentation, for example. This
! should also let MFC be built with CMake directly, without invoking mfc.sh.
 
! For pre-process.
# 8 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
 
! For moving immersed boundaries in simulation
# 12 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
# 6 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 167 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 167 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 126 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 197 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 211 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 236 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 247 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 249 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 260 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 310 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 320 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 356 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 366 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 373 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 379 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 385 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 391 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 397 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 403 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
! New line at end of file is required for FYPP
# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 167 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
 
# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
! GPU parallel region (scalar reductions, maxval/minval)
# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! GPU parallel loop over threads (most common GPU macro)
# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Required closing for GPU_PARALLEL_LOOP
# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Mark routine for device compilation
# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Declare device-resident data
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Inner loop within a GPU parallel region
# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Scoped GPU data region
# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Host code with device pointers (for MPI with GPU buffers)
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Allocate device memory (unscoped)
# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Free device memory
# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Atomic operation on device
# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! End atomic capture block
# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Copy data between host and device
# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Synchronization barrier
# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Import GPU library module (openacc or omp_lib)
# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for AMD compiler
# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-Cray compilers
# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for Cray compiler
# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-NVIDIA compilers
# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
! example see misc/nvidia_uvm/bind.sh.
# 55 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 75 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 83 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 107 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 123 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 148 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 161 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 7 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp" 2
 
!> @brief Global parameters for the computational domain, fluid properties, and simulation algorithm configuration
module m_global_parameters
 
#ifdef MFC_MPI
    use mpi  !< message passing interface (mpi) module
#endif
 
    use m_derived_types
    use m_helper_basic
    ! Shared state: generated_decls, generated_case_opt_decls, sys_size, eqn_idx, chemistry, shear_*
    use m_global_parameters_common
    ! $:USE_GPU_MODULE()
 
    implicit none
 
    real(wp) :: wall_time = 0
    real(wp) :: wall_time_avg = 0
 
    ! Logistics
    integer :: num_procs  !< Number of processors
    ! Computational Domain Parameters
    integer :: proc_rank  !< Rank of the local processor
 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(num_procs, proc_rank) 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (num_procs, proc_rank) 
# 30 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Max and min number of cells in a direction of each combination of x-,y-, and z-
    type(cell_num_bounds) :: cells_bounds
 
    !> @name Global number of cells in each direction
    !> @{
    integer :: m_glb, n_glb, p_glb
    !> @}
 
    !> @name Cylindrical coordinates (either axisymmetric or full 3D)
    !> @{
    integer :: grid_geometry
    !> @}
 
# 44 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 44 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(grid_geometry) 
# 44 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 44 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (grid_geometry) 
# 44 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Cell-boundary (CB) locations in the x-, y- and z-directions, respectively
    !> @{
    real(wp), target, allocatable, dimension(:) :: x_cb, y_cb, z_cb
    type(bounds_info), dimension(3)             :: glb_bounds
    !> @}
 
    !> @name Cell-center (CC) locations in the x-, y- and z-directions, respectively
    !> @{
    real(wp), target, allocatable, dimension(:) :: x_cc, y_cc, z_cc
    !> @}
    !> @name Cell-width distributions in the x-, y- and z-directions, respectively
    !> @{
    real(wp), target, allocatable, dimension(:) :: dx, dy, dz
    !> @}
 
 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(x_cb, y_cb, z_cb, x_cc, y_cc, z_cc, dx, dy, dz) 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (x_cb, y_cb, z_cb, x_cc, y_cc, z_cc, dx, dy, dz) 
# 61 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! dt, m, n, p, cfl_target: GPU-declared via generated_decls.fpp (registered params)
 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(glb_bounds) 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (glb_bounds) 
# 64 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: cfl_dt
    ! Simulation Algorithm Parameters generated_case_opt_decls.fpp: now in m_global_parameters_common
 
    !> Hypoelastic NC velocity-coupling mode; exactly one value, derived from riemann_solver + hypo_hll_interface_rhs.
    integer, parameter :: hypo_nc_mode_none = 0         !< no hypoelastic NC velocity coupling
    integer, parameter :: hypo_nc_mode_finite_diff = 1  !< velocity gradients by finite difference (HLL without interface RHS)
    !> interface-velocity export for the velocity-gradient tensor (HLL Method 2, HLLC)
    integer, parameter :: hypo_nc_mode_interface = 2
    integer, parameter :: hypo_nc_mode_dual_pass = 3  !< anchored dual-pass HLLD; all NC terms stay in the Riemann flux
    integer            :: hypo_nc_mode
    !> NC volume-fraction advection export mode; exactly one value, derived from riemann_solver + hll_u_interface.
    integer, parameter :: adv_src_mode_unset = 0        !< not yet derived
    integer, parameter :: adv_src_mode_alpha_iface = 1  !< flux_src exports per-fluid interface alpha
    integer, parameter :: adv_src_mode_vel_iface = 2    !< flux_src exports shared face-normal interface velocity
    integer, parameter :: adv_src_mode_none = 3         !< flux_src exports no NC advection quantity
    integer            :: adv_src_mode
    logical            :: use_nc_iface_vel              !< nc_iface_vel exports interface velocities needed outside flux_src
    ! chemistry: in m_global_parameters_common
    logical                :: shear_stress  !< Shear stresses
    logical                :: bulk_stress   !< Bulk stresses
    logical                :: bodyforces
    real(wp), dimension(3) :: accel_bf
 
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(accel_bf) 
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (accel_bf) 
# 88 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    ! $:GPU_DECLARE(create='[k_x,w_x,p_x,g_x,k_y,w_y,p_y,g_y,k_z,w_z,p_z,g_z]')
 
    !> Source fields for the spatially supported body force. `spatial_bf` and
    !> `bf_spatial_support` are auto-generated in generated_decls.fpp.
    real(wp), allocatable, dimension(:,:,:) :: spbf_source_x
    real(wp), allocatable, dimension(:,:,:) :: spbf_source_y
 
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(spbf_source_x, spbf_source_y) 
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (spbf_source_x, spbf_source_y) 
# 95 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! Synthetic turbulence (scalars auto-generated in generated_decls.fpp; their
    ! GPU_DECLARE lines live in m_global_parameters_common)
    integer, dimension(num_synth_shells_max)     :: synth_n_waves_per_shell
    real(wp), dimension(num_synth_shells_max)    :: synth_k_shell, synth_amp_shell
    real(wp), dimension(num_turb_sources_max, 3) :: turb_pos, synth_l
 
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(synth_n_waves_per_shell, synth_k_shell, synth_amp_shell) 
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (synth_n_waves_per_shell, synth_k_shell, synth_amp_shell) 
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(turb_pos, synth_L) 
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (turb_pos, synth_L) 
# 103 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: cpu_start, cpu_end, cpu_rate
 
 
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(shear_stress, bulk_stress) 
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (shear_stress, bulk_stress) 
# 107 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(hypo_nc_mode) 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (hypo_nc_mode) 
# 108 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical               :: bc_io
    logical, dimension(3) :: periodic_bc
    !> @name Boundary conditions (BC) in the x-, y- and z-directions, respectively
    !> @{
    type(int_bounds_info) :: bc_x, bc_y, bc_z
    type(bc_xyz_info)     :: bc
    !> @}
    !> @name Original boundary conditions preserved for immersed boundary code
    !> (bc_x/y/z get overwritten with MPI neighbor ranks during decomposition)
    !> @{
    type(int_bounds_info) :: ib_bc_x, ib_bc_y, ib_bc_z
    !> @}
#if defined(MFC_OpenACC)
 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc_x%vb1, bc_x%vb2, bc_x%vb3, bc_x%ve1, bc_x%ve2, bc_x%ve3) 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc_x%vb1, bc_x%vb2, bc_x%vb3, bc_x%ve1, bc_x%ve2, bc_x%ve3) 
# 123 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc_y%vb1, bc_y%vb2, bc_y%vb3, bc_y%ve1, bc_y%ve2, bc_y%ve3) 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc_y%vb1, bc_y%vb2, bc_y%vb3, bc_y%ve1, bc_y%ve2, bc_y%ve3) 
# 124 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc_z%vb1, bc_z%vb2, bc_z%vb3, bc_z%ve1, bc_z%ve2, bc_z%ve3) 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc_z%vb1, bc_z%vb2, bc_z%vb3, bc_z%ve1, bc_z%ve2, bc_z%ve3) 
# 125 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ib_bc_x%beg, ib_bc_y%beg, ib_bc_z%beg) 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ib_bc_x%beg, ib_bc_y%beg, ib_bc_z%beg) 
# 126 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#elif defined(MFC_OpenMP)
 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc_x, bc_y, bc_z) 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc_x, bc_y, bc_z) 
# 128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ib_bc_x, ib_bc_y, ib_bc_z) 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ib_bc_x, ib_bc_y, ib_bc_z) 
# 129 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#endif
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc) 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc) 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    type(bounds_info) :: neighbor_domain_x, neighbor_domain_y, neighbor_domain_z
    integer           :: num_gbl_ibs, num_local_ibs
 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(neighbor_domain_x, neighbor_domain_y, neighbor_domain_z, num_gbl_ibs) 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (neighbor_domain_x, neighbor_domain_y, neighbor_domain_z, num_gbl_ibs) 
# 134 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! proc_coords, start_idx, mpiiofs, mpi_info_int: in m_global_parameters_common
    ! down_sample: GPU-declared via generated_decls.fpp (registered param)
 
    !> @name MPI domain-decomposition state for Lagrangian-bubble exchange (#1290)
    !> @{
    type(bounds_info), allocatable, dimension(:) :: pcomm_coords    !< Local rank physical domain bounds
    type(int_bounds_info), dimension(3)          :: nidx            !< Neighbor index offsets per direction
    integer, allocatable, dimension(:,:,:)       :: neighbor_ranks  !< MPI ranks of neighbors
 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(pcomm_coords) 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (pcomm_coords) 
# 144 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
    type(mpi_io_var), public                      :: mpi_io_data
    type(mpi_io_ib_var), public                   :: mpi_io_ib_data
    type(mpi_io_airfoil_ib_var), public           :: mpi_io_airfoil_ib_data
    type(mpi_io_levelset_var), public             :: mpi_io_levelset_data
    type(mpi_io_levelset_norm_var), public        :: mpi_io_levelsetnorm_data
    real(wp), allocatable, dimension(:,:), public :: mpi_io_data_lag_bubbles
 
    ! sys_size and eqn_idx: in m_global_parameters_common (GPU_DECLARE there too)
    type(qbmm_idx_info) :: qbmm_idx  !< QBMM moment index mappings (allocatable; GPU-managed separately).
 
    ! Cell Indices for the (local) interior points (O-m, O-n, 0-p). Stands for "InDices With INTerior".
    type(int_bounds_info) :: idwint(1:3)
 
# 158 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 158 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(idwint) 
# 158 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 158 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (idwint) 
# 158 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! Cell Indices for the entire (local) domain. In simulation and post_process, this includes the buffer region. idwbuff and
    ! idwint are the same otherwise. Stands for "InDices With BUFFer".
    type(int_bounds_info) :: idwbuff(1:3)
 
# 163 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 163 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(idwbuff) 
# 163 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 163 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (idwbuff) 
# 163 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Herschel-Bulkley non-Newtonian viscosity: per-fluid flags and parameter arrays.
    !> @{
    logical                             :: any_non_newtonian  !< .true. if any fluid is non-Newtonian
    logical, allocatable, dimension(:)  :: is_non_newtonian   !< per-fluid NN flag
    real(wp), allocatable, dimension(:) :: hb_tau0, hb_k, hb_nn, hb_m_arr
    real(wp), allocatable, dimension(:) :: hb_mu_min, hb_mu_max
    real(wp), allocatable, dimension(:) :: fluid_inv_re       !< per-fluid Newtonian inverse-Re
    !> @}
 
 
# 174 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 174 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(any_non_newtonian, is_non_newtonian, hb_tau0, hb_K, hb_nn, hb_m_arr, hb_mu_min, hb_mu_max, fluid_inv_re) 
# 174 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 174 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (any_non_newtonian, is_non_newtonian, hb_tau0, hb_K, hb_nn, hb_m_arr, hb_mu_min, hb_mu_max, fluid_inv_re) 
# 174 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! WENO averaging flag: use arithmetic mean or unaltered WENO-reconstructed cell-boundary values
    !> @{
    real(wp) :: wa_flg
    !> @}
 
 
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(wa_flg) 
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (wa_flg) 
# 181 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name The coordinate direction indexes and flags (flg), respectively, for which the configurations will be determined with
    !! respect to a working direction and that will be used to isolate the contributions, in that direction, in the dimensionally
    !! split system of equations. Declared here rather than in m_global_parameters_common so the hot dimensionally-split kernels
    !! (Riemann solvers) read them from their own module: use-associating them from common costs ~18 kB/work-item of register spill
    !! on AMD OpenMP offload. Common code takes the mapping as explicit arguments instead.
    !> @{
    integer, dimension(3)  :: dir_idx
    real(wp), dimension(3) :: dir_flg
    integer, dimension(3)  :: dir_idx_tau  !< (nn, nt, nt2) stress indices for wave speeds and momentum flux
    integer, dimension(6)  :: stress_perm  !< Full tensor permutation: local basis -> physical storage index
    !> @}
 
 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(dir_idx, dir_flg, dir_idx_tau, stress_perm) 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (dir_idx, dir_flg, dir_idx_tau, stress_perm) 
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: buff_size  !< Number of ghost cells for boundary condition storage
 
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(buff_size) 
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (buff_size) 
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices: in m_global_parameters_common
 
    ! END: Simulation Algorithm Parameters
 
    ! Fluids Physical Parameters fluid_pp, bub_pp: auto-generated in generated_decls.fpp
 
    integer :: fd_number  !< Finite-difference half-stencil size: MAX(1, fd_order/2)
 
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(fd_number) 
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (fd_number) 
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Centered finite-difference coefficients in x-, y- and z-coordinate directions
    !> @{
    real(wp), allocatable, dimension(:,:) :: fd_coeff_x
    real(wp), allocatable, dimension(:,:) :: fd_coeff_y
    real(wp), allocatable, dimension(:,:) :: fd_coeff_z
    !> @}
 
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(fd_coeff_x, fd_coeff_y, fd_coeff_z) 
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (fd_coeff_x, fd_coeff_y, fd_coeff_z) 
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! probe: auto-generated in generated_decls.fpp
 
    !> @name Reference density and pressure for Tait EOS
    !> @{
    !> @name Immersed Boundaries
    !> patch_ib, ib_airfoil, stl_models, particle_cloud: auto-generated in generated_decls.fpp
    !> @{
    integer, dimension(num_local_ibs_max) :: local_ib_patch_ids  !< lookup table of IBs in the local compute domain
    integer, allocatable, dimension(:,:,:) :: ib_neighbor_ranks  !< MPI ranks of neighborhood domains, indexed (-N:N,-N:N,-N:N)
    type(ib_airfoil_grid), dimension(num_ib_airfoils_max) :: ib_airfoil_grids  !< Per-airfoil computed surface grids
 
 
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ib_airfoil_grids) 
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ib_airfoil_grids) 
# 228 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    !> @name Bubble modeling
    !> @{
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: nb
# 238 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
    real(wp) :: eu  !< Euler number
 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(Eu) 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (Eu) 
# 240 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp), dimension(:), allocatable :: weight  !< Simpson quadrature weights
    real(wp), dimension(:), allocatable :: r0      !< Bubble sizes
 
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(weight, R0) 
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (weight, R0) 
# 244 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp), allocatable, dimension(:,:,:) :: ptil  !< Pressure modification
 
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ptil) 
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ptil) 
# 247 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer, parameter :: nmom = 6  !< Number of carried moments per R0 location
    integer            :: nmomsp    !< Number of moments required by ensemble-averaging
    integer            :: nmomtot   !< Total number of carried moments moments/transport equations
 
# 252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(nmomsp, nmomtot) 
# 252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (nmomsp, nmomtot) 
# 252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(nb) 
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (nb) 
# 255 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 257 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
    type(scalar_field), allocatable, dimension(:)     :: mom_sp
    type(scalar_field), allocatable, dimension(:,:,:) :: mom_3d
 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(mom_sp, mom_3d) 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (mom_sp, mom_3d) 
# 260 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    ! chem_params: auto-generated in generated_decls.fpp
 
    !> @name Physical bubble parameters (see Ando 2010, Preston 2007)
    !> @{
    real(wp) :: phi_vg, phi_gv, pe_c, tw, k_vl, k_gl
 
# 268 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 268 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(phi_vg, phi_gv, Pe_c, Tw, k_vl, k_gl) 
# 268 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 268 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (phi_vg, phi_gv, Pe_c, Tw, k_vl, k_gl) 
# 268 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp), dimension(:), allocatable :: pb0, mass_g0, mass_v0, pe_t, k_v, k_g
    real(wp), dimension(:), allocatable :: re_trans_t, re_trans_c, im_trans_t, im_trans_c, omegan
 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(pb0, mass_g0, mass_v0, Pe_T, k_v, k_g) 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (pb0, mass_g0, mass_v0, Pe_T, k_v, k_g) 
# 272 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 273 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 273 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(Re_trans_T, Re_trans_c, Im_trans_T, Im_trans_c, omegaN) 
# 273 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 273 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (Re_trans_T, Re_trans_c, Im_trans_T, Im_trans_c, omegaN) 
# 273 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp) :: gam, gam_m
 
# 276 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 276 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(gam, gam_m) 
# 276 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 276 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (gam, gam_m) 
# 276 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp) :: p0ref, rho0ref, t0ref, ss, pv, vd, mu_l, mu_v, mu_g, gam_v, gam_g, m_v, m_g, cp_v, cp_g, r_v, r_g
 
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(p0ref, rho0ref, T0ref, ss, pv, vd, mu_l, mu_v, mu_g, gam_v, gam_g, M_v, M_g, cp_v, cp_g, R_v, R_g) 
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (p0ref, rho0ref, T0ref, ss, pv, vd, mu_l, mu_v, mu_g, gam_v, gam_g, M_v, M_g, cp_v, cp_g, R_v, R_g) 
# 279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    ! acoustic: auto-generated in generated_decls.fpp
 
    !> @name Surface tension parameters
    !> @{
    !> @}
 
    real(wp)                                    :: mytime     !< Current simulation time
    real(wp)                                    :: finaltime  !< Final simulation time
    type(pres_field), allocatable, dimension(:) :: pb_ts
    type(pres_field), allocatable, dimension(:) :: mv_ts
 
 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(mytime, pb_ts, mv_ts) 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (mytime, pb_ts, mv_ts) 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name lagrangian subgrid bubble parameters
    !> lag_params: auto-generated in generated_decls.fpp
    !> @{!
    ! lag_params (decl + GPU_DECLARE) auto-generated in generated_decls.fpp; bubbles_lagrange GPU-declared in
    ! m_global_parameters_common
    integer :: n_el_bubs_loc, n_el_bubs_glb  !< Number of Lagrangian bubbles (local and global)
    logical :: moving_lag_bubbles
    logical :: lag_pressure_force
    logical :: lag_gravity_force
    integer :: lag_vel_model, lag_drag_model
 
# 305 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 305 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(n_el_bubs_loc, n_el_bubs_glb) 
# 305 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 305 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (n_el_bubs_loc, n_el_bubs_glb) 
# 305 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 306 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 306 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(moving_lag_bubbles, lag_vel_model, lag_drag_model) 
# 306 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 306 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (moving_lag_bubbles, lag_vel_model, lag_drag_model) 
# 306 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 307 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 307 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(lag_pressure_force, lag_gravity_force) 
# 307 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 307 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (lag_pressure_force, lag_gravity_force) 
# 307 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    !> @name Continuum damage model parameters
    !> @{!
    !> @}
 
    !> @name MHD Hyperbolic cleaning parameters
    !> @{!
    !> @}
 
contains
 
    !> Assigns default values to the user inputs before reading them in. This enables for an easier consistency check of these
    !! parameters once they are read from the input file.
    impure subroutine s_assign_default_values_to_user_inputs
 
        integer :: i, j  !< Generic loop iterator
 
        ! Shared defaults (case_dir, m/n/p, cyl_coord, cfl flags, model_eqns, BC blocks,
        ! recon/weno/muscl/num_fluids/igr/mhd/relativity under case-opt guard, Tait EOS, bubble flags,
        ! IB flags, parallel I/O flags, fft_wrt)
 
        call s_assign_common_defaults
 
        ! Boundary conditions (bc_x/y/z are per-target declarations, not visible in common)
        bc_x%beg = dflt_int; bc_x%end = dflt_int
        bc_y%beg = dflt_int; bc_y%end = dflt_int
        bc_z%beg = dflt_int; bc_z%end = dflt_int
 
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb1 = 0._wp
                bc_x%ve1 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb2 = 0._wp
                bc_x%ve2 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb3 = 0._wp
                bc_x%ve3 = 0._wp
# 342 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb1 = 0._wp
                bc_y%ve1 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb2 = 0._wp
                bc_y%ve2 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb3 = 0._wp
                bc_y%ve3 = 0._wp
# 342 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb1 = 0._wp
                bc_z%ve1 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb2 = 0._wp
                bc_z%ve2 = 0._wp
# 339 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb3 = 0._wp
                bc_z%ve3 = 0._wp
# 342 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_x%isothermal_in = .false.
            bc_x%isothermal_out = .false.
            bc_x%Twall_in = dflt_real
            bc_x%Twall_out = dflt_real
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_y%isothermal_in = .false.
            bc_y%isothermal_out = .false.
            bc_y%Twall_in = dflt_real
            bc_y%Twall_out = dflt_real
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_z%isothermal_in = .false.
            bc_z%isothermal_out = .false.
            bc_z%Twall_in = dflt_real
            bc_z%Twall_out = dflt_real
# 350 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        call s_update_cell_bounds(cells_bounds, m, n, p)
 
        ! Logistics (sim-specific)
        run_time_info = .false.
        t_step_old = dflt_int
 
        ! Computational domain parameters (sim-specific)
        dt = dflt_real
        cfl_dt = .false.
        cfl_target = dflt_real
 
        t_step_stop = dflt_int
        t_step_save = dflt_int
        t_step_print = 1
 
        t_stop = dflt_real
        t_save = dflt_real
 
        ! NVIDIA UVM options
        nv_uvm_out_of_core = .false.
        nv_uvm_igr_temps_on_gpu = 3  ! => jac, jac_rhs, and jac_old on GPU (default)
        nv_uvm_pref_gpu = .false.
 
        ! Simulation algorithm parameters (sim-specific)
        mpp_lim = .false.
        time_stepper = dflt_int
        muscl_eps = dflt_real
        weno_eps = dflt_real
        teno_ct = dflt_real
        mp_weno = .false.
        weno_avg = .false.
        weno_re_flux = .false.
        riemann_hypo_adc = .false.
        adc_kappa = 1.0_wp
        hll_u_interface = .false.
        hypo_hll_interface_rhs = .false.
        hypo_nc_mode = hypo_nc_mode_none
        adv_src_mode = adv_src_mode_unset
        use_nc_iface_vel = .false.
        low_mach = 0
        wave_speeds = dflt_int
        null_weights = .false.
        precision = 2
        palpha_eps = dflt_real
        ptgalpha_eps = dflt_real
        int_comp = 0
        ic_eps = dflt_ic_eps
        ic_beta = dflt_ic_beta
        rdma_mpi = .false.
        shear_stress = .false.
        bulk_stress = .false.
        any_non_newtonian = .false.
        num_igr_iters = dflt_num_igr_iters
        num_igr_warm_start_iters = dflt_num_igr_warm_start_iters
        alf_factor = dflt_alf_factor
 
# 408 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            mapped_weno = .false.
            wenoz = .false.
            teno = .false.
            wenoz_q = dflt_real
            igr_order = dflt_int
            igr_pres_lim = .false.
            igr_iter_solver = 1
# 416 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        chem_params%diffusion = .false.
        chem_params%reactions = .false.
        chem_params%gamma_method = 1
        chem_params%transport_model = 1
 
        chem_params%reaction_substeps = 0
        chem_params%adap_substeps = .false.
        chem_params%reaction_substeps_max = 0
 
        num_bc_patches = 0
        bc_io = .false.
        periodic_bc = .false.
 
        ! bc_x/y/z (incl. vb/ve loop) already defaulted above; glb_bounds is #1290's grid-derived global extent
        glb_bounds(1)%beg = dflt_real; glb_bounds(1)%end = dflt_real
        glb_bounds(2)%beg = dflt_real; glb_bounds(2)%end = dflt_real
        glb_bounds(3)%beg = dflt_real; glb_bounds(3)%end = dflt_real
 
        bf_spatial_support = .false.
        spatial_bf%amp = 0._wp
        spatial_bf%x_centroid = 0._wp
        spatial_bf%y_centroid = 0._wp
        spatial_bf%conv_vel = 0._wp
        spatial_bf%sigma = 0._wp
        do i = 1, 8
            spatial_bf%freq(i) = 0._wp
            spatial_bf%phase(i) = 0._wp
        end do
 
        ! Fluids physical parameters (sim-specific; Re(:) and G=0._wp differ from post)
        do i = 1, num_fluids_max
            fluid_pp(i)%gamma = dflt_real
            fluid_pp(i)%pi_inf = dflt_real
            fluid_pp(i)%cv = 0._wp
            fluid_pp(i)%qv = 0._wp
            fluid_pp(i)%qvp = 0._wp
            fluid_pp(i)%Re(:) = dflt_real
            fluid_pp(i)%G = 0._wp
            fluid_pp(i)%non_newtonian = .false.
            fluid_pp(i)%K = dflt_real
            fluid_pp(i)%nn = dflt_real
            fluid_pp(i)%tau0 = 0._wp
            fluid_pp(i)%hb_m = dflt_real
            fluid_pp(i)%mu_min = dflt_real
            fluid_pp(i)%mu_max = dflt_real
            fluid_pp(i)%mu_bulk = dflt_real
        end do
 
        ! Subgrid bubble parameters (bub_pp struct + scalar companions; scalar companions are
        ! per-target manual declarations not in m_global_parameters_common scope)
        bub_pp%R0ref = dflt_real; r0ref = dflt_real
        bub_pp%p0ref = dflt_real; p0ref = dflt_real
        bub_pp%rho0ref = dflt_real; rho0ref = dflt_real
        bub_pp%T0ref = dflt_real; t0ref = dflt_real
        bub_pp%ss = dflt_real; ss = dflt_real
        bub_pp%pv = dflt_real; pv = dflt_real
        bub_pp%vd = dflt_real; vd = dflt_real
        bub_pp%mu_l = dflt_real; mu_l = dflt_real
        bub_pp%mu_v = dflt_real; mu_v = dflt_real
        bub_pp%mu_g = dflt_real; mu_g = dflt_real
        bub_pp%gam_v = dflt_real; gam_v = dflt_real
        bub_pp%gam_g = dflt_real; gam_g = dflt_real
        bub_pp%M_v = dflt_real; m_v = dflt_real
        bub_pp%M_g = dflt_real; m_g = dflt_real
        bub_pp%k_v = dflt_real
        bub_pp%k_g = dflt_real
        bub_pp%cp_v = dflt_real; cp_v = dflt_real
        bub_pp%cp_g = dflt_real; cp_g = dflt_real
        bub_pp%R_v = dflt_real; r_v = dflt_real
        bub_pp%R_g = dflt_real; r_g = dflt_real
 
        ! Immersed Boundaries (sim-specific extras)
        ib_neighborhood_radius = 0
        collision_model = 0
        coefficient_of_restitution = dflt_real
        collision_time = dflt_real
        ib_coefficient_of_friction = dflt_real
        ib_state_wrt = .false.
        many_ib_patch_parallelism = .false.
 
        ! Bubble modeling (sim-specific)
        bubble_model = 1
        polytropic = .true.
        thermal = dflt_int
 
# 503 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            nb = 1
            muscl_lim = dflt_int
# 506 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        adv_n = .false.
        adap_dt = .false.
        adap_dt_tol = dflt_adap_dt_tol
        adap_dt_max_iters = dflt_adap_dt_max_iters
 
        pi_fac = 1._wp
 
        eu = dflt_real
        ca = dflt_real
        re_inv = dflt_real
        web = dflt_real
 
        ! Acoustic source
        acoustic_source = .false.
        num_source = dflt_int
 
        bodyforces = .false.
        bf_x = .false.; bf_y = .false.; bf_z = .false.
        !> amplitude, frequency, and phase shift sinusoid in each direction
# 527 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_x = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_x = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_x = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_x = dflt_real
# 530 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 527 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_y = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_y = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_y = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_y = dflt_real
# 530 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 527 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_z = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_z = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_z = dflt_real
# 528 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_z = dflt_real
# 530 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 531 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        synthetic_turbulence = .false.
        synth_seed = 1234
        synth_n_shells = dflt_int
        num_turbulent_sources = 0
        synth_u_inf = dflt_real
        synth_n_waves_per_shell = 0
        synth_k_shell = dflt_real
        synth_amp_shell = dflt_real
        turb_pos = dflt_real
        synth_l = dflt_real
 
        do j = 1, num_probes_max
            acoustic(j)%pulse = dflt_int
            acoustic(j)%support = dflt_int
            acoustic(j)%dipole = .false.
            do i = 1, 3
                acoustic(j)%loc(i) = dflt_real
            end do
            acoustic(j)%mag = dflt_real
            acoustic(j)%length = dflt_real
            acoustic(j)%height = dflt_real
            acoustic(j)%wavelength = dflt_real
            acoustic(j)%frequency = dflt_real
            acoustic(j)%gauss_sigma_dist = dflt_real
            acoustic(j)%gauss_sigma_time = dflt_real
            acoustic(j)%npulse = dflt_real
            acoustic(j)%dir = dflt_real
            acoustic(j)%delay = dflt_real
            acoustic(j)%foc_length = dflt_real
            acoustic(j)%aperture = dflt_real
            acoustic(j)%element_spacing_angle = dflt_real
            acoustic(j)%element_polygon_ratio = dflt_real
            acoustic(j)%rotate_angle = dflt_real
            acoustic(j)%num_elements = dflt_int
            acoustic(j)%element_on = dflt_int
            acoustic(j)%bb_num_freq = dflt_int
            acoustic(j)%bb_lowest_freq = dflt_real
            acoustic(j)%bb_bandwidth = dflt_real
        end do
 
        fd_order = dflt_int
        probe_wrt = .false.
        num_probes = dflt_int
 
        do i = 1, num_probes_max
            probe(i)%x = dflt_real
            probe(i)%y = dflt_real
            probe(i)%z = dflt_real
        end do
 
        ! GRCBC flags
# 584 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_x%grcbc_in = .false.
            bc_x%grcbc_out = .false.
            bc_x%grcbc_vel_out = .false.
# 584 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_y%grcbc_in = .false.
            bc_y%grcbc_out = .false.
            bc_y%grcbc_vel_out = .false.
# 584 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            bc_z%grcbc_in = .false.
            bc_z%grcbc_out = .false.
            bc_z%grcbc_vel_out = .false.
# 588 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        ! Lagrangian subgrid bubble model
        lag_params%solver_approach = dflt_int
        lag_params%cluster_type = dflt_int
        lag_params%pressure_corrector = .false.
        lag_params%smooth_type = dflt_int
        lag_params%heatTransfer_model = .false.
        lag_params%massTransfer_model = .false.
        lag_params%write_bubbles = .false.
        lag_params%write_bubbles_stats = .false.
        lag_params%write_void_evol = .false.
        lag_params%nBubs_glb = dflt_int
        lag_params%vel_model = dflt_int
        lag_params%drag_model = dflt_int
        lag_params%pressure_force = .true.
        lag_params%gravity_force = .false.
        lag_params%kahan_summation = .true.
        lag_params%epsilonb = 1._wp
        lag_params%charwidth = dflt_real
        lag_params%charNz = dflt_int
        lag_params%valmaxvoid = dflt_real
        lag_params%input_path = 'input/lag_bubbles.dat'
        moving_lag_bubbles = .false.
        lag_vel_model = dflt_int
 
        ! Continuum damage model
        tau_star = dflt_real
        cont_damage_s = dflt_real
        alpha_bar = dflt_real
 
        ! MHD (sim-specific extras beyond common Bx0)
        hyper_cleaning_speed = dflt_real
        hyper_cleaning_tau = dflt_real
 
        do i = 1, num_ib_airfoils_max
            ib_airfoil(i)%c = dflt_real
            ib_airfoil(i)%p = dflt_real
            ib_airfoil(i)%t = dflt_real
            ib_airfoil(i)%m = dflt_real
            ib_airfoil_grids(i)%Np = 0
        end do
 
        num_particle_clouds = 0
        do i = 1, num_particle_clouds_max
            particle_cloud(i)%x_centroid = 0._wp
            particle_cloud(i)%y_centroid = 0._wp
            particle_cloud(i)%z_centroid = 0._wp
            particle_cloud(i)%length_x = dflt_real
            particle_cloud(i)%length_y = dflt_real
            particle_cloud(i)%length_z = dflt_real
            particle_cloud(i)%num_particles = 0
            particle_cloud(i)%radius = dflt_real
            particle_cloud(i)%mass = dflt_real
            particle_cloud(i)%min_spacing = 0._wp
            particle_cloud(i)%shell_inner_radius = dflt_real
            particle_cloud(i)%shell_outer_radius = dflt_real
            particle_cloud(i)%moving_ibm = 0
            particle_cloud(i)%seed = 0
            particle_cloud(i)%cloud_geometry = 1
            particle_cloud(i)%packing_method = dflt_int
            particle_cloud(i)%periodic = 0
        end do
 
        do i = 1, num_ib_patches_max_namelist
            patch_ib(i)%gbl_patch_id = i
            patch_ib(i)%geometry = dflt_int
            patch_ib(i)%x_centroid = 0._wp
            patch_ib(i)%y_centroid = 0._wp
            patch_ib(i)%z_centroid = 0._wp
            patch_ib(i)%length_x = dflt_real
            patch_ib(i)%length_y = dflt_real
            patch_ib(i)%length_z = dflt_real
            patch_ib(i)%radius = dflt_real
            patch_ib(i)%airfoil_id = 0
            patch_ib(i)%model_id = 0
            patch_ib(i)%slip = .false.
            patch_ib(i)%v_blow = 0._wp
            patch_ib(i)%inj_species = 0
            patch_ib(i)%burn_rate_exp = 0._wp
            patch_ib(i)%burn_rate_pref = 0._wp
 
            ! Variables to handle moving immersed boundaries, defaulting to no movement
            patch_ib(i)%moving_ibm = 0
            patch_ib(i)%vel(:) = 0._wp
            patch_ib(i)%angles(:) = 0._wp
            patch_ib(i)%angular_vel(:) = 0._wp
            patch_ib(i)%mass = dflt_real
            patch_ib(i)%moment = dflt_real
            patch_ib(i)%centroid_offset(:) = 0._wp
 
            ! sets values of a rotation matrix which can be used when calculating rotations
            patch_ib(i)%rotation_matrix = 0._wp
            patch_ib(i)%rotation_matrix(1, 1) = 1._wp
            patch_ib(i)%rotation_matrix(2, 2) = 1._wp
            patch_ib(i)%rotation_matrix(3, 3) = 1._wp
            patch_ib(i)%rotation_matrix_inverse = patch_ib(i)%rotation_matrix
        end do
 
        num_stl_models = 0
 
        do i = 1, num_stl_models_max
            stl_models(i)%model_filepath(:) = dflt_char
            stl_models(i)%model_translate(:) = 0._wp
            stl_models(i)%model_scale(:) = 1._wp
            stl_models(i)%model_threshold = ray_tracing_threshold
        end do
 
    end subroutine s_assign_default_values_to_user_inputs
 
    !> Initialize the global parameters module
    impure subroutine s_initialize_global_parameters_module
 
        integer :: i, j, k
        integer :: fac
 
# 704 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            ! Determining the degree of the WENO polynomials
            if (recon_type == recon_type_weno) then
                weno_polyn = (weno_order - 1)/2
                if (teno) then
                    weno_num_stencils = weno_order - 3
                else
                    weno_num_stencils = weno_polyn
                end if
            else if (recon_type == recon_type_muscl) then
                muscl_polyn = muscl_order
            end if
 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(weno_polyn, muscl_polyn) 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(weno_polyn, muscl_polyn) 
# 715 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 716 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 716 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(weno_num_stencils) 
# 716 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 716 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(weno_num_stencils) 
# 716 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 717 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 717 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(nb) 
# 717 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 717 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(nb) 
# 717 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 718 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 718 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(num_dims, num_vels, num_fluids) 
# 718 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 718 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(num_dims, num_vels, num_fluids) 
# 718 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 719 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 719 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(igr, igr_order, igr_iter_solver) 
# 719 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 719 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(igr, igr_order, igr_iter_solver) 
# 719 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 721 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        ! muscl_eps: use per-limiter defaults when user did not set it
        if (f_is_default(muscl_eps)) then
            if (muscl_lim <= 2) then
                muscl_eps = 1e-9_wp  ! minmod, MC
            else
                muscl_eps = 1e-6_wp  ! Van Albada, Van Leer, SUPERBEE
            end if
        end if
 
        ! Initialize counts: viscous fluids, surface-tension interfaces, curvature interfaces
        re_size = 0
        re_size_max = 0
 
        ! Populate eqn_idx, sys_size, shear_* (shared logic)
        call s_initialize_eqn_idx(nmom, nb, six_eqn_alf_is_advected=.true.)
 
        ! sim-only: GPU update for shear state after s_initialize_eqn_idx populated it
        if (model_eqns == model_eqns_5eq .or. model_eqns == model_eqns_6eq) then
            if (hypoelasticity) then
 
# 741 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 741 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 741 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 741 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 741 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
            end if
        end if
 
        ! Per-target (sim): nmomsp/nmomtot for qbmm, qbmm_idx alloc/fill, gam, Re_idx
        if (model_eqns == model_eqns_5eq .and. bubbles_euler) then
            if (qbmm) then
                nmomsp = 4  ! number of special moments
                if (nnode == 4) nmomtot = nmom*nb
            end if
 
#ifdef MFC_DEBUG
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:752: ', '@:ALLOCATE(qbmm_idx%rs(nb), qbmm_idx%vs(nb))'
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (qbmm_idx%rs(nb), qbmm_idx%vs(nb))
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(qbmm_idx%rs, qbmm_idx%vs) 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:qbmm_idx%rs, qbmm_idx%vs) 
# 752 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:753: ', '@:ALLOCATE(qbmm_idx%ps(nb), qbmm_idx%ms(nb))'
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (qbmm_idx%ps(nb), qbmm_idx%ms(nb))
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(qbmm_idx%ps, qbmm_idx%ms) 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:qbmm_idx%ps, qbmm_idx%ms) 
# 753 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
            gam = bub_pp%gam_g
 
            if (qbmm) then
#ifdef MFC_DEBUG
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:758: ', '@:ALLOCATE(qbmm_idx%moms(nb, nmom))'
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (qbmm_idx%moms(nb, nmom))
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(qbmm_idx%moms) 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:qbmm_idx%moms) 
# 758 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
                do i = 1, nb
                    do j = 1, nmom
                        qbmm_idx%moms(i, j) = eqn_idx%bub%beg + (j - 1) + (i - 1)*nmom
                    end do
                    qbmm_idx%rs(i) = qbmm_idx%moms(i, 2)
                    qbmm_idx%vs(i) = qbmm_idx%moms(i, 3)
                end do
            else
                do i = 1, nb
                    if (.not. polytropic) then
                        fac = 4
                    else
                        fac = 2
                    end if
 
                    qbmm_idx%rs(i) = eqn_idx%bub%beg + (i - 1)*fac
                    qbmm_idx%vs(i) = qbmm_idx%rs(i) + 1
 
                    if (.not. polytropic) then
                        qbmm_idx%ps(i) = qbmm_idx%vs(i) + 1
                        qbmm_idx%ms(i) = qbmm_idx%ps(i) + 1
                    end if
                end do
            end if
        end if
 
        ! sim-only: Re_idx (non-gamma-law models only)
        if (model_eqns /= model_eqns_gamma_law) then
            ! Count fluids with non-negligible viscous effects (Re > 0)
            do i = 1, num_fluids
                if (fluid_pp(i)%Re(1) > 0) re_size(1) = re_size(1) + 1
                if (fluid_pp(i)%Re(2) > 0) re_size(2) = re_size(2) + 1
            end do
 
            if (re_size(1) > 0._wp) shear_stress = .true.
            if (re_size(2) > 0._wp) bulk_stress = .true.
 
            re_size_max = maxval(re_size)
 
 
# 798 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 798 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(Re_size, Re_size_max, shear_stress, bulk_stress) 
# 798 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 798 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(Re_size, Re_size_max, shear_stress, bulk_stress) 
# 798 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
            ! Bookkeeping the indexes of any viscous fluids
            if (viscous) then
#ifdef MFC_DEBUG
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:802: ', '@:ALLOCATE(Re_idx(1:2, 1:Re_size_max))'
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (re_idx(1:2, 1:re_size_max))
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(Re_idx) 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:Re_idx) 
# 802 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
                k = 0
                do i = 1, num_fluids
                    if (fluid_pp(i)%Re(1) > 0) then
                        k = k + 1; re_idx(1, k) = i
                    end if
                end do
 
                k = 0
                do i = 1, num_fluids
                    if (fluid_pp(i)%Re(2) > 0) then
                        k = k + 1; re_idx(2, k) = i
                    end if
                end do
            end if
        end if
 
        ! Herschel-Bulkley non-Newtonian viscosity: gather per-fluid parameters into device arrays
#ifdef MFC_DEBUG
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:821: ', '@:ALLOCATE(is_non_newtonian(1:num_fluids))'
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (is_non_newtonian(1:num_fluids))
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(is_non_newtonian) 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:is_non_newtonian) 
# 821 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:822: ', '@:ALLOCATE(hb_tau0(1:num_fluids), hb_K(1:num_fluids), hb_nn(1:num_fluids), hb_m_arr(1:num_fluids))'
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (hb_tau0(1:num_fluids), hb_k(1:num_fluids), hb_nn(1:num_fluids), hb_m_arr(1:num_fluids))
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(hb_tau0, hb_K, hb_nn, hb_m_arr) 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:hb_tau0, hb_K, hb_nn, hb_m_arr) 
# 822 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:823: ', '@:ALLOCATE(hb_mu_min(1:num_fluids), hb_mu_max(1:num_fluids))'
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (hb_mu_min(1:num_fluids), hb_mu_max(1:num_fluids))
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(hb_mu_min, hb_mu_max) 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:hb_mu_min, hb_mu_max) 
# 823 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:824: ', '@:ALLOCATE(fluid_inv_re(1:num_fluids))'
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (fluid_inv_re(1:num_fluids))
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(fluid_inv_re) 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:fluid_inv_re) 
# 824 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        any_non_newtonian = .false.
        do i = 1, num_fluids
            is_non_newtonian(i) = fluid_pp(i)%non_newtonian
            if (is_non_newtonian(i)) any_non_newtonian = .true.
            hb_tau0(i) = fluid_pp(i)%tau0
            hb_k(i) = fluid_pp(i)%K
            hb_nn(i) = fluid_pp(i)%nn
            hb_m_arr(i) = fluid_pp(i)%hb_m
            hb_mu_min(i) = fluid_pp(i)%mu_min
            hb_mu_max(i) = fluid_pp(i)%mu_max
            if (fluid_pp(i)%Re(1) > 0._wp) then
                fluid_inv_re(i) = 1._wp/fluid_pp(i)%Re(1)
            else
                fluid_inv_re(i) = 0._wp
            end if
        end do
 
# 842 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 842 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(any_non_newtonian, is_non_newtonian, hb_tau0, hb_K, hb_nn, hb_m_arr, hb_mu_min, hb_mu_max, fluid_inv_re) 
# 842 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 842 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(any_non_newtonian, is_non_newtonian, hb_tau0, hb_K, hb_nn, hb_m_arr, hb_mu_min, hb_mu_max, fluid_inv_re) 
# 842 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (bubbles_euler .and. qbmm .and. .not. polytropic) then
            allocate (mpi_io_data%view(1:sys_size + 2*nb*nnode))
            allocate (mpi_io_data%var(1:sys_size + 2*nb*nnode))
        else if (bubbles_lagrange) then
            allocate (mpi_io_data%view(1:sys_size + 1))
            allocate (mpi_io_data%var(1:sys_size + 1))
        else
            allocate (mpi_io_data%view(1:sys_size))
            allocate (mpi_io_data%var(1:sys_size))
        end if
 
        if (.not. down_sample) then
            do i = 1, sys_size
                allocate (mpi_io_data%var(i)%sf(0:m,0:n,0:p))
                mpi_io_data%var(i)%sf => null()
            end do
        end if
        if (bubbles_euler .and. qbmm .and. .not. polytropic) then
            do i = sys_size + 1, sys_size + 2*nb*nnode
                allocate (mpi_io_data%var(i)%sf(0:m,0:n,0:p))
                mpi_io_data%var(i)%sf => null()
            end do
        else if (bubbles_lagrange) then
            do i = 1, sys_size + 1
                allocate (mpi_io_data%var(i)%sf(0:m,0:n,0:p))
                mpi_io_data%var(i)%sf => null()
            end do
        end if
 
        ! Configure WENO averaging flag (arithmetic mean vs. unaltered values)
        wa_flg = 0._wp; if (weno_avg) wa_flg = 1._wp
 
# 875 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 875 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wa_flg) 
# 875 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 875 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wa_flg) 
# 875 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        ! Resort to default WENO-JS if no other WENO scheme is selected
# 879 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            wenojs = .not. (mapped_weno .or. wenoz .or. teno)
# 881 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        if (ib) allocate (mpi_io_ib_data%var%sf(0:m,0:n,0:p))
 
        if (hypoelasticity .or. mhd .or. probe_wrt .or. ib .or. bubbles_lagrange) then
            fd_number = max(1, fd_order/2)
        end if
 
        hypo_nc_mode = hypo_nc_mode_none
        if (hypoelasticity) then
            if (riemann_solver == 1) then
                if (hypo_hll_interface_rhs) then
                    hypo_nc_mode = hypo_nc_mode_interface
                else
                    hypo_nc_mode = hypo_nc_mode_finite_diff
                end if
            else if (riemann_solver == 2) then
                hypo_nc_mode = hypo_nc_mode_interface
            else if (riemann_solver == 4) then
                hypo_nc_mode = hypo_nc_mode_dual_pass
            end if
        end if
 
        ! flux_src: choose exactly one export mode (adv_src_mode) for the NC volume fraction advection term.
        if (riemann_solver == 1 .and. .not. hll_u_interface) then
            ! HLL Method 1 (alpha-interface): flux_src(adv_idx%beg:adv_idx%end) carries interface alpha_k per fluid.
            adv_src_mode = adv_src_mode_alpha_iface
        else if ((riemann_solver == 1 .and. hll_u_interface) .or. riemann_solver == 2 .or. riemann_solver == 3 &
                 & .or. riemann_solver == 5) then
            ! HLLC, HLL Method 2 (u-interface), exact, LF: flux_src(adv_idx%beg) carries one shared face-normal velocity.
            adv_src_mode = adv_src_mode_vel_iface
        else if (riemann_solver == 4) then
            ! MHD HLLD: single species, no volume fraction to advect. Hypo HLLD: the dual-pass keeps all NC terms in the flux.
            adv_src_mode = adv_src_mode_none
        end if
 
        call s_configure_coordinate_bounds(recon_type, weno_polyn, muscl_polyn, igr_order, buff_size, idwint, idwbuff, viscous, &
                                           & bubbles_lagrange, m, n, p, num_dims, igr, ib, fd_number)
 
# 918 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 918 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(idwint, idwbuff) 
# 918 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 918 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(idwint, idwbuff) 
# 918 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        ! Configuring Coordinate Direction Indexes
        if (bubbles_euler) then
#ifdef MFC_DEBUG
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:922: ', '@:ALLOCATE(ptil( idwbuff(1)%beg:idwbuff(1)%end, idwbuff(2)%beg:idwbuff(2)%end, idwbuff(3)%beg:idwbuff(3)%end))'
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (ptil( idwbuff(1)%beg:idwbuff(1)%end, idwbuff(2)%beg:idwbuff(2)%end, idwbuff(3)%beg:idwbuff(3)%end))
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(ptil) 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:ptil) 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
        end if
 
 
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(fd_order, fd_number) 
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(fd_order, fd_number) 
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (cyl_coord .neqv. .true.) then  ! Cartesian grid
            grid_geometry = 1
        else if (cyl_coord .and. p == 0) then  ! Axisymmetric cylindrical grid
            grid_geometry = 2
        else
            grid_geometry = 3
        end if
 
        ! nc_iface_vel: use_nc_iface_vel enables a second export channel. Use it when the Riemann solver must expose interface
        ! velocities beyond what flux_src already provides:
        !
        ! 1. adv_src_mode_alpha_iface + alt_soundspeed: face-normal velocity only, for the KdivU correction (flux_src already
        ! carries alpha in this mode) 2. hypo_nc_mode_interface: all components for the hypoelastic velocity-gradient tensor
        ! 3. hypo_nc_mode_dual_pass + axisym: anchored radial face traces for the cylindrical completion (both velocity
        ! components are exported per face; the completion consumes the radial one from each pass)
        use_nc_iface_vel = hypo_nc_mode == hypo_nc_mode_interface .or. (hypo_nc_mode == hypo_nc_mode_dual_pass &
            & .and. grid_geometry == 2) .or. (adv_src_mode == adv_src_mode_alpha_iface .and. alt_soundspeed)
 
 
# 945 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 945 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(sys_size, buff_size, eqn_idx, adv_n, adap_dt, pi_fac, adap_dt_tol, adap_dt_max_iters) 
# 945 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 945 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(sys_size, buff_size, eqn_idx, adv_n, adap_dt, pi_fac, adap_dt_tol, adap_dt_max_iters) 
# 945 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 946 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 946 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(cfl_target, m, n, p) 
# 946 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 946 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(cfl_target, m, n, p) 
# 946 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 948 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 948 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(alt_soundspeed, acoustic_source, num_source) 
# 948 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 948 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(alt_soundspeed, acoustic_source, num_source) 
# 948 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 949 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 949 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(dt, sys_size, buff_size, eqn_idx, mpp_lim, bubbles_euler, hypoelasticity, alt_soundspeed, avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, low_Mach) 
# 949 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 949 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(dt, sys_size, buff_size, eqn_idx, mpp_lim, bubbles_euler, hypoelasticity, alt_soundspeed, avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, low_Mach) 
# 949 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(riemann_hypo_ADC, ADC_kappa, hll_u_interface, hypo_hll_interface_rhs, hypo_nc_mode) 
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(riemann_hypo_ADC, ADC_kappa, hll_u_interface, hypo_hll_interface_rhs, hypo_nc_mode) 
# 951 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 953 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 953 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(Bx0) 
# 953 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 953 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(Bx0) 
# 953 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 955 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 955 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(chem_params) 
# 955 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 955 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(chem_params) 
# 955 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 957 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 957 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(rburn) 
# 957 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 957 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(rburn) 
# 957 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 959 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 959 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(cont_damage, tau_star, cont_damage_s, alpha_bar) 
# 959 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 959 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(cont_damage, tau_star, cont_damage_s, alpha_bar) 
# 959 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 961 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 961 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(hyper_cleaning, hyper_cleaning_speed, hyper_cleaning_tau) 
# 961 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 961 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(hyper_cleaning, hyper_cleaning_speed, hyper_cleaning_tau) 
# 961 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wenojs, mapped_weno, wenoz, teno) 
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wenojs, mapped_weno, wenoz, teno) 
# 964 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wenoz_q) 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wenoz_q) 
# 965 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 966 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 966 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(mhd, relativity) 
# 966 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 966 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(mhd, relativity) 
# 966 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(muscl_order, muscl_lim) 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(muscl_order, muscl_lim) 
# 967 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 968 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 968 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(igr, igr_order) 
# 968 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 968 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(igr, igr_order) 
# 968 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 969 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 969 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(num_fluids, num_dims, viscous, num_vels, nb, muscl_lim) 
# 969 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 969 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(num_fluids, num_dims, viscous, num_vels, nb, muscl_lim) 
# 969 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 971 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
 
# 972 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 972 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(int_comp, ic_eps, ic_beta) 
# 972 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 972 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(int_comp, ic_eps, ic_beta) 
# 972 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 973 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 973 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(muscl_eps) 
# 973 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 973 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(muscl_eps) 
# 973 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 974 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 974 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(dir_idx, dir_flg, dir_idx_tau, stress_perm) 
# 974 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 974 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(dir_idx, dir_flg, dir_idx_tau, stress_perm) 
# 974 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 976 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 976 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 976 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 976 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 976 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (synthetic_turbulence) then
 
# 979 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 979 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(synthetic_turbulence, num_turbulent_sources) 
# 979 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 979 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(synthetic_turbulence, num_turbulent_sources) 
# 979 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 980 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 980 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(synth_U_inf, synth_n_waves_per_shell, synth_k_shell, synth_amp_shell) 
# 980 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 980 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(synth_U_inf, synth_n_waves_per_shell, synth_k_shell, synth_amp_shell) 
# 980 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 981 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 981 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(turb_pos, synth_L) 
# 981 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 981 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(turb_pos, synth_L) 
# 981 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
        end if
 
        ! Allocating grid variables for the x-, y- and z-directions
#ifdef MFC_DEBUG
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:985: ', '@:ALLOCATE(x_cb(-1 - buff_size:m + buff_size))'
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (x_cb(-1 - buff_size:m + buff_size))
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(x_cb) 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:x_cb) 
# 985 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:986: ', '@:ALLOCATE(x_cc(-buff_size:m + buff_size))'
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (x_cc(-buff_size:m + buff_size))
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(x_cc) 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:x_cc) 
# 986 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:987: ', '@:ALLOCATE(dx(-buff_size:m + buff_size))'
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dx(-buff_size:m + buff_size))
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dx) 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dx) 
# 987 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving x_cb to GPU => ", SHAPE(x_cb) set preferred location GPU
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cb), sizeof(x_cb), cudamemadvisesetpreferredlocation, 0)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cb), sizeof(x_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(x_cb), sizeof(x_cb), 0, 0)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 988 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving x_cc to GPU => ", SHAPE(x_cc) set preferred location GPU
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cc), sizeof(x_cc), cudamemadvisesetpreferredlocation, 0)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cc), sizeof(x_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(x_cc), sizeof(x_cc), 0, 0)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 989 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving dx to GPU => ", SHAPE(dx) set preferred location GPU
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dx), sizeof(dx), cudamemadvisesetpreferredlocation, 0)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dx), sizeof(dx), cudamemadvisesetaccessedby, cudacpudeviceid)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dx), sizeof(dx), 0, 0)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 990 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (n == 0) return
#ifdef MFC_DEBUG
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:993: ', '@:ALLOCATE(y_cb(-1 - buff_size:n + buff_size))'
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (y_cb(-1 - buff_size:n + buff_size))
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(y_cb) 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:y_cb) 
# 993 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:994: ', '@:ALLOCATE(y_cc(-buff_size:n + buff_size))'
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (y_cc(-buff_size:n + buff_size))
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(y_cc) 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:y_cc) 
# 994 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:995: ', '@:ALLOCATE(dy(-buff_size:n + buff_size))'
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dy(-buff_size:n + buff_size))
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dy) 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dy) 
# 995 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving y_cb to GPU => ", SHAPE(y_cb) set preferred location GPU
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cb), sizeof(y_cb), cudamemadvisesetpreferredlocation, 0)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cb), sizeof(y_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(y_cb), sizeof(y_cb), 0, 0)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 996 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving y_cc to GPU => ", SHAPE(y_cc) set preferred location GPU
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cc), sizeof(y_cc), cudamemadvisesetpreferredlocation, 0)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cc), sizeof(y_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(y_cc), sizeof(y_cc), 0, 0)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 997 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving dy to GPU => ", SHAPE(dy) set preferred location GPU
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dy), sizeof(dy), cudamemadvisesetpreferredlocation, 0)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dy), sizeof(dy), cudamemadvisesetaccessedby, cudacpudeviceid)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dy), sizeof(dy), 0, 0)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 998 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (p == 0) return
#ifdef MFC_DEBUG
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1001: ', '@:ALLOCATE(z_cb(-1 - buff_size:p + buff_size))'
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (z_cb(-1 - buff_size:p + buff_size))
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(z_cb) 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:z_cb) 
# 1001 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1002: ', '@:ALLOCATE(z_cc(-buff_size:p + buff_size))'
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (z_cc(-buff_size:p + buff_size))
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(z_cc) 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:z_cc) 
# 1002 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1003: ', '@:ALLOCATE(dz(-buff_size:p + buff_size))'
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dz(-buff_size:p + buff_size))
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dz) 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dz) 
# 1003 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving z_cb to GPU => ", SHAPE(z_cb) set preferred location GPU
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cb), sizeof(z_cb), cudamemadvisesetpreferredlocation, 0)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cb), sizeof(z_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(z_cb), sizeof(z_cb), 0, 0)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1004 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving z_cc to GPU => ", SHAPE(z_cc) set preferred location GPU
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cc), sizeof(z_cc), cudamemadvisesetpreferredlocation, 0)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cc), sizeof(z_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(z_cc), sizeof(z_cc), 0, 0)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1005 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! NVIDIA CUDA Fortran 25.3+: uses submodules (cuda_runtime_api, gpu_reductions, sort) See
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        ! https://docs.nvidia.com/hpc-sdk/compilers/cuda-fortran-prog-guide/index.html#fortran-host-modules
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! print*, "Moving dz to GPU => ", SHAPE(dz) set preferred location GPU
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dz), sizeof(dz), cudamemadvisesetpreferredlocation, 0)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dz), sizeof(dz), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dz), sizeof(dz), 0, 0)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    ! write(*,*) cudaGetErrorString(istat)
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1006 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    end subroutine s_initialize_global_parameters_module
 
    !> Initializes parallel infrastructure
    impure subroutine s_initialize_parallel_io
 
        ! proc_coords/start_idx/mpiiofs/mpi_info_int setup moved into the shared routine
        call s_initialize_parallel_io_common
 
        ! #1290: per-rank physical comm-domain bounds for Lagrangian-bubble exchange
#ifdef MFC_DEBUG
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1017: ', '@:ALLOCATE(pcomm_coords(1:num_dims))'
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (pcomm_coords(1:num_dims))
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(pcomm_coords) 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:pcomm_coords) 
# 1017 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    end subroutine s_initialize_parallel_io
 
    !> Module deallocation and/or disassociation procedures
    impure subroutine s_finalize_global_parameters_module
 
        integer :: i
 
        ! Deallocating the variables bookkeeping the indexes of any viscous fluids and any pairs of fluids whose interfaces
        ! supported effects of surface tension
 
        if (viscous) then
#ifdef MFC_DEBUG
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1030: ', '@:DEALLOCATE(Re_idx)'
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(Re_idx) 
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:Re_idx) 
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1030 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (re_idx)
        end if
 
        ! Herschel-Bulkley non-Newtonian viscosity arrays (always allocated)
#ifdef MFC_DEBUG
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1034: ', '@:DEALLOCATE(is_non_newtonian)'
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(is_non_newtonian) 
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:is_non_newtonian) 
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1034 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (is_non_newtonian)
#ifdef MFC_DEBUG
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1035: ', '@:DEALLOCATE(hb_tau0, hb_K, hb_nn, hb_m_arr)'
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(hb_tau0, hb_K, hb_nn, hb_m_arr) 
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:hb_tau0, hb_K, hb_nn, hb_m_arr) 
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1035 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (hb_tau0, hb_k, hb_nn, hb_m_arr)
#ifdef MFC_DEBUG
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1036: ', '@:DEALLOCATE(hb_mu_min, hb_mu_max)'
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(hb_mu_min, hb_mu_max) 
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:hb_mu_min, hb_mu_max) 
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1036 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (hb_mu_min, hb_mu_max)
#ifdef MFC_DEBUG
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1037: ', '@:DEALLOCATE(fluid_inv_re)'
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(fluid_inv_re) 
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:fluid_inv_re) 
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1037 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (fluid_inv_re)
 
        if (bubbles_euler) then
#ifdef MFC_DEBUG
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1040: ', '@:DEALLOCATE(ptil)'
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(ptil) 
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:ptil) 
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1040 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (ptil)
#ifdef MFC_DEBUG
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1041: ', '@:DEALLOCATE(qbmm_idx%rs, qbmm_idx%vs, qbmm_idx%ps, qbmm_idx%ms)'
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(qbmm_idx%rs, qbmm_idx%vs, qbmm_idx%ps, qbmm_idx%ms) 
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:qbmm_idx%rs, qbmm_idx%vs, qbmm_idx%ps, qbmm_idx%ms) 
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1041 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (qbmm_idx%rs, qbmm_idx%vs, qbmm_idx%ps, qbmm_idx%ms)
            if (qbmm) then
#ifdef MFC_DEBUG
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1043: ', '@:DEALLOCATE(qbmm_idx%moms)'
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(qbmm_idx%moms) 
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:qbmm_idx%moms) 
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1043 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (qbmm_idx%moms)
            end if
        end if
 
#ifdef MFC_DEBUG
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1047: ', '@:DEALLOCATE(pcomm_coords)'
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(pcomm_coords) 
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:pcomm_coords) 
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1047 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (pcomm_coords)
 
        ! Shared: deallocate proc_coords and start_idx
        call s_finalize_global_parameters_common
 
        if (parallel_io) then
            if (bubbles_lagrange) then
                do i = 1, sys_size + 1
                    mpi_io_data%var(i)%sf => null()
                end do
            else
                do i = 1, sys_size
                    mpi_io_data%var(i)%sf => null()
                end do
            end if
 
            deallocate (mpi_io_data%var)
            deallocate (mpi_io_data%view)
        end if
 
        if (ib) mpi_io_ib_data%var%sf => null()
 
        ! Deallocating grid variables for the x-, y- and z-directions
#ifdef MFC_DEBUG
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1070: ', '@:DEALLOCATE(x_cb, x_cc, dx)'
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(x_cb, x_cc, dx) 
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:x_cb, x_cc, dx) 
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1070 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (x_cb, x_cc, dx)
 
        if (n == 0) return
#ifdef MFC_DEBUG
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1073: ', '@:DEALLOCATE(y_cb, y_cc, dy)'
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(y_cb, y_cc, dy) 
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:y_cb, y_cc, dy) 
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1073 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (y_cb, y_cc, dy)
 
        if (p == 0) return
#ifdef MFC_DEBUG
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1076: ', '@:DEALLOCATE(z_cb, z_cc, dz)'
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(z_cb, z_cc, dz) 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:z_cb, z_cc, dz) 
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1076 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (z_cb, z_cc, dz)
 
        if (allocated(neighbor_ranks)) then
#ifdef MFC_DEBUG
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1079: ', '@:DEALLOCATE(neighbor_ranks)'
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(neighbor_ranks) 
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:neighbor_ranks) 
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1079 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (neighbor_ranks)
        end if
 
    end subroutine s_finalize_global_parameters_module
 
end module m_global_parameters