Documentation of compare_eof_regions

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

Cross-Reference Information

Listing of script compare_eof_regions

clear
cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load NP2_LP10_detrend_L1-7_yr101-550.mat; pcsnp = pcs;
load NH_LP10_detrend_L1-7_yr101-550.mat; pcsnh = pcs;
load SP2_LP10_detrend_L1-7_yr101-550.mat; pcssp = pcs;
load SH_LP10_detrend_L1-7_yr101-550.mat; pcssh = pcs;
load EQ_LP10_detrend_L1-7_yr101-550.mat; pcseq = pcs;
tim = 101:550;
lims = [106 312 -63 63];
temp = getnc('temp', lims, 1, tim);
[lat, lon] = getll('temp', lims);
lagind = -5;
 regnp = regress_eof(temp, pcsnp, lagind);
 regsp = regress_eof(temp, pcssp, lagind);
 regsh = regress_eof(temp, pcssh, lagind);
 regnh = regress_eof(temp, pcsnh, lagind);
clev = [-5:.025:5];
subplot(2,2,1);
  gcont(regnp, clev);
  title('NP [110 300 -15 60]');
  dc2(regnp);
subplot(2,2,3);
  gcont(-1*regsp, clev);
  title('SP [110 300 -60 15]');
  dc2(regnp);
subplot(2,2,2);
  gcont(regnh, clev);
  title('NH [110 300 0 60]');
  dc2(regnp);
subplot(2,2,4);
  gcont(-1*regsh, clev);
  title('SH [110 300 -60 0]');
  dc2(regnp);