Symptom is when trying to run merge_scans: ``` snrsel = np.where( (np.abs(prop[:, snrcol]) > snrmin) & (np.abs(prop[:, snrcol]) < snrmax) ) ``` The prop is an empty array of 1 dimension. Should be 2d empty?