Documentation of eof_slp

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

Cross-Reference Information

Listing of script eof_slp

clear
for biff = 1:3;
csirod
load ML_ANN_slp.mat
back
lims = [105 295 -60 -20];
[xk, yk] = keep_var(lims, lon, lat);
slp = slp(:, yk, xk);
lat = lat(yk); lon = lon(xk);
slp = detrend(slp);
if biff == 2;
  [b, a] = butter(9, 2/10);
  slp = slp - filtfilt(b, a, slp);
elseif biff == 3;
  [b, a] = butter(9, 2/9);
  slp = filtfilt(b, a, slp);
end
slp = cosweight(slp, lat);
[ntim, nlat, nlon] = size(slp);
slp = reshape(slp, ntim, nlat*nlon);
if biff == 1;
  hello('RAW EOF');
  [rlam, rlds, rpcs, rper] = eof_dan(slp);
elseif biff == 2;
  hello('HP EOF');
  [hlam, hlds, hpcs, hper] = eof_dan(slp);
elseif biff == 3;
  hello('LP EOF');
  [llam, llds, lpcs, lper] = eof_dan(slp);
end
end
data
save ML_SLP_eof_shem_-60to-20.mat ...
     rlam rlds rpcs rper ...
     hlam hlds hpcs hper ...
     llam llds lpcs lper lat lon lims
back