function [data] = select_wn(displayAxis)
% Select a BT region from the AIRS BT plot
% Input: axis of displayed data
% Output: data structure indicated wavenumber limit

displayAxis;

% Do the rrbox command
set(gcf,'Pointer','Crosshair');
waitforbuttonpress;
Pt1 = get(displayAxis,'currentpoint');
rbbox([get(gcf,'currentpoint') 0 0],get(gcf,'currentpoint'));
Pt2 = get(displayAxis,'currentpoint');
set(gcf,'Pointer','arrow');

% Get selected rectangle on the project map
rect = [min([Pt1(1,1:2);Pt2(1,1:2)]) abs(Pt2(1,1:2)-Pt1(1,1:2))];
data.rect = rect;
data.x_limit = [rect(1) rect(1) + rect(3)];
data.y_limit = [rect(2) rect(2) + rect(4)];
