From 9e38e1cd7889c7e7b14c9a5fc9f215a864d7c821 Mon Sep 17 00:00:00 2001 From: Stephen Neilson <36410751+s-neilson@users.noreply.github.com> Date: Tue, 14 Jul 2026 14:22:11 +1000 Subject: [PATCH 1/2] In read_dumpfiles, removed the use of the pandas.unique as it no longer supports using a Python list as an input in versions of pandas 3.0.0 or higher . It is now changed to use a inbuilt Python set instead; the output dumpfile list is also now guaranteed to always be in the same alphanumerical order in case this is useful in the future. --- src/utils/rdumpfiles.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/src/utils/rdumpfiles.py b/src/utils/rdumpfiles.py index 88aec92..9de891e 100644 --- a/src/utils/rdumpfiles.py +++ b/src/utils/rdumpfiles.py @@ -2,7 +2,6 @@ # -*- coding: utf-8 -*- import os -from pandas import unique def read_dumpfiles(files_prefix='binary_', path='./',evy_files=-1): ''' @@ -28,7 +27,8 @@ def read_dumpfiles(files_prefix='binary_', path='./',evy_files=-1): openlist = open('listfiles.txt') list_dumpfiles = openlist.read().splitlines() list_dumpfiles = [namefile.rstrip('.divv') for namefile in list_dumpfiles] - list_dumpfiles = list(unique(list_dumpfiles)) + list_dumpfiles = list(set(list_dumpfiles)) + list_dumpfiles.sort() openlist.close() os.system('rm listfiles.txt') if (evy_files == -1) or (evy_files == None): @@ -36,4 +36,4 @@ def read_dumpfiles(files_prefix='binary_', path='./',evy_files=-1): elif evy_files > 0: return list_dumpfiles[::evy_files] else: - print("Invalid option for evy_files.") \ No newline at end of file + print("Invalid option for evy_files.") From b6beee8903f340ba55d74443b14d470ec9378354 Mon Sep 17 00:00:00 2001 From: Stephen Neilson <36410751+s-neilson@users.noreply.github.com> Date: Thu, 16 Jul 2026 14:57:04 +1000 Subject: [PATCH 2/2] For newly unbound mass plots, all sink separations, not just the first, will now be plotted. --- src/newly_unbound.py | 6 +++++- 1 file changed, 5 insertions(+), 1 deletion(-) diff --git a/src/newly_unbound.py b/src/newly_unbound.py index dfd919b..45ca605 100644 --- a/src/newly_unbound.py +++ b/src/newly_unbound.py @@ -91,7 +91,11 @@ def unbound_mech_kipp(dumpfile_list): #Plot declaration ax = plt.gca() - ax.plot(ph_data['time']*yr, np.log10(ph_data['sep. 1'])) + + #Finds the total separation columns in the separation .ev file and plots them. + for columnName in ph_data.keys(): + if(columnName[0:5] == "sep. "): + ax.plot(ph_data['time']*yr, np.log10(ph_data[columnName])) #Histogram format nbins = min([int(len(dump_list)/1.05),500])