function diff_modis1b_scatter_callback

handles = guidata(gcf);
s = get(gcbo,'Label');
if strcmp(s, 'BT11-BT12')
  band1 = 31;
  band2 = 32;
elseif strcmp(s, 'BT8.6-BT11')
  band1 = 29;
  band2 = 31;
elseif strcmp(s, 'BT11-BT6.7')
  band1 = 31;
  band2 = 27;
elseif strcmp(s, 'BT3.9-BT3.7')
  band1 = 21;
  band2 = 20;
elseif strcmp(s, 'BT11-BT3.7')
  band1 = 31;
  band2 = 20;
elseif strcmp(s, 'BT12-BT4')
  band1 = 32;
  band2 = 23;
elseif strcmp(s, 'BT13.7-BT14')
  band1 = 34;
  band2 = 35;
elseif strcmp(s, 'BT11-BT3.9')
  band1 = 31;
  band2 = 21;
end

handles.modis.select_name = 'MODIS Level1b';
handles.modis.select_var = s;
  
% Read MODIS level1b data
x_data = read_mod1b(handles.modis.filename.level1b, band1);
y_data = read_mod1b(handles.modis.filename.level1b, band2);

switch handles.modis.axis_variable
  case 'x'
    handles.modis.x_label = [s ' (K) '];
    handles.modis.x_data = x_data - y_data;
    handles.modis.x_band = -1;
  case 'y'
    handles.modis.y_label = [s ' (K) '];
    handles.modis.y_data = x_data - y_data;
    handles.modis.y_band = -1;
end

handles = get_modis_scatter(handles);
scatter_plot(handles)

% Update handles structure
guidata(gcbo, handles);
