Introduced flux subtraction to time series so that flux median is 0.001
i.e. the mean of the Channel 4 flux from Jan 2010 to mid 2014.

Also replace time and flux variables after removal of artifacts with
clean_time and clean_flux.
Began writing two functions: lyra_event_list() to produce a LYRA
event list with just start and end times as inputs; and
get_lyra_timeseries() to create LYRA timeseries in a recarray,
downloading necessary files.
Wrote first draft of lyra_event_list().

Deleted get_lyra_timeseries.  Can use LYRALightCurve instead.

find_lyra_events():
Added "Recovery" to list of artifacts to be removed.

Altered way flare end condition is searched for.  Before peak was defined
as first time derivative went negative after start time.  Now this defines
the time from which end time begins being search for.  Each subsquent point
is checked to see if is satisfies the end condition where the peak flux is
defined as the maximum flux between the start time and the time being
currently check.  This has improved how end times are found and stopped some
"double flares" being found.

Reversed order of flux and time in definition of function.

_check_datetime():
Added facility to convert pandas time series into array of datetimes.

Instead of changing imput time, a copy of this variable is made, new_time,
and this is altered.

extract_combined_lytaf():
Changed default lytaf_path to be defined from hom directory.
Updated local LYRA annotation files.
Fixed bug where code crashed if flare start time was found but
flux did not decrease before end of time series.
Changed RISE_FACTOR threshold to 1.01 from 1.005.
Added capability to extract_combined_lytaf() to output a csv file
through a new csvfile kwarg.
Included start, peak and end fluxes in output recarray of
find_lyra_events().  Also removed calculation of flare_indices.
Edited doctrings.  Also changed start and end time variable names
in extract_combined_lytaf().