Global Index (short | long) | Local contents | Local Index (short | long)
[temp, lat, lon, depth, middepth] = getheat(lims, lev, tim);
This script calls | |
---|---|
clear varn = 'temp'; cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data load LP10_detrend_L1-7_yr101-1000.mat lev = 4:7; tim = 101:1000; lims = [105 300 0 75]; lags = -8:5; nfrm = lags; temp = getnc(varn, lims, 1, tim); [lat, lon, depth] = getll(varn, lims); treg = regress_eof(temp, pcs, nfrm); % Plot the various phases default_global; FRAME = [105 295 5 65]; cint = 1.5; clev = [-15:cint:-cint cint:cint:15]; for i = 1:7; subplot(7,2,2*i-1); gcont(-1e0*treg(i,:,:), clev); dc2(treg(i,:,:)); % dc ylabel(['Lag = ' num2str(lags(i))]); subplot(7,2,2*i); gcont(-1e0*treg(i+7,:,:), clev); dc2(treg(i+7,:,:)); % dc ylabel(['Lag = ' num2str(lags(i+7))]); end for i = 13:14 subplot(7,2,i); cint2 = 1e0*cint; % xlabel(['Contour Interval = ' num2str(cint2) ' K std^-^1']); % xlabel(['Contour Interval = ' num2str(cint2) 'x10^7 J m^-^2 std^-^1']); xlabel(['Contour Interval = ' num2str(cint2) 'x10^-^3 N m^-^2 std^-^1']); end for i = 1:2; subplot(7,2,i); title(['LP10: < HC80-270, PC1 > ' num2str(round(per(1))) ... '% Var Expl.']); title('LP10: < 12.5m Temp, PC1 >'); title('LP10: < TAUX, PC1 >'); end cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap4/Plots %print -dps2 HC270m_PC1_ML.ps