Documentation of pdo_calc

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

Help text

  look at SST regressions, to determine what modes 1 and 2 are

Cross-Reference Information

Listing of script pdo_calc

clear
cd ~/matlab/CSIRO/Thesis/Data
load RAW_detrend_L1-7_EOF_yr101-1000.mat; pc1 = pcs;
load PDO_detrend_L1_EOF_yr101-1000.mat; pc2 = pcs;
pc1 = getnc('temp', [180 270 -6 6], 1, tim);
pc1 = squeeze(mean2(mean2(shiftdim(pc1, 1))));
pc1 = detrend(pc1);

lims = [100 300 -60 60];
tim = 101:1000;
var = getnc('temp', lims, 1, tim);

reg1 = regress_eof(var, -pc2, 0);
reg2 = regress_eof(var, pc2(:, 2), 0);

[lat, lon, depth, lm] = getll('temp', lims);
default_global; FRAME = lims;

figure(1); clf; fo;
subplot(2, 1, 1);
  gcont(reg1, 0.05);
  dc2(lm);
  title('PC1 [100 255 20 60] of NP SST');
  ylabel([num2str(round(per(1))) '% Variance Explained']);
  xlabel('Contour Interval:  0.05 K std^-^1');
subplot(2, 1, 2);
  gcont(reg2, 0.05);
  dc2(lm);
  title('PC2 [100 255 20 60] of NP SST');
  title('PC2 [100 255 20 60] of NP SST');
  ylabel([num2str(round(per(2))) '% Variance Explained']);
  xlabel('Contour Interval:  0.05 K std^-^1');

%  look at lagged regression between pc1 and pc2

lags = -15:15;
for i = 1:length(lags);
  a(i) = corr(pc2(:, 1), pc2(:, 2), lags(i));
end
for i = 1:length(lags);
  b(i) = corr(pc1(:, 1), -pc2(:, 1), lags(i));
end
for i = 1:length(lags);
  c(i) = corr(pc1(:, 1), pc2(:, 2), lags(i));
end

figure(2); fo;
subplot(3, 1, 1);
  bar(lags, a);
  axis([min(lags-1) max(lags+1) -0.2 0.5]);
  grid on;
  set(gca, 'YTick', -1:.1:1, 'XTick', -20:2:20);
  title('Lagged Correlation between PC1 and PC2');
  xlabel('PC1 leads PC2                   PC2 leads PC1');
subplot(3, 1, 2);
  bar(lags, b);
  axis([min(lags-1) max(lags+1) -0.2 0.5]);
  grid on;
  set(gca, 'YTick', -1:.1:1, 'XTick', -20:2:20);
  title('Lagged Correlation between CT and PC1');
  xlabel('CT leads PC1                    PC1 leads CT');
subplot(3, 1, 3);
  bar(lags, c);
  axis([min(lags-1) max(lags+1) -0.2 0.5]);
  grid on;
  set(gca, 'YTick', -1:.1:1, 'XTick', -20:2:20);
  title('Lagged Correlation between CT and PC2');
  xlabel('CT leads PC2                    PC2 leads CT');

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap5