From 9817d40d4629ad17f8de9f380b467656ca8df932 Mon Sep 17 00:00:00 2001 From: Claire Kopenhafer Date: Tue, 17 May 2022 20:07:27 -0400 Subject: [PATCH] Testing tcool dist with disk --- save_tcool_mass_dist.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/save_tcool_mass_dist.py b/save_tcool_mass_dist.py index 34f623e..f2ace00 100644 --- a/save_tcool_mass_dist.py +++ b/save_tcool_mass_dist.py @@ -15,7 +15,7 @@ for my_storage, ds in datasets.piter(dynamic=False, storage=storage): sph = ds.sphere('c',(206,'kpc')) dsk = ds.disk([0.5,0.5,0.5], [0,0,1], (20,'kpc'), (1.3, 'kpc')) - cgm = sph - dsk + cgm = sph# - dsk cgm.set_field_parameter('center', ds.arr([0.5,0.5,0.5], 'code_length')) prof = yt.create_profile(cgm, "cooling_time", "cell_mass", @@ -31,6 +31,6 @@ for i in range(len(datasets)): data_arr[:,i+1] = storage[i] - np.savetxt("tcool_mass_dist_CGM.txt", data_arr, + np.savetxt("tcool_mass_dist_CGM-disk.txt", data_arr, header="tcool [Gyr] DD0000 mass [Msun] DD0001 mass ... etc")