function map_WHOTS_2022(lon,lat,LatLonLim,dstr)
%{
  Makes a map of the ship track using the m_map package.
  Inputs: 
	lon,        longitude (dec degrees, 0-360)
	lat,        latitude (dec degrees)
	LatLonLim   Map limits (dec degrees)
    dstr        date string for plot title
%}

disp('plotting track map');
width = (LatLonLim(2)-LatLonLim(1))*350;
height = (LatLonLim(4)-LatLonLim(3))*350;
maxh = max(width,height);
 if maxh>700 && height > width
     width = width*700/height;
     height = 700;
 elseif maxh > 700 && width > height
     height = height *700/width;
     width = 700;
 end

%from NGDC coastline extractor
%load([Fluxroot, 'flux\m_map\US_westcoast.dat'])
figure('position',[300,300,width,height]);
title(['Cruise track, WHOTS, ',dstr]);
m_proj('mercator','lon',[LatLonLim(1) LatLonLim(2)],'lat',[LatLonLim(3) LatLonLim(4)]);
% m_line(US_westcoast(:,1),US_westcoast(:,2),'linewidth',2);
% patch(US_westcoast(:,1),US_westcoast(:,2),[.7 .7 .7],'edgecolor','none');
%m_grid;
m_grid('box','fancy','tickdir','in');
m_gshhs_f('patch',[.5 .5 .5]);

%m_coast('linewidth',2,'color',0.5+[0 0 0]);
%m_coast('patch',[.7 .7 .7],'edgecolor','none');

% % m_line(-85,-20,'marker','x','markersize',8,'color','b');
% % m_line(-85,-20,'marker','o','markersize',8,'color','b');
% % m_line(-80,-27,'marker','square','markersize',7,'color','k');
% % m_line(-75,-20,'marker','diamond','markersize',7,'color','k');

m_line(lon,lat,'color','r','linestyle','-','linewidth',2);
%'marker','.','markersize',6);

end
