%
% clear previous plot
clg
%
% number of t and z values
nt=21;
nz=31;
%
% t grid
t=(0:0.05:1)';
%
% z grid - zl = 30
zl=30;
z30=(0:1:zl)';
%
% z grid - zl = 90
zl=90;
z90=(0:3:zl)';
%
% dependent variable
%
% zl = 30
load h5a.out;
for j=1:nt;
for i=1:nz;
    u30(i,j)=h5a((j-1)*nz+i,2);
end;
end;
%
% zl = 90
load h5b.out;
for j=1:nt;
for i=1:nz;
    u90(i,j)=h5b((j-1)*nz+i,2);
end;
end;
%
% plot u(0,z,t) vs t and z 
%
% zl = 30
figure(1);
mesh(t,z30,u30);
view(-135,45);
hold on;
%
% label three axes
xlabel(' t');
ylabel(' z');
zlabel(' u(0,z,t)');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV1');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV2');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV3');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV4');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV5');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV6');
%title(' Graetz problem with constant wall temperature, zl = 30, DERV7');
 title(' Graetz problem with constant wall temperature, zl = 30, DERV8');
%
% print output
% print h5p1a.eps -deps
% print h5p2a.eps -deps
% print h5p3a.eps -deps
% print h5p4a.eps -deps
% print h5p5a.eps -deps
% print h5p6a.eps -deps
% print h5p7a.eps -deps
  print h5p8a.eps -deps
%
% zl = 90
figure(2);
mesh(t,z90,u90);
view(-135,45);
hold on;
%
% label three axes
xlabel(' t');
ylabel(' z');
zlabel(' u(0,z,t)');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV1');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV2');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV3');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV4');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV5');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV6');
%title(' Graetz problem with constant wall temperature, zl = 90, DERV7');
 title(' Graetz problem with constant wall temperature, zl = 90, DERV8');
%
% print output
% print h5p1b.eps -deps
% print h5p2b.eps -deps
% print h5p3b.eps -deps
% print h5p4b.eps -deps
% print h5p5b.eps -deps
% print h5p6b.eps -deps
% print h5p7b.eps -deps
  print h5p8b.eps -deps
%
% send eps file to ihd09 Postscript printer
%!op  -s 1 -d cs1 -q ihd109 h5p1a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p1b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p2a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p2b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p3a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p3b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p4a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p4b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p5a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p5b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p6a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p6b.eps
%!op  -s 1 -d cs1 -q ihd109 h5p7a.eps
%!op  -s 1 -d cs1 -q ihd109 h5p7b.eps
 !op  -s 1 -d cs1 -q ihd109 h5p8a.eps
 !op  -s 1 -d cs1 -q ihd109 h5p8b.eps