Global Index (short | long) | Local contents | Local Index (short | long)
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/EARTH/atm/netcdf
| This script calls | |
|---|---|
clean
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/HTRAE/atm/netcdf
varn = 'PSL';
for i = 9:12;
ind = i-8;
if i < 10;
filin = ['0010-0' num2str(i) '.nc'];
else
filin = ['0010-' num2str(i) '.nc'];
end
disp(filin)
nc = netcdf(filin, 'nowrite');
tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
nc = close(nc);
end
for j = 11:19;
for i = 1:12;
ind = (j-11)*12+4+i;
if i < 10;
filin = ['00' num2str(j) '-0' num2str(i) '.nc'];
else
filin = ['00' num2str(j) '-' num2str(i) '.nc'];
end
disp(filin)
nc = netcdf(filin, 'nowrite');
tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
nc = close(nc);
end
end
for i = 1:8;
ind = 112+i;
if i < 10;
filin = ['0020-0' num2str(i) '.nc'];
else
filin = ['0020-' num2str(i) '.nc'];
end
disp(filin)
nc = netcdf(filin, 'nowrite');
tem(ind,:,:,:) = nc{varn, 1}(:,:,:,:);
nc = close(nc);
end
nc = netcdf(filin, 'nowrite');
lat = nc{'lat'}(:);
lon = nc{'lon'}(:);
lev = nc{'lev'}(:);
nc = close(nc);
% GET EOFS OF EARTH
psl = tem;
[ntim, nlat, nlon] = size(psl);
lims = [0 360 20 90];
timind = sort([3:12:ntim 4:12:ntim 5:12:ntim ...
6:12:ntim 7:12:ntim 8:12:ntim]);
[xk, yk] = keep_var(lims, lon, lat);
[psl, clim] = annave(psl);
psl = psl(timind,yk,xk);
lat = lat(yk); lon = lon(xk);
psl = cosweight(psl, lat);
[ntim, nlat, nlon] = size(psl);
psl = reshape(psl, ntim, nlat*nlon);
[lam, lds, pcs, per] = eof_dan(psl);
psl = tem;
[psl, clim] = annave(psl);
psl = psl(timind,:,:);
[reg1a, c1a] = regress_eof(psl, pcs(:,1), 0);
[reg2a, c2a] = regress_eof(psl, pcs(:,2), 0);
[reg3a, c3a] = regress_eof(psl, pcs(:,3), 0);
[reg4a, c4a] = regress_eof(psl, pcs(:,4), 0);
[reg5a, c5a] = regress_eof(psl, pcs(:,5), 0);
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
%save ANNULAR_EARTH.mat reg1a reg2a reg3a reg4a reg5a ...
% c1a c2a c3a c4a c5a
% GET EOFS OF HTRAE
psl = tem;
[ntim, nlat, nlon] = size(psl);
lims = [0 360 -90 -20];
timind = sort([9:12:ntim 10:12:ntim 11:12:ntim ...
12:12:ntim 1:12:ntim 2:12:ntim]);
[xk, yk] = keep_var(lims, lon, lat);
[psl, clim] = annave(psl);
psl = psl(timind,yk,xk);
lat = lat(yk); lon = lon(xk);
psl = cosweight(psl, lat);
[ntim, nlat, nlon] = size(psl);
psl = reshape(psl, ntim, nlat*nlon);
[lam, lds, pcs, per] = eof_dan(psl);
psl = tem;
[psl, clim] = annave(psl);
psl = psl(timind,:,:);
[reg1b, c1b] = regress_eof(psl, pcs(:,1), 0);
[reg2b, c2b] = regress_eof(psl, pcs(:,2), 0);
[reg3b, c3b] = regress_eof(psl, pcs(:,3), 0);
[reg4b, c4b] = regress_eof(psl, pcs(:,4), 0);
[reg5b, c5b] = regress_eof(psl, pcs(:,5), 0);
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
%save ANNULAR_HTRAE.mat reg1b reg2b reg3b reg4b reg5b ...
% c1b c2b c3b c4b c5b
cd /home/disk/tao/dvimont/matlab/CCM/Htrea/data
load ANNULAR_HTRAE.mat
load ANNULAR_EARTH.mat
figure(1); fo(1); clf;
cd /home/disk/hayes2/dvimont/ccm/ccm3.6/run/sun/HTRAE/HTRAE/atm/netcdf
filin = '0010-09.nc';
nc = netcdf(filin, 'nowrite');
lat = nc{'lat'}(:);
lon = nc{'lon'}(:);
lev = nc{'lev'}(:);
nc = close(nc);
default_global;
FRAME = [0 360 15 90];
subplot(3,2,1); cla
mcont3(0.01*reg1a, 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('EARTH EOF1');
subplot(3,2,2); cla
mcont3(0.01*rot90(rot90(squeeze(reg1b))), 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('HTRAE EOF1');
subplot(3,2,3); cla
mcont3(0.01*reg2a, 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('EARTH EOF2');
subplot(3,2,4); cla
mcont3(0.01*rot90(rot90(squeeze(reg2b))), 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('HTRAE EOF2');
subplot(3,2,5); cla
mcont3(0.01*reg3a, 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('EARTH EOF3');
subplot(3,2,6); cla
mcont3(0.01*rot90(rot90(squeeze(reg3b))), 0.5, 'stereo', [90 270]);
landshade
axis_limits
title('HTRAE EOF3');
cd ~
print -dps2 tem.ps