phunk¶
A python package to fit photometric phase curves of asteroids.[1]
import phunk
Observe the phase curve of your favourite asteroids.
# Observations of (20) Massalia from Gehrels 1956
phase = [0.57, 1.09, 3.20, 10.99, 14.69, 20.42]
mag = [6.555, 6.646, 6.793, 7.130, 7.210, 7.414]
epoch = [35193, 35194, 35198, 35214, 35223, 35242]
Load the observations into phunk.
pc = phunk.PhaseCurve(phase=phase, mag=mag, epoch=epoch, target="Massalia")
Fit your favourite model(s) to the observations.
pc.fit(["HG", "HG1G2", "sHG1G2"])
pc.plot()
Footnotes