Commit 69320935 authored by Médéric Boquien's avatar Médéric Boquien
Browse files

Do not attempt to save the SFH if we only fit the dust.

parent b7f66542
......@@ -38,7 +38,8 @@ def save_sed_to_fits(sed, prefix, norm=1.):
unit="W/nm")
table.write("{}_best_model.fits".format(prefix))
table = Table(meta=info)
table["time"] = Column(np.arange(sed.sfh.size), unit="Myr")
table["SFR"] = Column(norm * sed.sfh, unit="Msun/yr")
table.write("{}_SFH.fits".format(prefix))
if sed.sfh is not None:
table = Table(meta=info)
table["time"] = Column(np.arange(sed.sfh.size), unit="Myr")
table["SFR"] = Column(norm * sed.sfh, unit="Msun/yr")
table.write("{}_SFH.fits".format(prefix))
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment