Documentation of sst_tau_hc_fig

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

Help text

ct = getnc('temp', [180 270 -6 6], 1, tim);
ct = squeeze(mean(mean(ct, 2), 3));
[b, a] = butter(6, 2/10);
pcsl = filtfilt(b, a, ct);
pcsh = ct - pcsl;

Cross-Reference Information

Listing of script sst_tau_hc_fig

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Thesis/Data
load LP10_detrend_L1-7_yr101-1000.mat; pcsl = -1*pcs;
load HP10_detrend_L1-7_EOF_yr101-1000.mat; pcsh = -1*pcs;
lims = [100 300 -65 65];
tim = 101:1000;
lev = 1;
sst = getnc('temp', lims, lev, tim);
[lat, lon] = getll('temp', lims);

[treg, tcoef] = regress_eof(sst, pcsh, 0);
tem = zeros(size(tcoef));
tem(tcoef.^2 >= 0.25) = NaN;

cint = 0.05;  clev = [-1:cint:-cint cint:cint:1];
default_global; FRAME = [105 290 -60 60];
subplot(3,2,1);
  gcont(treg, clev);
  dc2(treg);
  shade_ccoef(tem);
  title('HP10:  12.5m Temp');
  xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);

[treg, tcoef] = regress_eof(sst, pcsl, 0);
tem = zeros(size(ccoef));
tem(tcoef.^2 >= 0.25) = NaN;

subplot(3,2,2);
  gcont(treg, clev);
  dc2(treg);
  shade_ccoef(tem);
  title('LP10:  12.5m Temp');
  xlabel(['Contour Interval:  ' num2str(cint) ' K std^-^1']);



[taux, tauy] = getnc('taux', 'tauy', lims, 1, tim);

[xreg, xcoef] = regress_eof(taux, pcsh, 0);
[yreg, ycoef] = regress_eof(tauy, pcsh, 0);
[lat, lon] = getll('taux', lims);
tem = zeros(size(xcoef));
tem(xcoef.^2 >= 0.25) = NaN;

[latt, lont] = getll('temp', lims);


default_global; FRAME = [105 290 -60 60];
subplot(3,2,3);
  gquiv(xreg, yreg, -.05, 2, 'dynes cm^-^1 std^-^1');
  XAX = lont; YAX = latt;
  dc2(treg);
  XAX = lon; YAX = lat;
  shade_ccoef(tem);
  title('HP10:  Wind Stress');

[xreg, xcoef] = regress_eof(taux, pcsl, 0);
[yreg, ycoef] = regress_eof(tauy, pcsl, 0);
[lat, lon] = getll('taux', lims);
tem = zeros(size(xcoef));
tem(xcoef.^2 >= 0.25) = NaN;

default_global; FRAME = [105 290 -60 60];
subplot(3,2,4);
  gquiv(xreg, yreg, -.05, 2, 'dynes cm^-^1 std^-^1');
  XAX = lont; YAX = latt;
  dc2(treg);
  XAX = lon; YAX = lat;
  shade_ccoef(tem);
  title('LP10:  Wind Stress');



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

[hreg, ccoef] = regress_eof(hc, pcsh, 0);
tem = zeros(size(ccoef));
tem(ccoef.^2 >= 0.25) = NaN;

cint = 0.3;  clev = [-9:cint:-cint cint:cint:9];
default_global; FRAME = [105 290 -60 60];
subplot(3,2,5);
  gcont(1e-8*hreg, clev);
  XAX = lont; YAX = latt;
  dc2(treg);
  XAX = lon; YAX = lat;
  shade_ccoef(tem);
  title('HP10:  80-270m Heat Content');
  xlabel(['Contour Interval:  ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']);

[hreg, ccoef] = regress_eof(hc, pcsl, 0);
tem = zeros(size(ccoef));
tem(ccoef.^2 >= 0.25) = NaN;

subplot(3,2,6);
  gcont(1e-8*hreg, clev);
  XAX = lont; YAX = latt;
  dc2(treg);
  XAX = lon; YAX = lat;
  shade_ccoef(tem);
  title('LP10:  80-270m Heat Content');
  xlabel(['Contour Interval:  ' num2str(100*cint) 'x10^6 J m^-^2 std^-^1']);

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