Documentation of gu_philander5

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

Help text

lag1 = -15:-10; lag2 = -9:-4; lp = length(lag1);

Cross-Reference Information

Listing of script gu_philander5

clear
loadpcs;
lags = [-3:8];
for biff = 2:2;
if biff == 1; lims = [180 210 -50 50];
elseif biff == 2; lims = [165 195 -50 50];
elseif biff == 3; lims = [150 180 -50 50];
elseif biff == 4; lims = [130 150 -50 50]; end
temp1 = getnc('temp', lims, 1:10, 101:1000);
[lat, lon, depth, lm] = getll('temp', lims);
temp1 = squeeze(mean2(shiftdim(temp1, 3)));
[treg, tcoef] = regress_eof(temp1, lpcs, lags);
cint = 0.025;
lag1 = -3:2; lag2 = 3:8; lp = length(lag1);

figure(biff); clf; fo(1);
for i = 1:lp;
  ind = find(lags == lag1(i));
  sptalk(6,2,2*i-1);
  tem = squeeze(treg(ind,:,:));
  clev = sort([0:-cint:min(min(tem)) cint:cint:(max(max(tem))) -100 100]);
  pncont(lat, -1*depth, tem, clev, 0, 'k');
  axis([-50 50 -500 0]);
  set(gca, 'XTick', -45:15:45);
end
for i = 1:lp;
  ind = find(lags == lag2(i));
  sptalk(6,2,2*i);
  tem = squeeze(treg(ind,:,:));
  clev = sort([0:-cint:min(min(tem)) cint:cint:(max(max(tem))) -100 100]);
  pncont(lat, -1*depth, tem, clev, 0, 'k');
  axis([-50 50 -500 0]);
  set(gca, 'XTick', -45:15:45);
end

sptalk(6,2,1); title([num2str(lims)]);

end

%  Look at temp tend terms

clear

cd ~/matlab/CSIRO/Thesis/Data
load gu_philand.mat
[latu, lonu, depthu] = getll('u', lims);
[latw, lonw, depthw] = getll('wl', lims);

order = ['ubtp'; 'uptb'; 'vbtp'; 'vptb'; 'wbtp'; 'wptb'];
llord = ['uuw'];

lims = [180 210 -45 45];
[xku, yku] = keep_var(lims, lonu, latu);
[xkw, ykw] = keep_var(lims, lonw, latw);

temp = getnc('temp', lims, 1:10, [min(tim-1) tim]);
ntim = length(tim)+1;
dtdt = (temp(3:ntim,:,:,:) - temp(1:(ntim-2),:,:,:)) ./ (2*3600*24*365);
dtdt = detrend(dtdt);

loadpcs; pcs = lpcs;
[dtreg, dtcoef] = regress_eof(dtdt, -1*pcs(1:(ntim-2),1), lags);
[treg, tcoef] = regress_eof(temp(2:ntim,:,:,:), -1*pcs, lags);

[latt, lont, deptht] = getll('temp', lims);
[latu, lonu, depthu] = getll('u', lims);
[latw, lonw, depthw] = getll('wl', lims);

ubtp2 = squeeze(mean2(shiftdim(ubtp(:,:,yku,xku), 3)));
vbtp2 = squeeze(mean2(shiftdim(vbtp(:,:,yku,xku), 3)));
wbtp2 = squeeze(mean2(shiftdim(wbtp(:,:,yku,xku), 3)));
uptb2 = squeeze(mean2(shiftdim(uptb(:,:,yku,xku), 3)));
vptb2 = squeeze(mean2(shiftdim(vptb(:,:,yku,xku), 3)));
wptb2 = squeeze(mean2(shiftdim(wptb(:,:,yku,xku), 3)));

lplot = [-10:2:-2]; lp = length(lplot);

for j = 1:3;
  figure(j); fo(1);
  eval(['tem1 = ' order(2*j-1,:) ';']);
  eval(['tem2 = ' order(2*j,:) ';']);
  eval(['lat = lat' llord(j) '; depth = depth' llord(j) ';']);
  for i = 1:5;
    lind = find(lags = lagplot(i));
    sptalk(5,2,2*i-1);
    tem = 1e9*squeeze(tem1(lind,:,:));
contour(

NOT FINISHED