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.
# 9 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
 
! For moving immersed boundaries in simulation
# 14 "/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"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 86 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 96 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 101 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 106 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 116 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 121 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 207 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 232 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 243 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 245 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 283 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 293 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 303 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 312 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 329 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 339 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 346 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 358 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 364 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 370 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 376 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
! New line at end of file is required for FYPP
# 3 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
# 1 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 1
# 1 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 1
# 2 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 3 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 4 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 8 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 9 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 10 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 46 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 98 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp" 2
 
# 7 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 17 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 22 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 27 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 32 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 42 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 47 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 52 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 57 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 62 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 73 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 78 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 83 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 103 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 131 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 160 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 175 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 192 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 213 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 241 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 256 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 266 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 275 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 291 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 301 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 308 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
# 21 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 37 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 50 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 76 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 91 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 102 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 115 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 143 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 154 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 165 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 176 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 187 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 198 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 208 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 214 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 220 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 226 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 232 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 234 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 235 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution:
! This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI rank.
! That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0.
! For an example see misc/nvidia_uvm/bind.sh.
# 63 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 88 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 111 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 153 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 159 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 167 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
! New line at end of file is required for FYPP
# 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        !< definitions of the derived types
 
    use m_helper_basic         !< functions to compare floating point numbers
 
    ! $:USE_GPU_MODULE()
 
    implicit none
 
    real(wp) :: wall_time = 0
    real(wp) :: wall_time_avg = 0
 
    ! Logistics
    integer :: num_procs             !< Number of processors
    character(LEN=path_len) :: case_dir              !< Case folder location
    logical :: run_time_info         !< Run-time output flag
    integer :: t_step_old            !< Existing IC/grid folder
 
    ! Computational Domain Parameters
    integer :: proc_rank !< Rank of the local processor
 
    !> @name Number of cells in the x-, y- and z-directions, respectively
    !> @{
    integer :: m, n, p
    !> @}
 
    !> @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)
    !> @{
    logical :: cyl_coord
    integer :: grid_geometry
    !> @}
 
# 53 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 53 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(cyl_coord, grid_geometry) 
# 53 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 53 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (cyl_coord, grid_geometry) 
# 53 "/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
    !> @}
 
    !> @name Cell-center (CC) locations in the x-, y- and z-directions, respectively
    !> @{
 
    real(wp), target, allocatable, dimension(:) :: x_cc, y_cc, z_cc
    !> @}
    !type(bounds_info) :: x_domain, y_domain, z_domain !<
    !! Locations of the domain bounds in the x-, y- and z-coordinate directions
    !> @name Cell-width distributions in the x-, y- and z-directions, respectively
    !> @{
 
    real(wp), target, allocatable, dimension(:) :: dx, dy, dz
    !> @}
 
    real(wp) :: dt !< Size of the time-step
 
 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 76 "/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, dt, m, n, p) 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 76 "/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, dt, m, n, p) 
# 76 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Starting time-step iteration, stopping time-step iteration and the number
    !! of time-step iterations between successive solution backups, respectively
    !> @{
    integer :: t_step_start, t_step_stop, t_step_save
    !> @}
 
    !> @name Starting time, stopping time, and time between backups, simulation time,
    !! and prescribed cfl respectively
    !> @{
    real(wp) :: t_stop, t_save, cfl_target
    integer :: n_start
    !> @}
 
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(cfl_target) 
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (cfl_target) 
# 90 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: cfl_adap_dt, cfl_const_dt, cfl_dt
 
    integer :: t_step_print !< Number of time-steps between printouts
 
    ! Simulation Algorithm Parameters
    integer :: model_eqns     !< Multicomponent flow model
# 102 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: num_dims       !< Number of spatial dimensions
        integer :: num_vels       !< Number of velocity components (different from num_dims for mhd)
# 105 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    logical :: mpp_lim        !< Mixture physical parameters (MPP) limits
    integer :: time_stepper   !< Time-stepper algorithm
    logical :: prim_vars_wrt
 
# 131 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: recon_type     !< Reconstruction Type
        integer :: weno_polyn     !< Degree of the WENO polynomials (polyn)
        integer :: muscl_polyn    !< Degree of the MUSCL polynomials (polyn)i
        integer :: weno_order     !< Order of the WENO reconstruction
        integer :: muscl_order    !< Order of the MUSCL reconstruction
        integer :: weno_num_stencils    !< Number of stencils for WENO reconstruction (only different from weno_polyn for TENO(>5))
        integer :: muscl_lim      !< MUSCL Limiter
        integer :: num_fluids     !< number of fluids in the simulation
        logical :: wenojs         !< WENO-JS (default)
        logical :: mapped_weno    !< WENO-M (WENO with mapping of nonlinear weights)
        logical :: wenoz          !< WENO-Z
        logical :: teno           !< TENO (Targeted ENO)
        real(wp) :: wenoz_q       !< Power constant for WENO-Z
        logical :: mhd            !< Magnetohydrodynamics
        logical :: relativity     !< Relativity (only for MHD)
        integer :: igr_iter_solver!< IGR elliptic solver
        integer :: igr_order      !< Reconstruction order for IGR
        logical :: igr            !< Use information geometric regularization
        logical :: igr_pres_lim   !< Limit to positive pressures for IGR
        logical :: viscous        !< Viscous effects
# 152 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
    !> @name Variables for our of core IGR computation on NVIDIA
    !> @{
    logical :: nv_uvm_out_of_core ! Enable out-of-core storage of q_cons_ts(2) in timestepping (default FALSE)
    integer :: nv_uvm_igr_temps_on_gpu ! 0 => jac, jac_rhs, and jac_old on CPU
    ! 1 => jac on GPU, jac_rhs and jac_old on CPU
    ! 2 => jac and jac_rhs on GPU, jac_old on CPU
    ! 3 => jac, jac_rhs, and jac_old on GPU (default)
    logical :: nv_uvm_pref_gpu ! Enable explicit gpu memory hints (default FALSE)
    !> @}
 
    real(wp) :: weno_eps       !< Binding for the WENO nonlinear weights
    real(wp) :: teno_ct        !< Smoothness threshold for TENO
    logical :: mp_weno        !< Monotonicity preserving (MP) WENO
    logical :: weno_avg       ! Average left/right cell-boundary states
    logical :: weno_re_flux   !< WENO reconstruct velocity gradients for viscous stress tensor
    integer :: riemann_solver !< Riemann solver algorithm
    integer :: low_mach       !< Low Mach number fix to HLLC Riemann solver
    integer :: wave_speeds    !< Wave speeds estimation method
    integer :: avg_state      !< Average state evaluation method
    logical :: alt_soundspeed !< Alternate mixture sound speed
    logical :: null_weights    !< Null undesired WENO weights
    logical :: mixture_err     !< Mixture properties correction
    logical :: hypoelasticity  !< hypoelasticity modeling
    logical :: hyperelasticity !< hyperelasticity modeling
    logical :: int_comp        !< THINC interface compression
    real(wp) :: ic_eps         !< THINC Epsilon to compress on surface cells
    real(wp) :: ic_beta        !< THINC Sharpness Parameter
    integer :: hyper_model     !< hyperelasticity solver algorithm
    logical :: elasticity      !< elasticity modeling, true for hyper or hypo
    logical, parameter :: chemistry = .false. !< Chemistry modeling
    logical :: shear_stress  !< Shear stresses
    logical :: bulk_stress   !< Bulk stresses
    logical :: cont_damage   !< Continuum damage modeling
    logical :: hyper_cleaning !< Hyperbolic cleaning for MHD for divB=0
    integer :: num_igr_iters !< number of iterations for elliptic solve
    integer :: num_igr_warm_start_iters !< number of warm start iterations for elliptic solve
    real(wp) :: alf_factor  !< alpha factor for IGR
 
    logical :: bodyforces
    logical :: bf_x, bf_y, bf_z !< body force toggle in three directions
    !< amplitude, frequency, and phase shift sinusoid in each direction
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: k_x
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: p_x
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: w_x
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: g_x
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: k_y
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: p_y
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: w_y
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: g_y
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 195 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: k_z
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: p_z
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: w_z
# 196 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            real(wp) :: g_z
# 198 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 199 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    real(wp), dimension(3) :: accel_bf
 
# 200 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 200 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(accel_bf) 
# 200 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 200 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (accel_bf) 
# 200 "/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]')
 
    integer :: cpu_start, cpu_end, cpu_rate
 
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(num_dims, num_vels, weno_polyn, weno_order) 
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (num_dims, num_vels, weno_polyn, weno_order) 
# 206 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 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(weno_num_stencils, num_fluids, wenojs) 
# 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 (weno_num_stencils, num_fluids, wenojs) 
# 207 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(mapped_weno, wenoz, teno, wenoz_q, mhd, relativity) 
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (mapped_weno, wenoz, teno, wenoz_q, mhd, relativity) 
# 208 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(igr_iter_solver, igr_order, viscous, igr_pres_lim, igr) 
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (igr_iter_solver, igr_order, viscous, igr_pres_lim, igr) 
# 209 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(recon_type, muscl_order, muscl_polyn, muscl_lim) 
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (recon_type, muscl_order, muscl_polyn, muscl_lim) 
# 210 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 212 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(mpp_lim, model_eqns, mixture_err, alt_soundspeed) 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (mpp_lim, model_eqns, mixture_err, alt_soundspeed) 
# 213 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(avg_state, mp_weno, weno_eps, teno_CT, hypoelasticity) 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (avg_state, mp_weno, weno_eps, teno_CT, hypoelasticity) 
# 214 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 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(hyperelasticity, hyper_model, elasticity, low_Mach) 
# 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 (hyperelasticity, hyper_model, elasticity, low_Mach) 
# 215 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 216 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 216 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(shear_stress, bulk_stress, cont_damage, hyper_cleaning) 
# 216 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 216 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (shear_stress, bulk_stress, cont_damage, hyper_cleaning) 
# 216 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: relax          !< activate phase change
    integer :: relax_model    !< Relaxation model
    real(wp) :: palpha_eps     !< trigger parameter for the p relaxation procedure, phase change model
    real(wp) :: ptgalpha_eps   !< trigger parameter for the pTg relaxation procedure, phase change model
 
 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (relax, relax_model, palpha_eps, ptgalpha_eps) 
# 223 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: num_bc_patches
    logical :: bc_io
    !> @name Boundary conditions (BC) in the x-, y- and z-directions, respectively
    !> @{
    type(int_bounds_info) :: bc_x, bc_y, bc_z
    !> @}
#if defined(MFC_OpenACC)
 
# 232 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 232 "/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) 
# 232 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 232 "/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) 
# 232 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 233 "/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) 
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 233 "/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) 
# 233 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 234 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 234 "/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) 
# 234 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 234 "/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) 
# 234 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#elif defined(MFC_OpenMP)
 
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bc_x, bc_y, bc_z) 
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bc_x, bc_y, bc_z) 
# 236 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#endif
    type(bounds_info) :: x_domain, y_domain, z_domain
    real(wp) :: x_a, y_a, z_a
    real(wp) :: x_b, y_b, z_b
 
    logical :: parallel_io !< Format of the data files
    logical :: file_per_process !< shared file or not when using parallel io
    integer :: precision !< Precision of output files
    logical :: down_sample !< down sample the output files
 
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(down_sample) 
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (down_sample) 
# 246 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer, allocatable, dimension(:) :: proc_coords !<
    !! Processor coordinates in MPI_CART_COMM
 
    integer, allocatable, dimension(:) :: start_idx !<
    !! Starting cell-center index of local processor in global grid
 
    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
 
    !> @name MPI info for parallel IO with Lustre file systems
    !> @{
    character(LEN=name_len) :: mpiiofs
    integer :: mpi_info_int
    !> @}
 
    !> @name Annotations of the structure of the state and flux vectors in terms of the
    !! size and the configuration of the system of equations to which they belong
    !> @{
    integer :: sys_size                                !< Number of unknowns in system of eqns.
    type(int_bounds_info) :: cont_idx                  !< Indexes of first & last continuity eqns.
    type(int_bounds_info) :: mom_idx                   !< Indexes of first & last momentum eqns.
    integer :: e_idx                                   !< Index of energy equation
    integer :: n_idx                                   !< Index of number density
    type(int_bounds_info) :: adv_idx                   !< Indexes of first & last advection eqns.
    type(int_bounds_info) :: internalenergies_idx      !< Indexes of first & last internal energy eqns.
    type(bub_bounds_info) :: bub_idx                   !< Indexes of first & last bubble variable eqns.
    integer :: alf_idx                                 !< Index of void fraction
    integer :: gamma_idx                               !< Index of specific heat ratio func. eqn.
    integer :: pi_inf_idx                              !< Index of liquid stiffness func. eqn.
    type(int_bounds_info) :: b_idx                     !< Indexes of first and last magnetic field eqns.
    type(int_bounds_info) :: stress_idx                !< Indexes of first and last shear stress eqns.
    type(int_bounds_info) :: xi_idx                    !< Indexes of first and last reference map eqns.
    integer :: b_size                                  !< Number of elements in the symmetric b tensor, plus one
    integer :: tensor_size                             !< Number of elements in the full tensor plus one
    type(int_bounds_info) :: species_idx               !< Indexes of first & last concentration eqns.
    integer :: c_idx                                   !< Index of color function
    integer :: damage_idx                              !< Index of damage state variable (D) for continuum damage model
    integer :: psi_idx                                 !< Index of hyperbolic cleaning state variable for MHD
    !> @}
 
# 291 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 291 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(sys_size, E_idx, n_idx, bub_idx, alf_idx, gamma_idx) 
# 291 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 291 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (sys_size, E_idx, n_idx, bub_idx, alf_idx, gamma_idx) 
# 291 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 292 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 292 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(pi_inf_idx, B_idx, stress_idx, xi_idx, b_size) 
# 292 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 292 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (pi_inf_idx, B_idx, stress_idx, xi_idx, b_size) 
# 292 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 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(tensor_size, species_idx, c_idx) 
# 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 (tensor_size, species_idx, c_idx) 
# 293 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! 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)
 
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(idwint) 
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 298 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (idwint) 
# 298 "/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)
 
# 304 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 304 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(idwbuff) 
# 304 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 304 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (idwbuff) 
# 304 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name The number of fluids, along with their identifying indexes, respectively,
    !! for which viscous effects, e.g. the shear and/or the volume Reynolds (Re)
    !! numbers, will be non-negligible.
    !> @{
    integer, dimension(2) :: re_size
    integer :: re_size_max
    integer, allocatable, dimension(:, :) :: re_idx
    !> @}
 
 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(Re_size, Re_size_max, Re_idx) 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (Re_size, Re_size_max, Re_idx) 
# 315 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! The WENO average (WA) flag regulates whether the calculation of any cell-
    ! average spatial derivatives is carried out in each cell by utilizing the
    ! arithmetic mean of the left and right, WENO-reconstructed, cell-boundary
    ! values or simply, the unaltered left and right, WENO-reconstructed, cell-
    ! boundary values.
    !> @{
    real(wp) :: wa_flg
    !> @}
 
 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(wa_flg) 
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (wa_flg) 
# 326 "/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.
    !> @{
    integer, dimension(3) :: dir_idx
    real(wp), dimension(3) :: dir_flg
    integer, dimension(3) :: dir_idx_tau !!used for hypoelasticity=true
    !> @}
 
 
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(dir_idx, dir_flg, dir_idx_tau) 
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (dir_idx, dir_flg, dir_idx_tau) 
# 338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: buff_size !<
    !! The number of cells that are necessary to be able to store enough boundary
    !! conditions data to march the solution in the physical computational domain
    !! to the next time-step.
 
 
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(buff_size) 
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (buff_size) 
# 345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: shear_num !! Number of shear stress components
    integer, dimension(3) :: shear_indices !<
    !! Indices of the stress components that represent shear stress
    integer :: shear_bc_flip_num !<
    !! Number of shear stress components to reflect for boundary conditions
    integer, dimension(3, 2) :: shear_bc_flip_indices !<
    !! Indices of shear stress components to reflect for boundary conditions.
    !! Size: (1:3, 1:shear_BC_flip_num) for (x/y/z, [indices])
 
 
# 356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    ! END: Simulation Algorithm Parameters
 
    ! Fluids Physical Parameters
 
    type(physical_parameters), dimension(num_fluids_max) :: fluid_pp !<
    !! Database of the physical parameters of each of the fluids that is present
    !! in the flow. These include the stiffened gas equation of state parameters,
    !! and the Reynolds numbers.
 
    ! Subgrid Bubble Parameters
    type(subgrid_bubble_physical_parameters) :: bub_pp
 
    integer :: fd_order !<
    !! The order of the finite-difference (fd) approximations of the first-order
    !! derivatives that need to be evaluated when the CoM or flow probe data
    !! files are to be written at each time step
 
    integer :: fd_number !<
    !! The finite-difference number is given by MAX(1, fd_order/2). Essentially,
    !! it is a measure of the half-size of the finite-difference stencil for the
    !! selected order of accuracy.
 
# 379 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 379 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(fd_order, fd_number) 
# 379 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 379 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (fd_order, fd_number) 
# 379 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: probe_wrt
    logical :: integral_wrt
    integer :: num_probes
    integer :: num_integrals
    type(vec3_dt), dimension(num_probes_max) :: probe
    type(integral_parameters), dimension(num_probes_max) :: integral
 
    !> @name Reference density and pressure for Tait EOS
    !> @{
    real(wp) :: rhoref, pref
    !> @}
 
# 392 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 392 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(rhoref, pref) 
# 392 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 392 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (rhoref, pref) 
# 392 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Immersed Boundaries
    !> @{
    logical :: ib
    integer :: num_ibs
 
    type(ib_patch_parameters), dimension(num_patches_max) :: patch_ib
    type(vec3_dt), allocatable, dimension(:) :: airfoil_grid_u, airfoil_grid_l
    integer :: np
    !! Database of the immersed boundary patch parameters for each of the
    !! patches employed in the configuration of the initial condition. Note that
    !! the maximum allowable number of patches, num_patches_max, may be changed
    !! in the module m_derived_types.f90.
 
 
# 407 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 407 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ib, num_ibs, patch_ib, Np, airfoil_grid_u, airfoil_grid_l) 
# 407 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 407 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ib, num_ibs, patch_ib, Np, airfoil_grid_u, airfoil_grid_l) 
# 407 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    !> @name Bubble modeling
    !> @{
# 415 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: nb       !< Number of eq. bubble sizes
# 417 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
    real(wp) :: eu       !< Euler number
    real(wp) :: ca       !< Cavitation number
    real(wp) :: web      !< Weber number
    real(wp) :: re_inv   !< Inverse Reynolds number
 
# 422 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 422 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(Eu, Ca, Web, Re_inv) 
# 422 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 422 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (Eu, Ca, Web, Re_inv) 
# 422 "/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
 
# 426 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 426 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(weight, R0) 
# 426 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 426 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (weight, R0) 
# 426 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: bubbles_euler      !< Bubbles euler on/off
    logical :: polytropic   !< Polytropic  switch
    logical :: polydisperse !< Polydisperse bubbles
 
# 431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bubbles_euler, polytropic, polydisperse) 
# 431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bubbles_euler, polytropic, polydisperse) 
# 431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: adv_n        !< Solve the number density equation and compute alpha from number density
    logical :: adap_dt      !< Adaptive step size control
    real(wp) :: adap_dt_tol !< Tolerance to control adaptive step size
    integer :: adap_dt_max_iters !< Maximum number of iterations
 
# 437 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 437 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(adv_n, adap_dt, adap_dt_tol, adap_dt_max_iters) 
# 437 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 437 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (adv_n, adap_dt, adap_dt_tol, adap_dt_max_iters) 
# 437 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    integer :: bubble_model !< Gilmore or Keller--Miksis bubble model
    integer :: thermal      !< Thermal behavior. 1 = adiabatic, 2 = isotherm, 3 = transfer
 
# 441 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 441 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bubble_model, thermal) 
# 441 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 441 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bubble_model, thermal) 
# 441 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp), allocatable, dimension(:, :, :) :: ptil  !< Pressure modification
 
    real(wp) :: poly_sigma  !< log normal sigma for polydisperse PDF
 
# 446 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 446 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(ptil, poly_sigma) 
# 446 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 446 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (ptil, poly_sigma) 
# 446 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: qbmm      !< Quadrature moment method
    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
 
    real(wp) :: pi_fac   !< Factor for artificial pi_inf
 
# 454 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 454 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(qbmm, nmomsp, nmomtot, pi_fac) 
# 454 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 454 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (qbmm, nmomsp, nmomtot, pi_fac) 
# 454 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(nb) 
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (nb) 
# 457 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 459 "/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
 
# 462 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 462 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(mom_sp, mom_3d) 
# 462 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 462 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (mom_sp, mom_3d) 
# 462 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @}
 
    type(chemistry_parameters) :: chem_params
 
# 467 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 467 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(chem_params) 
# 467 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 467 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (chem_params) 
# 467 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Physical bubble parameters (see Ando 2010, Preston 2007)
    !> @{
    real(wp) :: phi_vg, phi_gv, pe_c, tw, k_vl, k_gl
 
# 472 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 472 "/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) 
# 472 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 472 "/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) 
# 472 "/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
 
# 476 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 476 "/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) 
# 476 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 476 "/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) 
# 476 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 477 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 477 "/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) 
# 477 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 477 "/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) 
# 477 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp) :: gam, gam_m
 
# 480 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 480 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(gam, gam_m) 
# 480 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 480 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (gam, gam_m) 
# 480 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp) :: r0ref, 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
 
# 484 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 484 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(R0ref, 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) 
# 484 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 484 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (R0ref, 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) 
# 484 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 486 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    !> @}
 
    !> @name Acoustic acoustic_source parameters
    !> @{
    logical :: acoustic_source !< Acoustic source switch
    type(acoustic_parameters), dimension(num_probes_max) :: acoustic !< Acoustic source parameters
    integer :: num_source !< Number of acoustic sources
    !> @}
 
# 494 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 494 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(acoustic_source, acoustic, num_source) 
# 494 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 494 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (acoustic_source, acoustic, num_source) 
# 494 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name Surface tension parameters
    !> @{
 
    real(wp) :: sigma
    logical :: surface_tension
 
# 501 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 501 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(sigma, surface_tension) 
# 501 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 501 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (sigma, surface_tension) 
# 501 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    integer :: momxb, momxe
    integer :: advxb, advxe
    integer :: contxb, contxe
    integer :: intxb, intxe
    integer :: bubxb, bubxe
    integer :: strxb, strxe
    integer :: chemxb, chemxe
    integer :: xibeg, xiend
 
# 512 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 512 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(momxb, momxe, advxb, advxe, contxb, contxe) 
# 512 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 512 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (momxb, momxe, advxb, advxe, contxb, contxe) 
# 512 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 513 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 513 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(intxb, intxe, bubxb, bubxe) 
# 513 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 513 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (intxb, intxe, bubxb, bubxe) 
# 513 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 514 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 514 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(strxb, strxe, chemxb, chemxe) 
# 514 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 514 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (strxb, strxe, chemxb, chemxe) 
# 514 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 515 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 515 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(xibeg, xiend) 
# 515 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 515 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (xibeg, xiend) 
# 515 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp), allocatable, dimension(:) :: gammas, gs_min, pi_infs, ps_inf, cvs, qvs, qvps
 
# 518 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 518 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(gammas, gs_min, pi_infs, ps_inf, cvs, qvs, qvps) 
# 518 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 518 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (gammas, gs_min, pi_infs, ps_inf, cvs, qvs, qvps) 
# 518 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    real(wp) :: mytime       !< Current simulation time
    real(wp) :: finaltime    !< Final simulation time
 
    logical :: rdma_mpi
 
    type(pres_field), allocatable, dimension(:) :: pb_ts
 
    type(pres_field), allocatable, dimension(:) :: mv_ts
 
 
# 529 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 529 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(pb_ts, mv_ts) 
# 529 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 529 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (pb_ts, mv_ts) 
# 529 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    !> @name lagrangian subgrid bubble parameters
    !> @{!
    logical :: bubbles_lagrange                         !< Lagrangian subgrid bubble model switch
    type(bubbles_lagrange_parameters) :: lag_params     !< Lagrange bubbles' parameters
 
# 535 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 535 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(bubbles_lagrange, lag_params) 
# 535 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 535 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (bubbles_lagrange, lag_params) 
# 535 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    real(wp) :: bx0 !< Constant magnetic field in the x-direction (1D)
 
# 539 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 539 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(Bx0) 
# 539 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 539 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (Bx0) 
# 539 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
    logical :: fft_wrt
    logical :: dummy  !< AMDFlang workaround: keep a dummy logical to avoid a compiler case-optimization bug when a parameter+GPU-kernel conditional is false
 
    !> @name Continuum damage model parameters
    !> @{!
    real(wp) :: tau_star        !< Stress threshold for continuum damage modeling
    real(wp) :: cont_damage_s   !< Exponent s for continuum damage modeling
    real(wp) :: alpha_bar       !< Damage rate factor for continuum damage modeling
 
# 549 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 549 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(tau_star, cont_damage_s, alpha_bar) 
# 549 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 549 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (tau_star, cont_damage_s, alpha_bar) 
# 549 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
    !> @name MHD Hyperbolic cleaning parameters
    !> @{!
    real(wp) :: hyper_cleaning_speed    !< Hyperbolic cleaning wave speed (c_h)
    real(wp) :: hyper_cleaning_tau      !< Hyperbolic cleaning tau
 
# 556 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 556 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc declare create(hyper_cleaning_speed, hyper_cleaning_tau) 
# 556 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 556 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp declare target (hyper_cleaning_speed, hyper_cleaning_tau) 
# 556 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
    !> @}
 
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
 
        ! Logistics
        case_dir = '.'
        run_time_info = .false.
        t_step_old = dflt_int
 
        ! Computational domain parameters
        m = dflt_int; n = 0; p = 0
        call s_update_cell_bounds(cells_bounds, m, n, p)
 
        cyl_coord = .false.
 
        dt = dflt_real
 
        cfl_adap_dt = .false.
        cfl_const_dt = .false.
        cfl_dt = .false.
        cfl_target = dflt_real
 
        t_step_start = dflt_int
        t_step_stop = dflt_int
        t_step_save = dflt_int
        t_step_print = 1
 
        n_start = dflt_int
        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
        model_eqns = dflt_int
        mpp_lim = .false.
        time_stepper = dflt_int
        weno_eps = dflt_real
        teno_ct = dflt_real
        mp_weno = .false.
        weno_avg = .false.
        weno_re_flux = .false.
        riemann_solver = dflt_int
        low_mach = 0
        wave_speeds = dflt_int
        avg_state = dflt_int
        alt_soundspeed = .false.
        null_weights = .false.
        mixture_err = .false.
        parallel_io = .false.
        file_per_process = .false.
        precision = 2
        down_sample = .false.
        relax = .false.
        relax_model = dflt_int
        palpha_eps = dflt_real
        ptgalpha_eps = dflt_real
        hypoelasticity = .false.
        hyperelasticity = .false.
        int_comp = .false.
        ic_eps = dflt_ic_eps
        ic_beta = dflt_ic_beta
        elasticity = .false.
        hyper_model = dflt_int
        b_size = dflt_int
        tensor_size = dflt_int
        rdma_mpi = .false.
        shear_stress = .false.
        bulk_stress = .false.
        cont_damage = .false.
        hyper_cleaning = .false.
        num_igr_iters = dflt_num_igr_iters
        num_igr_warm_start_iters = dflt_num_igr_warm_start_iters
        alf_factor = dflt_alf_factor
 
# 643 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            mapped_weno = .false.
            wenoz = .false.
            teno = .false.
            wenoz_q = dflt_real
            igr = .false.
            igr_order = dflt_int
            igr_pres_lim = .false.
            viscous = .false.
            igr_iter_solver = 1
# 653 "/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
 
        num_bc_patches = 0
        bc_io = .false.
 
        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
 
# 667 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb1 = 0._wp
                bc_x%ve1 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb2 = 0._wp
                bc_x%ve2 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_x%vb3 = 0._wp
                bc_x%ve3 = 0._wp
# 671 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 667 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb1 = 0._wp
                bc_y%ve1 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb2 = 0._wp
                bc_y%ve2 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_y%vb3 = 0._wp
                bc_y%ve3 = 0._wp
# 671 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 667 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb1 = 0._wp
                bc_z%ve1 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb2 = 0._wp
                bc_z%ve2 = 0._wp
# 668 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                bc_z%vb3 = 0._wp
                bc_z%ve3 = 0._wp
# 671 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 672 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        x_domain%beg = dflt_real; x_domain%end = dflt_real
        y_domain%beg = dflt_real; y_domain%end = dflt_real
        z_domain%beg = dflt_real; z_domain%end = dflt_real
 
        ! Fluids physical parameters
        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
        end do
 
        ! Subgrid bubble parameters
        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
 
        ! Tait EOS
        rhoref = dflt_real
        pref = dflt_real
 
        ! Immersed Boundaries
        ib = .false.
        num_ibs = dflt_int
 
        ! Bubble modeling
        bubbles_euler = .false.
        bubble_model = 1
        polytropic = .true.
        polydisperse = .false.
        thermal = dflt_int
        r0ref = dflt_real
 
# 727 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            nb = 1
            recon_type = weno_type
            weno_order = dflt_int
            muscl_order = dflt_int
            muscl_lim = dflt_int
            num_fluids = dflt_int
# 734 "/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
 
        ! User inputs for qbmm for simulation code
        qbmm = .false.
 
        eu = dflt_real
        ca = dflt_real
        re_inv = dflt_real
        web = dflt_real
        poly_sigma = dflt_real
 
        ! Acoustic source
        acoustic_source = .false.
        num_source = dflt_int
 
        ! Surface tension
        sigma = dflt_real
        surface_tension = .false.
 
        bodyforces = .false.
        bf_x = .false.; bf_y = .false.; bf_z = .false.
        !< amplitude, frequency, and phase shift sinusoid in each direction
# 763 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_x = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_x = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_x = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_x = dflt_real
# 766 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 763 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_y = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_y = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_y = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_y = dflt_real
# 766 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 763 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                k_z = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                p_z = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                w_z = dflt_real
# 764 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                g_z = dflt_real
# 766 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
# 767 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        fft_wrt = .false.
        dummy = .false.
 
        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.
        integral_wrt = .false.
        num_probes = dflt_int
        num_integrals = 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
 
        do i = 1, num_probes_max
            integral(i)%xmin = dflt_real
            integral(i)%xmax = dflt_real
            integral(i)%ymin = dflt_real
            integral(i)%ymax = dflt_real
            integral(i)%ymin = dflt_real
            integral(i)%ymax = dflt_real
        end do
 
        ! GRCBC flags
# 823 "/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.
# 823 "/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.
# 823 "/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.
# 827 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        ! Lagrangian subgrid bubble model
        bubbles_lagrange = .false.
        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%nBubs_glb = dflt_int
        lag_params%epsilonb = 1._wp
        lag_params%charwidth = dflt_real
        lag_params%valmaxvoid = dflt_real
 
        ! Continuum damage model
        tau_star = dflt_real
        cont_damage_s = dflt_real
        alpha_bar = dflt_real
 
        ! MHD
        bx0 = dflt_real
        hyper_cleaning_speed = dflt_real
        hyper_cleaning_tau = dflt_real
 
# 854 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            mhd = .false.
            relativity = .false.
# 857 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        do i = 1, num_patches_max
            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)%theta = dflt_real
            patch_ib(i)%c = dflt_real
            patch_ib(i)%t = dflt_real
            patch_ib(i)%m = dflt_real
            patch_ib(i)%p = dflt_real
            patch_ib(i)%slip = .false.
 
            ! Proper default values for translating STL models
            patch_ib(i)%model_scale(:) = 1._wp
            patch_ib(i)%model_translate(:) = 0._wp
            patch_ib(i)%model_rotate(:) = 0._wp
            patch_ib(i)%model_filepath(:) = dflt_char
            patch_ib(i)%model_spc = num_ray
            patch_ib(i)%model_threshold = ray_tracing_threshold
 
            ! Variables to handle moving imersed 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
 
    end subroutine s_assign_default_values_to_user_inputs
 
    !>  The computation of parameters, the allocation of memory,
        !!      the association of pointers and/or the execution of any
        !!      other procedures that are necessary to setup the module.
    impure subroutine s_initialize_global_parameters_module
 
        integer :: i, j, k
        integer :: fac
 
# 910 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            ! Determining the degree of the WENO polynomials
            if (recon_type == weno_type) then
                weno_polyn = (weno_order - 1)/2
                if (teno) then
                    weno_num_stencils = weno_order - 3
                else
                    weno_num_stencils = weno_polyn
                end if
            elseif (recon_type == muscl_type) then
                muscl_polyn = muscl_order
            end if
 
# 921 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 921 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(weno_polyn, muscl_polyn) 
# 921 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 921 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(weno_polyn, muscl_polyn) 
# 921 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 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 update device(weno_num_stencils) 
# 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 update to(weno_num_stencils) 
# 922 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 923 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 923 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(nb) 
# 923 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 923 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(nb) 
# 923 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 924 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 924 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(num_dims, num_vels, num_fluids) 
# 924 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 924 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(num_dims, num_vels, num_fluids) 
# 924 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 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(igr, igr_order, igr_iter_solver) 
# 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(igr, igr_order, igr_iter_solver) 
# 925 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 927 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        ! Initializing the number of fluids for which viscous effects will
        ! be non-negligible, the number of distinctive material interfaces
        ! for which surface tension will be important and also, the number
        ! of fluids for which the physical and geometric curvatures of the
        ! interfaces will be computed
        re_size = 0
        re_size_max = 0
 
        ! Gamma/Pi_inf Model
        if (model_eqns == 1) then
 
            ! Annotating structure of the state and flux vectors belonging
            ! to the system of equations defined by the selected number of
            ! spatial dimensions and the gamma/pi_inf model
            cont_idx%beg = 1
            cont_idx%end = cont_idx%beg
            mom_idx%beg = cont_idx%end + 1
            mom_idx%end = cont_idx%end + num_vels
            e_idx = mom_idx%end + 1
            adv_idx%beg = e_idx + 1
            adv_idx%end = adv_idx%beg + 1
            gamma_idx = adv_idx%beg
            pi_inf_idx = adv_idx%end
            sys_size = adv_idx%end
 
            ! Volume Fraction Model
        else
 
            ! Annotating structure of the state and flux vectors belonging
            ! to the system of equations defined by the selected number of
            ! spatial dimensions and the volume fraction model
            if (model_eqns == 2) then
                cont_idx%beg = 1
                cont_idx%end = num_fluids
                mom_idx%beg = cont_idx%end + 1
                mom_idx%end = cont_idx%end + num_vels
                e_idx = mom_idx%end + 1
 
                if (igr) then
                    ! Volume fractions are stored in the indices immediately following
                    ! the energy equation. IGR tracks a total of (N-1) volume fractions
                    ! for N fluids, hence the "-1" in adv_idx%end. If num_fluids = 1
                    ! then adv_idx%end < adv_idx%beg, which skips all loops over the
                    ! volume fractions since there is no volume fraction to track
                    adv_idx%beg = e_idx + 1 ! Alpha for fluid 1
                    adv_idx%end = e_idx + num_fluids - 1
                else
                    ! Volume fractions are stored in the indices immediately following
                    ! the energy equation. WENO/MUSCL + Riemann tracks a total of (N)
                    ! volume fractions for N fluids, hence the lack of  "-1" in adv_idx%end
                    adv_idx%beg = e_idx + 1
                    adv_idx%end = e_idx + num_fluids
                end if
 
                sys_size = adv_idx%end
 
                if (bubbles_euler) then
                    alf_idx = adv_idx%end
                else
                    alf_idx = 1
                end if
 
                if (bubbles_euler) then
                    bub_idx%beg = sys_size + 1
                    if (qbmm) then
                        nmomsp = 4 !number of special moments
                        if (nnode == 4) then
                            ! nmom = 6 : It is already a parameter
                            nmomtot = nmom*nb
                        end if
                        bub_idx%end = adv_idx%end + nb*nmom
                    else
                        if (.not. polytropic) then
                            bub_idx%end = sys_size + 4*nb
                        else
                            bub_idx%end = sys_size + 2*nb
                        end if
                    end if
                    sys_size = bub_idx%end
 
                    if (adv_n) then
                        n_idx = bub_idx%end + 1
                        sys_size = n_idx
                    end if
 
#ifdef MFC_DEBUG
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1013: ', '@:ALLOCATE(bub_idx%rs(nb), bub_idx%vs(nb))'
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (bub_idx%rs(nb), bub_idx%vs(nb))
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(bub_idx%rs, bub_idx%vs) 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:bub_idx%rs, bub_idx%vs) 
# 1013 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1014: ', '@:ALLOCATE(bub_idx%ps(nb), bub_idx%ms(nb))'
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (bub_idx%ps(nb), bub_idx%ms(nb))
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(bub_idx%ps, bub_idx%ms) 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:bub_idx%ps, bub_idx%ms) 
# 1014 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
                    gam = bub_pp%gam_g
 
                    if (qbmm) then
#ifdef MFC_DEBUG
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1019: ', '@:ALLOCATE(bub_idx%moms(nb, nmom))'
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (bub_idx%moms(nb, nmom))
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(bub_idx%moms) 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:bub_idx%moms) 
# 1019 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
                        do i = 1, nb
                            do j = 1, nmom
                                bub_idx%moms(i, j) = bub_idx%beg + (j - 1) + (i - 1)*nmom
                            end do
                            bub_idx%rs(i) = bub_idx%moms(i, 2)
                            bub_idx%vs(i) = bub_idx%moms(i, 3)
                        end do
 
                    else
                        do i = 1, nb
                            if (.not. polytropic) then
                                fac = 4
                            else
                                fac = 2
                            end if
 
                            bub_idx%rs(i) = bub_idx%beg + (i - 1)*fac
                            bub_idx%vs(i) = bub_idx%rs(i) + 1
 
                            if (.not. polytropic) then
                                bub_idx%ps(i) = bub_idx%vs(i) + 1
                                bub_idx%ms(i) = bub_idx%ps(i) + 1
                            end if
                        end do
                    end if
                end if
 
                if (mhd) then
                    b_idx%beg = sys_size + 1
                    if (n == 0) then
                        b_idx%end = sys_size + 2 ! 1D: By, Bz
                    else
                        b_idx%end = sys_size + 3 ! 2D/3D: Bx, By, Bz
                    end if
                    sys_size = b_idx%end
                end if
 
            else if (model_eqns == 3) then
                cont_idx%beg = 1
                cont_idx%end = num_fluids
                mom_idx%beg = cont_idx%end + 1
                mom_idx%end = cont_idx%end + num_vels
                e_idx = mom_idx%end + 1
                adv_idx%beg = e_idx + 1
                adv_idx%end = e_idx + num_fluids
                alf_idx = adv_idx%end
                internalenergies_idx%beg = adv_idx%end + 1
                internalenergies_idx%end = adv_idx%end + num_fluids
                sys_size = internalenergies_idx%end
 
            else if (model_eqns == 4) then
                cont_idx%beg = 1 ! one continuity equation
                cont_idx%end = 1 !num_fluids
                mom_idx%beg = cont_idx%end + 1 ! one momentum equation in each direction
                mom_idx%end = cont_idx%end + num_vels
                e_idx = mom_idx%end + 1 ! one energy equation
                adv_idx%beg = e_idx + 1
                adv_idx%end = adv_idx%beg !one volume advection equation
                alf_idx = adv_idx%end
                sys_size = adv_idx%end
 
                if (bubbles_euler) then
                    bub_idx%beg = sys_size + 1
                    bub_idx%end = sys_size + 2*nb
                    if (.not. polytropic) then
                        bub_idx%end = sys_size + 4*nb
                    end if
                    sys_size = bub_idx%end
 
#ifdef MFC_DEBUG
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1089: ', '@:ALLOCATE(bub_idx%rs(nb), bub_idx%vs(nb))'
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (bub_idx%rs(nb), bub_idx%vs(nb))
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(bub_idx%rs, bub_idx%vs) 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:bub_idx%rs, bub_idx%vs) 
# 1089 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1090: ', '@:ALLOCATE(bub_idx%ps(nb), bub_idx%ms(nb))'
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (bub_idx%ps(nb), bub_idx%ms(nb))
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(bub_idx%ps, bub_idx%ms) 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:bub_idx%ps, bub_idx%ms) 
# 1090 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
                    do i = 1, nb
                        if (polytropic) then
                            fac = 2
                        else
                            fac = 4
                        end if
 
                        bub_idx%rs(i) = bub_idx%beg + (i - 1)*fac
                        bub_idx%vs(i) = bub_idx%rs(i) + 1
 
                        if (.not. polytropic) then
                            bub_idx%ps(i) = bub_idx%vs(i) + 1
                            bub_idx%ms(i) = bub_idx%ps(i) + 1
                        end if
                    end do
                end if
            end if
 
            ! Determining the number of fluids for which the shear and the
            ! volume Reynolds numbers, e.g. viscous effects, are important
            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)
 
 
# 1122 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1122 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(Re_size, Re_size_max, shear_stress, bulk_stress) 
# 1122 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1122 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(Re_size, Re_size_max, shear_stress, bulk_stress) 
# 1122 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
            ! Bookkeeping the indexes of any viscous fluids and any pairs of
            ! fluids whose interface will support effects of surface tension
            if (viscous) then
 
#ifdef MFC_DEBUG
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1128: ', '@:ALLOCATE(Re_idx(1:2, 1:Re_size_max))'
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (re_idx(1:2, 1:re_size_max))
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(Re_idx) 
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1128 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:Re_idx) 
# 1128 "/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
 
        if (model_eqns == 2 .or. model_eqns == 3) then
 
            if (hypoelasticity .or. hyperelasticity) then
                elasticity = .true.
                stress_idx%beg = sys_size + 1
                stress_idx%end = sys_size + (num_dims*(num_dims + 1))/2
                if (cyl_coord) stress_idx%end = stress_idx%end + 1
                ! number of stresses is 1 in 1D, 3 in 2D, 4 in 2D-Axisym, 6 in 3D
                sys_size = stress_idx%end
 
                ! shear stress index is 2 for 2D and 2,4,5 for 3D
                if (num_dims == 1) then
                    shear_num = 0
                else if (num_dims == 2) then
                    shear_num = 1
                    shear_indices(1) = stress_idx%beg - 1 + 2
                    shear_bc_flip_num = 1
                    shear_bc_flip_indices(1:2, 1) = shear_indices(1)
                    ! Both x-dir and y-dir: flip tau_xy only
                else if (num_dims == 3) then
                    shear_num = 3
                    shear_indices(1:3) = stress_idx%beg - 1 + (/2, 4, 5/)
                    shear_bc_flip_num = 2
                    shear_bc_flip_indices(1, 1:2) = shear_indices((/1, 2/))
                    shear_bc_flip_indices(2, 1:2) = shear_indices((/1, 3/))
                    shear_bc_flip_indices(3, 1:2) = shear_indices((/2, 3/))
                    ! x-dir: flip tau_xy and tau_xz
                    ! y-dir: flip tau_xy and tau_yz
                    ! z-dir: flip tau_xz and tau_yz
                end if
 
# 1178 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1178 "/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) 
# 1178 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1178 "/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) 
# 1178 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
            end if
 
            if (hyperelasticity) then
                ! number of entries in the symmetric btensor plus the jacobian
                b_size = (num_dims*(num_dims + 1))/2 + 1
                ! storing the jacobian in the last entry
                tensor_size = num_dims**2 + 1
                xi_idx%beg = sys_size + 1
                xi_idx%end = sys_size + num_dims
                ! adding three more equations for the \xi field and the elastic energy
                sys_size = xi_idx%end + 1
            end if
 
            if (surface_tension) then
                c_idx = sys_size + 1
                sys_size = c_idx
            end if
 
            if (cont_damage) then
                damage_idx = sys_size + 1
                sys_size = damage_idx
            end if
 
            if (hyper_cleaning) then
                psi_idx = sys_size + 1
                sys_size = psi_idx
            end if
 
        end if
 
        ! END: Volume Fraction Model
 
        if (chemistry) then
            species_idx%beg = sys_size + 1
            species_idx%end = sys_size + num_species
            sys_size = species_idx%end
        end if
 
        if (bubbles_euler .and. qbmm .and. .not. polytropic) then
            allocate (mpi_io_data%view(1:sys_size + 2*nb*4))
            allocate (mpi_io_data%var(1:sys_size + 2*nb*4))
        elseif (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*4
                allocate (mpi_io_data%var(i)%sf(0:m, 0:n, 0:p))
                mpi_io_data%var(i)%sf => null()
            end do
        elseif (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
 
        ! Configuring the WENO average flag that will be used to regulate
        ! whether any spatial derivatives are to computed in each cell by
        ! using the arithmetic mean of left and right, WENO-reconstructed,
        ! cell-boundary values or otherwise, the unaltered left and right,
        ! WENO-reconstructed, cell-boundary values
        wa_flg = 0._wp; if (weno_avg) wa_flg = 1._wp
 
# 1252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wa_flg) 
# 1252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wa_flg) 
# 1252 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        ! Resort to default WENO-JS if no other WENO scheme is selected
# 1256 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            wenojs = .not. (mapped_weno .or. wenoz .or. teno)
# 1258 "/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))
        np = 0
 
        if (elasticity) then
            fd_number = max(1, fd_order/2)
        end if
 
        if (mhd) then ! TODO merge with above; waiting for hyperelasticity PR
            fd_number = max(1, fd_order/2)
        end if
 
        if (probe_wrt) then
            fd_number = max(1, fd_order/2)
        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)
 
# 1279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(idwint, idwbuff) 
# 1279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(idwint, idwbuff) 
# 1279 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        ! Configuring Coordinate Direction Indexes
        if (bubbles_euler) then
#ifdef MFC_DEBUG
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1283: ', '@:ALLOCATE(ptil( idwbuff(1)%beg:idwbuff(1)%end, idwbuff(2)%beg:idwbuff(2)%end, idwbuff(3)%beg:idwbuff(3)%end))'
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1283 "/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))
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(ptil) 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:ptil) 
# 1283 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1287 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
 
 
# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(fd_order, fd_number) 
# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(fd_order, fd_number) 
# 1289 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (cyl_coord .neqv. .true.) then ! Cartesian grid
            grid_geometry = 1
        elseif (cyl_coord .and. p == 0) then ! Axisymmetric cylindrical grid
            grid_geometry = 2
        else ! Fully 3D cylindrical grid
            grid_geometry = 3
        end if
 
        momxb = mom_idx%beg
        momxe = mom_idx%end
        advxb = adv_idx%beg
        advxe = adv_idx%end
        contxb = cont_idx%beg
        contxe = cont_idx%end
        bubxb = bub_idx%beg
        bubxe = bub_idx%end
        strxb = stress_idx%beg
        strxe = stress_idx%end
        intxb = internalenergies_idx%beg
        intxe = internalenergies_idx%end
        xibeg = xi_idx%beg
        xiend = xi_idx%end
        chemxb = species_idx%beg
        chemxe = species_idx%end
 
 
# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(momxb, momxe, advxb, advxe, contxb, contxe, bubxb, bubxe, intxb, intxe, sys_size, buff_size, E_idx, alf_idx, n_idx, adv_n, adap_dt, pi_fac, strxb, strxe, chemxb, chemxe, c_idx, adap_dt_tol, adap_dt_max_iters) 
# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(momxb, momxe, advxb, advxe, contxb, contxe, bubxb, bubxe, intxb, intxe, sys_size, buff_size, E_idx, alf_idx, n_idx, adv_n, adap_dt, pi_fac, strxb, strxe, chemxb, chemxe, c_idx, adap_dt_tol, adap_dt_max_iters) 
# 1316 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(b_size, xibeg, xiend, tensor_size) 
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(b_size, xibeg, xiend, tensor_size) 
# 1320 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(species_idx) 
# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(species_idx) 
# 1322 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1323 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1323 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(cfl_target, m, n, p) 
# 1323 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1323 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(cfl_target, m, n, p) 
# 1323 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1325 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1325 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(alt_soundspeed, acoustic_source, num_source) 
# 1325 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1325 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(alt_soundspeed, acoustic_source, num_source) 
# 1325 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(dt, sys_size, buff_size, pref, rhoref, gamma_idx, pi_inf_idx, E_idx, alf_idx, stress_idx, mpp_lim, bubbles_euler, hypoelasticity, alt_soundspeed, avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, hyperelasticity, hyper_model, elasticity, xi_idx, B_idx, low_Mach) 
# 1326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(dt, sys_size, buff_size, pref, rhoref, gamma_idx, pi_inf_idx, E_idx, alf_idx, stress_idx, mpp_lim, bubbles_euler, hypoelasticity, alt_soundspeed, avg_state, model_eqns, mixture_err, grid_geometry, cyl_coord, mp_weno, weno_eps, teno_CT, hyperelasticity, hyper_model, elasticity, xi_idx, B_idx, low_Mach) 
# 1326 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1333 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
 
# 1334 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1334 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(Bx0) 
# 1334 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1334 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(Bx0) 
# 1334 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(chem_params) 
# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(chem_params) 
# 1336 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(cont_damage, tau_star, cont_damage_s, alpha_bar) 
# 1338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(cont_damage, tau_star, cont_damage_s, alpha_bar) 
# 1338 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1340 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1340 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(hyper_cleaning, hyper_cleaning_speed, hyper_cleaning_tau) 
# 1340 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1340 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(hyper_cleaning, hyper_cleaning_speed, hyper_cleaning_tau) 
# 1340 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wenojs, mapped_weno, wenoz, teno) 
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wenojs, mapped_weno, wenoz, teno) 
# 1343 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1344 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1344 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(wenoz_q) 
# 1344 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1344 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(wenoz_q) 
# 1344 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(mhd, relativity) 
# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(mhd, relativity) 
# 1345 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1346 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1346 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(muscl_order, muscl_lim) 
# 1346 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1346 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(muscl_order, muscl_lim) 
# 1346 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1347 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1347 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(igr, igr_order) 
# 1347 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1347 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(igr, igr_order) 
# 1347 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
# 1348 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1348 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(num_fluids, num_dims, viscous, num_vels, nb, muscl_lim) 
# 1348 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1348 "/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) 
# 1348 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1350 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
 
# 1351 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1351 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(dir_idx, dir_flg, dir_idx_tau) 
# 1351 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1351 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(dir_idx, dir_flg, dir_idx_tau) 
# 1351 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
 
# 1353 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1353 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc update device(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 1353 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1353 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target update to(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 1353 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        ! Allocating grid variables for the x-, y- and z-directions
#ifdef MFC_DEBUG
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1356: ', '@:ALLOCATE(x_cb(-1 - buff_size:m + buff_size))'
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (x_cb(-1 - buff_size:m + buff_size))
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(x_cb) 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:x_cb) 
# 1356 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1357: ', '@:ALLOCATE(x_cc(-buff_size:m + buff_size))'
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (x_cc(-buff_size:m + buff_size))
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(x_cc) 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:x_cc) 
# 1357 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1358: ', '@:ALLOCATE(dx(-buff_size:m + buff_size))'
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dx(-buff_size:m + buff_size))
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dx) 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dx) 
# 1358 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1359 "/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
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving x_cb to GPU => ", SHAPE(x_cb)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cb), sizeof(x_cb), cudamemadvisesetpreferredlocation, 0)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cb), sizeof(x_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(x_cb), sizeof(x_cb), 0, 0)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1359 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1360 "/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
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving x_cc to GPU => ", SHAPE(x_cc)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cc), sizeof(x_cc), cudamemadvisesetpreferredlocation, 0)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(x_cc), sizeof(x_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(x_cc), sizeof(x_cc), 0, 0)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1360 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1361 "/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
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving dx to GPU => ", SHAPE(dx)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dx), sizeof(dx), cudamemadvisesetpreferredlocation, 0)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dx), sizeof(dx), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dx), sizeof(dx), 0, 0)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1361 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (n == 0) return; 
#ifdef MFC_DEBUG
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1364: ', '@:ALLOCATE(y_cb(-1 - buff_size:n + buff_size))'
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (y_cb(-1 - buff_size:n + buff_size))
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(y_cb) 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:y_cb) 
# 1364 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1365: ', '@:ALLOCATE(y_cc(-buff_size:n + buff_size))'
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (y_cc(-buff_size:n + buff_size))
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(y_cc) 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:y_cc) 
# 1365 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1366: ', '@:ALLOCATE(dy(-buff_size:n + buff_size))'
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dy(-buff_size:n + buff_size))
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dy) 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dy) 
# 1366 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1367 "/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
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving y_cb to GPU => ", SHAPE(y_cb)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cb), sizeof(y_cb), cudamemadvisesetpreferredlocation, 0)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cb), sizeof(y_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(y_cb), sizeof(y_cb), 0, 0)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1367 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1368 "/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
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving y_cc to GPU => ", SHAPE(y_cc)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cc), sizeof(y_cc), cudamemadvisesetpreferredlocation, 0)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(y_cc), sizeof(y_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(y_cc), sizeof(y_cc), 0, 0)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1368 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1369 "/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
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving dy to GPU => ", SHAPE(dy)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dy), sizeof(dy), cudamemadvisesetpreferredlocation, 0)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dy), sizeof(dy), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dy), sizeof(dy), 0, 0)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1369 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
 
        if (p == 0) return; 
#ifdef MFC_DEBUG
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1372: ', '@:ALLOCATE(z_cb(-1 - buff_size:p + buff_size))'
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (z_cb(-1 - buff_size:p + buff_size))
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(z_cb) 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:z_cb) 
# 1372 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1373: ', '@:ALLOCATE(z_cc(-buff_size:p + buff_size))'
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (z_cc(-buff_size:p + buff_size))
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(z_cc) 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:z_cc) 
# 1373 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_DEBUG
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1374: ', '@:ALLOCATE(dz(-buff_size:p + buff_size))'
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    allocate (dz(-buff_size:p + buff_size))
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc enter data create(dz) 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target enter data map(always,alloc:dz) 
# 1374 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1375 "/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
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving z_cb to GPU => ", SHAPE(z_cb)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cb), sizeof(z_cb), cudamemadvisesetpreferredlocation, 0)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cb), sizeof(z_cb), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(z_cb), sizeof(z_cb), 0, 0)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1375 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1376 "/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
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving z_cc to GPU => ", SHAPE(z_cc)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cc), sizeof(z_cc), cudamemadvisesetpreferredlocation, 0)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(z_cc), sizeof(z_cc), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(z_cc), sizeof(z_cc), 0, 0)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1376 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
#ifdef MFC_SIMULATION
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#ifdef __NVCOMPILER_GPU_UNIFIED_MEM
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! Beginning in the 25.3 release, the structure of the cudafor module has been changed slightly.
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The module now includes, or “uses” 3 submodules: cuda_runtime_api, gpu_reductions, and sort.
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! The cudafor functionality has not changed. But for new users, or users who have needed to
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! work-around name conflicts in the module, it may be better to use cuda_runtime_api to expose
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
! interfaces to the CUDA runtime calls described in Chapter 4 of this guide.
# 1377 "/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
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if __NVCOMPILER_MAJOR__ < 25 || (__NVCOMPILER_MAJOR__ == 25 && __NVCOMPILER_MINOR__ < 3)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cudafor, gpu_sum => sum, gpu_maxval => maxval, gpu_minval => minval
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#else
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use cuda_runtime_api
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        integer :: istat
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        if (nv_uvm_pref_gpu) then
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                !print*, "Moving dz to GPU => ", SHAPE(dz)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set preferred location GPU
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dz), sizeof(dz), cudamemadvisesetpreferredlocation, 0)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! set accessed by CPU
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemadvise(c_devloc(dz), sizeof(dz), cudamemadvisesetaccessedby, cudacpudeviceid)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                ! prefetch to GPU - physically populate memory pages
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                istat = cudamemprefetchasync(c_devloc(dz), sizeof(dz), 0, 0)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                if (istat /= cudasuccess) then
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    write (*, "('Error code: ',I0, ': ')") istat
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                    !write(*,*) cudaGetErrorString(istat)
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
                end if
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        end if
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1377 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1377 "/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
 
#ifdef MFC_MPI
        integer :: ierr !< Generic flag used to identify and report MPI errors
#endif
 
# 1389 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
            num_dims = 1 + min(1, n) + min(1, p)
 
            if (mhd) then
                num_vels = 3
            else
                num_vels = num_dims
            end if
# 1397 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
        allocate (proc_coords(1:num_dims))
 
        if (parallel_io .neqv. .true.) return
 
#ifdef MFC_MPI
 
        ! Option for Lustre file system (Darter/Comet/Stampede)
        write (mpiiofs, '(A)') '/lustre_'
        mpiiofs = trim(mpiiofs)
 
        call mpi_info_create(mpi_info_int, ierr)
        call mpi_info_set(mpi_info_int, 'romio_ds_write', 'disable', ierr)
 
        ! Option for UNIX file system (Hooke/Thomson)
        ! WRITE(mpiiofs, '(A)') '/ufs_'
        ! mpiiofs = TRIM(mpiiofs)
        ! mpi_info_int = MPI_INFO_NULL
 
        allocate (start_idx(1:num_dims))
 
#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
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1431: ', '@:DEALLOCATE(Re_idx)'
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(Re_idx) 
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:Re_idx) 
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1431 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (re_idx)
        end if
 
        deallocate (proc_coords)
        if (parallel_io) then
            deallocate (start_idx)
 
            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
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1455: ', '@:DEALLOCATE(x_cb, x_cc, dx)'
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(x_cb, x_cc, dx) 
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:x_cb, x_cc, dx) 
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1455 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (x_cb, x_cc, dx)
 
        if (n == 0) return; 
#ifdef MFC_DEBUG
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1458: ', '@:DEALLOCATE(y_cb, y_cc, dy)'
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(y_cb, y_cc, dy) 
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:y_cb, y_cc, dy) 
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1458 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (y_cb, y_cc, dy)
 
        if (p == 0) return; 
#ifdef MFC_DEBUG
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    block
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        use iso_fortran_env, only: output_unit
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        print *, 'm_global_parameters.fpp:1461: ', '@:DEALLOCATE(z_cb, z_cc, dz)'
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
        call flush (output_unit)
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    end block
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
 
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#if defined(MFC_OpenACC)
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$acc exit data delete(z_cb, z_cc, dz) 
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#elif defined(MFC_OpenMP)
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
!$omp target exit data map(release:z_cb, z_cc, dz) 
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
#endif
# 1461 "/home/runner/work/MFC/MFC/src/simulation/m_global_parameters.fpp"
    deallocate (z_cb, z_cc, dz)
 
    end subroutine s_finalize_global_parameters_module
 
end module m_global_parameters