Documentation of regress_NP_eof

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

Help text

lims = [110 265 20 65]; %  North Pacific Region

Cross-Reference Information

Listing of script regress_NP_eof

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load NP_LP10_T1_EOF_yr101-550.mat;  tim = 101:550;
load NP_LP10_T1_EOF_yr551-1000.mat; tim = 551:1000;
load EQ_LP10_T1_EOF_yr101-1000.mat; tim = 101:1000;
lims = [110 300 -20 20]; %  EQ Region

temp = getnc('temp', lims, 1, tim);
[lat, lon, depth] = getll('temp', lims);

treg = regress_eof(temp, pcs(:,1), 0);
treg2 = regress_eof(temp, pcs(:,2), 0);

default_global; FRAME = lims;

figure(3); clf; figure_orient;

cint = 0.1; clev = [-1:cint:-cint cint:cint:1];



subplot(3,1,1);
  gcont(treg, clev);
  dc2(treg);
  axis([110 260 20 60]);
  ylabel('RAW EOF1');
  title('North Pacific EOF 1 and 2 for years 101-1000');
 
subplot(3,1,2);
  gcont(treg2, clev);
  dc2(treg2);
  axis([110 260 20 60]);
  ylabel('RAW EOF2');

subplot(6,1,5);
  plot(tim, pcs(:,1), '-k');
  axis([min(tim)-1 max(tim)+1 -3 3]);
  grid on;
  ylabel(['PC1 - ' num2str(round(per(1))) '%']);

subplot(6,1,6);
  plot(tim, pcs(:,2), '-k');
  axis([min(tim)-1 max(tim)+1 -3 3]);
  grid on;
  ylabel(['PC2 - ' num2str(round(per(2))) '%']);

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Mid_Lats/Plot_NP


load EQ_LP10_T1_EOF_yr101-1000.mat; pcseq = pcs;
load NP_LP10_T1_EOF_yr101-1000.mat;

lags = -30:30; nlag = length(lags);
figure(4); figure_orient;

for i = 1:nlag;
  a(i) = corr(pcseq(:,1), pcs(:,1), lags(i));
end
subplot(3,1,1);
  bar(lags, a);
  axis([min(lags) max(lags) -1 1]);
  grid on
  title('PC\_EQ, PC1\_NP');

for i = 1:nlag;
  a(i) = corr(pcseq(:,1), pcs(:,2), lags(i));
end
figure(4);% figure_landscape(2);
subplot(3,1,2);
  bar(lags, a);
  axis([min(lags) max(lags) -1 1]);
  grid on
  title('PC\_EQ, PC2\_NP');

for i = 1:nlag;
  a(i) = corr(pcs(:,1), pcs(:,2), lags(i));
end
figure(4);% figure_landscape(2);
subplot(3,1,3);
  bar(lags, a);
  axis([min(lags) max(lags) -1 1]);
  grid on
  title('PC1\_NP, PC2\_NP');

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/Mid_Lats/Plot_NP