Documentation of heat_content_clim

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

Help text

cwat = 4.218e3;  %  heat capacity of liquid water, J/(kg K)
rhowat = 1e3;    %  density of liquid water, kg/m^3
tem = 1e-9*cwat * rhowat * squeeze(mean(heat));

Cross-Reference Information

Listing of script heat_content_clim

clear
cd /home/disk/tao/dvimont/matlab/CSIRO/Heat
lev = 1:3;
tim = 101:550;
ctlim = [-0.1 360 -90 90];
[heat, lat, lon, depth, middepth] = getheat(lev, tim, ctlim);
tem = 1e-8*squeeze(mean(heat));

default_global
figure(1);
  sp(1);
  gcont(tem, [-10:5:100]);
  dc
  title('0:80m Heat Content, years 101-550');
  xlabel('Contour Interval:  5x10^8 J m^-^2');

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

lev = 1:7;
tim = 101:550;
ctlim = [-0.1 360 -90 90];

[heat, lat, lon, depth, middepth] = getheat(lev, tim, ctlim);
%cwat = 4.218e3;  %  heat capacity of liquid water, J/(kg K)
%rhowat = 1e3;    %  density of liquid water, kg/m^3
%tem = 1e-9*cwat * rhowat * squeeze(mean(heat));
tem = 1e-8*squeeze(mean(heat)) - tem;

default_global
figure(1);
  sp(2);
  gcont(tem, [-1000:10:1000]);
  dc
  title('80:270m Heat Content, years 101-550');
  xlabel('Contour Interval:  1x10^9 J m^-^2');

cd /home/disk/tao/dvimont/matlab/CSIRO/Plots_Clim
print -dps2 heat_L1-3_L4-7_yr101-550.ps