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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
! New line at end of file is required for FYPP
# 2 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp" 2
 
# 4 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 6 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 7 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 8 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 20 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 43 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 48 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 58 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 63 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 68 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 81 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 95 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 100 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 110 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 115 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 120 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 125 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 130 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 135 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 165 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 219 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 233 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 258 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 269 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 271 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
# 282 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 332 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 342 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 352 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 361 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 378 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 388 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 395 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 401 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 407 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 413 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 419 "/home/runner/work/MFC/MFC/src/common/include/omp_macros.fpp"
 
# 425 "/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"
 
# 12 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
# 1 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp" 1
! Default for Fypp runs outside the toolchain (e.g. documentation). Searched last, so the
! build/include/<target>/thermochem.fpp the toolchain writes from m_thermochem takes precedence.
# 4 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 5 "/home/runner/work/MFC/MFC/src/common/include/defaults/thermochem.fpp"
# 13 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp" 2
 
# 16 "/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"
 
# 24 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 53 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 65 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 75 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 105 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 117 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 127 "/home/runner/work/MFC/MFC/src/common/include/shared_parallel_macros.fpp"
 
# 174 "/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.
# 52 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Allocate and create GPU device memory
# 72 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Free GPU device memory and deallocate
# 80 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for vector fields
# 104 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for scalar fields
# 120 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
! Cray-specific GPU pointer setup for acoustic source spatials
# 145 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 151 "/home/runner/work/MFC/MFC/src/common/include/macros.fpp"
 
# 158 "/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: num_fluids_max, model_eqns_gamma_law, model_eqns_5eq, model_eqns_6eq, recon_type_weno, &
        & recon_type_muscl, name_len, dflt_int, dflt_real, eos_stiffened_gas, eos_ideal_gas
 
    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                 :: alt_soundspeed
integer                 :: avg_state
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
character(LEN=path_len) :: file_extension
logical                 :: file_per_process
character(LEN=path_len) :: files_dir
real(wp), dimension(num_fluids_max) :: fluid_rho
logical                 :: hyper_cleaning
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
logical                 :: mixture_err
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_particle_clouds
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
integer                 :: precision
real(wp)                :: ptgalpha_eps
logical                 :: qbmm
logical                 :: reactive_burn
integer                 :: recon_type
logical                 :: relativity
logical                 :: relax
integer                 :: relax_model
real(wp)                :: rhorv
integer                 :: riemann_solver
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(reactive_burn_parameters)      :: rburn !< Condensed-phase reactive-burn (programmed detonation) parameters
type(particle_cloud_parameters), dimension(num_particle_clouds_max) :: particle_cloud !< Particle bed specifications
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
!
integer                 :: num_dims
integer                 :: num_vels
integer                 :: weno_polyn
integer                 :: muscl_polyn
logical                 :: any_state_dependent_eos
# 29 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp" 2
 
    !> @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
    !> @}
 
    !> @name Chemistry modeling (Fypp compile-time constant; same value in all targets)
    !> @{
    logical, parameter :: chemistry = .false.
# 40 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
    !> @}
 
    !> @name Hypoelastic shear stress state (identical across all three executables)
    !> @{
    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)
    !> @}
 
    !> @name Material properties derived from fluid_pp
    !> @{ One declaration is shared by all executables and initialized by m_eos after the case parameters have been read.
    !> gammas is the stored form 1/(gamma - 1), not the ratio of specific heats; isentrope_n and isentrope_B are the same EOS
    !! written as p + B = const*rho**n.
    real(wp), allocatable, dimension(:) :: gammas, isentrope_n, pi_infs, isentrope_b, cvs, qvs, qvps
 
# 55 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 55 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(gammas, isentrope_n, pi_infs, isentrope_B, cvs, qvs, qvps) 
# 55 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 55 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (gammas, isentrope_n, pi_infs, isentrope_B, cvs, qvs, qvps) 
# 55 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
# 56 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
    !> Per-fluid EOS selector and Mie-Gruneisen reference curve, resolved once at init like the arrays above
    integer, allocatable, dimension(:) :: eoss
    !> Per-fluid EOS coefficients, whatever the family; see type eos_coefficients.
    type(eos_coefficients), dimension(num_fluids_max) :: eos_coeffs
    !> any_state_dependent_eos is declared with the case-optimization block above: a parameter when the case is baked in, so the
    !! compiler drops the whole state-dependent chain from kernels that never need it.
 
# 62 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 62 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(eoss, eos_coeffs) 
# 62 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 62 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (eoss, eos_coeffs) 
# 62 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
# 63 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
    !> Fourier heat conduction: true when any fluid sets k_therm > 0. Derived, never read from the namelist.
    logical                             :: heat_conduction
    real(wp), allocatable, dimension(:) :: fluid_k_therm
 
# 66 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#if defined(MFC_OpenACC)
# 66 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$acc declare create(heat_conduction, fluid_k_therm) 
# 66 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#elif defined(MFC_OpenMP)
# 66 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
!$omp declare target (heat_conduction, fluid_k_therm) 
# 66 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
#endif
# 67 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
# 68 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
 
# 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(any_state_dependent_eos) 
# 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 (any_state_dependent_eos) 
# 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"
# 70 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
    !> @}
 
    !> @name Fluids participating in shear and bulk viscosity
    !> @{
    integer, dimension(2)                :: re_size = 0
    integer                              :: re_size_max = 0
    integer, allocatable, dimension(:,:) :: re_idx
 
# 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(Re_size, Re_size_max, Re_idx) 
# 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 (Re_size, Re_size_max, Re_idx) 
# 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"
    !> @}
 
    !> Minimum cell widths. These are distinct from the per-cell width arrays named dx, dy, and dz in simulation and post-process.
    real(wp) :: dx_min, dy_min, dz_min
 
 
# 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(sys_size, eqn_idx) 
# 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 (sys_size, eqn_idx) 
# 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"
 
# 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(shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 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 (shear_num, shear_indices, shear_BC_flip_num, shear_BC_flip_indices) 
# 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"
 
    !> @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
    integer                            :: num_procs_x = 1  !< Number of MPI ranks in x-direction
    integer                            :: num_procs_y = 1  !< Number of MPI ranks in y-direction
    integer                            :: num_procs_z = 1  !< Number of MPI ranks in z-direction
    !> @}
 
    !> @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 and sys_size) from the namelist parameters. This is the shared skeleton: it covers
    !! the model_eqns dispatch, all eqn_idx field assignments, and the hypoelastic/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, six_eqn_alf_is_advected)
 
        integer, intent(in) :: nmom_in
        integer, intent(in) :: nb_in
        logical, intent(in) :: six_eqn_alf_is_advected
 
        ! 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
            if (six_eqn_alf_is_advected) eqn_idx%alf = eqn_idx%adv%end
            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
        end if
 
        if (model_eqns == model_eqns_5eq .or. model_eqns == model_eqns_6eq) then
            if (hypoelasticity) then
                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. Readers test the whole array
                ! rather than the first shear_num entries, so unused slots must not be garbage.
                shear_indices = 0
                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 (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
 
        ! Resolved here, not with the other fluid properties, because the MPI halo buffers are sized before
        ! m_variables_conversion runs and conduction adds temperature to the exchange.
        heat_conduction = any(fluid_pp(:)%k_therm > 0._wp)
 
    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
 
        ! Under case-optimization, num_dims and num_vels are compile-time parameters; skip assignment.
# 279 "/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
# 287 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.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_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.
        cont_damage = .false.
        hyper_cleaning = .false.
 
        ! Condensed-phase reactive burn
        reactive_burn = .false.
        rburn%k = dflt_real
        rburn%pign = dflt_real
        rburn%pref = dflt_real
        rburn%n = dflt_real
        rburn%ta = 0._wp
        rburn%substeps = 0
 
        ! 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.
# 364 "/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.
            igr_order = dflt_int
            mhd = .false.
            relativity = .false.
            viscous = .false.
# 374 "/home/runner/work/MFC/MFC/src/common/m_global_parameters_common.fpp"
 
        ! Not a case-optimization parameter, so it stays a runtime variable in every build.
        riemann_solver = dflt_int
 
        ! Tait EOS
 
        ! 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.
 
        ! Mixture conversion and sound-speed behavior
        avg_state = dflt_int
        alt_soundspeed = .false.
        mixture_err = .false.
        sigr = dflt_real
        dx_min = dflt_real
        dy_min = dflt_real
        dz_min = dflt_real
 
    end subroutine s_assign_common_defaults
 
end module m_global_parameters_common