Newer
Older
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
function plot_ClimatologyB(hMainFig, hPlotAxes)
% Function to plot mean climatology and standard deviation
%
% Input
% -----
% hTsgGUI ............ Handel to the main user interface
% hPlotAxes .......... Handels to the 3 graphic axes
%
% Output
% ------
% none
%
% $Id$
% Read surface climatology (annual, seasonal or monthly)
% ------------------------------------------------------
read_Climatology(hMainFig);
% Get data after read_Climatology because it could be change tsg.levitus.type
% -------------------------------------------------------------------------
tsg = getappdata( hMainFig, 'tsg_data' );
% if reading error, tsg.levitus.data is empty, no action
% ------------------------------------------------------
if isempty(tsg.levitus.data)
return
end
% get last selected climatology structure
% ---------------------------------------
s = get(findobj('Tag', 'TAG_UIMENU_CLIMATO_MAIN'), 'Userdata');
% select time dimension for climatology form saved structure s
% ------------------------------------------------------------
time_dim = s.time;
% round positive latitude and Longitude toward zero
% -------------------------------------------------
ind = find(tsg.LATX > 0);
lat(ind) = fix(tsg.LATX(ind)) + 0.5;
ind = find(tsg.LONX > 0);
lon(ind) = fix(tsg.LONX(ind)) + 0.5;
% rounds negative latitude and Longitudeto the nearest lowest integers
% ---------------------------------------------------------------------
ind = find(tsg.LATX <= 0);
lat(ind) = floor(tsg.LATX(ind)) + 0.5;
ind = find(tsg.LONX <= 0);
lon(ind) = floor(tsg.LONX(ind)) + 0.5;
% Calculates differences between adjacent elements of X.
% 0 if adajacent latitude or longitude are equal
% - 1 or -1 otherwise
% ------------------------------------------------------------
lat_diff = [diff( lat )'; 0];
lon_diff = [diff( lon )'; 0];
% Select latitude and longitude
% -----------------------------
ind = find(abs(lat_diff) == 1 | abs(lon_diff == 1));
lat2 = lat( ind );
lon2 = lon( ind );
dayd = tsg.DAYD( ind );
% Get Climatology
% LATX(80) = -0.5 et LATX(81) = 0.5
% LONX(180) = -0.5 et LONX(181) = 0.5
% ----------------
para = {'SSTP'; 'SSPS'};
for i = 1:2
mean = zeros(size(ind));
std = zeros(size(ind));
for ii=1:length(ind)
ilat = find(tsg.levitus.data.WOA01_LATX == lat2(ii));
ilon = find(tsg.levitus.data.WOA01_LONX == lon2(ii));
mean(ii) = tsg.levitus.data.(['WOA01_MEAN_' para{i}])(time_dim,1,ilat,ilon);
std(ii) = tsg.levitus.data.(['WOA01_STD_' para{i}])(time_dim,1,ilat,ilon);
end
% get handle of axe when SSPS is plotted
% --------------------------------------
hdls = findobj('-regexp','Tag', ['TAG_PLOT\d_LINE_' para{i}]);
% get handle of axe when SSTP is plotted
% --------------------------------------
if strcmp( para{i}, 'SSTP' )
hdls = findobj('-regexp','Tag', 'TAG_PLOT\d_LINE_SSJT');
if isempty(hdls)
hdls = findobj('-regexp','Tag', 'TAG_PLOT\d_LINE_SSTP');
end
end
if ~isempty(hdls)
hdl_current_axes = get(hdls(1),'parent');
axes( hdl_current_axes );
% Plot mean salinity climatology
% ------------------------------
line(dayd, mean, ...
'Tag', ['TAG_LINE_CLIMATO_MEAN_' para{i}], 'Linestyle', '-', 'Color','k');
% Plot with 3 standard deviation
% ------------------------------
line(dayd, mean + 3 * std, ...
'Tag', ['TAG_LINE_CLIMATO_STDDEV_PLUS_' para{i}], 'Linestyle', '-', 'Color','r');
line(dayd, mean - 3 * std, ...
'Tag', ['TAG_LINE_CLIMATO_STDDEV_MINUS_' para{i}], 'Linestyle', '-', 'Color','r');
end
end
% save tsg structure
% ------------------
setappdata( hMainFig, 'tsg_data', tsg );
end