def timeseq(time, covering, refpt):
    """Extract time sequence from a covering"""
    ts = time.take(covering[refpt][3], 0)
    print("time sequence min = %f, max = %f"%(min(ts),max(ts)))
    return ts