% so if there's a map of our glider transects along line 67 (can be just 
% the transect ones, not the ones only in monterey bay), current line 90, 
% deployment from North Carolina, and LOCNESS

%% load data
load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\21603401\QC\21603401_QC.mat","s")
g21603401 = s;

load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\24403401\QC\24403401_QC.mat","s")
g24403401 = s;

load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25820301\real_time_binned\25820301sat.mat","s")
g25820301 = s;

load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25420901\real_time_binned\25420901sat.mat","s")
g25420901 = s;

load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25720901\real_time_binned\25720901sat.mat","s")
g25720901 = s;

load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25821001\real_time_binned\25821001sat.mat","s")
g25821001 = s;

%%
GliderMap = figure;
clf
Maptitle = 'Glider pH Profiles';

%% Create a geographic axes
ax = geoaxes('Basemap','satellite',Position=[0, 0, 1, .95]);
ax.Title.String = Maptitle;
ax.Title.Visible = "on";
%% Hold the current plot
hold(ax, 'on');
%% Plot Spray g21603401 Location History
SprayLatHist = g21603401.lat;
SprayLonHist = g21603401.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Plot Spray g24403401 Location History
SprayLatHist = g24403401.lat;
SprayLonHist = g24403401.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Plot Spray g25820301 Location History
SprayLatHist = g25820301.lat;
SprayLonHist = g25820301.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Plot Spray g25420901 Location History
SprayLatHist = g25420901.lat;
SprayLonHist = g25420901.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Plot Spray g25720901 Location History
SprayLatHist = g25720901.lat;
SprayLonHist = g25720901.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Plot Spray g25821001 Location History
SprayLatHist = g25821001.lat;
SprayLonHist = g25821001.lon;
geoplot(ax, SprayLatHist, SprayLonHist, '.', 'Color', '#D3D3D3', 'MarkerSize', 10);
%% Add a legend
% legend(ax, 'Location', 'northwest');

