Documentation of plot_regressions

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

Help text

  Old stuff

Cross-Reference Information

Listing of script plot_regressions

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);