Documentation of wbar_dtpdz_hc0-80

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

Help text

  Load PCS:

Cross-Reference Information

Listing of script wbar_dtpdz_hc0-80

clear

cd /home/disk/hayes2/dvimont/csiro/matlab_data/Heat_Content
load LP10_detrend_L1-7_yr101-1000.mat
pcs = -1*pcs;

lev = 1:3;
lims = [150 290 -3 3]; 
nfrm = -15:15;
tim = 101:1000;

cd ~/matlab/CSIRO/Heat/Old_routines
   [wtem, lat_out, lon_out] = ...
           get_wbar_dtprimedz(pcs, lims, nfrm, tim, lev);

   heat = getheat(4:7, tim, lims);
   heat = regress_eof(heat, pcs, nfrm);

   heat1 = getheat(1:3, tim, lims);
   heat1 = regress_eof(heat1, pcs, nfrm);

wtem = squeeze(mean(wtem, 2));
heat = squeeze(mean(heat, 2));
heat1 = squeeze(mean(heat1, 2));

xlab = num2str(['150E'; '180 '; '150W'; '120W'; '90W ']);

figure(2);% figure_landscape;
subplot(1,3,1);
  pncont(lon_out, nfrm, 1e-8*heat1, [-5:.125:5], 'k');
  axis([150 270 -10 10]);
  set(gca, 'XTick', 150:30:270, 'XTickLabel', xlab);
  ylabel('Lag in Years')
  title('0-80m Heat Content');

subplot(1,3,2);
  pncont(lon_out, nfrm, -wtem, [-5:0.5:5], 'k');
  axis([150 270 -10 10]);
  set(gca, 'XTick', 150:30:270, 'XTickLabel', xlab);
  title('0-80m Wbar dT''/dz');

subplot(1,3,3);
  pncont(lon_out, nfrm, 1e-8*heat, [-1:.2:2], 'k');
  axis([150 270 -10 10]);
  set(gca, 'XTick', 150:30:270, 'XTickLabel', xlab);
  title('80-270m Heat Content');

cd /home/disk/tao/dvimont/matlab/CSIRO/Heat/GRL_Figs



%  Color plot

figure(2);% figure_landscape;
subplot(1,3,2);
  contourf(lon_out, nfrm, -wtem, [-2.2:.2:2]);
  hold on;
    [c,h] = contour(lon_out, nfrm, wtem+50, [50 50], 'k');
  hold off;
  set(h, 'linewidth', 2);
  h = colorbar;
  caxis([-2.2 2]);

subplot(1,3,3);
  contourf(lon_out, nfrm, 1e-8*heat, [-.6:.1:1.2]);
  hold on
    [c,h] = contour(lon_out, nfrm, heat, [0 0], 'k');
  hold off
  set(h, 'linewidth', 2);
  h = colorbar;
  caxis([-.6 1.2]);

subplot(1,3,1);
  contourf(lon_out, nfrm, 1e-8*heat1, [-.6:.1:1.2]);
  hold on
    [c,h] = contour(lon_out, nfrm, heat1, [0 0], 'k');
  hold off
  set(h, 'linewidth', 2);
  h = colorbar;
  caxis([-.6 1.2]);