From b8deac0f852209403d4b789d09931f6c5ccc121f Mon Sep 17 00:00:00 2001 From: Ben Prather Date: Fri, 4 Oct 2024 08:12:23 -0600 Subject: [PATCH] Excising the polar fluxes has proven unstable, disable by default --- pars/bondi/bondi_b_vertical.par | 4 +++- pars/tori_3d/mad.par | 5 +++-- pars/tori_3d/sane.par | 5 +++-- 3 files changed, 9 insertions(+), 5 deletions(-) diff --git a/pars/bondi/bondi_b_vertical.par b/pars/bondi/bondi_b_vertical.par index 31380300..6e165083 100644 --- a/pars/bondi/bondi_b_vertical.par +++ b/pars/bondi/bondi_b_vertical.par @@ -55,7 +55,9 @@ on = true # Transmitting boundaries, since there will be material at the pole inner_x2 = transmitting outer_x2 = transmitting -excise_polar_flux = true +# This might additionally help, but the option is unstable +#excise_polar_flux = true + # Allow in material from boundary check_inflow_outer_x1 = false diff --git a/pars/tori_3d/mad.par b/pars/tori_3d/mad.par index 47579314..809217c9 100644 --- a/pars/tori_3d/mad.par +++ b/pars/tori_3d/mad.par @@ -66,8 +66,9 @@ two_sync = false inner_x2 = transmitting outer_x2 = transmitting -# This reduces the polar "wake," but may shrink timestep -excise_polar_flux = true +# This reduces the polar "wake," but shrinks the timestep, +# and has caused instabilities. YMMV. +#excise_polar_flux = true frame = drift diff --git a/pars/tori_3d/sane.par b/pars/tori_3d/sane.par index b1f05cbb..eccdf7a4 100644 --- a/pars/tori_3d/sane.par +++ b/pars/tori_3d/sane.par @@ -64,8 +64,9 @@ two_sync = false inner_x2 = transmitting outer_x2 = transmitting -# This reduces the polar "wake," but may shrink timestep -excise_polar_flux = true +# This reduces the polar "wake," but shrinks the timestep, +# and has caused instabilities. YMMV. +#excise_polar_flux = true # Now with FOFC it might work to set "frame = normal",