Documentation of hc_eof

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

Help text

save heat_pcs.mat rlam rlds rpcs rper lims lev lat lon ...
save heat_pcs_eq.mat rlam rlds rpcs rper lims lev lat lon ...
save heat_pcs_shem.mat rlam rlds rpcs rper lims lev lat lon ...

Cross-Reference Information

Listing of script hc_eof

clear
lims = [110 285 20 70];
tim = 101:1000;
lev = 1:7;
[heat, lat, lon, depth, middepth] = getheat(lims, lev, tim, 'temp');
for hplp = 0:2;
  if hplp == 0;
    heat2 = heat;
  elseif hplp == 1;
    [b, a] = butter(9, 2/10);
    heat2 = heat - filtfilt(b, a, heat);
  elseif hplp == 2;
    [b, a] = butter(9, 2/9);
    heat2 = filtfilt(b, a, heat);
  end
  hplp
  [heat2, clim] = remove_mean(heat2);
  heat2 = detrend(heat2);
  kp = find(~isnan(clim));
  heat2 = cosweight(heat2, lat);
  heat2 = heat2(:, kp);
  if hplp == 0;
    [rlam, rlds, rpcs, rper] = eof_dan(heat2);
  elseif hplp == 1; 
    [hlam, hlds, hpcs, hper] = eof_dan(heat2);
  elseif hplp == 2; 
    [llam, llds, lpcs, lper] = eof_dan(heat2);
  end
end
cd ~/matlab/CSIRO/New_calcs/data
save heat_pcs_nhem.mat rlam rlds rpcs rper lims lev lat lon ...
                  hlam hlds hpcs hper ...
                  llam llds lpcs lper