%%%%%%%%%%%%%%%%%%%%%%%%%%
% Function to plot g(r)s %
%%%%%%%%%%%%%%%%%%%%%%%%%%
function gr()
global gr T D;

if exist('gr')==0
 load gr.out;
end;

ss=input('T = ','s');
if isstr(ss) T = ss; end;
ss=input('Density = ','s');
if isstr(ss) D = ss; end;

minx = min(gr(:,2));
maxx = max(gr(:,2));
clf;

%%%%%%%%%%%%%%%%
% PLOT G_OO(R) %
%%%%%%%%%%%%%%%%
plot(gr(:,2),gr(:,3));
xlabel('Angstrom');
title(['Goo(r) at T= ' T 'K  Density = ' D]);
CHG_AXIS(minx,maxx);
PRINT;

%%%%%%%%%%%%%%%%
% PLOT G_OH(R) %
%%%%%%%%%%%%%%%%
plot(gr(:,2),gr(:,4));
xlabel('Angstrom');
title(['Goh(r) at T= ' T 'K  Density = ' D]);
CHG_AXIS(minx,maxx);
PRINT;

%%%%%%%%%%%%%%%%
% PLOT G_OO(R) %
%%%%%%%%%%%%%%%%
plot(gr(:,2),gr(:,5));
xlabel('Angstrom');
title(['Ghh(r) at T= ' T 'K  Density = ' D]);
CHG_AXIS(minx,maxx);
PRINT;

end;
