Global Index (short | long) | Local contents | Local Index (short | long)
Old stuff
This script calls | |
---|---|
clear cd /home/disk/tao/dvimont/matlab/CSIRO/New_calcs/data load heat_pcs_nhem.mat pcs = rpcs; ct = getct; ct = detrend(ct); [b, a] = butter(9, 2/6); [b2, a2] = butter(9, 2/12); pcs = filtfilt(b, a, ct) - filtfilt(b2, a2, ct); tim = 101:1000; lims = [90 310 -80 80]; [var, lat, lon, depth, middepth] = getheat(lims, 4:7, tim, 'temp'); var = getnc('temp', lims, 1, tim); [lat, lon, depth, lm] = getll('temp', lims); var2 = getflx('psl', lims, tim); [lat2, lon2] = getll('temp', lims); var3 = getnc('taux', lims, 1, tim); [lat3, lon3] = getll('taux', lims); var4 = getnc('heat', lims, 1, tim); [lat4, lon4] = getll('heat', lims); var4 = getnc('temp', lims, 5, tim); [lat4, lon4] = getll('heat', lims); lag = -3:1:0; nlag = length(lag); reg1 = regress_eof(var, pcs(:,1), lag); reg2 = regress_eof(var2, pcs(:,1), lag); reg3 = regress_eof(var3, pcs(:,1), lag); reg4 = regress_eof(var4, pcs(:,1), lag); default_global; FRAME = [95 305 -75 75]; figure(1); fo(1); clf; cint = 0.05; cint2 = 0.1;; for i = 1:nlag; subplot(4,2,2*i-1); dg2(lat, lon); gcont(reg1(i,:,:), cint); dc2(lm, 0.3, 1000); subplot(4,2,2*i); dg2(lat2, lon2); gcont(reg2(i,:,:), cint2); dg2(lat, lon); dc2(lm, 0.3, -1000); end figure(2); fo(1); clf; cint = 0.01; cint2 = 0.05;; for i = 1:nlag; subplot(4,2,2*i-1); dg2(lat3, lon3); gcont(reg3(i,:,:), cint); dg2(lat, lon); dc2(lm, 0.3, 1000); subplot(4,2,2*i); dg2(lat, lon); gcont(reg4(i,:,:), cint2); dg2(lat, lon); dc2(lm, 0.3, -1000); end subplot(3,2,2); gcont(-reg2, cint); dc2(lm); pcs = hpcs; reg1 = regress_eof(var, pcs(:,1), lag); reg2 = regress_eof(var, pcs(:,2), lag); subplot(3,2,3); gcont(reg1, cint); dc2(lm); subplot(3,2,4); gcont(-reg2, cint); dc2(lm); pcs = lpcs; reg1 = regress_eof(var, pcs(:,1), lag); reg2 = regress_eof(var, pcs(:,2), lag); subplot(3,2,5); gcont(-reg1, cint); dc2(lm); subplot(3,2,6); gcont(-reg2, cint); dc2(lm);