diff --git a/.github/scripts/check_coverage_map_health.py b/.github/scripts/check_coverage_map_health.py index f4dddc9838..f8307b5244 100644 --- a/.github/scripts/check_coverage_map_health.py +++ b/.github/scripts/check_coverage_map_health.py @@ -1,11 +1,10 @@ """Fail loudly if the committed coverage map is stale or under-covers. Used by coverage-health.yml.""" import datetime -import subprocess import sys from pathlib import Path sys.path.insert(0, str(Path(__file__).resolve().parents[2] / "toolchain")) -from mfc.test.coverage import COVERAGE_MAP_PATH, load_map, map_health # noqa: E402 +from mfc.test.coverage import COVERAGE_MAP_PATH, _git, load_map, map_health # noqa: E402 from mfc.test.cases import list_cases # noqa: E402 (returns the current test list) MAX_AGE_DAYS = 10 @@ -40,7 +39,7 @@ def verified_sha(cwd=None): caller must read that as undeterminable and fall back to the wall-clock age rule, not as a failure -- an absent ref is not evidence of a broken refresh. """ - rev = subprocess.run(["git", "rev-parse", "--verify", "--quiet", f"{VERIFIED_REF}^{{commit}}"], capture_output=True, text=True, check=False, cwd=cwd) + rev = _git(["rev-parse", "--verify", "--quiet", f"{VERIFIED_REF}^{{commit}}"], cwd) return rev.stdout.strip() or None @@ -53,10 +52,10 @@ def verified_after_last_change(git_sha, cwd=None): """ if not git_sha: return None - last = subprocess.run(["git", "log", "-1", "--format=%H", "--", *COVERAGE_RELEVANT_PATHS], capture_output=True, text=True, check=False, cwd=cwd) + last = _git(["log", "-1", "--format=%H", "--", *COVERAGE_RELEVANT_PATHS], cwd) if last.returncode != 0 or not last.stdout.strip(): return None # shallow clone or no such commit -> fall back to the age rule - ancestor = subprocess.run(["git", "merge-base", "--is-ancestor", last.stdout.strip(), git_sha], capture_output=True, check=False, cwd=cwd) + ancestor = _git(["merge-base", "--is-ancestor", last.stdout.strip(), git_sha], cwd) return {0: True, 1: False}.get(ancestor.returncode) # anything else -> None (unknown sha, shallow history) diff --git a/examples/1D_advection_convergence/case.py b/examples/1D_advection_convergence/case.py index 8da9db09d0..d661dc8870 100644 --- a/examples/1D_advection_convergence/case.py +++ b/examples/1D_advection_convergence/case.py @@ -95,10 +95,8 @@ "patch_icpp(1)%alpha(2)": "0.5 - 0.2 * sin(2.0 * pi * x / lx)", "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, **scheme_params, } ) diff --git a/examples/1D_brio_wu/case.py b/examples/1D_brio_wu/case.py index f94ed1350b..b8639f924b 100644 --- a/examples/1D_brio_wu/case.py +++ b/examples/1D_brio_wu/case.py @@ -74,7 +74,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (2.0e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_brio_wu_hlld/case.py b/examples/1D_brio_wu_hlld/case.py index 1dc2fbd2db..da26c45d17 100644 --- a/examples/1D_brio_wu_hlld/case.py +++ b/examples/1D_brio_wu_hlld/case.py @@ -75,7 +75,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (2.0e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_brio_wu_rmhd/case.py b/examples/1D_brio_wu_rmhd/case.py index 73020d18da..989ddb842d 100644 --- a/examples/1D_brio_wu_rmhd/case.py +++ b/examples/1D_brio_wu_rmhd/case.py @@ -76,7 +76,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (2.0e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_convergence/case.py b/examples/1D_convergence/case.py index 5be7a4f13c..140c3ed3ea 100755 --- a/examples/1D_convergence/case.py +++ b/examples/1D_convergence/case.py @@ -74,10 +74,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_dai_woodward/case.py b/examples/1D_dai_woodward/case.py index 8b0ae4e6fa..acab1b119d 100644 --- a/examples/1D_dai_woodward/case.py +++ b/examples/1D_dai_woodward/case.py @@ -83,7 +83,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / ((5.0 / 3.0) - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, }, ) ) diff --git a/examples/1D_dai_woodward_hlld/case.py b/examples/1D_dai_woodward_hlld/case.py index 1b51af0f8a..aa068f3bb2 100644 --- a/examples/1D_dai_woodward_hlld/case.py +++ b/examples/1D_dai_woodward_hlld/case.py @@ -84,7 +84,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / ((5.0 / 3.0) - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, }, ) ) diff --git a/examples/1D_euler_convergence/case.py b/examples/1D_euler_convergence/case.py index eacb56661a..7d84d611f0 100644 --- a/examples/1D_euler_convergence/case.py +++ b/examples/1D_euler_convergence/case.py @@ -100,7 +100,6 @@ "patch_icpp(1)%alpha(1)": 1.0, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, **scheme_params, } ) diff --git a/examples/1D_hypo_2materials/case.py b/examples/1D_hypo_2materials/case.py index 86a81365fc..a6920ce363 100755 --- a/examples/1D_hypo_2materials/case.py +++ b/examples/1D_hypo_2materials/case.py @@ -83,7 +83,6 @@ "fluid_pp(1)%G": 1.0e09, "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": 0.0, } ) diff --git a/examples/1D_kapilashocktube/case.py b/examples/1D_kapilashocktube/case.py index aa43b5f36a..00f35288fd 100755 --- a/examples/1D_kapilashocktube/case.py +++ b/examples/1D_kapilashocktube/case.py @@ -69,7 +69,6 @@ "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/1D_laxshocktube/case.py b/examples/1D_laxshocktube/case.py index 4d5ab64138..8612c099f8 100644 --- a/examples/1D_laxshocktube/case.py +++ b/examples/1D_laxshocktube/case.py @@ -68,7 +68,6 @@ # air: an ideal gas, so it carries no stiffness and pi_inf is not read "fluid_pp(1)%eos": "ideal_gas", "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_mhd_smooth_alfven_wave/case.py b/examples/1D_mhd_smooth_alfven_wave/case.py index ab19a23ebe..52bae03830 100644 --- a/examples/1D_mhd_smooth_alfven_wave/case.py +++ b/examples/1D_mhd_smooth_alfven_wave/case.py @@ -62,7 +62,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_analytical/case.py b/examples/1D_shuosher_analytical/case.py index 36c360e8c1..77470f16b8 100644 --- a/examples/1D_shuosher_analytical/case.py +++ b/examples/1D_shuosher_analytical/case.py @@ -69,7 +69,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_old/case.py b/examples/1D_shuosher_old/case.py index b07243a48f..be0d226f06 100644 --- a/examples/1D_shuosher_old/case.py +++ b/examples/1D_shuosher_old/case.py @@ -70,7 +70,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_teno5/case.py b/examples/1D_shuosher_teno5/case.py index 3c2e9fe75a..b2bc388d4b 100644 --- a/examples/1D_shuosher_teno5/case.py +++ b/examples/1D_shuosher_teno5/case.py @@ -72,7 +72,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_teno7/case.py b/examples/1D_shuosher_teno7/case.py index 1eb9b88ca5..3fdc2d03e9 100644 --- a/examples/1D_shuosher_teno7/case.py +++ b/examples/1D_shuosher_teno7/case.py @@ -72,7 +72,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_wenojs5/case.py b/examples/1D_shuosher_wenojs5/case.py index 9ca624bdbf..f833c6592d 100644 --- a/examples/1D_shuosher_wenojs5/case.py +++ b/examples/1D_shuosher_wenojs5/case.py @@ -71,7 +71,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_wenom5/case.py b/examples/1D_shuosher_wenom5/case.py index 1fca5b958c..af8f76c642 100644 --- a/examples/1D_shuosher_wenom5/case.py +++ b/examples/1D_shuosher_wenom5/case.py @@ -71,7 +71,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_shuosher_wenoz5/case.py b/examples/1D_shuosher_wenoz5/case.py index 990ba033e2..323806701a 100644 --- a/examples/1D_shuosher_wenoz5/case.py +++ b/examples/1D_shuosher_wenoz5/case.py @@ -71,7 +71,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_sod_convergence/case.py b/examples/1D_sod_convergence/case.py index ee03382297..43021676a0 100644 --- a/examples/1D_sod_convergence/case.py +++ b/examples/1D_sod_convergence/case.py @@ -99,7 +99,6 @@ "patch_icpp(2)%alpha(1)": 1.0, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, **scheme_params, } ) diff --git a/examples/1D_sodshocktube/case.py b/examples/1D_sodshocktube/case.py index 6c33032b6b..02b85761f3 100755 --- a/examples/1D_sodshocktube/case.py +++ b/examples/1D_sodshocktube/case.py @@ -70,7 +70,6 @@ # air: an ideal gas, so it carries no stiffness and pi_inf is not read "fluid_pp(1)%eos": "ideal_gas", "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_sodshocktube_muscl/case.py b/examples/1D_sodshocktube_muscl/case.py index 30387944c8..e4e6724d72 100755 --- a/examples/1D_sodshocktube_muscl/case.py +++ b/examples/1D_sodshocktube_muscl/case.py @@ -65,7 +65,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_titarevtorro/case.py b/examples/1D_titarevtorro/case.py index 29cbeff6da..2d9bc2f256 100644 --- a/examples/1D_titarevtorro/case.py +++ b/examples/1D_titarevtorro/case.py @@ -68,7 +68,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_titarevtorro_analytical/case.py b/examples/1D_titarevtorro_analytical/case.py index fef323f011..e743e979f4 100644 --- a/examples/1D_titarevtorro_analytical/case.py +++ b/examples/1D_titarevtorro_analytical/case.py @@ -67,7 +67,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/1D_vacuum/case.py b/examples/1D_vacuum/case.py index f9e688d9fd..1cf87170f5 100644 --- a/examples/1D_vacuum/case.py +++ b/examples/1D_vacuum/case.py @@ -69,7 +69,6 @@ "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/1D_vacuum_restart/case.py b/examples/1D_vacuum_restart/case.py index 067b9fbc0d..2a2d4de0cc 100644 --- a/examples/1D_vacuum_restart/case.py +++ b/examples/1D_vacuum_restart/case.py @@ -69,7 +69,6 @@ "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_GreshoVortex/case.py b/examples/2D_GreshoVortex/case.py index 0b9240504f..1d730f6ab4 100644 --- a/examples/2D_GreshoVortex/case.py +++ b/examples/2D_GreshoVortex/case.py @@ -81,7 +81,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_IGR_2fluid/case.py b/examples/2D_IGR_2fluid/case.py index b4967a0567..a0451aeef9 100644 --- a/examples/2D_IGR_2fluid/case.py +++ b/examples/2D_IGR_2fluid/case.py @@ -69,11 +69,9 @@ # Fluid Parameters (Gas) "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(1)%Re(1)": 1e5, "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1e5, # Ambient pressure "patch_icpp(1)%geometry": 3, diff --git a/examples/2D_IGR_triple_point/case.py b/examples/2D_IGR_triple_point/case.py index 6efbb9cf67..ff78d73e99 100755 --- a/examples/2D_IGR_triple_point/case.py +++ b/examples/2D_IGR_triple_point/case.py @@ -96,10 +96,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0 / (1.5 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_TaylorGreenVortex/case.py b/examples/2D_TaylorGreenVortex/case.py index e82cf7089a..231a30028e 100644 --- a/examples/2D_TaylorGreenVortex/case.py +++ b/examples/2D_TaylorGreenVortex/case.py @@ -71,7 +71,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, # Shear viscosity of STD air "fluid_pp(1)%Re(1)": 1 / Mu, } diff --git a/examples/2D_acoustic_broadband/case.py b/examples/2D_acoustic_broadband/case.py index b61aa6e86f..ef7e4db43f 100644 --- a/examples/2D_acoustic_broadband/case.py +++ b/examples/2D_acoustic_broadband/case.py @@ -78,7 +78,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/2D_acoustic_pulse/case.py b/examples/2D_acoustic_pulse/case.py index 86239faf09..f0abf8007c 100644 --- a/examples/2D_acoustic_pulse/case.py +++ b/examples/2D_acoustic_pulse/case.py @@ -98,7 +98,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_acoustic_pulse_analytical/case.py b/examples/2D_acoustic_pulse_analytical/case.py index 2a513a1e07..676c9f2972 100644 --- a/examples/2D_acoustic_pulse_analytical/case.py +++ b/examples/2D_acoustic_pulse_analytical/case.py @@ -97,7 +97,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_advection/case.py b/examples/2D_advection/case.py index 351e4e64e8..3af58db8a7 100644 --- a/examples/2D_advection/case.py +++ b/examples/2D_advection/case.py @@ -81,7 +81,6 @@ "fluid_pp(1)%pi_inf": 2.35e00 * 1.0e09 / (2.35e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_advection_convergence/case.py b/examples/2D_advection_convergence/case.py index 8fffaad603..75ecfcdf09 100644 --- a/examples/2D_advection_convergence/case.py +++ b/examples/2D_advection_convergence/case.py @@ -111,7 +111,6 @@ "patch_icpp(1)%alpha(1)": 1.0, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, **scheme_params, } ) diff --git a/examples/2D_advection_muscl/case.py b/examples/2D_advection_muscl/case.py index 5cedea6273..e53b83cce0 100644 --- a/examples/2D_advection_muscl/case.py +++ b/examples/2D_advection_muscl/case.py @@ -79,7 +79,6 @@ "fluid_pp(1)%pi_inf": 2.35e00 * 1.0e09 / (2.35e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_axisym_shockbubble/case.py b/examples/2D_axisym_shockbubble/case.py index 0ad63dbb10..9d148dd601 100644 --- a/examples/2D_axisym_shockbubble/case.py +++ b/examples/2D_axisym_shockbubble/case.py @@ -112,10 +112,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / (1.6666e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_backward_facing_step/case.py b/examples/2D_backward_facing_step/case.py index f7ac307bf1..574d631b09 100644 --- a/examples/2D_backward_facing_step/case.py +++ b/examples/2D_backward_facing_step/case.py @@ -82,7 +82,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0 / (gam_a - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1 / mu, }, indent=4, diff --git a/examples/2D_bingham_poiseuille_nn/case.py b/examples/2D_bingham_poiseuille_nn/case.py index f457c30c0f..61b243d804 100644 --- a/examples/2D_bingham_poiseuille_nn/case.py +++ b/examples/2D_bingham_poiseuille_nn/case.py @@ -119,7 +119,6 @@ # Fluids Physical Parameters: single Bingham (HB with n = 1, tau0 > 0) fluid "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.0 / K, "fluid_pp(1)%non_newtonian": "T", "fluid_pp(1)%K": K, diff --git a/examples/2D_forward_facing_step/case.py b/examples/2D_forward_facing_step/case.py index df12918d41..015b643f93 100644 --- a/examples/2D_forward_facing_step/case.py +++ b/examples/2D_forward_facing_step/case.py @@ -79,7 +79,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0 / (gam_a - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "viscous": "T", "fluid_pp(1)%Re(1)": 1 / mu, }, diff --git a/examples/2D_hardcoded_ic/case.py b/examples/2D_hardcoded_ic/case.py index 7626817edc..af51250d2d 100644 --- a/examples/2D_hardcoded_ic/case.py +++ b/examples/2D_hardcoded_ic/case.py @@ -74,7 +74,6 @@ "fluid_pp(1)%pi_inf": 2.35e00 * 1.0e09 / (2.35e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_herschel_bulkley_poiseuille_nn/case.py b/examples/2D_herschel_bulkley_poiseuille_nn/case.py index 96ac84ba8a..b0bb7f3368 100644 --- a/examples/2D_herschel_bulkley_poiseuille_nn/case.py +++ b/examples/2D_herschel_bulkley_poiseuille_nn/case.py @@ -126,7 +126,6 @@ # Fluids Physical Parameters: single Herschel-Bulkley (n != 1, tau0 > 0) fluid "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.0 / K, "fluid_pp(1)%non_newtonian": "T", "fluid_pp(1)%K": K, diff --git a/examples/2D_hypo_shear_contact/case.py b/examples/2D_hypo_shear_contact/case.py index 0593b205db..5c70959dd5 100644 --- a/examples/2D_hypo_shear_contact/case.py +++ b/examples/2D_hypo_shear_contact/case.py @@ -83,11 +83,9 @@ "fd_order": 4, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%G": G, "fluid_pp(2)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": G, # Left/background state. "patch_icpp(1)%geometry": 3, diff --git a/examples/2D_ibm/case.py b/examples/2D_ibm/case.py index 70b224d7b1..f2ddf16fdd 100644 --- a/examples/2D_ibm/case.py +++ b/examples/2D_ibm/case.py @@ -91,7 +91,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_ibm_airfoil/case.py b/examples/2D_ibm_airfoil/case.py index 06c16caf45..9ccbd71299 100644 --- a/examples/2D_ibm_airfoil/case.py +++ b/examples/2D_ibm_airfoil/case.py @@ -101,10 +101,8 @@ # Use the same stiffness as the air bubble "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50 (Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(2)%gamma": 1.0e00 / (gam_b - 1.0e00), # 2.50 (Not 1.40) "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0, } ) ) diff --git a/examples/2D_ibm_cfl_dt/case.py b/examples/2D_ibm_cfl_dt/case.py index c4fd269877..3db9e70f07 100644 --- a/examples/2D_ibm_cfl_dt/case.py +++ b/examples/2D_ibm_cfl_dt/case.py @@ -94,7 +94,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 250000, } ) diff --git a/examples/2D_ibm_ellipse/case.py b/examples/2D_ibm_ellipse/case.py index 20e46f04d1..bf8c0e03e9 100644 --- a/examples/2D_ibm_ellipse/case.py +++ b/examples/2D_ibm_ellipse/case.py @@ -92,7 +92,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_ibm_multiphase/case.py b/examples/2D_ibm_multiphase/case.py index 6006606da7..221aff57e0 100644 --- a/examples/2D_ibm_multiphase/case.py +++ b/examples/2D_ibm_multiphase/case.py @@ -91,10 +91,8 @@ # Specify 2 fluids "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(2)%gamma": 1.0e00 / (gam_b - 1.0e00), # 2.50(Not 1.40) "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0, } ) ) diff --git a/examples/2D_ibm_poiseuille_nn/case.py b/examples/2D_ibm_poiseuille_nn/case.py index 5e6f69e586..72fc408349 100644 --- a/examples/2D_ibm_poiseuille_nn/case.py +++ b/examples/2D_ibm_poiseuille_nn/case.py @@ -131,7 +131,6 @@ # Fluids Physical Parameters (EOS shared by all modes) "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } if MODE == "powerlaw": diff --git a/examples/2D_ibm_stl_MFCCharacter/case.py b/examples/2D_ibm_stl_MFCCharacter/case.py index c404fbf820..120cb02729 100644 --- a/examples/2D_ibm_stl_MFCCharacter/case.py +++ b/examples/2D_ibm_stl_MFCCharacter/case.py @@ -80,7 +80,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 10000, } ) diff --git a/examples/2D_ibm_stl_test/case.py b/examples/2D_ibm_stl_test/case.py index 0230d5f36e..0518549d02 100644 --- a/examples/2D_ibm_stl_test/case.py +++ b/examples/2D_ibm_stl_test/case.py @@ -81,7 +81,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/2D_ibm_stl_wedge/case.py b/examples/2D_ibm_stl_wedge/case.py index bd2da045d2..a0e90fb147 100644 --- a/examples/2D_ibm_stl_wedge/case.py +++ b/examples/2D_ibm_stl_wedge/case.py @@ -81,7 +81,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 7535533.2, } ) diff --git a/examples/2D_ibm_viscous_drag_over_cylinder/case.py b/examples/2D_ibm_viscous_drag_over_cylinder/case.py index 91fabc921c..de54bf6999 100644 --- a/examples/2D_ibm_viscous_drag_over_cylinder/case.py +++ b/examples/2D_ibm_viscous_drag_over_cylinder/case.py @@ -78,7 +78,6 @@ # --- Fluid properties (calorically perfect air) --- "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), # MFC: 1/(gamma-1) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": Re_p / (d_cyl * U_inf), # --- Boundary conditions --- "bc_x%beg": -7, diff --git a/examples/2D_icpp_stl_circle/case.py b/examples/2D_icpp_stl_circle/case.py index afca9089d8..a68406be40 100644 --- a/examples/2D_icpp_stl_circle/case.py +++ b/examples/2D_icpp_stl_circle/case.py @@ -88,7 +88,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/2D_isentropicvortex/case.py b/examples/2D_isentropicvortex/case.py index 9ca22b4e0c..378a4c62b5 100644 --- a/examples/2D_isentropicvortex/case.py +++ b/examples/2D_isentropicvortex/case.py @@ -110,7 +110,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_isentropicvortex_analytical/case.py b/examples/2D_isentropicvortex_analytical/case.py index 3c2b768bc9..792c642c6a 100644 --- a/examples/2D_isentropicvortex_analytical/case.py +++ b/examples/2D_isentropicvortex_analytical/case.py @@ -110,7 +110,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_jet/case.py b/examples/2D_jet/case.py index 93e8ea9484..6cf73d72f3 100644 --- a/examples/2D_jet/case.py +++ b/examples/2D_jet/case.py @@ -109,10 +109,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_kelvin_helmholtz/case.py b/examples/2D_kelvin_helmholtz/case.py index 861490890f..2de4acaf5c 100644 --- a/examples/2D_kelvin_helmholtz/case.py +++ b/examples/2D_kelvin_helmholtz/case.py @@ -81,10 +81,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(2)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_lagrange_bubblescreen/case.py b/examples/2D_lagrange_bubblescreen/case.py index 3771df419a..adb3c7c199 100644 --- a/examples/2D_lagrange_bubblescreen/case.py +++ b/examples/2D_lagrange_bubblescreen/case.py @@ -242,7 +242,6 @@ def generate_bubble_cloud(): # Bubble gas state "fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)), }, indent=4, diff --git a/examples/2D_lagrange_in_crossflow/case.py b/examples/2D_lagrange_in_crossflow/case.py index 4682e3c631..8d15478de1 100644 --- a/examples/2D_lagrange_in_crossflow/case.py +++ b/examples/2D_lagrange_in_crossflow/case.py @@ -182,7 +182,6 @@ def write_lag_bubbles_file(): # Bubble gas state "fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)), } ) diff --git a/examples/2D_lagrange_rising_bubble/case.py b/examples/2D_lagrange_rising_bubble/case.py index f2c134add1..0273e63a3f 100644 --- a/examples/2D_lagrange_rising_bubble/case.py +++ b/examples/2D_lagrange_rising_bubble/case.py @@ -185,7 +185,6 @@ # Bubble gas state "fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)), }, indent=4, diff --git a/examples/2D_laplace_pressure_jump/case.py b/examples/2D_laplace_pressure_jump/case.py index b6a8f8ab7c..46a8fd000f 100644 --- a/examples/2D_laplace_pressure_jump/case.py +++ b/examples/2D_laplace_pressure_jump/case.py @@ -78,7 +78,6 @@ # Fluid Parameters (Gas) "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, # Air Patch "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0, diff --git a/examples/2D_lid_driven_cavity/case.py b/examples/2D_lid_driven_cavity/case.py index ae9a957dbb..9a4ae1eee1 100644 --- a/examples/2D_lid_driven_cavity/case.py +++ b/examples/2D_lid_driven_cavity/case.py @@ -66,11 +66,9 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1e4, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%Re(1)": 1e4, } ) diff --git a/examples/2D_lid_driven_cavity_nn/case.py b/examples/2D_lid_driven_cavity_nn/case.py index fdbdd8c7f8..a6e27a3a5e 100644 --- a/examples/2D_lid_driven_cavity_nn/case.py +++ b/examples/2D_lid_driven_cavity_nn/case.py @@ -90,7 +90,6 @@ # mu_eff at each cell via the Papanastasiou-regularized formula. "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.0 / K, "fluid_pp(1)%non_newtonian": "T", "fluid_pp(1)%tau0": 0.0, @@ -101,7 +100,6 @@ "fluid_pp(1)%hb_m": 1000.0, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%Re(1)": 1.0 / K, "fluid_pp(2)%non_newtonian": "T", "fluid_pp(2)%tau0": 0.0, diff --git a/examples/2D_mhd_magnetic_vortex/case.py b/examples/2D_mhd_magnetic_vortex/case.py index 3118767de2..1cd2e4d2c3 100644 --- a/examples/2D_mhd_magnetic_vortex/case.py +++ b/examples/2D_mhd_magnetic_vortex/case.py @@ -72,7 +72,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_mhd_rotor/case.py b/examples/2D_mhd_rotor/case.py index 011e0241ca..5aa001e522 100644 --- a/examples/2D_mhd_rotor/case.py +++ b/examples/2D_mhd_rotor/case.py @@ -79,7 +79,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_mibm_cylinder_in_cross_flow/case.py b/examples/2D_mibm_cylinder_in_cross_flow/case.py index 62ccc83fec..0c984be6e5 100644 --- a/examples/2D_mibm_cylinder_in_cross_flow/case.py +++ b/examples/2D_mibm_cylinder_in_cross_flow/case.py @@ -100,7 +100,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_mibm_particle_cloud/case.py b/examples/2D_mibm_particle_cloud/case.py index a503eb260c..d9e73a0221 100644 --- a/examples/2D_mibm_particle_cloud/case.py +++ b/examples/2D_mibm_particle_cloud/case.py @@ -133,7 +133,6 @@ # Fluid properties: air "fluid_pp(1)%gamma": 1.0 / (gam_a - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_mibm_shock_cylinder/case.py b/examples/2D_mibm_shock_cylinder/case.py index a7bdc86177..4b904332c3 100644 --- a/examples/2D_mibm_shock_cylinder/case.py +++ b/examples/2D_mibm_shock_cylinder/case.py @@ -134,7 +134,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_moving_lag_bubs/case.py b/examples/2D_moving_lag_bubs/case.py index 09710c2ae4..adaf4e5c04 100644 --- a/examples/2D_moving_lag_bubs/case.py +++ b/examples/2D_moving_lag_bubs/case.py @@ -125,7 +125,6 @@ # Fluid Parameters Gas "fluid_pp(2)%gamma": 1.0 / (gamma_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, # Bubble parameters "bub_pp%R0ref": 1.0, "bub_pp%p0ref": 1.0, diff --git a/examples/2D_orszag_tang/case.py b/examples/2D_orszag_tang/case.py index 6f34d3f298..123c9018b4 100644 --- a/examples/2D_orszag_tang/case.py +++ b/examples/2D_orszag_tang/case.py @@ -71,7 +71,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_orszag_tang_hyper_cleaning/case.py b/examples/2D_orszag_tang_hyper_cleaning/case.py index c5440ad6bc..a85a5c0219 100644 --- a/examples/2D_orszag_tang_hyper_cleaning/case.py +++ b/examples/2D_orszag_tang_hyper_cleaning/case.py @@ -75,7 +75,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (5.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_patch_modal_shape/case.py b/examples/2D_patch_modal_shape/case.py index e28164c5b5..ab726ccd19 100644 --- a/examples/2D_patch_modal_shape/case.py +++ b/examples/2D_patch_modal_shape/case.py @@ -106,10 +106,8 @@ "bc_y%pres_out": p_inf, "fluid_pp(1)%gamma": 1.0 / (gam - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0 / (gam - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, # Acoustic source (similar to 2D_acoustic_support5) "acoustic_source": "T", "num_source": 1, diff --git a/examples/2D_patch_modal_shape_exp/case.py b/examples/2D_patch_modal_shape_exp/case.py index 14b00ef120..5795a7eeb2 100644 --- a/examples/2D_patch_modal_shape_exp/case.py +++ b/examples/2D_patch_modal_shape_exp/case.py @@ -108,10 +108,8 @@ "bc_y%pres_out": p_inf, "fluid_pp(1)%gamma": 1.0 / (gam - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0 / (gam - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, # Acoustic source (similar to 2D_acoustic_support5) "acoustic_source": "T", "num_source": 1, diff --git a/examples/2D_poiseuille_nn/case.py b/examples/2D_poiseuille_nn/case.py index 252089772c..4b93eb85e1 100644 --- a/examples/2D_poiseuille_nn/case.py +++ b/examples/2D_poiseuille_nn/case.py @@ -124,7 +124,6 @@ # Fluids Physical Parameters: single power-law (HB with tau0 = 0) fluid "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.0 / K, "fluid_pp(1)%non_newtonian": "T", "fluid_pp(1)%K": K, diff --git a/examples/2D_poiseuille_thickening_nn/case.py b/examples/2D_poiseuille_thickening_nn/case.py index f459f8f6e8..223b9cd140 100644 --- a/examples/2D_poiseuille_thickening_nn/case.py +++ b/examples/2D_poiseuille_thickening_nn/case.py @@ -118,7 +118,6 @@ # Fluids Physical Parameters: single power-law (HB with tau0 = 0) fluid "fluid_pp(1)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.0 / K, "fluid_pp(1)%non_newtonian": "T", "fluid_pp(1)%K": K, diff --git a/examples/2D_rayleigh_taylor/case.py b/examples/2D_rayleigh_taylor/case.py index c37a280c79..da6d1ec57c 100644 --- a/examples/2D_rayleigh_taylor/case.py +++ b/examples/2D_rayleigh_taylor/case.py @@ -73,12 +73,10 @@ # Fluid Parameters (Heavy Gas) "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(1)%Re(1)": 1 / 0.0219, # Fluid Parameters (Light Gas) "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1 / 0.0073, # Body Forces "bf_y": "T", diff --git a/examples/2D_richtmyer_meshkov/case.py b/examples/2D_richtmyer_meshkov/case.py index b9bfdad0bd..f06b5c7973 100644 --- a/examples/2D_richtmyer_meshkov/case.py +++ b/examples/2D_richtmyer_meshkov/case.py @@ -86,10 +86,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.093 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(2)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "viscous": "T", "fluid_pp(1)%Re(1)": 1 / mu, "fluid_pp(2)%Re(1)": 1 / mu, diff --git a/examples/2D_riemann_test/case.py b/examples/2D_riemann_test/case.py index b3af4e5578..d6e5f566a3 100644 --- a/examples/2D_riemann_test/case.py +++ b/examples/2D_riemann_test/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_riemann_test_muscl/case.py b/examples/2D_riemann_test_muscl/case.py index 79c1a24fd0..b386ba7080 100644 --- a/examples/2D_riemann_test_muscl/case.py +++ b/examples/2D_riemann_test_muscl/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_shock_cloud_rmhd/case.py b/examples/2D_shock_cloud_rmhd/case.py index 26445bad16..52636a009f 100644 --- a/examples/2D_shock_cloud_rmhd/case.py +++ b/examples/2D_shock_cloud_rmhd/case.py @@ -102,7 +102,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (4.0 / 3.0 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_shockbubble/case.py b/examples/2D_shockbubble/case.py index 00c8894ae4..d808469158 100644 --- a/examples/2D_shockbubble/case.py +++ b/examples/2D_shockbubble/case.py @@ -107,10 +107,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / (1.6666e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_shockdroplet/case.py b/examples/2D_shockdroplet/case.py index cf778a3bd2..95160ab9b3 100755 --- a/examples/2D_shockdroplet/case.py +++ b/examples/2D_shockdroplet/case.py @@ -159,7 +159,6 @@ "fluid_pp(1)%pi_inf": pi_w * gam_w / (gam_w - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_shockdroplet_muscl/case.py b/examples/2D_shockdroplet_muscl/case.py index 67767cbe16..61b904ec85 100755 --- a/examples/2D_shockdroplet_muscl/case.py +++ b/examples/2D_shockdroplet_muscl/case.py @@ -119,7 +119,6 @@ "fluid_pp(1)%pi_inf": pi_w * gam_w / (gam_w - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/2D_synthetic_turbulence/case.py b/examples/2D_synthetic_turbulence/case.py index 9f3a9803e5..a52bccb460 100644 --- a/examples/2D_synthetic_turbulence/case.py +++ b/examples/2D_synthetic_turbulence/case.py @@ -122,10 +122,8 @@ # Use the same stiffness as the air bubble "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50 (Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(2)%gamma": 1.0e00 / (gam_b - 1.0e00), # 2.50 (Not 1.40) "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0, # -- Synthetic turbulence -- "synthetic_turbulence": "T", "synth_seed": 42, diff --git a/examples/2D_triple_point/case.py b/examples/2D_triple_point/case.py index 65da52f2a9..96ee052daf 100755 --- a/examples/2D_triple_point/case.py +++ b/examples/2D_triple_point/case.py @@ -99,10 +99,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0 / (1.5 - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_tumbling_rectangle/case.py b/examples/2D_tumbling_rectangle/case.py index 0eddad239e..a28a55bb3c 100644 --- a/examples/2D_tumbling_rectangle/case.py +++ b/examples/2D_tumbling_rectangle/case.py @@ -100,7 +100,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/2D_viscous_shock_tube/case.py b/examples/2D_viscous_shock_tube/case.py index ead361d30e..a839c8b4e4 100644 --- a/examples/2D_viscous_shock_tube/case.py +++ b/examples/2D_viscous_shock_tube/case.py @@ -82,11 +82,9 @@ "fluid_pp(1)%gamma": 1.0e00 / ((7.0 / 5.0) - 1.0e00), "fluid_pp(1)%Re(1)": 1 / mu, "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / ((7.0 / 5.0) - 1.0e00), "fluid_pp(2)%Re(1)": 1 / mu, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_zero_circ_vortex/case.py b/examples/2D_zero_circ_vortex/case.py index 87974f48e8..a891751c3d 100644 --- a/examples/2D_zero_circ_vortex/case.py +++ b/examples/2D_zero_circ_vortex/case.py @@ -103,7 +103,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/2D_zero_circ_vortex_analytical/case.py b/examples/2D_zero_circ_vortex_analytical/case.py index 81ac639b21..c56cbe74f5 100644 --- a/examples/2D_zero_circ_vortex_analytical/case.py +++ b/examples/2D_zero_circ_vortex_analytical/case.py @@ -104,7 +104,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/3D_IGR_33jet/case.py b/examples/3D_IGR_33jet/case.py index e30ffd75b1..8845d839b0 100644 --- a/examples/3D_IGR_33jet/case.py +++ b/examples/3D_IGR_33jet/case.py @@ -111,7 +111,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "viscous": "T", "fluid_pp(1)%Re(1)": 5e5, }, diff --git a/examples/3D_IGR_TaylorGreenVortex/case.py b/examples/3D_IGR_TaylorGreenVortex/case.py index fec431e35a..5ca2b90990 100644 --- a/examples/3D_IGR_TaylorGreenVortex/case.py +++ b/examples/3D_IGR_TaylorGreenVortex/case.py @@ -91,7 +91,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 1 / mu, } ) diff --git a/examples/3D_IGR_TaylorGreenVortex_nvidia/case.py b/examples/3D_IGR_TaylorGreenVortex_nvidia/case.py index e2dd18f929..72838e8a1e 100644 --- a/examples/3D_IGR_TaylorGreenVortex_nvidia/case.py +++ b/examples/3D_IGR_TaylorGreenVortex_nvidia/case.py @@ -96,7 +96,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 1 / mu, # NVIDIA UVM Options "nv_uvm_out_of_core": "T", diff --git a/examples/3D_IGR_jet/case.py b/examples/3D_IGR_jet/case.py index e94a654e88..7fabd9248c 100644 --- a/examples/3D_IGR_jet/case.py +++ b/examples/3D_IGR_jet/case.py @@ -118,10 +118,8 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "viscous": "T", "fluid_pp(1)%Re(1)": 5e4, "fluid_pp(2)%Re(1)": 5e4, diff --git a/examples/3D_IGR_jet_1fluid/case.py b/examples/3D_IGR_jet_1fluid/case.py index 06931c6d5c..a495f81f10 100644 --- a/examples/3D_IGR_jet_1fluid/case.py +++ b/examples/3D_IGR_jet_1fluid/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "viscous": "T", "fluid_pp(1)%Re(1)": 5e4, }, diff --git a/examples/3D_TaylorGreenVortex/case.py b/examples/3D_TaylorGreenVortex/case.py index a46dcf670f..4896230e94 100644 --- a/examples/3D_TaylorGreenVortex/case.py +++ b/examples/3D_TaylorGreenVortex/case.py @@ -95,7 +95,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 1 / mu, } ) diff --git a/examples/3D_TaylorGreenVortex_analytical/case.py b/examples/3D_TaylorGreenVortex_analytical/case.py index 9339f0e7fa..2d6bcf65ca 100644 --- a/examples/3D_TaylorGreenVortex_analytical/case.py +++ b/examples/3D_TaylorGreenVortex_analytical/case.py @@ -94,7 +94,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 1 / mu, } ) diff --git a/examples/3D_advection_convergence/case.py b/examples/3D_advection_convergence/case.py index 768f3e6727..c983d5aeb8 100644 --- a/examples/3D_advection_convergence/case.py +++ b/examples/3D_advection_convergence/case.py @@ -113,7 +113,6 @@ "patch_icpp(1)%alpha(1)": 1.0, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, **scheme_params, } ) diff --git a/examples/3D_brio_wu/case.py b/examples/3D_brio_wu/case.py index 136ad35b44..c254842388 100644 --- a/examples/3D_brio_wu/case.py +++ b/examples/3D_brio_wu/case.py @@ -90,7 +90,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (2.0e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, } ) ) diff --git a/examples/3D_ibm_bowshock/case.py b/examples/3D_ibm_bowshock/case.py index 656964d71c..78814d5f5c 100644 --- a/examples/3D_ibm_bowshock/case.py +++ b/examples/3D_ibm_bowshock/case.py @@ -102,7 +102,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 7535533.2, } ) diff --git a/examples/3D_ibm_stl_ellipsoid/case.py b/examples/3D_ibm_stl_ellipsoid/case.py index 40932f1a05..85d2da0e35 100644 --- a/examples/3D_ibm_stl_ellipsoid/case.py +++ b/examples/3D_ibm_stl_ellipsoid/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 7535533.2, } ) diff --git a/examples/3D_ibm_stl_pyramid/case.py b/examples/3D_ibm_stl_pyramid/case.py index b8d3cd56a1..94dc1cb431 100644 --- a/examples/3D_ibm_stl_pyramid/case.py +++ b/examples/3D_ibm_stl_pyramid/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 7535533.2, } ) diff --git a/examples/3D_ibm_stl_test/case.py b/examples/3D_ibm_stl_test/case.py index 0c28eb63a7..83fae2efe6 100644 --- a/examples/3D_ibm_stl_test/case.py +++ b/examples/3D_ibm_stl_test/case.py @@ -93,7 +93,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/3D_icpp_stl_cube/case.py b/examples/3D_icpp_stl_cube/case.py index c4cea9a090..1a72d2fab3 100644 --- a/examples/3D_icpp_stl_cube/case.py +++ b/examples/3D_icpp_stl_cube/case.py @@ -97,7 +97,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/3D_lagrange_bubblescreen/case.py b/examples/3D_lagrange_bubblescreen/case.py index db0c8a6835..d2fb323b75 100644 --- a/examples/3D_lagrange_bubblescreen/case.py +++ b/examples/3D_lagrange_bubblescreen/case.py @@ -230,7 +230,6 @@ def generate_bubble_cloud(): # Bubble gas state "fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)), }, indent=4, diff --git a/examples/3D_lagrange_shbubcollapse/case.py b/examples/3D_lagrange_shbubcollapse/case.py index 0ee96f8ef7..8d637ee3f3 100644 --- a/examples/3D_lagrange_shbubcollapse/case.py +++ b/examples/3D_lagrange_shbubcollapse/case.py @@ -184,7 +184,6 @@ # Bubble gas state "fluid_pp(2)%gamma": 1.0 / (gam_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1.0 / (mu_g / (rho0 * c0 * x0)), } ) diff --git a/examples/3D_mibm_periodic_collision/case.py b/examples/3D_mibm_periodic_collision/case.py index 8e0e854366..e6f499c6c5 100644 --- a/examples/3D_mibm_periodic_collision/case.py +++ b/examples/3D_mibm_periodic_collision/case.py @@ -152,7 +152,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 1.0 / mu, "collision_model": 1, # soft-sphere collision model "ib_coefficient_of_friction": 0.1, diff --git a/examples/3D_mibm_sphere_head_on_collision/case.py b/examples/3D_mibm_sphere_head_on_collision/case.py index a14b0815c8..86a70bb001 100644 --- a/examples/3D_mibm_sphere_head_on_collision/case.py +++ b/examples/3D_mibm_sphere_head_on_collision/case.py @@ -125,7 +125,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, } ) ) diff --git a/examples/3D_moving_lag_particles/case.py b/examples/3D_moving_lag_particles/case.py index b185e8fd10..1377266187 100644 --- a/examples/3D_moving_lag_particles/case.py +++ b/examples/3D_moving_lag_particles/case.py @@ -120,7 +120,6 @@ # Fluid Parameters Gas "fluid_pp(2)%gamma": 1.0 / (gamma_g - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, # Bubble parameters "bub_pp%R0ref": 1.0, "bub_pp%p0ref": 1.0, diff --git a/examples/3D_patch_spherical_harmonic/case.py b/examples/3D_patch_spherical_harmonic/case.py index 3efa3a1e2b..3db6c2370c 100644 --- a/examples/3D_patch_spherical_harmonic/case.py +++ b/examples/3D_patch_spherical_harmonic/case.py @@ -126,10 +126,8 @@ "patch_icpp(2)%alter_patch(1)": "T", "fluid_pp(1)%gamma": 1.0 / (gam - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(2)%gamma": 1.0 / (gam - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, # Acoustic source (similar to 3D_acoustic_support7) "acoustic_source": "T", "num_source": 1, diff --git a/examples/3D_performance_test/case.py b/examples/3D_performance_test/case.py index ffad84ee66..6749ab35e9 100644 --- a/examples/3D_performance_test/case.py +++ b/examples/3D_performance_test/case.py @@ -103,7 +103,6 @@ "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/3D_rayleigh_taylor/case.py b/examples/3D_rayleigh_taylor/case.py index 3c273651ec..489aea8675 100644 --- a/examples/3D_rayleigh_taylor/case.py +++ b/examples/3D_rayleigh_taylor/case.py @@ -80,12 +80,10 @@ # Fluid Parameters (Heavy Gas) "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(1)%Re(1)": 1 / 0.0219, # Fluid Parameters (Light Gas) "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1 / 0.0073, # Body Forces "bf_y": "T", diff --git a/examples/3D_rayleigh_taylor_muscl/case.py b/examples/3D_rayleigh_taylor_muscl/case.py index f8e7324dba..4161ec6c33 100644 --- a/examples/3D_rayleigh_taylor_muscl/case.py +++ b/examples/3D_rayleigh_taylor_muscl/case.py @@ -78,12 +78,10 @@ # Fluid Parameters (Heavy Gas) "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0e00, "fluid_pp(1)%Re(1)": 1 / 0.0219, # Fluid Parameters (Light Gas) "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%Re(1)": 1 / 0.0073, # Body Forces "bf_y": "T", diff --git a/examples/3D_recovering_sphere/case.py b/examples/3D_recovering_sphere/case.py index ada09cf9f4..6a964adc85 100644 --- a/examples/3D_recovering_sphere/case.py +++ b/examples/3D_recovering_sphere/case.py @@ -86,7 +86,6 @@ # Fluid Parameters (Gas) "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, # Air Patch "patch_icpp(1)%geometry": 9, "patch_icpp(1)%x_centroid": 0, diff --git a/examples/3D_rotating_sphere/case.py b/examples/3D_rotating_sphere/case.py index 141e10f198..18b790331a 100644 --- a/examples/3D_rotating_sphere/case.py +++ b/examples/3D_rotating_sphere/case.py @@ -105,7 +105,6 @@ # Fluids Physical Parameters "fluid_pp(1)%gamma": 1.0e00 / (gam_a - 1.0e00), # 2.50(Not 1.40) "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "fluid_pp(1)%Re(1)": 2500000, } ) diff --git a/examples/3D_sphbubcollapse/case.py b/examples/3D_sphbubcollapse/case.py index b57230f2f3..3d6e0df42c 100644 --- a/examples/3D_sphbubcollapse/case.py +++ b/examples/3D_sphbubcollapse/case.py @@ -129,7 +129,6 @@ "fluid_pp(1)%pi_inf": 4.4e00 * 6.0e08 / (4.4e00 - 1.0e00), "fluid_pp(2)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, } ) ) diff --git a/examples/3D_turb_mixing/case.py b/examples/3D_turb_mixing/case.py index 81f0e44769..a617c76d26 100644 --- a/examples/3D_turb_mixing/case.py +++ b/examples/3D_turb_mixing/case.py @@ -118,7 +118,6 @@ # Surrounding liquid "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": Re0, } ) diff --git a/toolchain/mfc/case_validator.py b/toolchain/mfc/case_validator.py index 0dbe510ce4..9d73512455 100644 --- a/toolchain/mfc/case_validator.py +++ b/toolchain/mfc/case_validator.py @@ -44,8 +44,7 @@ "category": "Thermodynamic Constraints", "math": r"\Pi_\infty = 0 \;\; \text{for an ideal gas}", "explanation": ( - "An ideal gas is the stiffened-gas equation of state with no stiffness, so selecting it and also supplying a " - "nonzero pi_inf is contradictory. The selector determines the stiffness, not the input." + "An ideal gas is the stiffened-gas equation of state with no stiffness, so a case that selects it may not " "set pi_inf at all: the selector determines the stiffness, not the input." ), "references": ["Wilfong26"], }, @@ -958,10 +957,9 @@ def check_eos_selector(self): if eos is None: continue self.prohibit(eos not in eos_values, f"fluid_pp({i})%eos must be 'stiffened_gas' or 'ideal_gas'") - pi_inf = self.get(f"fluid_pp({i})%pi_inf") self.prohibit( - eos == eos_ideal_gas and pi_inf is not None and pi_inf != 0, - f"fluid_pp({i})%eos = 'ideal_gas' requires fluid_pp({i})%pi_inf = 0 (an ideal gas has no stiffness)", + eos == eos_ideal_gas and self.get(f"fluid_pp({i})%pi_inf") is not None, + f"fluid_pp({i})%eos = 'ideal_gas' has no stiffness; do not set fluid_pp({i})%pi_inf", ) def check_stiffened_eos(self): diff --git a/toolchain/mfc/params_tests/test_eos_selector.py b/toolchain/mfc/params_tests/test_eos_selector.py index 418481dbc8..02554b556f 100644 --- a/toolchain/mfc/params_tests/test_eos_selector.py +++ b/toolchain/mfc/params_tests/test_eos_selector.py @@ -30,8 +30,10 @@ def test_stiffened_gas_accepts_stiffness(): assert _errors(**{"fluid_pp(1)%eos": EOS["stiffened_gas"], "fluid_pp(1)%pi_inf": 1.0e5}) == [] -def test_ideal_gas_accepts_zero_stiffness(): - assert _errors(**{"fluid_pp(1)%eos": EOS["ideal_gas"], "fluid_pp(1)%pi_inf": 0.0}) == [] +def test_ideal_gas_owns_no_stiffness_parameter(): + assert _errors(**{"fluid_pp(1)%eos": EOS["ideal_gas"], "fluid_pp(1)%pi_inf": None}) == [] + errors = _errors(**{"fluid_pp(1)%eos": EOS["ideal_gas"], "fluid_pp(1)%pi_inf": 0.0}) + assert len(errors) == 1 and "do not set fluid_pp(1)%pi_inf" in errors[0] def test_ideal_gas_rejects_nonzero_stiffness(): diff --git a/toolchain/mfc/test/case.py b/toolchain/mfc/test/case.py index 14dfc4ef94..2d9f241b3b 100644 --- a/toolchain/mfc/test/case.py +++ b/toolchain/mfc/test/case.py @@ -101,7 +101,6 @@ def get_post_process_mods(case_params: dict) -> dict: "patch_icpp(3)%alpha(1)": 1.0, "fluid_pp(1)%gamma": 1.0e00 / (1.4 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%cv": 0.0, "fluid_pp(1)%qv": 0.0, "fluid_pp(1)%qvp": 0.0, diff --git a/toolchain/mfc/test/cases.py b/toolchain/mfc/test/cases.py index d86ceea636..34a8c4bc5a 100644 --- a/toolchain/mfc/test/cases.py +++ b/toolchain/mfc/test/cases.py @@ -531,9 +531,7 @@ def alter_muscl(): stack.pop() def alter_eos(): - # BASE_CFG already sets fluid_pp(1)%pi_inf = 0, so the fluid is an ideal gas either way and - # naming it one must not change a single digit. That is the point: it exercises the selector - # end to end without moving the physics. + # BASE_CFG's fluid is already an ideal gas, so this re-selection must not change a single digit. cases.append(define_case_d(stack, "eos=ideal_gas", {"fluid_pp(1)%eos": "ideal_gas"})) def alter_riemann_solvers(num_fluids): @@ -647,7 +645,6 @@ def alter_num_fluids(dimInfo): { "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "patch_icpp(1)%alpha_rho(1)": 0.81, "patch_icpp(1)%alpha(1)": 0.9, "patch_icpp(1)%alpha_rho(2)": 0.19, @@ -804,7 +801,6 @@ def alter_2d(): "cyl_coord": "T", "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "patch_icpp(1)%alpha_rho(1)": 0.81, "patch_icpp(1)%alpha(1)": 0.9, "patch_icpp(1)%alpha_rho(2)": 0.19, @@ -883,7 +879,6 @@ def alter_3d(): "num_fluids": 2, "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "patch_icpp(1)%alpha_rho(1)": 0.81, "patch_icpp(1)%alpha(1)": 0.9, "patch_icpp(1)%alpha_rho(2)": 0.19, @@ -1672,7 +1667,6 @@ def alter_mixlayer_perturb(dimInfo): "wenoz": "T", "fluid_pp(1)%gamma": 2.5, "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 1.6881644098979287, "viscous": "T", "patch_icpp(1)%geometry": 9, @@ -1743,7 +1737,6 @@ def alter_phasechange(dimInfo): "fluid_pp(1)%qvp": 0.0, "fluid_pp(2)%gamma": 2.3266, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0e00, "fluid_pp(2)%cv": 1040, "fluid_pp(2)%qv": 2030000, "fluid_pp(2)%qvp": -23400, @@ -1771,7 +1764,6 @@ def alter_phasechange(dimInfo): { "fluid_pp(3)%gamma": 2.4870, "fluid_pp(3)%eos": "ideal_gas", - "fluid_pp(3)%pi_inf": 0.0e00, "fluid_pp(3)%cv": 717.5, "fluid_pp(3)%qv": 0.0e00, "fluid_pp(3)%qvp": 0.0, @@ -1923,7 +1915,6 @@ def alter_lag_bubbles(dimInfo): "fluid_pp(1)%pi_inf": 3515.0, "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "patch_icpp(1)%alpha_rho(1)": 0.96, "patch_icpp(1)%alpha(1)": 4e-02, "patch_icpp(1)%alpha_rho(2)": 0.0, @@ -2464,11 +2455,9 @@ def apply_solver(case: dict, base_cfg: dict, solver_mods: dict, alt_soundspeed: "dt": 0.2 * dx / c_outer, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%G": G_solid, "fluid_pp(2)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": 0.0, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0.5, @@ -2537,11 +2526,9 @@ def make_shear_guard_cfg(regime): "dt": 0.2 * dx / c_outer, "fluid_pp(1)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%G": G, "fluid_pp(2)%gamma": 1.0 / (gamma - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": G, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0.5, @@ -2613,11 +2600,9 @@ def make_axisym_fallback_cfg(adc): "dt": 2.0e-3, "fluid_pp(1)%gamma": 2.5, "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%G": G, "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": G, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0.5, @@ -2696,7 +2681,6 @@ def make_cbc_uniform_cfg(alt_soundspeed): "fluid_pp(1)%G": 0.0, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": 0.0, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0.5, @@ -2779,7 +2763,6 @@ def make_cbc_uniform_cfg(alt_soundspeed): "fluid_pp(1)%G": 1.0e6, "fluid_pp(2)%gamma": 1.0 / (1.4 - 1.0), "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, "fluid_pp(2)%G": 1.0e6, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%x_centroid": 0.5, @@ -3093,7 +3076,6 @@ def chemistry_cases(): "viscous": "T", "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 100000, "patch_icpp(1)%geometry": 1, "patch_icpp(1)%hcid": 191, @@ -3147,7 +3129,6 @@ def chemistry_cases(): "viscous": "T", "fluid_pp(1)%gamma": 1.0e00 / (1.4e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0.0, "fluid_pp(1)%Re(1)": 100000, "patch_icpp(1)%geometry": 3, "patch_icpp(1)%hcid": 291, @@ -3201,7 +3182,6 @@ def chemistry_cases(): "patch_icpp(1)%alpha(1)": 1, "fluid_pp(1)%gamma": 1.0e00 / (1.9326e00 - 1.0e00), "fluid_pp(1)%eos": "ideal_gas", - "fluid_pp(1)%pi_inf": 0, "cantera_file": "h2o2.yaml", "t_step_start": 0, "t_step_stop": 50, @@ -3669,7 +3649,6 @@ def kernel_golden_tests(): "num_fluids": 2, "fluid_pp(2)%gamma": 2.5, "fluid_pp(2)%eos": "ideal_gas", - "fluid_pp(2)%pi_inf": 0.0, # Patch 1: fluid 1 background rectangle; length covers stretched extent (~1.39). # vel(1)=0.5 provides advection so MTHINC reconstruction affects the solution. "patch_icpp(1)%geometry": 3, diff --git a/toolchain/mfc/test/coverage.py b/toolchain/mfc/test/coverage.py index 26d664affb..e279fc4b93 100644 --- a/toolchain/mfc/test/coverage.py +++ b/toolchain/mfc/test/coverage.py @@ -216,8 +216,14 @@ def select_tests(cases, coverage_map, changed_files): return to_run, skipped, f"selected {len(to_run)}/{len(cases)} by coverage overlap" +def _env_without_git(): + # Git exports GIT_DIR and GIT_INDEX_FILE to hooks, and neither cwd nor `git -C` overrides them: + # under the pre-commit hook every call below would otherwise act on the committing repository. + return {k: v for k, v in os.environ.items() if not k.startswith("GIT_")} + + def _git(args, cwd, timeout=60): - return subprocess.run(["git", *args], capture_output=True, text=True, cwd=cwd, timeout=timeout, check=False) + return subprocess.run(["git", *args], capture_output=True, text=True, cwd=cwd, timeout=timeout, check=False, env=_env_without_git()) def _merge_base(cwd, branch): diff --git a/toolchain/mfc/test/test_coverage_unit.py b/toolchain/mfc/test/test_coverage_unit.py index 24f52ec1aa..5e9ad49e3c 100644 --- a/toolchain/mfc/test/test_coverage_unit.py +++ b/toolchain/mfc/test/test_coverage_unit.py @@ -6,7 +6,7 @@ from pathlib import Path from unittest.mock import patch -from mfc.test.coverage import canonicalize_param_paths, entries_equal, format_summary, get_changed_files, is_always_run_all, load_map, map_health, param_hash, save_map, select_tests +from mfc.test.coverage import _env_without_git, canonicalize_param_paths, entries_equal, format_summary, get_changed_files, is_always_run_all, load_map, map_health, param_hash, save_map, select_tests def test_param_hash_is_order_independent(): @@ -462,17 +462,6 @@ def test_health_fails_immediately_when_no_refresh_ran_since_last_source_change() CHANGED_SCRIPT = Path(__file__).resolve().parents[3] / ".github" / "scripts" / "coverage_map_changed.py" -def _env_without_git(): - """The environment minus every GIT_* variable. - - `git -C ` changes directory but does NOT override an inherited GIT_DIR or - GIT_INDEX_FILE. Git exports both when it runs a hook, and MFC's pre-commit hook runs - precheck, which runs this suite -- so without this scrub the commits below are made - against the real repository instead of the throwaway one. - """ - return {k: v for k, v in os.environ.items() if not k.startswith("GIT_")} - - def _repo_with_committed_map(d, entries): """A throwaway git repo whose HEAD holds `entries` as the coverage map.""" repo = Path(d)