plot_map.m

function plot_map(hTsgGUI, hPlotAxes)
% Function to plot the earth map and ship trackline
%
% Input
% -----
% hTsgGUI ............ Handel to the main user interface
% hPlotAxes .......... Handels to the graphic axes
%
% Output
% ------
%
% Library : 'M_MAP'
%
% TODOS
%
% note that a decimal degree notation is used, so that a longitude of
% 120?30'W is specified as -120.5
%
% PROBLEM NOT YET RESOLVED
%    Does M_MAP use 0-360? or -180? +180? covention ?
%
% $Id$

% Retrieve named application data
% -------------------------------
tsg = getappdata( hTsgGUI, 'tsg_data');

% save pointer and change it to watch
% -----------------------------------
pointer = get(hTsgGUI, 'pointer');
set(hTsgGUI, 'pointer', 'watch');
set(get(hPlotAxes(4), 'parent'), 'pointer', 'watch');

% Retrieve map border
% --------------------
border = tsg.preference.map_border;

% Retrieve map resolution 
% set 'checked' property to 'on' for current map resolution if user
% change it from menu option/preferences
% ------------------------------------------------------------
hdl = flipud(findobj( '-regexp', 'tag', 'TAG_UIMENU_MAP_RESOLUTION_'));
set(hdl, 'Checked', 'off');  % set all menu off
set(hdl(tsg.preference.map_resolution), 'checked', 'on');

% Get the Geographic limit of the TSG time series
% -----------------------------------------------
dateLim = get(hPlotAxes(1), 'Xlim');
ind = find( tsg.DAYD >= dateLim(1) & tsg.DAYD <= dateLim(2));

% m_grid need the use of double instead single
% --------------------------------------------
latx = double(tsg.LATX);
lonx = double(tsg.LONX);

if ~isempty( ind )
  
  latMin = min(latx(ind) );
  if latMin < -90
    latMin = -90;
  end
  latMax = max( latx(ind) );
  if latMax > 90
    latMax = 90;
  end
  lonMin = min( lonx(ind) );
  lonMax = max( lonx(ind) );
  
  % Oversize window due to the large frame
  %--------------------------------------
  latRange = (latMax-latMin);
  latMin = max(floor(latMin), -90);
  latMax = min(ceil(latMax), 90);
  
  lonRange = (lonMax-lonMin);
  lonMin = floor(lonMin);
  lonMax = ceil(lonMax);
  
  lonRange = (lonMax-lonMin);
  if lonRange>=360
    lonMin = -183.6;     %to account for fancy frame
    lonMax = 183.6;
    tsg.lonplus = 180;
    tsg.lonmod = 360;
  else
    tsg.lonplus = 0;
    tsg.lonmod = 0;
  end
  
  % Positionning the right axes (set map current axe)
  % -------------------------------------------------
  axes(hPlotAxes(4));
  
  % Use of Mercator projection
  % --------------------------
  m_proj('Mercator','lat',[latMin-border latMax+border],'long',[lonMin-border lonMax+border]);
  
  % plot climatology mean contour if available
  % ------------------------------------------
  hdl = findobj( '-regexp', 'tag', 'PUSHTOOL_CLIM');
  if strcmp(get(hdl, 'state'), 'on')
    climato = tsg.levitus.data.(['WOA_MEAN_' tsg.plot.sample]);
    % get last selected climatology structure
    s = get(findobj('Tag', 'TAG_UIMENU_CLIMATO_MAIN'), 'Userdata');
    level = tsg.preference.levitus_depth_value;
    % extract the matrix
    climato = squeeze(climato(s.time, level,:,:));
    % use local var for quick access
    LATX = tsg.levitus.data.WOA_LATX;
    LONX = tsg.levitus.data.WOA_LONX;
    
    % swap negative value (west) to 180-360
    if ~isempty(tsg.indcross)
      LONX(LONX < 0) = LONX(LONX < 0) + 360;
      [val, indval] = sort(LONX);
      LONX          = LONX(indval);
      climato       = climato(:,indval);
    end
    
    % find  indices for the selected area
    indLat = find(LATX > latMin-border & LATX <latMax+border);
    indLon = find(LONX > lonMin-border & LONX <lonMax+border);
    lon = LONX(indLon);
    lat = LATX(indLat);
    
    % plot climatogogy
    m_pcolor(lon, lat, climato(indLat, indLon));
    shading flat;
    %colormap(m_colmap('jet','step',10));
    colorbar;
  else
    colorbar('off');
  end
  
  % select map type and resolution
  % ------------------------------
  switch tsg.preference.map_resolution
    
    case 1
      % Low-resolution coast lines,
      % --------------------------
      if tsg.preference.map_patch_value == 2
        m_coast('patch',[.7 .7 .7], 'TAG', 'TAG_PLOT4_LINE_COAST');
      else
        s = load('world_c.mat');
      end
      
    case 2
      % Medium-resolution coast lines
      % -----------------------------
      if tsg.preference.map_patch_value  == 2
        m_gshhs_l('patch',[.7 .7 .7], 'TAG', 'TAG_PLOT4_LINE_COAST');
      else
        s = load('world_l.mat');
      end
      
      
    case 3
      % Intermediate-resolution coast lines
      % ----------------------------
      if tsg.preference.map_patch_value  == 2
        m_gshhs_i('patch',[.7 .7 .7], 'TAG', 'TAG_PLOT4_LINE_COAST');
      else
        s = load('world_i.mat');
      end
      
      
    case 4
      % High-resolution coast lines
      % ----------------------------
      if tsg.preference.map_patch_value  == 2
        m_gshhs_h('patch',[.7 .7 .7], 'TAG', 'TAG_PLOT4_LINE_COAST');
      else
        s = load('world_h.mat');
      end
      
      
    otherwise
      tsg.preference.map_resolution = 1;
      % Low-resolution coast lines
      % --------------------------
      if tsg.preference.map_patch_value  == 2
        m_coast('patch',[.7 .7 .7], 'TAG', 'TAG_PLOT4_LINE_COAST');
      else
        s = load('world_c.mat');
      end
      
  end
  
  
  % Make a grid on the map with fancy box
  % -------------------------------------
  m_grid('box', 'fancy', 'tickdir', 'in', 'TAG', 'TAG_PLOT4_LINE_GRID', ...
    'Fontsize', tsg.fontSize);
  
  
  % use 'world_x.mat' as tsg.preference.map_patch_string == 'off'
  % --------------------------------------------------------------
  if tsg.preference.map_patch_value == 1
    x = s.ncst(:,1);
    y = s.ncst(:,2);
    % swap negative value (west) to 180-360
    if ~isempty(tsg.indcross)
      x(x < 0) = x(x < 0) + 360;
      % remove line horizontal lines
      x(x == 0) = nan;
    end
    
    % plot coast
    m_line(x,y, 'color','k');
  end
  
  
  
  % Plot the ship track - Check if LAT-LON in deg-decimal
  % -----------------------------------------------------
  %m_line(lonx, latx, 'TAG', 'TAG_PLOT4_LINE_SHIPTRACK', 'marker','.','markersize',3,'color','b');
  %   m_line(lonx, latx, 'LineStyle', 'none', 'marker','*','markersize',2,'color','b');
  
  % Plot using QC
  % -------------
  
  % Get list of keys from hashtable tsg.qc.hash, defined inside
  % tsg_initialisation.m
  % -----------------------------------------------------------
  qc_list = keys(tsg.qc.hash);
  
  % Plot Sample/TSG differences on axe 2
  % iterate (loop) on each key store inside hastable
  % Iterate from the end of the cell array : fliplr
  % ------------------------------------------------
  for key = fliplr( qc_list )
    
    % get  some values in hashtable
    % -----------------------------
    %qcState = tsg.qc.hash.(key).state;
    qcCode  = tsg.qc.hash.(key).code;
    qcColor = tsg.qc.hash.(key).color;
    
    % plot tsg salinity sample with right code/color
    % ----------------------------------------------
    ind = find( tsg.SSPS_QC == qcCode & isnan(tsg.SSPS) == 0);
    if ~isempty( ind )
      
      % Create line on the map
      % ----------------------
      m_line( mod(lonx(ind) + tsg.lonplus, tsg.lonmod) - tsg.lonplus,...
        latx(ind), 'LineStyle', 'none', 'marker','*',...
        'markersize', 2, 'color', qcColor);
    end
  end
  
  % Write title
  % -----------
  title( ['Cruise: ' tsg.file.name], 'fontsize',tsg.fontSize +2, ...
    'fontweight', 'bold','interpreter','none');
  
  % Set tag for each line
  % ---------------------
  hLines = get(hPlotAxes(4), 'Children');
  for i=1:length(hLines)
    set(hLines(i), 'Tag', ['TAG_PLOT4_LINE_MAP' int2str(i)]);
  end
  
  % Save tsg structure
  % ------------------
  setappdata( hTsgGUI, 'tsg_data', tsg);
  
  % restore pointer
  % ---------------
  set(hTsgGUI, 'pointer', pointer);
  set(get(hPlotAxes(4), 'parent'), 'pointer', 'arrow');
  
end