ind2 = np.where((h_scale[0] <= tn2) & (tn2 <= h_scale[1]))[0]