Documentation of save_20deg_ceof

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

Help text

  Do complex eof:

Cross-Reference Information

Listing of script save_20deg_ceof

clear
for biff = [1 3 5 7];
ind = 1;
cd /home/disk/hayes2/dvimont/csiro/matlab_data
if ind == 1;
  load 20degiso_yr101-550_temp_noave.mat
  lag = -60;
  ttit = ['101 - 550'];
elseif ind == 2;
  load 20degiso_yr551-1000_temp_noave.mat
  lag = 90;
  ttit = ['551 - 1000'];
end
ctlim = [110 300 -60 60];
ntim = length(tim); nlat = length(lat); nlon = length(lon);
temp3 = NaN * ones(ntim, nlat*nlon);
temp3(:,kp20) = temp20;
temp3 = detrend(temp3);
temp3 = cosweight(temp3, lat);
kp = find(~isnan(mean(temp3)));
temp3 = temp3(:,kp);
if biff == 1;
  yr1 = 10; yr2 = 10; tit = ['HP ( < 10 yr ) '];
  scomm = '20deg_CEOF_HP10_yr101-550.mat';
elseif biff == 2;
  yr1 = 20; yr2 = 20; tit = ['HP ( < 20 yr ) ']; 
  scomm = '20deg_CEOF_HP20_yr101-550.mat';
elseif biff == 3; 
  yr1 = 10; yr2 = 60; tit = ['LP ( > 10 yr ) '];
  scomm = '20deg_CEOF_LP10_yr101-550.mat';
elseif biff == 4;
  yr1 = 20; yr2 = 60; tit = ['LP ( > 20 yr ) '];
  scomm = '20deg_CEOF_LP20_yr101-550.mat';
elseif biff == 5;
  yr1 = 8; yr2 = 20; tit = ['HP ( 8 yr ) '];
  scomm = '20deg_CEOF_HP8_yr101-550.mat';
elseif biff == 6;
  yr1 = 10; yr2 = 50; tit = ['BP ( 10-50 yr ) '];
  scomm = '20deg_CEOF_BP10-50_yr101-550.mat';
elseif biff == 7;
  tit = ['Unfiltered Data '];
  scomm = '20deg_CEOF_RAW_yr101-550.mat';
end
[b1, a1] = butter(6, 2/yr1); [b2, a2] = butter(6, 2/yr2);
if biff < 3 | biff == 5;
  temp3 = temp3 - filtfilt(b1, a1, temp3); disp('HP Filtered');
elseif biff > 2 & biff < 5;
  temp3 = filtfilt(b1, a1, temp3); disp('LP Filtered');
elseif biff ~= 7
  temp3 = filtfilt(b1, a1, temp3) - filtfilt(b2, a2, temp3); disp('BP Filtered');
end

[lam, lds, pcs, per] = complex_eof(temp3);

cd /home/disk/hayes2/dvimont/csiro/matlab_data
eval(['save ' scomm ' lam lds pcs per lat lon tim kp20 ; ']);

end