function map_FATIMA_2022(lon,lat,LatLonLim)
% Plots map of the ship track using the m_map package and
% gshhs high resolution coastline data
% Inputs: 
%   lon,        track longitude (dec degrees)
%   lat,        track latitude (dec degrees)
%   LatLonLim   Map limits: [lon_min,lon_max,lat_min,lat_max](dec degrees)
%   ddd         Day-of-Year

width = (LatLonLim(2)-LatLonLim(1))*200;
height = (LatLonLim(4)-LatLonLim(3))*200;
maxh = max(width,height);
if maxh>1000 && height>width
    width = width*1000/height;
    height = 1000;
elseif maxh>1000 && width>height
    height = height*1000/width;
    width = 1000;
end
    
% from NGDC coastline extractor
figure('position',[300,300,width,height]);
title(['Cruise Track, FATIMA 2022']);
m_proj('mercator','lon',[LatLonLim(1) LatLonLim(2)],'lat',[LatLonLim(3) LatLonLim(4)]);
m_grid('box','fancy','tickdir','in');
m_gshhs_i('patch',[.5 .5 .5]);
m_line(lon,lat,'color','r','linestyle','-','linewidth',2);
