load("full_timeseries.mat");
[ww2023_rm] = readncTHOON('THWW01-d03_20230830-20231020_rm_profiles.nc',5,1,0);
[ww2023_rc] = readncTHOON('THWW01-d02_20230725-20230825_rc_profiles.nc',5,1,0);
[ww2023_t] = readncTHOON('THWW01-d01_20230516-20230720_t_profiles.nc',5,1,0);
%%
figure(1)
clf

scatter(ww2023_t.Dissolved_Oxygen(:),ww2023_t.Sigma0(:),'.')
hold on
scatter(ww2023_rm.Dissolved_Oxygen(:),ww2023_rm.Sigma0(:),'.')
scatter(ww2023_rc.Dissolved_Oxygen(:),ww2023_rc.Sigma0(:),'.')

scatter(wwf.Dissolved_Oxygen(:),wwf.Sigma0(:),'.')
grid on
xlabel('doxy [\mumol/kg]')
ylabel('sigma [kg/m^3]')

legend( '2023 t','2023 rm','2023 rc','2024')

%%
figure(2)
clf

scatter(ww2023_t.Practical_Salinity(:),ww2023_t.Temperature(:),'.')
hold on
scatter(ww2023_rm.Practical_Salinity(:),ww2023_rm.Temperature(:),'.')
% scatter(ww2023_rc.Absolute_Salinity(:),ww2023_rc.Conservative_Temperature(:),'.')

scatter(wwf.Practical_Salinity(:),wwf.Temperature(:),'.','MarkerEdgeColor',[0.4940, 0.1840, 0.5560], 'MarkerFaceColor', [0.4940, 0.1840, 0.5560])
grid on
ylabel('Temperature [^oC]')
xlabel('Salinity [PSU]')

legend( '2023 t','2023 rm','2024')
%%
% figure(3)
% clf
% 
% scatter(ww2023_rc.Absolute_Salinity(:),ww2023_rc.Conservative_Temperature(:),'.')
% hold on
% scatter(wwf.(:),wwf.Temperature(:),'.','MarkerEdgeColor',[0.4940, 0.1840, 0.5560], 'MarkerFaceColor', [0.4940, 0.1840, 0.5560])
% grid on
% ylabel('Temperature [^oC]')
% xlabel('Salinity [PSU]')
% 
% legend( '2023 t','2023 rm','2024')