k1, = np.where((di_low < di) & (di < di_up))