Global Index (short | long) | Local contents | Local Index (short | long)
load gr_pcs.mat lam lds pcs per lims cthp gr grpcs = pcs; for i = 1:900; ind = 12*(i-1)+[1:12]; tem(i,:) = mean(grpcs(ind,:)); end pcs = tem;
This script calls | |
---|---|
clear cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data load LP9_detrend_L1-7_EOF_yr101-1000.mat cd ~/matlab/CSIRO/Heat/Old_routines nfrm = -12:2:12; lags = nfrm; tim = 101:1000; lev = 1:3; lims = [100 305 -35 35]; tit = 'LP9'; varn1 = 'taux'; varn2 = 'tauy'; [var1, lat, lon] = getheat(lims, 1:3, tim); %var1 = getnc(varn1, lims, 1, tim); %[lat, lon, depth] = getll(varn1, lims); [var2, lat2, lon2] = getheat(lims, 4:7, tim); %var2 = getnc(varn2, lims, 1, tim); %[lat2, lon2, depth2] = getll(varn2, lims); %lev = 4:7; %[var1] = getheat(lims, lev, tim, 'u'); %[var2, lat2, lon2] = getheat(lims, lev, tim, 'v'); reg1 = regress_eof(detrend(var1), detrend(pcs), nfrm); reg2 = regress_eof(detrend(var2), detrend(pcs), nfrm); % Plot the data default_global; XAX = lon; YAX = lat; FRAME = [105 300 -30 30]; num = 1 cint = 0.1; figure(num); figure_orient; clf; for i = 8:10; ind = i-7; pos = [.2 .92-ind*.1415 .2920 .1072]; subplot('position', pos); % subplot(6,2,2*i-2); gcont(-1e-8*reg1(i,:,:), cint); dc2(reg1(i,:,:)); ylabel(['Lag = ' num2str(lags(i))]); if ind == 1; % title([tit ' 0-80m HC']); end % xlabel(['Contour Interval: ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']); % xlabel(['Contour Interval: ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']); end end %subplot(6,2,1) % title([tit ' TAUX']); % title([tit ' 0-80m HC']); default_global; XAX = lon2; YAX = lat2; FRAME = [105 300 -30 30]; cint = 0.2; figure(num); figure_orient; for i = 8:10; ind = i-7; pos = [.558 .92-ind*.1415 .2920 .1072]; subplot('position', pos); % subplot(6,2,2*i); gcont(-1e-8*reg2(i,:,:), cint); dc2(reg2(i,:,:)); % ylabel(['Lag = ' num2str(lags(i))]); if ind == 1; % title([tit ' 80-270m HC']); end % xlabel(['Contour Interval: ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']); % xlabel(['Contour Interval: ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']); end end %subplot(6,2,2) % title([tit ' Zonal Wind Stress']); % title([tit ' 80-270m HC']); cd ~/Thesis/Chap4 cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap2/Plots %print -dps2 sfc_T_taux_PC1.ps % Plot HC80, HC80-270, TAUX on one figure figure(1); figure_landscape; default_global; XAX = lon; YAX = lat; FRAME = [110 299 -30 30]; default_global; XAX = lon2; YAX = lat2; FRAME = [110 299 -30 30]; cint = 1; clev = [-12:cint:-cint cint:cint:12]; for i = 1:5; ind = i+1; subplot(5,3,3*i); gcont(1e2*reg2(ind,:,:), clev); dc2(reg1(i,:,:)); ylabel(['Lag = ' num2str(lags(ind))]); if i == 5; xlabel(['Contour Interval: ' num2str(cint) 'x10^-^3 N m^-^2 std^-^1']); % xlabel(['Contour Interval: ' num2str(10*cint) 'x10^7 J m^-^2 std^-^1']); end end subplot(5,3,3) title([tit ' Zonal Wind Stress']); title([tit ' 80-270m HC']); title([tit ' 0-80m HC']); cd ~/Thesis/Chap3 % Plot 60s to 60n data lev = 1:3; [var1] = getheat(lims, lev, tim, 'wl'); [var2, lat2, lon2] = getheat(lims, lev, tim, 'wl'); reg1 = regress_eof(var1, pcs, nfrm); reg2 = regress_eof(var2, pcs, nfrm); get_global; FRAME = [105 300 -60 60]; lm = getnc('temp', lims, 1, 1); [latt, lont] = getll('temp', lims); lind = 0; ind = find(lags == lind); figure(1); figure_landscape; subplot(2,2,1); XAX = lon2; YAX = lat2; gquiv(reg1(ind,:,:), reg2(ind,:,:), -0.75, 1); XAX = lont; YAX = latt; dc2(lm); title('LP9: 0-80m < U'', V'' >'); subplot(2,2,2); XAX = lon2; YAX = lat2; gquiv(reg1(ind,:,:), reg2(ind,:,:), -0.25, 1); XAX = lont; YAX = latt; dc2(lm); title('LP9: 80-270m < U'', V'' >'); cint = 0.05; clev = [-5:cint:-cint cint:cint:5]; subplot(2,2,3); xlabel('Contour Interval: 5x10^-^6 cm s^-1 std^1') XAX = lon2; YAX = lat2; gcont(reg1(ind,:,:), clev); XAX = lont; YAX = latt; dc2(lm); title('LP9: 0-80m W'''); cint = 0.025; clev = [-5:cint:-cint cint:cint:5]; subplot(2,2,4); xlabel('Contour Interval: 2.5x10^-^6 cm s^-1 std^1') XAX = lon2; YAX = lat2; gcont(reg2(ind,:,:), clev); XAX = lont; YAX = latt; dc2(lm); title('LP9: 80-270m W'''); cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap3/Plots % print -dps2 U_V_W_zerolag_regression.ps