Merge pull request #32 from jeonggyukim/core_formation

Core formation
jeonggyukim · Aug 21, 2024 · 9d7df15 · 9d7df15
2 parents c37015f + 653799c
commit 9d7df15
Show file tree

Hide file tree

Showing 19 changed files with 3,554 additions and 1,767 deletions.
diff --git a/pyathena/core_formation/.gitignore b/pyathena/core_formation/.gitignore
@@ -0,0 +1,3 @@
+*.p
+visit_script.py
+visitlog.py
diff --git a/pyathena/core_formation/config.py b/pyathena/core_formation/config.py
@@ -1,5 +1,6 @@
-GRID_NUM_START = 30
-PLOT_PREFIX_TCOLL_CORES = "tcoll_cores"
+GRID_NUM_START = 1
+PLOT_PREFIX_CORE_EVOLUTION = "core_evolution"
 PLOT_PREFIX_SINK_HISTORY = "sink_history"
 PLOT_PREFIX_PDF_PSPEC = "PDF_Pspecs"
 PLOT_PREFIX_RHOC_EVOLUTION = "central_density_evolution"
+PLOT_PREFIX_MASS_RADIUS= "mass_radius"
diff --git a/pyathena/core_formation/do_tasks.py b/pyathena/core_formation/do_tasks.py
diff --git a/pyathena/core_formation/hst.py b/pyathena/core_formation/hst.py
@@ -20,10 +20,10 @@ def read_hst(self, savdir=None, force_override=False):
         # Time in code unit
         h['time'] = hst['time']
         # Time in unit of free-fall time
-        h['time_in_tff'] = h['time']/self.tff
+        h['time_in_tff'] = h['time']/self.tff0
         # Timestep
         h['dt'] = hst['dt']
-        h['dt_in_tff'] = h['dt']/self.tff
+        h['dt_in_tff'] = h['dt']/self.tff0
 
         # Total gas mass
         h['mass'] = hst['mass']