Global Index (short | long) | Local contents | Local Index (short | long)
Look at zonal cross sections
| This script calls | |
|---|---|
clear
lims = [105 310 -2 2];
lev = 1:10;
tim = 101:450;
[u, v] = getnc('u', 'v', lims, lev, tim);
w = getnc('wl', lims, 1:9, tim);
[lat, lon, depth] = getll('u', lims);
[latw, lonw, depthw] = getll('wl', lims);
u = squeeze(mean(u)); v = squeeze(mean(v)); w = squeeze(mean(w));
u = squeeze(mean(u, 2)); v = squeeze(mean(v, 2)); w = squeeze(mean(w, 2));
w = [zeros(1, size(w, 2)); w];
w(1,find(isnan(w(2,:)))) = NaN;
figure(1); figure_orient;
subplot(4,1,1);
pncont(lon, -1*depth, u, [-40:1:40], 'k');
dcv(lon, -1*depth, u);
ylabel('Mean U');
axis([105 299 -500 0]);
set(gca, 'XTick', 105:15:300, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 cm s^-^1'])
subplot(4,1,2);
pncont(lon, -1*depth, v, [-40:1:40], 'k');
dcv(lon, -1*depth, v);
ylabel('Mean V');
axis([105 299 -500 0]);
set(gca, 'XTick', 105:15:300, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 cm s^-^1'])
subplot(4,1,3);
pncont(lonw, -1*[0; depthw], w, [-40:0.5:40], 'k');
dcv(lonw, -1*[0; depthw], w);
ylabel('Mean W');
axis([105 299 -500 0]);
set(gca, 'XTick', 105:15:300, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 0.5 x 10^4 cm s^-^1'])
temp = getnc('temp', lims, lev, tim);
temp = squeeze(mean(temp));
temp = squeeze(mean(temp, 2));
[lat, lon, depth] = getll('temp', lims);
subplot(4,1,4);
pncont(lon, -1*depth, temp, [0:40], 1, 'k');
dcv(lonw, -1*depth, w);
ylabel('Mean Temperature');
axis([105 299 -500 0]);
set(gca, 'XTick', 105:15:300, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 Degree C'])
subplot(4,1,1);
title('Mean SH Midlatitude Velocity and Temperature (50s to 40s), years 101-550');
cd ~/matlab/CSIRO/Plots_Clim
print -dps2 ML45SH_vel_temp_yr101-550.ps
clear
lev = 1:10;
tim = 101:550;
%lims = [115 150 -60 60];
%lims = [150 195 -60 60];
%lims = [195 240 -60 60];
lims = [240 275 -60 60];
[u, v, temp] = getnc('u', 'v', 'temp', lims, lev, tim);
u = squeeze(mean(mean(u,4)));
v = squeeze(mean(mean(v,4)));
temp = squeeze(mean(mean(temp,4)));
[latu, lonu, depthu] = getll('u', lims);
[latt, lont, deptht] = getll('temp', lims);
w = getnc('wl', lims, 1:9, tim);
w = squeeze(mean(mean(w,4)));
w = [zeros(1, size(w, 2)); w];
[latw, lonw, depthw] = getll('wl', lims);
% Plot the data
figure(1); figure_orient;
subplot(4,1,1);
pncont(latu, -1*depthu, u, [-40:1:40], 'k');
dcv(latu, -1*depthu, u);
ylabel('Mean U');
axis([-60 60 -500 0]);
set(gca, 'XTick', -60:15:60, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 cm s^-^1'])
subplot(4,1,2);
pncont(latu, -1*depthu, v, [-40:1:40], 'k');
dcv(latu, -1*depthu, v);
ylabel('Mean V');
axis([-60 60 -500 0]);
set(gca, 'XTick', -60:15:60, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 cm s^-^1'])
subplot(4,1,3);
pncont(latw, -1*[0; depthw], w, [-40:0.5:40], 'k');
dcv(latw, -1*[0; depthw], w);
ylabel('Mean W');
axis([-60 60 -500 0]);
set(gca, 'XTick', -60:15:60, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 0.5 x 10^4 cm s^-^1'])
subplot(4,1,4);
pncont(latt, -1*deptht, temp, [0:40], 1, 'k');
dcv(latt, -1*deptht, temp);
ylabel('Mean Temperature');
axis([-60 60 -500 0]);
set(gca, 'XTick', -60:15:60, ...
'YTick', -500:125:0, 'YTickLabel', ['500'; ' '; '250'; ' '; ' 0']);
box on
xlabel(['Contour Interval: 1 Degree C'])
subplot(4,1,1);
title(['Mean Velocity and Temperature (' num2str(lims(1)) 'E to ' ...
num2str(lims(2)) 'E), years 101-550']);
cd ~/matlab/CSIRO/Plots_Clim
print -dps2 zonal_240to275_vel_temp_yr101-550.ps