% plt_map_xbt.m - plot sampling runs and injection

clear h

figure(1)
clf

Lat_min=30+15/60;
Lat_max=31+30/60;
Lon_min=-67.333;
Lon_max=-65.66;

m_proj('mercator','lat',[Lat_min Lat_max],'lon',[Lon_min Lon_max])

hold on


latlon= ...
    [31+15.21899/60  -65- 50.42285/60
     31+8.62891/60    -65- 56.49902/60
     31+0.51611 /60   -66- 7.95898/60
     30+54.23706/60   -66- 18.396/60
     30+47.77002/60   -66- 29.19092/60
     30+41.51294/60   -66- 39.604/60
     30+34.98096/60  -66- 50.45801/60
     30+29.46191/60   -67- 1.65479/60
     30+28.11108/60   -67- 12.09424/60
     30+38.86401/60  -67- 12.896/60
     30+44.88501/60  -67- 6.40479/60
     30+40.46899/60  -66- 59.38086/60
     30+31.37012/60  -66- 59.49707/60];
%%%%%%%%%%%%%%% XBT
for cast = 1:13
  
  m_plot(latlon(cast,2),latlon(cast,1),'.');
  m_text(latlon(cast,2),latlon(cast,1),[int2str(cast)]);

end

m_grid('box','on','tickdir','in','fontweight','bold'); 

% this makes the color reappear
set(findobj('tag','m_grid_color'),'facecolor','none')

xlabel('Longitude')
ylabel('Latitude')
title('EDDIES OC415-2 - XBT stations')

orient portrait
print -dpng xbt_map.png
