Testing tcool dist with disk

clairekope · May 18, 2022 · 9817d40 · 9817d40
1 parent 37e4ff6
commit 9817d40
Showing 1 changed file with 2 additions and 2 deletions.
diff --git 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")