remove fluid_rho

Author: sbryngelson

src/pre_process/m_global_parameters.fpp

@@ -78,16 +78,13 @@ module m_global_parameters
     type(int_bounds_info) :: mom_idx                    !< Indexes of first & last momentum eqns.
     integer :: E_idx                      !< Index of total energy equation
     type(int_bounds_info) :: adv_idx                    !< Indexes of first & last advection eqns.
-    integer :: gamma_idx                  !< Index of specific heat ratio func. eqn.
-    integer :: pi_inf_idx                 !< Index of liquid stiffness func. eqn.
 
     type(int_bounds_info) :: bc_x, bc_y !<
     !! Boundary conditions in the x-, y- and z-coordinate directions
 
     logical :: parallel_io !< Format of the data files
     integer :: precision !< Precision of output files
 
-    real(kind(0d0)), dimension(num_fluids_max) :: fluid_rho
 
     integer, allocatable, dimension(:) :: proc_coords !<
     !! Processor coordinates in MPI_CART_COMM
@@ -187,7 +184,6 @@ contains
 
         parallel_io = .false.
         precision = 2
-        fluid_rho = dflt_real
 
         ! Initial condition parameters
         num_patches = dflt_int

src/pre_process/m_mpi_proxy.fpp

@@ -148,7 +148,6 @@ contains
         #:for VAR in [ 'old_grid','old_ic','stretch_x','stretch_y', 'parallel_io'  ]
             call MPI_BCAST(${VAR}$, 1, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr)
         #:endfor
-        call MPI_BCAST(fluid_rho(1), num_fluids_max, MPI_LOGICAL, 0, MPI_COMM_WORLD, ierr)
 
 
         #:for VAR in [ 'x_domain%beg', 'x_domain%end', 'y_domain%beg',         &

src/pre_process/m_start_up.fpp

@@ -94,7 +94,6 @@ contains
             weno_order, bc_x, bc_y, num_patches, &
             patch_icpp, fluid_pp, &
             precision, parallel_io, &
-            fluid_rho, &
             loops_x, loops_y
 
         ! Inquiring the status of the pre_process.inp file

src/simulation/m_global_parameters.fpp

@@ -40,8 +40,8 @@ module m_global_parameters
     character(LEN=path_len) :: case_dir              !< Case folder location
     logical :: run_time_info         !< Run-time output flag
     integer :: t_step_old            !< Existing IC/grid folder
-    real(kind(0d0)), parameter :: small_alf = 1d-7 !< Small alf tolerance
     ! ==========================================================================
+
     ! Computational Domain Parameters ==========================================
     integer :: proc_rank !< Rank of the local processor
 
@@ -72,7 +72,7 @@ module m_global_parameters
 
     real(kind(0d0)) :: dt !< Size of the time-step
 
-!$acc declare create(x_cb, y_cb, x_cc, y_cc, dx, dy, dt, m, n)
+    !$acc declare create(x_cb, y_cb, x_cc, y_cc, dx, dy, dt, m, n)
 
     !> @name Starting time-step iteration, stopping time-step iteration and the number
     !! of time-step iterations between successive solution backups, respectively
@@ -112,7 +112,7 @@ module m_global_parameters
 
     !> @name Boundary conditions (BC) in the x-, y- and z-directions, respectively
     !> @{
-    type(int_bounds_info) :: bc_x, bc_y, bc_z
+    type(int_bounds_info) :: bc_x, bc_y
     type(int_bounds_info) :: bc_x_glb, bc_y_glb
     !> @}
 
@@ -258,7 +258,6 @@ contains
 
         bc_x%beg = dflt_int; bc_x%end = dflt_int
         bc_y%beg = dflt_int; bc_y%end = dflt_int
-        bc_z%beg = dflt_int; bc_z%end = dflt_int
 
         ! Fluids physical parameters
         do i = 1, num_fluids_max

toolchain/mfc/run/case_dicts.py

@@ -11,8 +11,7 @@
 
 
 PRE_PROCESS = COMMON + [
-    'old_grid', 'old_ic', 't_step_old',
-    'fluid_rho', 'num_patches'
+    'old_grid', 'old_ic', 't_step_old', 'num_patches'
 ]
 
 for cmp in ["x", "y"]:
@@ -26,8 +25,6 @@
     PRE_PROCESS.append(f"bc_{cmp}%end")
 
 for f_id in range(1, 10+1):
-    PRE_PROCESS.append(f'fluid_rho({f_id})')
-
     for attribute in ["gamma", "pi_inf"]:
         PRE_PROCESS.append(f"fluid_pp({f_id})%{attribute}")