Documentation of plot_lag_regressions

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

Help text

    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);

Cross-Reference Information

Listing of script plot_lag_regressions

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -10:-2;
sst = getnc('temp', lims, lev, tim);
[lat, lon] = getll('temp', lims);
taux = getnc('taux', lims, 1, tim);
[latx, lonx] = getll('taux', lims);
[treg, tcoef] = regress_eof(sst, pcs, lags);
[xreg, xcoef] = regress_eof(taux, pcs, lags);
default_global; FRAME = [105 295 -40 60];
ind = 4;
cint = 0.05; clev = [-2:cint:-cint cint:cint:2];
for i = 1:4;
  subplot(4,2,2*i-1);
    gcont(treg(i+ind,:,:), clev);
    dc2(treg(i,:,:));
    ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);

cint = 1;  clev = [-12:cint:-cint cint:cint:12];
for i = 1:4;
  subplot(4,2,2*i);
    XAX = lonx; YAX = latx;
    gcont(100*xreg(i+ind,:,:), clev);
    XAX = lon; YAX = lat;
    dc2(treg(i,:,:));
    XAX = lonx; YAX = latx;
%    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
    ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval:  ' num2str(1*cint) 'x10^-^3 N m^-^2 std^-^1']);
subplot(4,2,1);
  title('LP10:  12.5m TEMP');
subplot(4,2,2);
  title('LP10:  ZONAL WIND STRESS');


cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots


%  Heat content and ocean currents

clear

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;

lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -6:-3;

hc = getheat(lims, 4:7, tim);
[lat, lon] = getll('temp', lims);
u = getheat(lims, 4:7, tim, 'u');
v = getheat(lims, 4:7, tim, 'v');
[latv, lonv] = getll('v', lims);

[hreg, hcoef] = regress_eof(hc, pcs, lags);
[ureg, ucoef] = regress_eof(u, pcs, lags);
[vreg, vcoef] = regress_eof(v, pcs, lags);

default_global;  FRAME = [105 295 -40 60];

ind = 0;
cint = 0.15; clev = [-3:cint:-cint cint:cint:3];
for i = 1:4;
  subplot(4,2,2*i-1);
    gcont(1e-8*hreg(i+ind,:,:), clev);
    dc2(hreg(i,:,:));
%    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
    ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval:  ' num2str(cint) ' J m¯^2 std^-^1']);
yl = get(gca, 'YLim');

FRAME = [105 295 -30 60];
ind = 0;
for i = 1:4;
  subplot(4,2,2*i);
    XAX = lonv; YAX = latv;
    gquiv(ureg(i+ind,:,:), vreg(i+ind,:,:), -.15, 2, 'cm s^-^1 std^-^1');
    XAX = lon; YAX = lat;
    dc2(hreg(i,:,:));
    set(gca, 'YLim', yl);
%    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
    ylabel(['Lag = ' num2str(lags(i+ind))]);
end

subplot(4,2,1);
  title('LP10:  80-270m HEAT CONTENT');
subplot(4,2,2);
  title('LP10:  80-270m OCEAN VELOCITY');

cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots



%  Heat Flux and ocean currents

clear

cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcs = -1*pcs;

lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
lags = -6:-3;

hc = getnc('heat', lims, 1, tim);
[lat, lon] = getll('temp', lims);
u = getheat(lims, 1:3, tim, 'u');
v = getheat(lims, 1:3, tim, 'v');
[latv, lonv] = getll('v', lims);

[hreg, hcoef] = regress_eof(hc, pcs, lags);
[ureg, ucoef] = regress_eof(u, pcs, lags);
[vreg, vcoef] = regress_eof(v, pcs, lags);

default_global;  FRAME = [105 295 -40 60];

ind = 0;
cint = 0.5; clev = [-5:cint:-cint cint:cint:5];
for i = 1:4;
  subplot(4,2,2*i-1);
    gcont(hreg(i+ind,:,:), clev);
    dc2(hreg(i,:,:));
%    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
    ylabel(['Lag = ' num2str(lags(i+ind))]);
end
xlabel(['Contour Interval:  ' num2str(cint) ' W m^-^2 std^-^1']);
yl = get(gca, 'YLim');

FRAME = [105 295 -30 60];
ind = 0;
for i = 1:4;
  subplot(4,2,2*i);
    XAX = lonv; YAX = latv;
    gquiv(ureg(i+ind,:,:), vreg(i+ind,:,:), -.25, 2, 'cm s^-^1 std^-^1');
    XAX = lon; YAX = lat;
    dc2(hreg(i,:,:));
    set(gca, 'YLim', yl);
%    color_shade(squeeze(tcoef(i+ind,:,:).^2), 0.1, [.9 .9 .9]);
    ylabel(['Lag = ' num2str(lags(i+ind))]);
    if i ~= 4;  xlabel('');  end;
end

subplot(4,2,1);
  title('LP10:  NET HEAT FLUX (Pos. Down)');
subplot(4,2,2);
  title('LP10:  0-80m OCEAN VELOCITY');

cd /home/disk/tao/dvimont/matlab/CSIRO/IUGG_talk/Plots