Documentation of ct_gr_z500

Global Index (short | long) | Local contents | Local Index (short | long)

Help text

  get CTstar:

Cross-Reference Information

Listing of script ct_gr_z500

clear
cd /cdrom
filin = 'oi_dynam.nc';
nc = netcdf(filin, 'nowrite');
lev = nc{'lev'}(:);
lon = nc{'lon'}(:);
lat = nc{'lat'}(:);
ts  = nc{'TS'}(:);
z3  = nc{'Z3'}(:,8,:,:);
slp = nc{'PSL'}(:);
airt = nc{'AIRT'}(:,15,:,:);
nc = close(nc);

cd /home/disk/hayes2/dvimont/ccm/ccm3.6/data
load ct_gr.mat
whos

ct = ctstar(252:408);
ct = (ct-mean(ct))/std(ct-mean(ct));
gr = gr(252:408);
gr = (gr-mean(gr))/(std(gr-mean(gr)));

subplot(2,1,1)
     plot(1:157,ct)
     set(gca,'Xtick',2:12:157,'XTickLabel',82:98);
     axis([0 158 -3 3])
     grid
subplot(2,1,2)
     plot(1:157,gr)
     set(gca,'Xtick',2:12:157,'XTickLabel',82:98);
     axis([0 158 -3 3])
     grid

z3 = z3(1:157,:,:);
[ntim, nlat, nlon] = size(z3);
[z3, clim] = annave(z3(1:157,:,:));

ctpat = ct'*z3/(ntim);

ctpat = reshape(ctpat, nlat, nlon);

sd(1)
     mcont(ctpat, [-20:5:20])
sd(2)
     plot(1:157,ct)
     set(gca,'Xtick',2:12:157,'XTickLabel',82:98);
     axis([0 158 -3 3])
     grid

grpat = reshape(gr'*z3/(ntim),nlat,nlon);

sd(1)
     mcont(grpat, [-20:2:20])
sd(2)
     plot(1:157,ct)
     set(gca,'Xtick',2:12:157,'XTickLabel',82:98);
     axis([0 158 -3 3])
     grid

[slp, clim] = annave(slp);
slp = cosweight(slp, lat);
c = slp*slp'/(nlat*nlon);
[lam,pc,per]=eof(c);

p = pc(:,1:10) ./ (ones(ntim, 1) * std(pc(:,1:10)));
eof = p'*acosweight(slp, lat)./ntim;
eof = reshape(eof, 10, nlat, nlon);

num = 1;
sd(1)
     tem = squeeze(eof(num,:,:));
     mcont(tem,[-1000:100:1000],'stereo',[-90,270]);
sd(2)
     plot(1:199,p(:,num))
     set(gca,'Xtick',2:12:157,'XTickLabel',82:98);
     axis([0 200 -3 3])
     grid