py code to extract lightcurve from fits in 2013