Documentation of fig1

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

Help text

hc = getheat(lims, 1:7, tim);

Cross-Reference Information

Listing of script fig1

clear
loadpcs
lims = [-0.1 360 -90 90];
tim = 101:1000;
temp = getnc('temp', lims, 1, tim);
slp = getflx('psl', lims, tim);
taux = getnc('taux', lims, 1, tim);
tauy = getnc('tauy', lims, 1, tim);
tamp = sqrt(taux.^2 + tauy.^2);
prec = getflx('rnd', lims, tim);
temp = detrend(temp);
slp = detrend(slp);
taux = detrend(taux);
tauy = detrend(tauy);
tamp = detrend(tamp);
prec = detrend(prec);

[reg1, c1] = regress_eof(temp, rpcs, 0);
[reg2, c2] = regress_eof(slp, rpcs, 0);
[reg3, c3] = regress_eof(taux, rpcs, 0);
[reg4, c4] = regress_eof(tauy, rpcs, 0);
[reg5, c5] = regress_eof(prec, rpcs, 0);
[reg6, c6] = regress_eof(tamp, rpcs, 0);

[lat, lon, depth, lm] = getll('temp', lims);

default_global; FRAME = [0 360 -75 75];

figure(1); fo; clf
subplot(3,1,1);
  gcont(reg1, 0.05);
  dc2(lm, 0.3, 1);
  color_shade(squeeze(c1.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.05 K std^-^1');
  xbox = [110 110 300 300 110];
  ybox = [-60 60 60 -60 -60];
  zbox = [100 100 100 100 100];
  hold on;
    h = plot3(xbox, ybox, zbox, '-k');
  hold off
  set(h, 'linewidth', 2);
  title('< SST, PC1 (Unfiltered) >')

subplot(3,1,2);
  gcont(reg2, 0.1);
  dc2(lm, 0.3, -1000);
  color_shade(squeeze(c2.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.1 hPa std^-^1');
  title('< SLP, PC1 (Unfiltered) >')

[lat2, lon2] = getll('taux', lims);
subplot(3,1,3);  cla;
  XAX = lon2; YAX = lat2; tem1 = reg3; tem2 = reg4;
  gquiv(tem1, tem2, -.1, 1, 'dynes cm^-^1 std^-^1');
  XAX = lon; YAX = lat;
  dc2(lm, 0.3, -1000);
  color_shade(squeeze(reg5), .3, [1 0 0]);
  color_shade(squeeze(-reg5), .3, [0 0 1]);
  color_shade(squeeze(reg5), .1, [1 .6 .6]);
  color_shade(squeeze(-reg5), .1, [.6 .6 1]);
  xlabel(['Reference Arrow = 0.1 dynes cm^-^2 std^-^1;'...
        '  Shading => Precip Anom exceeds 0.15 mm day^-^1 std^-^1']);
  title(['< \fontsize{12}\tau\fontsize{10}, PC1 (Unfiltered) > ;'...
         '  < PRECIP, PC1 (Unfiltered) >'])

cd ~/Thesis/Talk
print -dpsc2 SST_SLP_TAU_PC1_Unfiltered_CSIRO.ps


%  Do the same for HPPC1

[reg1, c1] = regress_eof(temp, hpcs, 0);
[reg2, c2] = regress_eof(slp, hpcs, 0);
[reg3, c3] = regress_eof(taux, hpcs, 0);
[reg4, c4] = regress_eof(tauy, hpcs, 0);
[reg5, c5] = regress_eof(prec, hpcs, 0);

default_global; FRAME = [0 360 -75 75];

figure(1); fo; clf;
subplot(3,1,1);
  gcont(reg1, 0.05);
  dc2(lm, 0.3, 1);
  color_shade(squeeze(c1.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.05 K std^-^1');
  title('< SST, HPPC1 >')

subplot(3,1,2);
  gcont(reg2, 0.1);
  dc2(lm, 0.3, -200);
  color_shade(squeeze(c2.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.1 hPa std^-^1');
  title('< SLP, HPPC1 >')

[lat2, lon2] = getll('taux', lims);
subplot(3,1,3);  
  XAX = lon2; YAX = lat2;
  gquiv(reg3, reg4, -.1, 1, 'dynes cm^-^1 std^-^1');
  XAX = lon; YAX = lat;
  dc2(lm, 0.3, -1000);
  color_shade(squeeze(reg5), .1, [1 .6 .6]);
  color_shade(squeeze(-reg5), .1, [.6 .6 1]);
  color_shade(squeeze(reg5), .3, [1 0 0]);
  color_shade(squeeze(-reg5), .3, [0 0 1]);
  xlabel(['Reference Arrow = 0.1 dynes cm^-^2 std^-^1;'...
        '  Shading => Precip Anom exceeds 0.1 mm day^-^1 std^-^1']);
  title('< \fontsize{12}\tau\fontsize{10}, HPPC1 > ;  < PRECIP, HPPC1 >')

cd ~/Thesis/Talk
print -dpsc2 SST_SLP_TAU_HPPC1_CSIRO.ps

%  Low pass filtered

[reg1, c1] = regress_eof(temp, lpcs, 0);
[reg2, c2] = regress_eof(slp, lpcs, 0);
[reg3, c3] = regress_eof(taux, lpcs, 0);
[reg4, c4] = regress_eof(tauy, lpcs, 0);
[reg5, c5] = regress_eof(prec, lpcs, 0);

default_global; FRAME = [0 360 -75 75];

figure(1); fo;
subplot(3,1,1);
  gcont(reg1, 0.05);
  dc2(lm, 0.3, 1);
  color_shade(squeeze(c1.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.05 K std^-^1');
  title('< SST, LPPC1 >')

subplot(3,1,2);
  gcont(reg2, 0.1);
  dc2(lm, 0.3, -100);
  color_shade(squeeze(c2.^2), 0.25, .7*[1 1 1]);
  xlabel('Contour Interval:  0.1 hPa std^-^1');
  title('< SLP, LPPC1 >')

[lat2, lon2] = getll('taux', lims);
subplot(3,1,3);  
  XAX = lon2; YAX = lat2;
  gquiv(reg3, reg4, -.1, 1, 'dynes cm^-^1 std^-^1');
  XAX = lon; YAX = lat;
  dc2(lm, 0.3, -1000);
  color_shade(squeeze(reg5), .1, [1 .6 .6]);
  color_shade(squeeze(-reg5), .1, [.6 .6 1]);
  color_shade(squeeze(reg5), .3, [1 0 0]);
  color_shade(squeeze(-reg5), .3, [0 0 1]);
  xlabel(['Reference Arrow = 0.1 dynes cm^-^2 std^-^1;'...
        '  Shading => Precip Anom exceeds 0.1 mm day^-^1 std^-^1']);
  title('< \fontsize{12}\tau\fontsize{10}, LPPC1 > ;  < PRECIP, LPPC1 >')

cd ~/Thesis/Talk
print -dpsc2 SST_SLP_TAU_LPPC1_CSIRO.ps