Skip to content

Commit 4ffa1e1

Browse files
hsosikhsosik
committed
Gordon's transect segments start times for TN368_transects.m
1 parent 9695a6c commit 4ffa1e1

File tree

2 files changed

+344
-3
lines changed

2 files changed

+344
-3
lines changed

‎Tools/btlmat2table.m

+10-3
Original file line numberDiff line numberDiff line change
@@ -9,9 +9,10 @@
99
%file2load = 'C:\work\SPIROPA\ar29_bottle_data_Apr_2019.mat';
1010
%file2load = 'C:\work\SPIROPA\RB1904\rb1904_bottle_data_Jul_2019.mat';
1111
%file2load = '\\sosiknas1\Lab_data\SPIROPA\20180414_AR29\fromOlga\ar29_bottle_data_Apr_2019.mat';
12-
file2load = '\\sosiknas1\Lab_data\SPIROPA\20180503_RB1904\fromOlga\rb1904_bottle_data_Apr_2020.mat';
13-
file2load = '\\sosiknas1\Lab_data\SPIROPA\20190705_TN368\fromOlga\tn368_bottle_data_Apr_2020.mat';
14-
12+
%file2load = '\\sosiknas1\Lab_data\SPIROPA\20180503_RB1904\fromOlga\rb1904_bottle_data_Apr_2020.mat';
13+
%file2load = '\\sosiknas1\Lab_data\SPIROPA\20190705_TN368\fromOlga\tn368_bottle_data_Apr_2020.mat';
14+
file2load = 'C:\work\SPIROPA\TN368\tn368_bottle_data_Dec_2021.mat';
15+
CTDlist = readtable("C:\work\SPIROPA\TN368\tn368_ctd_list.txt");
1516
if ~exist(file2load,'file')
1617
[FileName,PathName] = uigetfile('*.mat','Select BTL mat file');
1718
file2load = fullfilename(PathName, FileName);
@@ -64,6 +65,12 @@
6465
%Olga's longitudes are missing the negative sign
6566
BTL.Longitude_decimalDeg = -1*BTL.Longitude_decimalDeg;
6667

68+
%fix messed up times in Olga's BTL files
69+
[ii,ia] = ismember(BTL.Cast, CTDlist.Var1);
70+
BTL.hour = CTDlist.Var5(ia); BTL.minute = CTDlist.Var6(ia); BTL.second = CTDlist.Var7(ia);
71+
%
72+
BTL.datetime = datetime(BTL.Year, BTL.Month, BTL.Day, BTL.hour, BTL.minute, BTL.second);
73+
6774
clear ii s t f n
6875

6976
[p,n,e] = fileparts(file2load);

‎transect_plotting/TN368_transects.m

+334
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,334 @@
1+
%S19 = load("\\sosiknas1\IFCB_products\SPIROPA\summary\summary_biovol_allHDF_min20_2019.mat");
2+
S19 = load("c:\work\IFCB_products\SPIROPA\summary\summary_biovol_allHDF_min20_2019.mat");
3+
L19 = load("c:\work\IFCB_products\SPIROPA\summary\summary_biovol_allHDF_min20_2019lists.mat");
4+
5+
TNuwind = find(strcmp(S19.meta_data.cruise, 'TN368') & strcmp(S19.meta_data.sample_type, 'underway') & ~(S19.meta_data.skip));
6+
7+
run("\\sosiknas1\Lab_data\SPIROPA\20190705_TN368\scripts\tn368_transect_startime.m")
8+
trnum = 1:length(tstime);
9+
S19.meta_data.transect(TNuwind) = floor(interp1(tstime,trnum, S19.mdate(TNuwind)));
10+
ii1 = (S19.meta_data.transect(TNuwind)==1);
11+
ii5_7 = (ismember(S19.meta_data.transect(TNuwind),[5 6 7]));
12+
ii8 = (S19.meta_data.transect(TNuwind)==8);
13+
ii3_7 = (ismember(S19.meta_data.transect(TNuwind),[3 4 5 6 7]));
14+
15+
%%
16+
class2use = S19.class2use;
17+
group_table = readtable('\\sosiknas1\training_sets\IFCB\config\IFCB_classlist_type.csv');
18+
[~,ia,ib] = intersect(group_table.CNN_classlist, class2use);
19+
diatom_ind = ib(find(group_table.Diatom(ia)));
20+
dino_ind = ib(find(group_table.Dinoflagellate(ia)));
21+
nano_ind = ib(find(group_table.Nano(ia) | group_table.flagellate(ia) | group_table.Coccolithophore(ia)));
22+
ciliate_ind = ib(find(group_table.Ciliate(ia)));
23+
artifact_ind = ib(find(group_table.IFCBArtifact(ia)));
24+
otherPhyto_ind = [strmatch('Pseudochattonella_farcimen', class2use); strmatch('Phaeocystis', class2use, 'exact')];
25+
phyto_ind = [diatom_ind; dino_ind; nano_ind; otherPhyto_ind];
26+
particle_ind = setdiff(1:length(class2use), artifact_ind);
27+
%%
28+
load \\sosiknas1\Lab_data\Attune\cruise_data\20190705_TN368\bead_calibrated\Attune_uw_match
29+
Attune_mdate = datenum(Attune_uw_match.StartDate_Atable);
30+
31+
%%
32+
phytoCconc = table;
33+
Twin = 12/60/24; %12 minutes as days
34+
for count = 1:length(TNuwind)
35+
temp = array2table(cat(1,L19.classFeaList{TNuwind(count),[phyto_ind]}), 'VariableNames', L19.classFeaList_variables);
36+
ind = (temp.ESD>=10 & temp.score<20);
37+
phytoCconc.IFCB_C10_20(count) = sum(temp.cellC(ind))./S19.meta_data.ml_analyzed(TNuwind(count))/1000;
38+
ind = (temp.ESD>=20 & temp.score<100);
39+
phytoCconc.IFCB_C20_100(count) = sum(temp.cellC(ind))./S19.meta_data.ml_analyzed(TNuwind(count))/1000;
40+
attune_ind = find(Attune_mdate>S19.mdate(TNuwind(count))-Twin & Attune_mdate<S19.mdate(TNuwind(count))+Twin);
41+
phytoCconc.FCM_C0to2(count) = sum(Attune_uw_match.carbon0to2(attune_ind))./sum(Attune_uw_match.VolAnalyzed_ml(attune_ind))/1000;
42+
phytoCconc.FCM_C2to10(count) = sum(Attune_uw_match.carbon2to5(attune_ind)+Attune_uw_match.carbon5to10(attune_ind))./sum(Attune_uw_match.VolAnalyzed_ml(attune_ind))/1000;
43+
phytoCconc.FCM_C10to20(count) = sum(Attune_uw_match.carbon10to20(attune_ind))./sum(Attune_uw_match.VolAnalyzed_ml(attune_ind))/1000;
44+
end
45+
phytoCarbonConcentration_underway.longitude = S19.meta_data.longitude(TNuwind);
46+
phytoCarbonConcentration_underway.latitude = S19.meta_data.latitude(TNuwind);
47+
phytoCarbonConcentration_underway.transect = S19.meta_data.transect(TNuwind);
48+
49+
phytoCarbonConcentration_underway = table;
50+
phytoCarbonConcentration_underway.datetime = datetime(S19.mdate(TNuwind), 'ConvertFrom', 'datenum');
51+
phytoCarbonConcentration_underway.latitude = S19.meta_data.latitude(TNuwind);
52+
phytoCarbonConcentration_underway.longitude = S19.meta_data.longitude(TNuwind);
53+
phytoCarbonConcentration_underway.transect = S19.meta_data.transect(TNuwind);
54+
phytoCarbonConcentration_underway.IFCBfilename = S19.meta_data.pid(TNuwind);
55+
phytoCarbonConcentration_underway.ESD_less_than2micron = phytoCconc.FCM_C0to2;
56+
phytoCarbonConcentration_underway.ESD_2to10micron = phytoCconc.FCM_C2to10;
57+
phytoCarbonConcentration_underway.ESD_10to20micron = phytoCconc.IFCB_C10_20;
58+
phytoCarbonConcentration_underway.ESD_20to100micron = phytoCconc.IFCB_C20_100;
59+
phytoCarbonConcentration_underway.ESD_less_than_100micron = phytoCconc.FCM_C0to2+phytoCconc.FCM_C2to10+phytoCconc.IFCB_C10_20+phytoCconc.IFCB_C20_100;
60+
61+
notes = {'Phytoplankton carbon concentration, units milligrams per cubic meter'; 'Less than 10 microns from Attune; Greater than 10 microns from IFCB'; 'Produced from TN368_transects.m'};
62+
save('\\sosiknas1\IFCB_products\SPIROPA\summary\TN368_phytoCarbonConcentration_underway', 'notes', 'phytoCarbonConcentration_underway')
63+
64+
65+
66+
67+
%%
68+
iA8 = (Attune_uw_match.transect==8);
69+
iA1 = (Attune_uw_match.transect==1);
70+
iA5_7 = (ismember(Attune_uw_match.transect,[5 6 7]));
71+
iA3_7 = (ismember(Attune_uw_match.transect,[3 4 5 6 7]));
72+
73+
%%
74+
figure, set(gcf, 'position', [350 50 550 580])
75+
th = tiledlayout(4,1);
76+
nexttile
77+
yyaxis left
78+
plot(Attune_uw_match.lon_flr(iA8), Attune_uw_match.flr(iA8), '-', 'linewidth', 2)
79+
ylabel('Chl fluor')
80+
yyaxis right
81+
plot(Attune_uw_match.lon_flr(iA8), Attune_uw_match.T(iA8), '-', 'linewidth', 2)
82+
ylabel('Temperature (\circC)')
83+
title(th, 'TN368 Transect 8')
84+
grid on
85+
86+
nexttile
87+
yyaxis left
88+
plot(Attune_uw_match.longitude_fullres(iA8), Attune_uw_match.(" Syn_count ")(iA8)./Attune_uw_match.VolAnalyzed_ml(iA8), '.-', 'linewidth', 2)
89+
ylabel('Syn (ml^{-1})')
90+
ylim([0 inf])
91+
yyaxis right
92+
plot(Attune_uw_match.longitude_fullres(iA8), Attune_uw_match.(" Euk_count ")(iA8)./Attune_uw_match.VolAnalyzed_ml(iA8), '.-', 'linewidth', 2)
93+
ylabel('Picoeukaryotes (ml^{-1})')
94+
grid on
95+
96+
nexttile
97+
yyaxis left
98+
plot(S19.meta_data.longitude(TNuwind(ii8)), sum(S19.classcount(TNuwind(ii8),nano_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii8)), '.-', 'linewidth', 2)
99+
ylabel('Misc. nanoplankton (ml^{-1})')
100+
grid on
101+
102+
nexttile
103+
yyaxis left
104+
plot(S19.meta_data.longitude(TNuwind(ii8)), sum(S19.classcount(TNuwind(ii8),diatom_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii8)), '.-', 'linewidth', 2)
105+
hold on
106+
ylabel('Diatoms (ml^{-1})')
107+
yyaxis right
108+
plot(S19.meta_data.longitude(TNuwind(ii8)), sum(S19.classcount(TNuwind(ii8),dino_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii8)), '.-', 'linewidth', 2)
109+
ylabel('Dinoflagellates (ml^{-1})')
110+
grid on
111+
112+
%%
113+
figure, set(gcf, 'position', [350 50 550 580])
114+
th = tiledlayout(4,1);
115+
nexttile
116+
yyaxis left
117+
plot(Attune_uw_match.lat_flr(iA1), Attune_uw_match.flr(iA1), '-', 'linewidth', 2)
118+
set(gca, 'xdir', 'rev')
119+
ylabel('Chl fluor')
120+
yyaxis right
121+
plot(Attune_uw_match.lat_flr(iA1), Attune_uw_match.T(iA1), '-', 'linewidth', 2)
122+
set(gca, 'xdir', 'rev')
123+
ylabel('Temperature (\circC)')
124+
title(th, 'TN368 Transect 1')
125+
grid on
126+
127+
nexttile
128+
yyaxis left
129+
plot(Attune_uw_match.latitude_fullres(iA1), Attune_uw_match.(" Syn_count ")(iA1)./Attune_uw_match.VolAnalyzed_ml(iA1), '.-', 'linewidth', 2)
130+
set(gca, 'xdir', 'rev')
131+
ylabel('Syn (ml^{-1})')
132+
ylim([0 inf])
133+
yyaxis right
134+
plot(Attune_uw_match.latitude_fullres(iA1), Attune_uw_match.(" Euk_count ")(iA1)./Attune_uw_match.VolAnalyzed_ml(iA1), '.-', 'linewidth', 2)
135+
ylabel('Picoeukaryotes (ml^{-1})')
136+
grid on
137+
138+
nexttile
139+
yyaxis left
140+
plot(S19.meta_data.latitude(TNuwind(ii1)), sum(S19.classcount(TNuwind(ii1),nano_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii1)), '.-', 'linewidth', 2)
141+
ylabel('Misc. nanoplankton (ml^{-1})')
142+
grid on
143+
set(gca, 'xdir', 'rev')
144+
145+
nexttile
146+
yyaxis left
147+
plot(S19.meta_data.latitude(TNuwind(ii1)), sum(S19.classcount(TNuwind(ii1),diatom_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii1)), '.-', 'linewidth', 2)
148+
hold on
149+
ylabel('Diatoms (ml^{-1})')
150+
yyaxis right
151+
plot(S19.meta_data.latitude(TNuwind(ii1)), sum(S19.classcount(TNuwind(ii1),dino_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii1)), '.-', 'linewidth', 2)
152+
ylabel('Dinoflagellates (ml^{-1})')
153+
grid on
154+
set(gca, 'xdir', 'rev')
155+
%%
156+
figure, set(gcf, 'position', [350 50 550 580])
157+
th = tiledlayout(2,4);
158+
nexttile
159+
scatter(Attune_uw_match.lon_flr(iA5_7), Attune_uw_match.lat_flr(iA5_7), 20, Attune_uw_match.T(iA5_7), 'filled')
160+
cb = colorbar('Location', 'northoutside'); title(cb,'Temperature (\circC)')
161+
grid on
162+
nexttile
163+
scatter(Attune_uw_match.lon_flr(iA5_7), Attune_uw_match.lat_flr(iA5_7), 20, Attune_uw_match.flr(iA5_7), 'filled')
164+
cb = colorbar('Location', 'northoutside'); title(cb, 'Chl fluor')
165+
title(th, 'TN368 Transects 5-7')
166+
grid on
167+
168+
nexttile
169+
scatter(Attune_uw_match.longitude_fullres(iA5_7),Attune_uw_match.latitude_fullres(iA5_7), 20, log10(Attune_uw_match.(" Syn_count ")(iA5_7)./Attune_uw_match.VolAnalyzed_ml(iA5_7)), 'filled')
170+
cb = colorbar('Location', 'northoutside'); title(cb, 'log10 Syn (ml^{-1})')
171+
grid on
172+
173+
nexttile
174+
scatter(Attune_uw_match.longitude_fullres(iA5_7), Attune_uw_match.latitude_fullres(iA5_7), 20, log10(Attune_uw_match.(" Euk_count ")(iA5_7)./Attune_uw_match.VolAnalyzed_ml(iA5_7)), 'filled')
175+
cb = colorbar('Location', 'northoutside'); title(cb,'log10 Picoeukaryotes (ml^{-1})')
176+
grid on
177+
178+
nexttile
179+
scatter(S19.meta_data.longitude(TNuwind(ii5_7)), S19.meta_data.latitude(TNuwind(ii5_7)), 20, log10(sum(S19.classcount(TNuwind(ii5_7),nano_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii5_7))), 'filled')
180+
cb = colorbar('Location', 'northoutside'); title(cb,'log10 Misc. nanoplankton (ml^{-1})')
181+
grid on
182+
183+
nexttile
184+
scatter(S19.meta_data.longitude(TNuwind(ii5_7)), S19.meta_data.latitude(TNuwind(ii5_7)), 20, sum(S19.classcount(TNuwind(ii5_7),dino_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii5_7)), 'filled')
185+
cb = colorbar('Location', 'northoutside'); title(cb,'log10 Dinoflagellates (ml^{-1})')
186+
grid on
187+
188+
nexttile
189+
scatter(S19.meta_data.longitude(TNuwind(ii5_7)), S19.meta_data.latitude(TNuwind(ii5_7)), 20, log10(sum(S19.classcount(TNuwind(ii5_7),diatom_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii5_7))), 'filled')
190+
cb = colorbar('Location', 'northoutside'); title(cb,'log10 Diatoms (ml^{-1})')
191+
grid on
192+
%%
193+
194+
figure, set(gcf, 'position', [350 50 550 580])
195+
th = tiledlayout(2,4);
196+
nexttile
197+
Z = (Attune_uw_match.biovolume0to2(iA3_7) + Attune_uw_match.biovolume2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
198+
scatter(Attune_uw_match.longitude_fullres(iA3_7),Attune_uw_match.latitude_fullres(iA3_7), 20, Z, 'filled')
199+
cb = colorbar('Location', 'northoutside'); title(cb, 'Phyto 1-5 \mum ESD (\mum^{3} ml^{-1})')
200+
grid on
201+
202+
nexttile
203+
Z = sum(S19.classbiovol(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
204+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, Z, 'filled')
205+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto ~5-100 \mum ESD (\mum^{3} ml^{-1})')
206+
grid on
207+
208+
nexttile
209+
Z1 = (Attune_uw_match.biovolume0to2(iA3_7) + Attune_uw_match.biovolume2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
210+
Z2 = sum(S19.classbiovol(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
211+
Z1int = interp1(datenum(Attune_uw_match.StartDate(iA3_7)), Z1, S19.mdate(TNuwind(ii4_7)));
212+
Z = Z1int+Z2;
213+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, Z, 'filled')
214+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto 1-100 \mum ESD (\mum^{3} ml^{-1})')
215+
grid on
216+
%%
217+
figure, set(gcf, 'position', [350 50 550 580])
218+
th = tiledlayout(2,4);
219+
nexttile
220+
Z = (Attune_uw_match.carbon0to2(iA3_7) + Attune_uw_match.carbon2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
221+
scatter(Attune_uw_match.longitude_fullres(iA3_7),Attune_uw_match.latitude_fullres(iA3_7), 20, Z/1000, 'filled')
222+
cb = colorbar('Location', 'northoutside'); title(cb, 'Phyto 1-5 \mum ESD (mgC m^{-3})')
223+
caxis([0 35])
224+
grid on
225+
226+
nexttile
227+
Z = sum(S19.classC(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
228+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, Z/1000, 'filled')
229+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto ~5-100 \mum ESD (mgC m^{-3})')
230+
grid on
231+
caxis([0 35])
232+
233+
nexttile
234+
Z1 = (Attune_uw_match.carbon0to2(iA3_7) + Attune_uw_match.carbon2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
235+
Z2 = sum(S19.classC(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
236+
Z1int = interp1(datenum(Attune_uw_match.StartDate(iA3_7)), Z1, S19.mdate(TNuwind(ii4_7)));
237+
Z = Z1int+Z2;
238+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, Z/1000, 'filled')
239+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto 1-100 \mum ESD (mgC m^{-3})')
240+
grid on
241+
caxis([0 70])
242+
243+
nexttile
244+
temp = gca;
245+
246+
nexttile
247+
Z = (Attune_uw_match.carbon0to2(iA3_7) + Attune_uw_match.carbon2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
248+
scatter(Attune_uw_match.longitude_fullres(iA3_7),Attune_uw_match.latitude_fullres(iA3_7), 20, log10(Z/1000), 'filled')
249+
cb = colorbar('Location', 'northoutside'); title(cb, 'Phyto 1-5 \mum ESD (log10 mgC m^{-3})')
250+
caxis(log10([5 35]))
251+
grid on
252+
253+
nexttile
254+
Z = sum(S19.classC(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
255+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, log10(Z/1000), 'filled')
256+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto ~5-100 \mum ESD (log10 mgC m^{-3})')
257+
grid on
258+
caxis(log10([5 35]))
259+
260+
nexttile
261+
Z1 = (Attune_uw_match.carbon0to2(iA3_7) + Attune_uw_match.carbon2to5(iA3_7))./Attune_uw_match.VolAnalyzed_ml(iA3_7);
262+
Z2 = sum(S19.classC(TNuwind(ii4_7),[diatom_ind; nano_ind; dino_ind]),2)./S19.meta_data.ml_analyzed(TNuwind(ii4_7));
263+
Z1int = interp1(datenum(Attune_uw_match.StartDate(iA3_7)), Z1, S19.mdate(TNuwind(ii4_7)));
264+
Z = Z1int+Z2;
265+
scatter(S19.meta_data.longitude(TNuwind(ii4_7)), S19.meta_data.latitude(TNuwind(ii4_7)), 20, log10(Z/1000), 'filled')
266+
cb = colorbar('Location', 'northoutside'); title(cb,'Phyto 1-100 \mum ESD (log10 mgC m^{-3})')
267+
grid on
268+
caxis(log10([10 70]))
269+
delete(temp)
270+
271+
%%
272+
castind = find(strcmp(S19.meta_data.cruise, 'TN368') & strcmp(S19.meta_data.sample_type, 'cast') & ~(S19.meta_data.skip));
273+
trnum = 1:length(tstime);
274+
S19.meta_data.transect(castind) = floor(interp1(tstime,trnum, S19.mdate(castind)));
275+
ic1 = (S19.meta_data.transect(castind)==1);
276+
ic5_7 = (ismember(S19.meta_data.transect(castind),[5 6 7]));
277+
ic8 = (S19.meta_data.transect(castind)==8);
278+
279+
castind = find(strcmp(S19.meta_data.cruise, 'TN368') & strcmp(S19.meta_data.sample_type, 'cast') & ~(S19.meta_data.skip));
280+
%castind = find(strcmp(S19.meta_data.cruise, 'TN368') & strcmp(S19.meta_data.sample_type, 'cast') & ~(S19.meta_data.skip) & S19.mdate < datenum('7-7-2019') );
281+
%%
282+
t1 = find(strcmp(S19.meta_data.cruise, 'TN368') & strcmp(S19.meta_data.sample_type, 'cast') & ~(S19.meta_data.skip) & S19.meta_data.transect==1 );
283+
figure
284+
tiledlayout(3,1)
285+
nexttile
286+
scatter( S19.meta_data.latitude(t1), S19.meta_data.depth(t1),40, log10(sum(S19.classcount(t1,diatom_ind),2)./S19.meta_data.ml_analyzed(t1)), 'filled')
287+
set(gca, 'ydir', 'rev', 'xdir', 'rev')
288+
nexttile
289+
scatter( S19.meta_data.latitude(t1), S19.meta_data.depth(t1),40, log10(sum(S19.classcount(t1,dino_ind),2)./S19.meta_data.ml_analyzed(t1)), 'filled')
290+
set(gca, 'ydir', 'rev', 'xdir', 'rev')
291+
nexttile
292+
scatter( S19.meta_data.latitude(t1), S19.meta_data.depth(t1),40, log10(sum(S19.classcount(t1,nano_ind),2)./S19.meta_data.ml_analyzed(t1)), 'filled')
293+
set(gca, 'ydir', 'rev', 'xdir', 'rev')
294+
colorbar
295+
296+
%%
297+
figure
298+
th = tiledlayout(2,4)
299+
[~,s] = sort(sum(S19.classcount(TNuwind(ii5_7),diatom_ind),1), 'descend');
300+
for ind = 1:8
301+
nexttile
302+
scatter(S19.meta_data.longitude(TNuwind(ii5_7)), S19.meta_data.latitude(TNuwind(ii5_7)), 20, log10(S19.classcount(TNuwind(ii5_7),diatom_ind(s(ind)))./S19.meta_data.ml_analyzed(TNuwind(ii5_7))), 'filled')
303+
title(class2use(diatom_ind(s(ind))), 'interpreter', 'none')
304+
end
305+
%%
306+
%%
307+
figure
308+
th = tiledlayout(2,4)
309+
[~,s] = sort(sum(S19.classcount(TNuwind(ii5_7),dino_ind),1), 'descend');
310+
for ind = 1:8
311+
nexttile
312+
scatter(S19.meta_data.longitude(TNuwind(ii5_7)), S19.meta_data.latitude(TNuwind(ii5_7)), 20, (S19.classcount(TNuwind(ii5_7),dino_ind(s(ind)))./S19.meta_data.ml_analyzed(TNuwind(ii5_7))), 'filled')
313+
title(class2use(dino_ind(s(ind))), 'interpreter', 'none')
314+
end
315+
%%
316+
317+
return
318+
319+
320+
%%
321+
figure
322+
[~,s] = sort(sum(S19.classC(TNuwind(ii8),diatom_ind),1), 'descend');
323+
figure
324+
tiledlayout(2,1)
325+
nexttile
326+
plot(S19.meta_data.longitude(TNuwind(ii8)), S19.classcount(TNuwind(ii8),diatom_ind(s(1:2)))./S19.meta_data.ml_analyzed(TNuwind(ii8)), '.-', 'linewidth', 2)
327+
hold on
328+
plot(S19.meta_data.longitude(TNuwind(ii8)), sum(S19.classcount(TNuwind(ii8),diatom_ind),2)./S19.meta_data.ml_analyzed(TNuwind(ii8)), '.-', 'linewidth', 2)
329+
330+
legend([class2use(diatom_ind(s(1:2))); 'all diatoms'])
331+
332+
cc = strmatch('Chaetoceros', class2use);
333+
334+
%%

0 commit comments

Comments
 (0)