Documentation of plot_eq_currents

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

Help text

tser = lpcs(:,1); lags = -8:2:6; biff = 2;

Cross-Reference Information

Listing of script plot_eq_currents

clear
loadpcs;
lims = [100 300 -35 35];
tim = 101:1000;
ut = getheat(lims, 1:3, tim, 'u');
vt = getheat(lims, 1:3, tim, 'v');
ub = getheat(lims, 4:7, tim, 'u');
vb = getheat(lims, 4:7, tim, 'v');
wt = getheat(lims, 1:3, tim, 'wl');
wb = getheat(lims, 4:7, tim, 'wl');
[lat, lon] = getll('u', lims);
[lat1, lon1, tem, lm] = getll('temp', lims);
[lat2, lon2] = getll('wl', lims);
tser = hpcs(:,1); lags = -1:1; biff = 1;
nlag = length(lags);

utr = regress_eof(ut, tser, lags);
vtr = regress_eof(vt, tser, lags);
ubr = regress_eof(ub, tser, lags);
vbr = regress_eof(vb, tser, lags);
wtr = regress_eof(wt, tser, lags);
wbr = regress_eof(wb, tser, lags);

default_global; FRAME = [105 240 -30 30];
figure(biff); fo(1); clf;

for i = 1:nlag
  XAX = lon; YAX = lat;
  subplot(3,2,2*i-1);
    gquiv(utr(i,:,:), vtr(i,:,:), -1, 2, '');
    hold on;
    XAX = lon2; YAX = lat2;
    gcont(wtr(i,:,:), 0.05);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
    ylabel(['Lag = ' num2str(lags(i))]);
  XAX = lon; YAX = lat;
  subplot(3,2,2*i-1);
    gquiv(ubr(i,:,:), vbr(i,:,:), -.5, 2, '');
    hold on;
    XAX = lon2; YAX = lat2;
    gcont(wbr(i,:,:), 0.025);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
end

figure(biff+2); fo(1); clf;

for i = 1:nlag
  XAX = lon2; YAX = lat2;
  subplot(6,2,2*i-1);
    hold on;
    gcont(wtr(i,:,:), 0.05);
%    XAX = lon1; YAX = lat1;
%    dc2(lm);
%    ylabel(['Lag = ' num2str(lags(i))]);
  XAX = lon2; YAX = lat2;
  subplot(6,2,2*i);
    gcont(wbr(i,:,:), 0.015);
    XAX = lon1; YAX = lat1;
    dc2(lm);
end



tser = hpcs(:,1); lags = -1:1; biff = 1;
nlag = length(lags);

ubr = regress_eof(ub, tser, lags);
vbr = regress_eof(vb, tser, lags);
wbr = regress_eof(wb, tser, lags);

tser = lpcs(:,1); lags = -4:3:2; biff = 1;
nlag = length(lags);

ubr2 = regress_eof(ub, tser, lags);
vbr2 = regress_eof(vb, tser, lags);
wbr2 = regress_eof(wb, tser, lags);

figure(1); fl; clf;
default_global; FRAME = [105 270 -30 30];

for i = 1:nlag;
  XAX = lon; YAX = lat;
  subplot(3,2,2*i-1);
%    gquiv(ubr(i,:,:), vbr(i,:,:), -0.5, 2, '');
%    hold on;
    XAX = lon2; YAX = lat2;
    gcont(wbr(i,:,:), 0.025);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
    ylabel(['Lag = ' num2str(lags(i))]);
end
for i = 1:nlag;
  XAX = lon; YAX = lat;
  subplot(3,2,2*i);
%    gquiv(ubr2(i,:,:), vbr2(i,:,:), -0.5, 2, '');
%    hold on;
    XAX = lon2; YAX = lat2;
    gcont(wbr2(i,:,:), 0.015);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
    ylabel(['Lag = ' num2str(lags(i))]);
end

subplot(3,2,1);
title('< 80-270m W, HPPC1 >');
subplot(3,2,2);
title('< 80-270m W, LPPC1 >');



cd ~/Thesis/Talk
print -dps2 HP10_LP9_W_eq_lags.ps



%  Look at wind stress

tx = getnc('taux', lims, 1, tim);
ty = getnc('tauy', lims, 1, tim);

tser = hpcs(:,1); lags = -1:1; biff = 1;
nlag = length(lags);

ubr = regress_eof(tx, tser, lags);
vbr = regress_eof(ty, tser, lags);
wbr = regress_eof(wb, tser, lags);

tser = lpcs(:,1); lags = -4:3:2; biff = 1;
nlag = length(lags);

ubr2 = regress_eof(tx, tser, lags);
vbr2 = regress_eof(ty, tser, lags);
wbr2 = regress_eof(wb, tser, lags);

figure(1); fl; clf;
default_global; FRAME = [105 270 -30 30];

for i = 1:nlag;
  XAX = lon; YAX = lat;
  subplot(3,2,2*i-1);
    gquiv(ubr(i,:,:), vbr(i,:,:), -0.05, 1, '');
    hold on;
    XAX = lon2; YAX = lat2;
%    gcont(wbr(i,:,:), 0.025);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
    ylabel(['Lag = ' num2str(lags(i))]);
end
for i = 1:nlag;
  XAX = lon; YAX = lat;
  subplot(3,2,2*i);
    gquiv(ubr2(i,:,:), vbr2(i,:,:), -0.05, 1, '');
    hold on;
    XAX = lon2; YAX = lat2;
%    gcont(wbr2(i,:,:), 0.015);
    hold off;
    XAX = lon1; YAX = lat1;
    dc2(lm);
    xlabel('');
    ylabel(['Lag = ' num2str(lags(i))]);
end