Add split_multiselected_dims = T in obs call. Add checks for dataset consistency