Documentation of HC_clim


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


Help text

 print -dps2 HC0-80_80-270.ps

Cross-Reference Information

This script calls

Listing of script HC_clim


clear
lims = [-0.1 360 -90 90];
tim = 101:1000;
h1 = getheat(lims, 1:3, tim);
h2 = getheat(lims, 4:7, tim);
h1 = squeeze(mean(h1));
h2 = squeeze(mean(h2));
[lat, lon, depth, lm] = getll('temp', lims);
default_global; FRAME = [0 360 -90 90];
figure(1); 
spthes(1);
  gcont(1e-9*h1, 1);
  dc2(lm, 0.55, 1);
  title('0-80m Heat Content');
  xlabel('Contour Interval:  1 \times 10^9 J m^-^2');
spthes(2);
  gcont(1e-9*h2, 1.5);
  dc2(lm, 0.55, 1);
  title('80-270m Heat Content');
  xlabel('Contour Interval:  1.5 \times 10^9 J m^-^2');
cd ~/Thesis/Clim

h2 = h1+h2;

h1 = getnc('temp', lims, 1, tim);

h1 = squeeze(mean(h1));

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

default_global; FRAME = [0 360 -90 90];
figure(2); fo(1); clf; 
spthes(1);
gcont(h1, 2);
  dc2(lm, 0.55, 1);
  title('SST');
  xlabel('Contour Interval:  2 K');

spthes(2);
  gcont(1e-9*h2, 2.5);
  dc2(lm, 0.55, 1);
  title('0-270m Heat Content');
  xlabel('Contour Interval:  2.5 \times 10^9 J m^-^2');

cd ~/Thesis/Clim
% print -dps2 SST_HC0-270.ps