function [newHandles] = plotModisImageRefine(handles)

caxis('auto')
visible=[1:19,26];

drawnow
figure(handles.modis_figure)
axes(handles.modis_axes)
cla
while 0
h = findobj(handles.modis_figure, 'Type', 'axes');
for n = 1:length(h)
  axes(h(n))
  cla
end
end

% Get map keywords
if strcmp(handles.modis.select_name, 'Classification Mask')
  keywords = {'mercator', 'colorbar', 'magentamap', 'cmaskcolormap', 'overwrite'};
  imageTitle = ['MODIS Classification Mask'];
elseif strcmp(handles.modis.select_name, 'Cloud Phase')
  keywords = {'mercator', 'colorbar', 'magentamap', 'cphasecolormap', 'overwrite'};
  imageTitle = ['MODIS Cloud Phase Mask'];
elseif strcmp(handles.modis.select_name, 'Total Precipitable Water')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap', 'overwrite'};
  imageTitle = ['MODIS Total Precipitable Water'];
elseif strcmp(handles.modis.select_name, 'MODIS Level1b') & handles.modis.band == -1
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap', 'overwrite'};
  imageTitle = ['MODIS Level1b ' handles.modis.select_var];
elseif strcmp(handles.modis.select_name, 'MODIS Level1b') & ...
       ismember(handles.modis.band,visible)
  keywords = {'mercator', 'magentamap', 'colorbar', 'overwrite'};
  imageTitle = ['MODIS Channel ' num2str(handles.modis.band) ' Reflectance'];
elseif strcmp(handles.modis.select_name, 'MODIS Level1b') & ...
       ~ismember(handles.modis.band,visible)
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap', 'overwrite'};
  imageTitle = ['MODIS Channel ' num2str(handles.modis.band) ' Brightness Temperature'];
elseif strcmp(handles.modis.select_name, 'MODIS Cloud Mask')
  keywords = {'mercator', 'magentamap', 'colorbar', 'cloudmaskcolormap', 'overwrite'};
  imageTitle = ['MODIS Cloud Mask'];
elseif strcmp(handles.modis.select_name, 'Cloud Top Pressure')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap', 'overwrite'};
  imageTitle = ['MODIS Cloud Top Pressure'];
elseif strcmp(handles.modis.select_name, 'Water Path')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['MODIS Water Path'];
elseif strcmp(handles.modis.select_name, 'Effective Particle Radius')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['MODIS Effective Particle Radius'];
elseif strcmp(handles.modis.select_name, 'Cloud Optical Thickness')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['MODIS Cloud Optical Thickness'];
elseif strcmp(handles.modis.select_name, 'Cloud Effective Emissivity')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['MODIS Cloud Effective Emissivity'];
elseif strcmp(handles.modis.select_name, 'Layered Cloud Mask')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['AIRS Layered Cloud Mask from MODIS'];
elseif strcmp(handles.modis.select_name, 'Layered Cloud Phase')
  keywords = {'mercator', 'magentamap', 'colorbar', 'jetcolormap'};
  imageTitle = ['AIRS Layered Cloud Phase from MODIS'];
end

for n = 1:length(handles.modis.refine)
  handles.modis_axes = plotSingleImageAxes(handles.modis.refine(n).sub_data, ...
      handles.modis.refine(n).sub_lat, handles.modis.refine(n).sub_lon, ...
      handles.modis_axes, keywords, handles.modis.sub_geo);
end

if strcmp(handles.modis.select_name, 'Total Precipitable Water')
  caxis([0 70])
  colorbar('horizontal')
end

title(imageTitle)

disp(sprintf('Plot %s', handles.modis.select_name))

newHandles = handles;
