ArgoSpainWebPage/createDataSetStatus_FunctionProfiles.m at master · PedroVelez/ArgoSpainWebPage · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
function createDataSetStatus_FunctionProfiles(OneFloatData,Limits)
%Colors
cl=parula;
color=linspace(1,length(cl),size(OneFloatData.sals,2));

fprintf('profiles, ')

%Postions
RepresentaOxigen=0;
hPTPosition=[0.14 0.05 0.30 0.42];
hPSPosition=[0.56 0.05 0.30 0.42];
if isfield(OneFloatData,'oxys')
    if length(OneFloatData.oxys)==length(OneFloatData.pres)
        RepresentaOxigen=1;
        hPTPosition=[0.05 0.05 0.27 0.42];
        hPSPosition=[0.37 0.05 0.27 0.42];
        hPOPosition=[0.69 0.05 0.27 0.42];
    end
end

%% Temperature
if exist('hPTPosition','var') ==1
    hT=axes('position',hPTPosition);
elseif RepresentaOxigen==0
    hT=subplot(1,2,1);set(hT,'clipping','on');
else
    hT=subplot(1,3,1);set(hT,'clipping','on');
end

%CTD profile in the case is was sampled during deployment.
if isfield(OneFloatData,'CTD0')==1
    hCTD0=plot(OneFloatData.CTD0.ptmp,-OneFloatData.CTD0.pres,'-','color',[0.75 0.75 0.75],'linewidth',3);hold on
end

%Fist profile
hfp=plot(OneFloatData.ptms(:,1),-OneFloatData.pres(:,1),'color',cl(ceil(color(1)),:),'linewidth',1.25);hold on

%Following profiles
for j=1:size(OneFloatData.tems,2)
    plot(OneFloatData.ptms(:,j),-OneFloatData.pres(:,j),'color',cl(ceil(color(j)),:),'linewidth',1.25);
end

%Last profile
plot(OneFloatData.ptms(:,end),-OneFloatData.pres(:,end),'color',cl(ceil(color(end)),:),'linewidth',3);grid on
hlp=plot(OneFloatData.ptms(:,end),-OneFloatData.pres(:,end),'color','k','linewidth',1);grid on

if isfield(OneFloatData,'CTD0')==1
    hl=legend([hCTD0,hfp,hlp],sprintf('Perfil CTD inicial'),sprintf('Primer perfil: %s',datestr(OneFloatData.julds(1),1)),sprintf('Último perfil: %s',datestr(OneFloatData.julds(end),1)),'Location','southeast');
else
    hl=legend(hlp,sprintf('Último perfil: %s',datestr(OneFloatData.julds(end),1)),'Location','southeast');
end
hl.Box='off';

colormap(parula)
if size(OneFloatData.tems,2)==1
    caxis([1 2])
else
    caxis([1 size(OneFloatData.tems,2)])
end

axis([Limits.minT Limits.maxT -Limits.maxP 0])
hT.FontSize=7;
hT.TitleFontWeight='normal';
hT.XColor=[0.5 0.5 0.5];
hT.YColor=[0.5 0.5 0.5];
hT.ClippingStyle='rectangle';
title('Potential temperature [ITS-90]','fontsize',8,'VerticalAlignment','baseline')

%% Salinity
if exist('hPSPosition','var') ==1
    hS=axes('position',hPSPosition);
elseif RepresentaOxigen==0
    hS=subplot(1,2,2);set(hS,'clipping','on');
else
    hS=subplot(1,3,2);set(hS,'clipping','on');
end

%CTD profile in the case is was sampled during deployment.
if isfield(OneFloatData,'CTD0')==1
    hCTD0=plot(OneFloatData.CTD0.salt,-OneFloatData.CTD0.pres,'-','color',[0.75 0.75 0.75],'linewidth',3);hold on
end
%Fist profile
hfp=plot(OneFloatData.sals(:,1),-OneFloatData.pres(:,1),'color',cl(ceil(color(1)),:),'linewidth',1.25);hold on

%Following profiles
for j=1:size(OneFloatData.tems,2)
    plot(OneFloatData.sals(:,j),-OneFloatData.pres(:,j),'color',cl(ceil(color(j)),:),'linewidth',1.25);
end

%Last profile
plot(OneFloatData.sals(:,end),-OneFloatData.pres(:,end),'color',cl(ceil(color(end)),:),'linewidth',3);grid on
hlp=plot(OneFloatData.sals(:,end),-OneFloatData.pres(:,end),'color','k','linewidth',1);grid on

if isfield(OneFloatData,'CTD0')==1
    hl=legend([hCTD0,hfp,hlp],sprintf('Perfil CTD inicial'),sprintf('Primer perfil: %s',datestr(OneFloatData.julds(1),1)),sprintf('Último perfil: %s',datestr(OneFloatData.julds(end),1)),'Location','southeast');
else
    hl=legend(hlp,sprintf('Último perfil: %s',datestr(OneFloatData.julds(end),1)),'Location','southeast');
end
hl.Box='off';

colormap(parula)
if size(OneFloatData.tems,2)==1
    caxis([1 2])
else
    caxis([1 size(OneFloatData.tems,2)])
end
axis([Limits.minS Limits.maxS -Limits.maxP 0])
hS.FontSize=7;
hS.TitleFontWeight='normal';
hS.XColor=[0.5 0.5 0.5];
hS.YColor=[0.5 0.5 0.5];
hS.ClippingStyle='rectangle';
title('Salinity [psu]','fontsize',8,'VerticalAlignment','baseline')

%% Oxygen
if RepresentaOxigen==1
    if exist('hPOPosition','var') ==1
        hO=axes('position',hPOPosition);
    elseif RepresentaOxigen==0
        hO=subplot(1,2,3);set(hO,'clipping','on');
    else
        hO=subplot(1,3,3);set(hO,'clipping','on');
    end

    %Fist profile
    plot(OneFloatData.oxys(:,1),-OneFloatData.pres(:,1),'color',cl(ceil(color(1)),:),'linewidth',1.25);hold on
    %Following profiles
    for j=1:size(tems,2)
        plot(OneFloatData.oxys(:,j),-OneFloatData.pres(:,j),'color',cl(ceil(color(j)),:),'linewidth',1.25);
    end
    %Last profiles
    plot(OneFloatData.oxys(:,end),-OneFloatData.pres(:,end),'color',cl(ceil(color(end)),:),'linewidth',3);grid on
    plot(OneFloatData.oxys(:,end),-OneFloatData.pres(:,end),'color','k','linewidth',1);grid on

    colormap(parula)
    caxis([1 size(OneFloatData.tems,2)])
    axis([Limits.minO Limits.maxO -Limits.maxP 0])
    hO.FontSize=7;
    hO.TitleFontWeight='normal';
    hO.XColor=[0.5 0.5 0.5];
    hO.YColor=[0.5 0.5 0.5];
    hO.ClippingStyle='rectangle';
    title('Dissolved oxygen [micromole/kg]','fontsize',8,'VerticalAlignment','baseline')
end