Documentation of eq_depth_terms

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

Cross-Reference Information

Listing of script eq_depth_terms

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP9_detrend_L1-7_EOF_yr101-1000.mat
lims = [101 299 4.5 15];
lags = -8:2:6;
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Chap5
[ubtp, uptb] = vert_ubtp_uptb(pcs, lims, lags);
[vbtp, wbtp] = vert_vdtdy_wdtdz(pcs, lims, lags);
[vptb, wptb] = vert_vptb_wptb(pcs, lims, lags);
[lat, lon, depth] = getll('u', lims);
var1 = 1e8*ubtp;
var2 = 1e8*uptb;
figure(1); fo(1); clf;
[lat, lon, depth] = getll('v', lims);
var1 = 1e8*vbtp;
var2 = 1e8*vptb;
figure(2); fo(1); clf;
[lat, lon, depth] = getll('wl', lims);
var1 = 1e8*wbtp;
var2 = 1e8*wptb;
figure(3); fo(1); clf;
cint = 0.1; 
clev1 = sort([-cint:-cint:min(min(min(min(var1)))) ...
	  cint:cint:max(max(max(max(var1))))]);
clev2 = sort([-cint:-cint:min(min(min(min(var2)))) ...
       	  cint:cint:max(max(max(max(var2))))]);
for i = 1:6;
  tem = squeeze(mean2(shiftdim(var1(i,:,:,:), 2)))';
  subplot(6,2,2*i-1);
    pncont(lon, -depth, tem, clev1, 0, 'k');
    axis([105 290 -300 0]);
    set(gca, 'YTick', -500:100:0, 'YTickLabel', 500:-100:0);
    ylabel(['Lag = ' num2str(lags(i))]);
  tem = squeeze(mean2(shiftdim(var2(i,:,:,:), 2)))';
  subplot(6,2,2*i);
    pncont(lon, -depth, tem, clev2, 0, 'k');
    axis([105 290 -300 0]);
    set(gca, 'YTick', -500:100:0, 'YTickLabel', 500:-100:0);
end