m_global_parameters_common.fpp.f90

Go to the documentation of this file.

# 1 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!>
!! @file
!! @brief Contains module m_global_parameters_common
 
# 1 "/home/runner/work/MFC/MFC/src/common/include/case.fpp" 1
! This file exists so that Fypp can be run without generating case.fpp files for
! each target. This is useful when generating documentation, for example. This
! should also let MFC be built with CMake directly, without invoking mfc.sh.
 
! For pre-process.
# 8 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
 
! For moving immersed boundaries in simulation
# 12 "/home/runner/work/MFC/MFC/src/common/include/case.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.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"
 
# 206 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 231 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 242 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 255 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 284 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 313 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 330 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 340 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 347 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 353 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 359 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 365 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 371 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 377 "/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"
 
# 193 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 215 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 244 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 259 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 278 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 294 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 304 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
 
# 311 "/home/runner/work/MFC/MFC/src/common/include/acc_macros.fpp"
! New line at end of file is required for FYPP
# 4 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp" 2
 
! GPU parallel region (scalar reductions, maxval/minval)
# 23 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! GPU parallel loop over threads (most common GPU macro)
# 43 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Required closing for GPU_PARALLEL_LOOP
# 55 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Mark routine for device compilation
# 112 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Declare device-resident data
# 130 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Inner loop within a GPU parallel region
# 145 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Scoped GPU data region
# 164 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Host code with device pointers (for MPI with GPU buffers)
# 193 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Allocate device memory (unscoped)
# 207 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Free device memory
# 219 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Atomic operation on device
# 231 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! End atomic capture block
# 242 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Copy data between host and device
# 254 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Synchronization barrier
# 266 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Import GPU library module (openacc or omp_lib)
# 275 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for AMD compiler
# 282 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-Cray compilers
# 289 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code only for Cray compiler
# 296 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
! Emit code for non-NVIDIA compilers
# 303 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
 
# 305 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
# 306 "/home/runner/work/MFC/MFC/src/common/include/parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp" 2
 
# 14 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Caution: This macro requires the use of a binding script to set CUDA_VISIBLE_DEVICES, such that we have one GPU device per MPI
! rank. That's because for both cudaMemAdvise (preferred location) and cudaMemPrefetchAsync we use location = device_id = 0. For an
! example see misc/nvidia_uvm/bind.sh. NVIDIA unified memory page placement hint
# 57 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 77 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 85 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 109 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 125 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 150 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 156 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 163 "/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/common/m_global_parameters_common.fpp" 2
 
!> @brief Shared global parameters and equation-index setup for all three executables. Each per-target m_global_parameters uses this
!! module (default-public, so all symbols re-export to downstream consumers via two-hop use-association).
module m_global_parameters_common
 
#ifdef MFC_MPI
    use mpi
#endif
 
    use m_derived_types
    use m_thermochem, only: num_species
    use m_constants, only: model_eqns_gamma_law, model_eqns_5eq, model_eqns_6eq, model_eqns_4eq, recon_type_weno, &
        & recon_type_muscl, name_len, dflt_int, dflt_real
 
    implicit none
 
    ! All namelist-bound scalar and array declarations (per-target, regenerated at build time by the ninja custom command)
# 1 "/home/runner/work/MFC/MFC/build/include/pre_process/generated_decls.fpp" 1
! AUTO-GENERATED - do not edit directly. Regenerate: cmake reconfigure
!
real(wp)                :: bx0
real(wp)                :: ca
real(wp)                :: r0ref
real(wp)                :: re_inv
real(wp)                :: web
real(wp)                :: a_x
real(wp)                :: a_y
real(wp)                :: a_z
logical                 :: adv_n
logical                 :: bubbles_euler
logical                 :: bubbles_lagrange
character(LEN=path_len) :: case_dir
logical                 :: cfl_adap_dt
logical                 :: cfl_const_dt
logical                 :: cont_damage
logical                 :: cyl_coord
integer                 :: dist_type
logical                 :: down_sample
logical                 :: elliptic_smoothing
integer                 :: elliptic_smoothing_iters
logical                 :: fft_wrt
logical                 :: file_per_process
real(wp), dimension(num_fluids_max) :: fluid_rho
logical                 :: hyper_cleaning
logical                 :: hyperelasticity
logical                 :: hypoelasticity
logical                 :: ib
logical                 :: igr
integer                 :: igr_order
integer                 :: loops_x
integer                 :: loops_y
integer                 :: loops_z
integer                 :: m
logical                 :: mhd
logical                 :: mixlayer_perturb
real(wp)                :: mixlayer_perturb_k0
integer                 :: mixlayer_perturb_nk
real(wp)                :: mixlayer_vel_coef
logical                 :: mixlayer_vel_profile
integer                 :: model_eqns
logical                 :: mpp_lim
integer                 :: muscl_order
integer                 :: n
integer                 :: n_start
integer                 :: n_start_old
integer                 :: nb
integer                 :: num_bc_patches
integer                 :: num_fluids
integer                 :: num_ibs
integer                 :: num_patches
integer                 :: num_stl_models
logical                 :: old_grid
logical                 :: old_ic
integer                 :: p
real(wp)                :: palpha_eps
logical                 :: parallel_io
logical                 :: perturb_flow
integer                 :: perturb_flow_fluid
real(wp)                :: perturb_flow_mag
logical                 :: perturb_sph
integer                 :: perturb_sph_fluid
real(wp)                :: pi_fac
real(wp)                :: poly_sigma
logical                 :: polydisperse
logical                 :: polytropic
logical                 :: pre_stress
integer                 :: precision
real(wp)                :: pref
real(wp)                :: ptgalpha_eps
logical                 :: qbmm
integer                 :: recon_type
logical                 :: relativity
logical                 :: relax
integer                 :: relax_model
real(wp)                :: rhorv
real(wp)                :: rhoref
real(wp)                :: sigr
real(wp)                :: sigv
real(wp)                :: sigma
logical                 :: simplex_perturb
logical                 :: stretch_x
logical                 :: stretch_y
logical                 :: stretch_z
logical                 :: surface_tension
integer                 :: t_step_old
integer                 :: t_step_start
integer                 :: thermal
logical                 :: viscous
integer                 :: weno_order
real(wp)                :: x_a
real(wp)                :: x_b
real(wp)                :: y_a
real(wp)                :: y_b
real(wp)                :: z_a
real(wp)                :: z_b
type(physical_parameters), dimension(num_fluids_max) :: fluid_pp !< Per-fluid stiffened-gas EOS parameters, Reynolds numbers, and shear modulus
type(subgrid_bubble_physical_parameters) :: bub_pp !< Subgrid bubble physical parameters
type(ic_patch_parameters), dimension(num_patches_max) :: patch_icpp !< IC patch parameters
type(bc_patch_parameters), dimension(num_bc_patches_max) :: patch_bc !< Boundary condition patch parameters
type(ib_patch_parameters), dimension(num_ib_patches_max_namelist) :: patch_ib !< Immersed boundary patch parameters
type(ib_airfoil_parameters), dimension(num_ib_airfoils_max) :: ib_airfoil !< Per-airfoil NACA user inputs
type(ib_stl_parameters), dimension(num_stl_models_max) :: stl_models !< Per-STL model parameters
type(simplex_noise_params)          :: simplex_params
# 25 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp" 2
 
    ! Case-optimization declarations: parameters (under MFC_CASE_OPTIMIZATION) or plain variables for
    ! num_dims, num_vels, weno_polyn, muscl_polyn, weno_num_stencils, wenojs, igr, etc.
# 1 "/home/runner/work/MFC/MFC/build/include/pre_process/generated_case_opt_decls.fpp" 1
! AUTO-GENERATED - do not edit directly. Regenerate: cmake reconfigure
!
! (no case-optimization declarations for this target)
# 29 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp" 2
 
    ! For pre_process and post_process: num_dims and num_vels are declared manually here
    ! (sim gets them from generated_case_opt_decls.fpp above)
#ifndef MFC_SIMULATION
    integer :: num_dims  !< Number of spatial dimensions
    integer :: num_vels  !< Number of velocity components (different from num_dims for mhd)
#endif
 
    ! For pre_process: weno_polyn and muscl_polyn are declared manually here
    ! (sim gets them from generated_case_opt_decls.fpp; post does not use them)
#ifdef MFC_PRE_PROCESS
    integer :: weno_polyn   !< Degree of the WENO polynomials
    integer :: muscl_polyn  !< Degree of the MUSCL polynomials
#endif
 
    !> @name Annotations of the structure of the state and flux vectors in terms of the size and configuration of the system of
    !! equations
    !> @{
    integer            :: sys_size     !< Number of unknowns in system of equations
    type(eqn_idx_info) :: eqn_idx      !< All conserved-variable equation index ranges and scalars
    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
    !> @}
 
    !> @name Chemistry modeling (Fypp compile-time constant; same value in all targets)
    !> @{
    logical, parameter :: chemistry = .false.
    !> @}
 
    !> @name Elasticity and shear stress state (identical across all three executables)
    !> @{
    logical                  :: elasticity             !< elasticity modeling, true for hyper or hypo
    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  !< Shear stress BC reflection indices (1:3, 1:shear_BC_flip_num)
    !> @}
 
#ifdef MFC_SIMULATION
 
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(sys_size, eqn_idx, b_size, tensor_size) 
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (sys_size, eqn_idx, b_size, tensor_size) 
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 69 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 69 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 69 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 69 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 69 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
    ! Device residency for namelist/case-opt state declared above via the generated
    ! includes: declare directives must live in the declaring module (Cray ftn rejects
    ! declare-target on use-associated names), so these moved here from simulation.
 
# 73 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 73 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(cyl_coord) 
# 73 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 73 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (cyl_coord) 
# 73 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 74 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 74 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(dt, m, n, p) 
# 74 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 74 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (dt, m, n, p) 
# 74 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 75 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 75 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(cfl_target) 
# 75 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 75 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (cfl_target) 
# 75 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 76 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 76 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(int_comp, ic_eps, ic_beta) 
# 76 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 76 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (int_comp, ic_eps, ic_beta) 
# 76 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 77 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 77 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(muscl_eps) 
# 77 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 77 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (muscl_eps) 
# 77 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 78 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 78 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(mpp_lim, model_eqns, mixture_err, alt_soundspeed) 
# 78 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 78 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (mpp_lim, model_eqns, mixture_err, alt_soundspeed) 
# 78 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 79 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 79 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(avg_state, mp_weno, weno_eps, teno_CT, hypoelasticity) 
# 79 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 79 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (avg_state, mp_weno, weno_eps, teno_CT, hypoelasticity) 
# 79 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 80 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 80 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(hyperelasticity, elasticity, low_Mach) 
# 80 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 80 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (hyperelasticity, elasticity, low_Mach) 
# 80 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 81 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 81 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(cont_damage, hyper_cleaning) 
# 81 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 81 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (cont_damage, hyper_cleaning) 
# 81 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 82 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 82 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(relax, relax_model, palpha_eps, ptgalpha_eps) 
# 82 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 82 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (relax, relax_model, palpha_eps, ptgalpha_eps) 
# 82 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 83 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 83 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(down_sample) 
# 83 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 83 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (down_sample) 
# 83 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 84 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 84 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(fd_order) 
# 84 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 84 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (fd_order) 
# 84 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 85 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 85 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(rhoref, pref) 
# 85 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 85 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (rhoref, pref) 
# 85 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 86 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 86 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(ib, num_ibs) 
# 86 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 86 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (ib, num_ibs) 
# 86 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 87 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 87 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(ib_coefficient_of_friction) 
# 87 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 87 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (ib_coefficient_of_friction) 
# 87 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 88 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 88 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(Ca, Web, Re_inv) 
# 88 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 88 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (Ca, Web, Re_inv) 
# 88 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 89 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 89 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(bubbles_euler, polytropic, polydisperse) 
# 89 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 89 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (bubbles_euler, polytropic, polydisperse) 
# 89 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 90 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 90 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(adv_n, adap_dt, adap_dt_tol, adap_dt_max_iters) 
# 90 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 90 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (adv_n, adap_dt, adap_dt_tol, adap_dt_max_iters) 
# 90 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 91 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 91 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(bubble_model, thermal) 
# 91 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 91 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (bubble_model, thermal) 
# 91 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 92 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 92 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(poly_sigma) 
# 92 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 92 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (poly_sigma) 
# 92 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 93 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 93 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(qbmm, pi_fac) 
# 93 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 93 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (qbmm, pi_fac) 
# 93 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 94 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 94 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(R0ref) 
# 94 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 94 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (R0ref) 
# 94 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 95 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 95 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(acoustic_source, num_source) 
# 95 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 95 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (acoustic_source, num_source) 
# 95 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 96 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 96 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(sigma, surface_tension) 
# 96 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 96 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (sigma, surface_tension) 
# 96 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 97 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 97 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(bubbles_lagrange) 
# 97 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 97 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (bubbles_lagrange) 
# 97 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 98 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 98 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(Bx0) 
# 98 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 98 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (Bx0) 
# 98 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 99 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 99 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(tau_star, cont_damage_s, alpha_bar) 
# 99 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 99 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (tau_star, cont_damage_s, alpha_bar) 
# 99 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 100 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 100 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(hyper_cleaning_speed, hyper_cleaning_tau) 
# 100 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 100 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (hyper_cleaning_speed, hyper_cleaning_tau) 
# 100 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
 
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(num_dims, num_vels, weno_polyn, weno_order) 
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (num_dims, num_vels, weno_polyn, weno_order) 
# 102 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 103 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 103 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(weno_num_stencils, num_fluids, wenojs) 
# 103 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 103 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (weno_num_stencils, num_fluids, wenojs) 
# 103 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 104 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 104 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(mapped_weno, wenoz, teno, wenoz_q, mhd, relativity) 
# 104 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 104 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (mapped_weno, wenoz, teno, wenoz_q, mhd, relativity) 
# 104 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 105 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 105 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(igr_iter_solver, igr_order, viscous, igr_pres_lim, igr) 
# 105 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 105 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (igr_iter_solver, igr_order, viscous, igr_pres_lim, igr) 
# 105 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
# 106 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 106 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(recon_type, muscl_order, muscl_polyn, muscl_lim) 
# 106 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 106 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (recon_type, muscl_order, muscl_polyn, muscl_lim) 
# 106 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
# 108 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
 
    !> @name Processor coordinates and parallel-IO addressing (identical declaration across all three targets)
    !> @{
    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
    !> @}
 
    !> @name MPI info for parallel IO with Lustre file systems (identical across all three targets)
    !> @{
    character(len=name_len) :: mpiiofs
    integer                 :: mpi_info_int
    !> @}
 
contains
 
    !> Initialize equation-index state (eqn_idx, sys_size, b_size, tensor_size) from the namelist parameters. This is the shared
    !! skeleton: it covers the model_eqns dispatch, all eqn_idx field assignments, and the elasticity/surface-tension/chemistry
    !! extensions.
    !!
    !! @param nmom_in  Number of carried moments per R0 location (per-target: pre/post pass an
    !!   integer variable; sim passes its integer parameter nmom = 6).  Used only in the 5eq
    !!   qbmm bubble-index calculation (eqn_idx%bub%end = eqn_idx%adv%end + nb_in*nmom_in).
    !!
    !! Per-target callers are responsible for the following after this call:
    !!   - qbmm_idx allocations and fills (diverge between pre vs sim/post)
    !!   - sim-only: gam = bub_pp%gam_g, nmomsp/nmomtot, Re_idx allocation, GPU_UPDATE calls
    !!   - post-only: beta_idx increment (bubbles_lagrange path), offset/grid allocations
    impure subroutine s_initialize_eqn_idx(nmom_in, nb_in)
 
        integer, intent(in) :: nmom_in
        integer, intent(in) :: nb_in
 
        ! Gamma/Pi_inf Model
 
        if (model_eqns == model_eqns_gamma_law) 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
            eqn_idx%cont%beg = 1
            eqn_idx%cont%end = eqn_idx%cont%beg
            eqn_idx%mom%beg = eqn_idx%cont%end + 1
            eqn_idx%mom%end = eqn_idx%cont%end + num_vels
            eqn_idx%E = eqn_idx%mom%end + 1
            eqn_idx%adv%beg = eqn_idx%E + 1
            eqn_idx%adv%end = eqn_idx%adv%beg + 1
            eqn_idx%gamma = eqn_idx%adv%beg
            eqn_idx%pi_inf = eqn_idx%adv%end
            sys_size = eqn_idx%adv%end
 
            ! Volume Fraction Model (5-equation model)
        else if (model_eqns == model_eqns_5eq) 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 volume fraction model
            eqn_idx%cont%beg = 1
            eqn_idx%cont%end = num_fluids
            eqn_idx%mom%beg = eqn_idx%cont%end + 1
            eqn_idx%mom%end = eqn_idx%cont%end + num_vels
            eqn_idx%E = eqn_idx%mom%end + 1
 
            if (igr) then
                ! IGR: volume fractions after energy (N-1 for N fluids; skipped when num_fluids=1)
                eqn_idx%adv%beg = eqn_idx%E + 1
                eqn_idx%adv%end = eqn_idx%E + num_fluids - 1
            else
                ! WENO/MUSCL + Riemann tracks a total of (N) volume fractions for N fluids
                eqn_idx%adv%beg = eqn_idx%E + 1
                eqn_idx%adv%end = eqn_idx%E + num_fluids
            end if
 
            sys_size = eqn_idx%adv%end
 
            if (bubbles_euler) then
                eqn_idx%alf = eqn_idx%adv%end
            else
                eqn_idx%alf = 1
            end if
 
            if (bubbles_euler) then
                eqn_idx%bub%beg = sys_size + 1
                if (qbmm) then
                    eqn_idx%bub%end = eqn_idx%adv%end + nb_in*nmom_in
                else
                    if (.not. polytropic) then
                        eqn_idx%bub%end = sys_size + 4*nb_in
                    else
                        eqn_idx%bub%end = sys_size + 2*nb_in
                    end if
                end if
                sys_size = eqn_idx%bub%end
 
                if (adv_n) then
                    eqn_idx%n = eqn_idx%bub%end + 1
                    sys_size = eqn_idx%n
                end if
            end if
 
            if (mhd) then
                eqn_idx%B%beg = sys_size + 1
                if (n == 0) then
                    eqn_idx%B%end = sys_size + 2  ! 1D: By, Bz
                else
                    eqn_idx%B%end = sys_size + 3  ! 2D/3D: Bx, By, Bz
                end if
                sys_size = eqn_idx%B%end
            end if
 
            ! Volume Fraction Model (6-equation model)
        else if (model_eqns == model_eqns_6eq) 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 volume fraction model
            eqn_idx%cont%beg = 1
            eqn_idx%cont%end = num_fluids
            eqn_idx%mom%beg = eqn_idx%cont%end + 1
            eqn_idx%mom%end = eqn_idx%cont%end + num_vels
            eqn_idx%E = eqn_idx%mom%end + 1
            eqn_idx%adv%beg = eqn_idx%E + 1
            eqn_idx%adv%end = eqn_idx%E + num_fluids
#ifdef MFC_SIMULATION
            eqn_idx%alf = eqn_idx%adv%end
#endif
            eqn_idx%int_en%beg = eqn_idx%adv%end + 1
            eqn_idx%int_en%end = eqn_idx%adv%end + num_fluids
            sys_size = eqn_idx%int_en%end
        else if (model_eqns == model_eqns_4eq) then
            ! 4-equation model with subgrid bubbles
            eqn_idx%cont%beg = 1
            eqn_idx%cont%end = 1
            eqn_idx%mom%beg = eqn_idx%cont%end + 1
            eqn_idx%mom%end = eqn_idx%cont%end + num_vels
            eqn_idx%E = eqn_idx%mom%end + 1
            eqn_idx%adv%beg = eqn_idx%E + 1
            eqn_idx%adv%end = eqn_idx%adv%beg
            eqn_idx%alf = eqn_idx%adv%end
            sys_size = eqn_idx%adv%end
 
            if (bubbles_euler) then
                eqn_idx%bub%beg = sys_size + 1
                eqn_idx%bub%end = sys_size + 2*nb_in
                if (.not. polytropic) then
                    eqn_idx%bub%end = sys_size + 4*nb_in
                end if
                sys_size = eqn_idx%bub%end
            end if
        end if
 
        if (model_eqns == model_eqns_5eq .or. model_eqns == model_eqns_6eq) then
            if (hypoelasticity .or. hyperelasticity) then
                elasticity = .true.
                eqn_idx%stress%beg = sys_size + 1
                eqn_idx%stress%end = sys_size + (num_dims*(num_dims + 1))/2
                if (cyl_coord) eqn_idx%stress%end = eqn_idx%stress%end + 1
                ! number of stresses is 1 in 1D, 3 in 2D, 4 in 2D-Axisym, 6 in 3D
                sys_size = eqn_idx%stress%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) = eqn_idx%stress%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) = eqn_idx%stress%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
            end if
 
            if (hyperelasticity) then
                ! number of entries in the symmetric b tensor plus the jacobian
                b_size = (num_dims*(num_dims + 1))/2 + 1
                tensor_size = num_dims**2 + 1
                eqn_idx%xi%beg = sys_size + 1
                eqn_idx%xi%end = sys_size + num_dims
                ! adding equations for the xi field and the elastic energy
                sys_size = eqn_idx%xi%end + 1
            end if
 
            if (surface_tension) then
                eqn_idx%c = sys_size + 1
                sys_size = eqn_idx%c
            end if
 
            if (cont_damage) then
                eqn_idx%damage = sys_size + 1
                sys_size = eqn_idx%damage
            end if
 
            if (hyper_cleaning) then
                eqn_idx%psi = sys_size + 1
                sys_size = eqn_idx%psi
            end if
        end if
 
        if (chemistry) then
            eqn_idx%species%beg = sys_size + 1
            eqn_idx%species%end = sys_size + num_species
            sys_size = eqn_idx%species%end
        end if
 
    end subroutine s_initialize_eqn_idx
 
    !> Configure MPI parallel I/O settings and allocate processor coordinate arrays. Shared across all three executables;
    !! num_dims/num_vels are computed here for pre/post unconditionally and for sim only when not case-optimized (in which case they
    !! are compile-time parameters). Callers must have already populated n and p (grid dimensions).
    impure subroutine s_initialize_parallel_io_common
 
#ifdef MFC_MPI
        integer :: ierr  !< Generic flag used to identify and report MPI errors
#endif
 
#ifdef MFC_SIMULATION
        ! Under case-optimization, num_dims and num_vels are compile-time parameters; skip assignment.
# 328 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
            num_dims = 1 + min(1, n) + min(1, p)
 
            if (mhd) then
                num_vels = 3
            else
                num_vels = num_dims
            end if
# 336 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#else
        num_dims = 1 + min(1, n) + min(1, p)
 
        if (mhd) then
            num_vels = 3
        else
            num_vels = num_dims
        end if
#endif
 
        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_common
 
    !> Shared finalize core: deallocate proc_coords and start_idx. Per-target finalize routines call this first, then handle their
    !! own extras (qbmm_idx, grid arrays, MPI_IO_DATA null/dealloc - those reference per-target typed variables and stay
    !! per-target).
    impure subroutine s_finalize_global_parameters_common
 
        deallocate (proc_coords)
 
#ifdef MFC_MPI
        if (parallel_io) then
            deallocate (start_idx)
        end if
#endif
 
    end subroutine s_finalize_global_parameters_common
 
    !> Assign default values to the user-input parameters that are shared across all three executables (pre_process, simulation,
    !! post_process). Per-target defaults (bc_io, cfl_dt, precision, nb, chem_params, bub_pp scalar companions, fluid_pp loop, patch
    !! arrays, output flags) remain in the per-target s_assign_default_values_to_user_inputs routines, which call this subroutine
    !! first, then call s_update_cell_bounds(cells_bounds, m, n, p) (cells_bounds is per-target), then apply their own assignments.
    impure subroutine s_assign_common_defaults
 
        ! Logistics
        case_dir = '.'
 
        ! Computational domain parameters (m/n/p set here; caller must call s_update_cell_bounds after)
        m = dflt_int; n = 0; p = 0
 
        cyl_coord = .false.
 
        ! CFL adaptive time-stepping flags
        cfl_adap_dt = .false.
        cfl_const_dt = .false.
 
        ! Time-stepping bookkeeping
        n_start = dflt_int
        t_step_start = dflt_int
 
        ! Simulation algorithm
        model_eqns = dflt_int
        relax = .false.
        relax_model = dflt_int
        hypoelasticity = .false.
        hyperelasticity = .false.
        elasticity = .false.
        b_size = dflt_int
        tensor_size = dflt_int
        cont_damage = .false.
        hyper_cleaning = .false.
 
        ! Case-optimization params: under case-opt these are compile-time constants in sim (skip assignment); in pre/post
        ! MFC_CASE_OPTIMIZATION is always False so the block always executes there.
# 417 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
            recon_type = recon_type_weno
            weno_order = dflt_int
            muscl_order = dflt_int
            num_fluids = dflt_int
            igr = .false.
            mhd = .false.
            relativity = .false.
# 425 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
 
        ! Tait EOS
        rhoref = dflt_real
        pref = dflt_real
 
        ! Bubble modeling flags and parameters
        r0ref = dflt_real
        bubbles_euler = .false.
        polydisperse = .false.
        poly_sigma = dflt_real
        qbmm = .false.
        surface_tension = .false.
        adv_n = .false.
        sigma = dflt_real
        bubbles_lagrange = .false.
 
        ! Immersed boundaries
        ib = .false.
        num_ibs = dflt_int
 
        ! MHD (background field)
        bx0 = dflt_real
 
        ! Output and I/O options
        parallel_io = .false.
        file_per_process = .false.
        down_sample = .false.
        fft_wrt = .false.
 
    end subroutine s_assign_common_defaults
 
end module m_global_parameters_common