%% Compare 212 and 203
s203 = load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25820301\25820301sat.mat");
s212 = load("\\atlas.shore.mbari.org\ProjectLibrary\901805_Coastal_Biogeochemical_Sensing\Spray_Data\25921201\25921201sat.mat");
dnum_203 = 4;
dnum_212 = 5;
fig = figure;
clf
tiledlayout(1,2)
ax1 = nexttile;
plot(s203.s.pHin(:,dnum_203)-s203.s.pHin_canb(:,dnum_203),s203.s.depth(:,dnum_203))
set(gca,'YDir','reverse')
title(['SN203 Dive ',num2str(dnum_203)])
xlabel('pHin-canB')
grid on;

ax2 = nexttile;
plot(s212.s.pHin(:,dnum_212)-s212.s.pHin_canb(:,dnum_212),s212.s.depth(:,dnum_212))
set(gca,'YDir','reverse')
title(['SN212 Dive ',num2str(dnum_212)])
xlabel('pHin-canB')
grid on;

linkaxes([ax1,ax2],'x')
%%
canBdelta203 = [s203.s.pHin(:,dnum_203)-s203.s.pHin_canb(:,dnum_203), s203.s.depth(:,dnum_203)];
canBdelta212 = [s212.s.pHin(:,dnum_212)-s212.s.pHin_canb(:,dnum_212), s212.s.depth(:,dnum_212)];


SN203_25_350 = 0.0859875644732888 - 0.0277217158180880
SN212_25_350 = -1.00202826433637 - -1.05354684093830