added sys_size under eqn_idx

Author: mohdsaid497566

src/common/m_boundary_common.fpp

@@ -25,7 +25,7 @@ contains
         ! conservative variables.
 
         type(scalar_field), intent(inout) :: q_T_sf
-        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_cons_vf
         type(int_bounds_info), dimension(1:3), intent(in) :: bounds
 
         integer :: x, y, z, eqn
@@ -62,7 +62,7 @@ contains
     subroutine s_compute_T_from_primitives(q_T_sf, q_prim_vf, bounds)
 
         type(scalar_field), intent(inout) :: q_T_sf
-        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_prim_vf
         type(int_bounds_info), dimension(1:3), intent(in) :: bounds
 
         integer :: x, y, z, i
@@ -87,9 +87,9 @@ contains
 
     subroutine s_compute_chemistry_reaction_flux(rhs_vf, q_cons_qp, q_T_sf, q_prim_qp, bounds)
 
-        type(scalar_field), dimension(sys_size), intent(inout) :: rhs_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: rhs_vf
         type(scalar_field), intent(inout) :: q_T_sf
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_qp, q_prim_qp
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_qp, q_prim_qp
         type(int_bounds_info), dimension(1:3), intent(in) :: bounds
 
         integer :: x, y, z

src/common/m_chemistry.fpp

@@ -140,6 +140,7 @@ module m_derived_types
     !! size and the configuration of the system of equations to which they belong
     !> @{
     type system_of_equations
+        integer :: sys_size                            !< Size of the system of equations
         type(int_bounds_info) :: cont                  !< Indexes of first & last continuity eqns.
         type(int_bounds_info) :: mom                   !< Indexes of first & last momentum eqns.
         integer :: E                                   !< Index of energy equation

src/common/m_derived_types.fpp

@@ -48,7 +48,7 @@ contains
         !!      some parameter. It replaces the procedure pointer, which CCE
         !!      is breaking on.
     impure subroutine s_relaxation_solver(q_cons_vf)
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_vf
         ! This is empty because in current master the procedure pointer
         ! was never assigned
         @:ASSERT(.false., "s_relaxation_solver called but it currently does nothing")
@@ -80,7 +80,7 @@ contains
         !!  @param q_cons_vf Cell-average conservative variables
     pure subroutine s_infinite_relaxation_k(q_cons_vf)
 
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_vf
         real(wp) :: pS, pSOV, pSSL !< equilibrium pressure for mixture, overheated vapor, and subcooled liquid
         real(wp) :: TS, TSOV, TSSL, TSatOV, TSatSL !< equilibrium temperature for mixture, overheated vapor, and subcooled liquid. Saturation Temperatures at overheated vapor and subcooled liquid
         real(wp) :: rhoe, dynE, rhos !< total internal energy, kinetic energy, and total entropy
@@ -293,7 +293,7 @@ contains
         integer, intent(in) :: j, k, l, MFL
         real(wp), intent(out) :: pS
         real(wp), dimension(num_fluids), intent(out) :: p_infpT
-        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_cons_vf
         real(wp), intent(in) :: rhoe
         real(wp), intent(out) :: TS
         real(wp) :: gp, gpp, hp, pO, mCP, mQ !< variables for the Newton Solver
@@ -397,7 +397,7 @@ contains
         real(wp), intent(inout) :: pS
         real(wp), dimension(num_fluids), intent(in) :: p_infpT
         real(wp), intent(in) :: rhoe
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_vf
         real(wp), intent(inout) :: TS
 
         real(wp), dimension(num_fluids) :: p_infpTg !< stiffness for the participating fluids for pTg-equilibrium
@@ -521,7 +521,7 @@ contains
         !> @name variables for the correction of the reacting partial densities
         !> @{
         real(wp), intent(out) :: MCT
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_vf
         real(wp), intent(inout) :: rM
         integer, intent(in) :: j, k, l
         !> @}
@@ -586,7 +586,7 @@ contains
         real(wp), dimension(2, 2), intent(out) :: Jac
         integer, intent(in) :: k, l
         real(wp), intent(in) :: mCPD, mCVGP, mCVGP2, pS
-        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_cons_vf
         real(wp), dimension(2, 2), intent(out) :: TJac
 
         real(wp) :: ml, mT, TS, dFdT, dTdm, dTdp ! mass of the reacting fluid, total reacting mass, and auxiliary variables
@@ -690,7 +690,7 @@ contains
 
         integer, intent(in) :: j, k, l
         real(wp), intent(in) :: mCPD, mCVGP, mQD
-        type(scalar_field), dimension(sys_size), intent(in) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_cons_vf
         real(wp), intent(in) :: pS, rhoe
         real(wp), dimension(2), intent(out) :: R2D
 

src/common/m_mpi_common.fpp

@@ -82,7 +82,7 @@ contains
     subroutine s_convert_to_mixture_variables(q_vf, i, j, k, &
                                               rho, gamma, pi_inf, qv, Re_K, G_K, G)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: q_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_vf
         integer, intent(in) :: i, j, k
         real(wp), intent(out), target :: rho, gamma, pi_inf, qv
         real(wp), optional, dimension(2), intent(out) :: Re_K
@@ -208,7 +208,7 @@ contains
     subroutine s_convert_mixture_to_mixture_variables(q_vf, i, j, k, &
                                                       rho, gamma, pi_inf, qv)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: q_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_vf
         integer, intent(in) :: i, j, k
 
         real(wp), intent(out), target :: rho
@@ -250,7 +250,7 @@ contains
     subroutine s_convert_species_to_mixture_variables_bubbles(q_vf, j, k, l, &
                                                               rho, gamma, pi_inf, qv, Re_K)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: q_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_vf
 
         integer, intent(in) :: j, k, l
 
@@ -372,7 +372,7 @@ contains
     subroutine s_convert_species_to_mixture_variables(q_vf, k, l, r, rho, &
                                                       gamma, pi_inf, qv, Re_K, G_K, G)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: q_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_vf
 
         integer, intent(in) :: k, l, r
 
@@ -735,7 +735,7 @@ contains
     !Initialize mv at the quadrature nodes based on the initialized moments and sigma
     pure subroutine s_initialize_mv(qK_cons_vf, mv)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: qK_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: qK_cons_vf
 
         real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: mv
 
@@ -767,7 +767,7 @@ contains
 
     !Initialize pb at the quadrature nodes using isothermal relations (Preston model)
     pure subroutine s_initialize_pb(qK_cons_vf, mv, pb)
-        type(scalar_field), dimension(sys_size), intent(in) :: qK_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: qK_cons_vf
 
         real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(in) :: mv
         real(wp), dimension(idwint(1)%beg:, idwint(2)%beg:, idwint(3)%beg:, 1:, 1:), intent(inout) :: pb
@@ -811,9 +811,9 @@ contains
                                                              qK_prim_vf, &
                                                              ibounds)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: qK_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: qK_cons_vf
         type(scalar_field), intent(inout) :: q_T_sf
-        type(scalar_field), dimension(sys_size), intent(inout) :: qK_prim_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: qK_prim_vf
         type(int_bounds_info), dimension(1:3), intent(in) :: ibounds
 
         real(wp), dimension(num_fluids) :: alpha_K, alpha_rho_K
@@ -1156,8 +1156,8 @@ contains
     impure subroutine s_convert_primitive_to_conservative_variables(q_prim_vf, &
                                                                     q_cons_vf)
 
-        type(scalar_field), dimension(sys_size), intent(in) :: q_prim_vf
-        type(scalar_field), dimension(sys_size), intent(inout) :: q_cons_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(in) :: q_prim_vf
+        type(scalar_field), dimension(eqn_idx%sys_size), intent(inout) :: q_cons_vf
 
         ! Density, specific heat ratio function, liquid stiffness function
         ! and dynamic pressure, as defined in the incompressible flow sense,