more off;
load struc.out;
clf;
name = {
't', 'ps'
'T', 'K'
'potential energy / atom', 'J/mol'
'internal energy / atom', 'J/mol'
'Hamiltonian / atom', 'J/mol'
'structure factor', ''
'Mean squared displacement', 'A^2'
'H_{1x}', 'A'
'H_{2y}', 'A'
'H_{3z}', 'A'
'H_{2z}', 'A'
'H_{1z}', 'A'
'H_{1y}', 'A'
'volume', 'A^3'
'\sigma_{xx}', 'MPa'
'\sigma_{yy}', 'MPa'
'\sigma_{zz}', 'MPa'
'\sigma_{yz}', 'MPa'
'\sigma_{xz}', 'MPa'
'\sigma_{xy}', 'MPa'
'internal pressure', 'MPa'
'internal von Mises shear stress', 'MPa'
'C_{xxxx}', 'MPa'
'C_{xxyy}', 'MPa'
'C_{xxzz}', 'MPa'
'C_{xxyz}', 'MPa'
'C_{xxxz}', 'MPa'
'C_{xxxy}', 'MPa'
'C_{yyyy}', 'MPa'
'C_{yyzz}', 'MPa'
'C_{yyyz}', 'MPa'
'C_{yyxz}', 'MPa'
'C_{yyxy}', 'MPa'
'C_{zzzz}', 'MPa'
'C_{zzyz}', 'MPa'
'C_{zzxz}', 'MPa'
'C_{zzxy}', 'MPa'
'C_{yzyz}', 'MPa'
'C_{yzxz}', 'MPa'
'C_{yzxy}', 'MPa'
'C_{xzxz}', 'MPa'
'C_{xzxy}', 'MPa'
'C_{xyxy}', 'MPa'
};

fprintf (1, '\n');
n = size (name,1);
for i = 1:n,
  if (size(name{i,1},2) > 21) buf = '\t';
  elseif (size(name{i,1},2) > 13) buf = '\t\t';
  elseif (size(name{i,1},2) > 1) buf = '\t\t\t';
  else buf = '\t\t\t\t';
  end;
  fprintf (1, ['%d. %s' buf], i, name{i,1});
  if (mod(i,2) == 0) fprintf (1, '\n'); end;
end;
fprintf (1, '\n');

c = 1;
[minx minxd] = min(struc(:,c));
[maxx maxxd] = max(struc(:,c));
centerd = round (maxxd * 2 / (1+sqrt(5)) + minxd * (sqrt(5)-1) /(1+sqrt(5))); 
ss = input (['Select quantities to plot with ' name{c,1} ...
      ' [' name{c,2} ']: '], 's');
eval(['d = [' ss '];']);

clear rcc;
for i=1:size(d,2)
  rc;
  set (plot(struc(:,c), struc(:,d(i)), 'o'), 'Color', rcc);
  if (size(d,2) > 1) 
    set (text(struc(centerd,c), ...
	struc(centerd,d(i))+0.035, ...
	[name{d(i),1} ' [' name{d(i),2} ']']), 'Color', rcc);
  end;
  hold on;
end;

if (size(d,2) == 1)
  ylabel([name{d(1),1} ' [' name{d(1),2} ']']);
end;
xlabel([name{c,1} ' [' name{c,2} ']']);

v = axis;
axis ([minx maxx v(3) v(4)]);
