import numpy as np


def zero_crossing(eod)
threshold = 0;
shift_eod = np.roll(eod, 1)
eod_times = time[(eod >= threshold) & (shift_eod < threshold)]
return eod_times