Global Index (short | long) | Local contents | Local Index (short | long)
Start with EARTH
| This script calls | |
|---|---|
clean
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/Data2
load u_atlev_T42.mat; var1a = var1; var2a = var2;
load u_atlev_T42_2.mat;
u1 = 0.5*(var1a+var1);
u2 = 0.5*(var2a+var2);
load v_atlev_T42.mat; var1a = var1; var2a = var2;
load v_atlev_T42_2.mat;
v1 = 0.5*(var1a+var1);
v2 = 0.5*(var2a+var2);
load omega_atlev_T42.mat; var1a = var1; var2a = var2;
load omega_atlev_T42_2.mat;
omega1 = 0.5*(var1a+var1);
omega2 = 0.5*(var2a+var2);
load airt_atlev_T42.mat; var1a = var1; var2a = var2;
load airt_atlev_T42_2.mat;
airt1 = 0.5*(var1a+var1);
airt2 = 0.5*(var2a+var2);
clear var1 var2 var1a var2a
[ntim, nlev, nlat, nlon] = size(airt1);
tim = [12 1 2];
tim = [6 7 8];
tim = [3 4 5];
tim = [9 10 11];
figure(1); fo(1); clf;
subplot(3,1,1);
tem1 = squeeze(mean(u1(tim,:,:,:))); % Seasonally average
tem1 = squeeze(mean2(shiftdim(tem1, 2))); % Zonally average
tem2 = squeeze(mean(v1(tim,:,:,:)));
tem2 = squeeze(mean2(shiftdim(tem2, 2)));
tem3 = squeeze(mean(omega1(tim,:,:,:)));
tem3 = squeeze(mean2(shiftdim(tem3, 2)));
tem4 = squeeze(mean(airt1(tim,:,:,:)));
tem4 = squeeze(mean2(shiftdim(tem4, 2)));
tem4 = tem4.*(((1000./newlev').^(2/7))*ones(1, nlat)); % POT
cint = 5; clev = [-20:cint:-cint 0 cint:cint:60];
tem = tem1;
tem(abs(tem)<cint) = NaN;
h3 = surface(lat-mean(diff(lat)), -.1*(newlev+25), tem);
zd = get(h3, 'ZData'); set(h3, 'ZData', -1*ones(size(zd)));
shading flat;
% [h1, c1] = pncont(lat, -.1*newlev, tem1, clev, 0, 'k');
hold on;
nthin = 2;
h2 = quiver2(thin(lat, nthin), -.1*thin(newlev, nthin), ...
thin(tem2, nthin), -100*thin(tem3, nthin), 2);
[h4, c4] = pncont(lat, -.1*newlev, tem4, 50, 0, 'k');
hold off;
set(h2, 'color', 'k', 'linewidth', 2);
caxis([-40 40]);
set(gca, 'YScale', 'log', 'YTick', [-100 -80 -60 -40 -25 -15]);
set(gca, 'YTicklabel', [1000 800 600 400 250 150], 'box', 'on');
axis([-90 90 -110 -10]);
colorbar2
subplot(3,1,2);
tem1a = squeeze(mean(u2(tim,:,:,:))); % Seasonally average
tem1a = -squeeze(mean2(shiftdim(tem1a, 2))); % Zonally average
tem2a = squeeze(mean(v2(tim,:,:,:)));
tem2a = squeeze(mean2(shiftdim(tem2a, 2)));
tem3a = squeeze(mean(omega2(tim,:,:,:)));
tem3a = squeeze(mean2(shiftdim(tem3a, 2)));
tem4a = squeeze(mean(airt2(tim,:,:,:)));
tem4a = squeeze(mean2(shiftdim(tem4a, 2)));
tem4a = tem4a.*(((1000./newlev').^(2/7))*ones(1, nlat)); % POT
cint = 5; clev = [-60:cint:-cint 0 cint:cint:60];
tem = tem1a;
tem(abs(tem)<cint) = NaN;
h3 = surface(lat-mean(diff(lat)), -.1*(newlev+25), tem);
zd = get(h3, 'ZData'); set(h3, 'ZData', -1*ones(size(zd)));
shading flat;
% [h1, c1] = pncont(lat, -.1*newlev, tem1, clev, 0, 'k');
hold on;
nthin = 2;
h2 = quiver2(thin(lat, nthin), -.1*thin(newlev, nthin), ...
thin(tem2a, nthin), -100*thin(tem3a, nthin), 2);
[h4, c4] = pncont(lat, -.1*newlev, tem4a, 50, 0, 'k');
hold off;
set(h2, 'color', 'k', 'linewidth', 2);
caxis([-40 40]);
set(gca, 'YScale', 'log', 'YTick', [-100 -80 -60 -40 -25 -15]);
set(gca, 'YTicklabel', [1000 800 600 400 250 150], 'box', 'on');
axis([-90 90 -110 -10]);
colorbar2
subplot(3,1,3);
cint = 0.5; clev = [-6:cint:-cint 0 cint:cint:6];
tem = -tem1a-tem1;
tem(abs(tem)<cint) = NaN;
h3 = surface(lat-mean(diff(lat)), -.1*(newlev+25), tem);
zd = get(h3, 'ZData'); set(h3, 'ZData', -1*ones(size(zd)));
shading flat;
% [h1, c1] = pncont(lat, -.1*newlev, tem1, clev, 0, 'k');
hold on;
nthin = 2;
h2 = quiver2(thin(lat, nthin), -.1*thin(newlev, nthin), ...
thin(tem2a-tem2, nthin), ...
-100*thin(tem3a-tem3, nthin), 2);
[h4, c4] = pncont(lat, -.1*newlev, tem4a-tem4, 50, 0, 'k');
hold off;
set(h2, 'color', 'k', 'linewidth', 2);
caxis([-3 3]);
set(gca, 'YScale', 'log', 'YTick', [-100 -80 -60 -40 -25 -15]);
set(gca, 'YTicklabel', [1000 800 600 400 250 150], 'box', 'on');
axis([-90 90 -110 -10]);
colorbar2