Documentation of gu_philander_20degiso

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

Help text

  Note:  everything in kp2 is also in kp1.  However, kp1 has two 
  extra points -- the following command finds the common points, 
  then we will eliminate the uncommon points.

Cross-Reference Information

Listing of script gu_philander_20degiso

clear
data
load LP9_detrend_L1-7_EOF_yr101-1000.mat 
cd /home/disk/hayes2/dvimont/csiro/matlab_data/20Deg_Isotherm
load 20degiso_yr101-550_temp_noave.mat 
temp1 = temp20; kp1 = kp20; lat1 = lat; lon1 = lon;
load 20degiso_yr551-1000_temp_noave.mat 
temp2 = temp20; kp2 = kp20; lat2 = lat; lon2 = lon;

common = find(ismember(kp1, kp2));
kp1 = kp1(common);
temp1 = temp1(:, common);

temp = [temp1; temp2];

%  Get regressions

lags = -12:12;

[reg1, c1] = regress_eof(temp, pcs, lags);

nlat = length(lat); nlon = length(lon);

default_global; FRAME = [100 300 -30 30];

for i = 1:12;
  subplot(6,2,i);
    tem = repmat(NaN, nlat, nlon);
    tem(kp2) = squeeze(reg1(i,:));
    gcont(tem, 0.025);
end