Documentation of lagged_regression

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

Help text

  Load the pcs:

Cross-Reference Information

Listing of script lagged_regression

clear
varn = 'temp';

biff2 = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
if biff2 == 1;
  load BP10-60yr_L1-7_EOF_yr2_corr.mat; tit = 'Bandpass Filtered Data ( 10 - 60 Years )';
elseif biff2 == 2;
  load HP8_L1-7_EOF.mat; tit = 'Highpass Filtered Data ( < 8 Years )';
elseif biff2 == 3;
  load RAW_L1-7_EOF.mat; tit = 'Unfiltered Data';
end
pcs = pcs;

nfrm = [-10:2:12]; nfrm = [-6:5];
tim = 551:1000; % 101:550;
lev = 1;
lims = [106 308 -65 65];

%temp = getnc(varn, lims, lev, tim);
cd ~/matlab/CSIRO/Heat
temp = (getheat(1:7, tim, lims));% - getheat(1:3, tim, lims))/1e8;
[temp, tem] = remove_mean(temp);
[b, a] = butter(6, 2/10); [b2, a2] = butter(6, 2/60);
temp = filtfilt(b, a, temp) - filtfilt(b2, a2, temp);

%[lat, lon, depth] = getll(varn, lims);

[treg, tcoef] = regress_eof(temp, pcs, nfrm);
score = tscore(450/10, 2.5);
tval = tcoef * sqrt(450/10 - 1) ./ sqrt(1 - tcoef.^2);

%  Plot the data

cint = 0.5;
figure(1); figure_orient(1);
clev = [-10:cint:-cint 0  cint:cint:10];
default_global;  FRAME = [110 299 -60 60];
for i = 1:6;
  tem = squeeze(1e-8*treg(i,:,:));
%  tem = squeeze(tcoef(i,:,:)).^2;
  tem2 = zeros(size(tem));
  tem2(find(abs(tval(i,:,:)) > score)) = NaN;
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
    shade_ccoef(tem2);
    title(['Lag:  ' num2str(nfrm(i))]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2']);
  end
end
subplot(3,2,3);
%  ylabel([tit ':  L' num2str(lev) ' Temp regressed on HC (L1-7) PC1']);
  ylabel([tit ':  L1-7 Heat Content regressed on HC (L1-7) PC1']);
%  ylabel([tit ':  Heat Flux regressed on HC (L1-7) PC1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_EOF
%eval(['print -dps2 ' varn '_L' num2str(lev) '_HC1-7_neglag.ps']);

figure(2); figure_orient(2);
default_global;  FRAME = [110 299 -60 60];
for i = 1:6;
  tem = squeeze(1e-8*treg(i+6,:,:));
%  tem = squeeze(tcoef(i+6,:,:)).^2;
  tem2 = zeros(size(tem));
  tem2(find(abs(tval(i+6,:,:)) > score)) = NaN;
  subplot(3,2,i);
    gcont(tem, clev);
    dc2(tem);
    shade_ccoef(tem2);
    title(['Lag:  ' num2str(nfrm(i+6))]);
  if i >= 5;
    xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);
  end
end
subplot(3,2,3);
  ylabel([tit ':  L' num2str(lev) ' Temp regressed on HC (L1-7) PC1']);
  ylabel([tit ':  Heat Flux regressed on HC (L1-7) PC1']);
  ylabel([tit ':  L1-7 Heat Content regressed on HC (L1-7) PC1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_EOF
%eval(['print -dps2 ' varn '_L' num2str(lev) '_HC1-7_poslag.ps']);




%  Wind stress, or other vector

lev = 5;
%varn1 = 'taux'; varn2 = 'tauy'; varn = 'stress';
varn1 = 'u'; varn2 = 'v'; varn = ['ocevel'];
[taux, tauy] = getnc(varn1, varn2, lims, lev, tim);
[taux, tem] = remove_mean(taux);
[tauy, tem] = remove_mean(tauy);
[b, a] = butter(6, 10/2);
taux = filtfilt(b, a, taux);
tauy = filtfilt(b, a, tauy);

[lat, lon] = getll(varn1, lims);

xreg = regress_eof(taux, pcs, nfrm);
yreg = regress_eof(tauy, pcs, nfrm);

%  Plot the data

figure(1); figure_orient(1);
default_global;  FRAME = [110 299 -60 60];

cint = -0.2;
figure(1); figure_orient;
for i = 1:6;
  tem1 = squeeze(xreg(i,:,:));
  tem2 = squeeze(yreg(i,:,:));
  subplot(3,2,i);
    gquiv(tem1, tem2, cint, 2, 'cm s^-^2');
    dc2(tem1);
    title(['Lag:  ' num2str(nfrm(i))]);
end
for i = 1:4;
  subplot(3,2,i);
  xlabel('');
end
subplot(3,2,3);
  ylabel([tit ':  L5 Ocean Velocity regressed on HC (L1-7) PC1']);
%  ylabel([tit ':  Wind Stress regressed on HC (L1-7) PC1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_EOF
eval(['print -dps2 ' varn '_L' num2str(lev) '_HC1-7_neglag.ps']);

figure(2); figure_orient;
for i = 7:12;
  tem1 = squeeze(xreg(i,:,:));
  tem2 = squeeze(yreg(i,:,:));
  subplot(3,2,i-6);
    gquiv(tem1, tem2, cint, 2, 'dynes cm^-^2');
    dc2(tem1);
    title(['Lag:  ' num2str(nfrm(i))]);
end
for i = 1:4;
  subplot(3,2,i);
  xlabel('');
end
subplot(3,2,3);
  ylabel([tit ':  Wind Stress regressed on HC (L1-7) PC1']);
  ylabel([tit ':  L5 Ocean Velocity regressed on HC (L1-7) PC1']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Lagged_regressions/Plot_EOF
eval(['print -dps2 ' varn '_L' num2str(lev) '_HC1-7_poslag.ps']);