Global Index (short | long) | Local contents | Local Index (short | long)
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 ...
This script calls | |
---|---|
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