torus problem in tutorial

AFD-Illinois · Oct 21, 2024 · 33aa664 · 33aa664
1 parent 609ba00
commit 33aa664
Show file tree

Hide file tree

Showing 4 changed files with 215 additions and 13 deletions.
diff --git a/scripts/tutorial/.ipynb_checkpoints/tutorial-checkpoint.ipynb b/scripts/tutorial/.ipynb_checkpoints/tutorial-checkpoint.ipynb
diff --git a/scripts/tutorial/analysis/plot_torus.py b/scripts/tutorial/analysis/plot_torus.py
@@ -71,8 +71,6 @@ def plot(dumpno):
 
     # load dump
     dump = pyharm.load_dump(os.path.join(dumpsdir, 'torus.out0.{:05d}.phdf'.format(dumpno)))
-
-    bsq = np.maximum(dump['bsq'], SMALL)
 
     x_slice   = xz_slice(np.squeeze(dump['x']), dump, patch_pole=True)
     z_slice   = xz_slice(np.squeeze(dump['z']), dump)

diff --git a/scripts/tutorial/parfiles/sane2d.par b/scripts/tutorial/parfiles/sane2d.par
@@ -66,7 +66,7 @@ bsq_over_rho_max = 100
 
 <parthenon/output0>
 file_type = hdf5
-dt = 10.0
+dt = 5.0
 single_precision_output = true
 # Fields beginning with a specifier here will be included (e.g., all prims.XXX)
 # Fields specified but not present are silently skipped

diff --git a/scripts/tutorial/tutorial.ipynb b/scripts/tutorial/tutorial.ipynb