Currently the code detects eclipse boundaries by finding the phases where the flux returns to one; however, this only works if the light curve has already been normalized to one before using eclipsebin. Add a function to test if the median value of the inputted light curve is near one, and if not, change the eclipse boundary detection value to the median value of the light curve.