ClusterMonteCarlo / CMC-COSMIC

14 stars 19 forks source link

Convert old fits to new HDF5 so you can test using the latest code compiled without fits #24

Closed scottcoughlin2014 closed 3 years ago

scottcoughlin2014 commented 3 years ago

Install latest COSMIC to have these features

pip install git+https://github.com/COSMIC-PopSynth/COSMIC.git

from cosmic.sample import InitialCMCTable
Singles, Binaries = InitialCMCTable.read("input.fits")
Binaries = Binaries.iloc[1:]
Singles = Singles.iloc[1:-1]
Singles = InitialCMCTable(Singles)
Binaries = InitialCMCTable(Binaries)
Singles.mass_of_cluster = 1.0
Singles.scaled_to_nbody_units = True
Singles.metallicity = 0.001
InitialCMCTable.write(Singles, Binaries, filename="input_test.hdf5")

carlrodriguez commented 3 years ago

Don't think this is needed anymore; can reopen if it is