% plt_map_ctd.m - plot sampling runs and injection

mat_path='..//process/';

clear h

figure(1)
clf

Lat_min=29.5;
Lat_max=30.5;
Lon_min=-70.5;
Lon_max=-69;


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

hold on

%%%%%%%%%%%%%%% CTD
for cast = [1:3 5:21]
  disp(int2str(cast))
  if cast <10
    fn = ['oc41540',int2str(cast)];
  elseif cast <100
    fn = ['oc4154',int2str(cast)];
  end

  %single casts
  eval(['load ',mat_path,fn])
  readhdr;
  m_plot(-lon,lat,'.');
  m_text(-lon,lat,[int2str(cast)]);
  if 0
    figure(2)
    plot(lon,lat,'.')
    if ~ishold, hold on; end
    figure(1)
  end
  
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-4 - CTD stations')

orient portrait
print -dpng ctd_map.png
