
M=read("thermoR134A.txt",-1,12);
T=M(:,1);
P=M(:,2)/100;
vliq=M(:,3); 
vgaz=M(:,4);
hliq=M(:,7);  
hgaz=M(:,9);
sliq=M(:,10);
sgaz=M(:,11);
			
plotPT=scf();

plot2d(T,P,style=5,logflag="nl");
xgrid(color("grey"));
xtitle("P (bar) - T (C)","T","P");

plotPV=scf();

plot2d(vliq,P,style=5,logflag="ll")
plot2d(vgaz,P,style=5)
xgrid(color("grey"));
xtitle("P (bar) - v (m^3/kg)","v","P");
			
plotPH=scf();

plot2d(hliq,P,style=5,logflag="nl");
plot2d(hgaz,P,style=5);
xgrid(color("grey"));
xtitle("P (bar)-H (kJ/kg)","H","P")
			
plotTS=scf();

plot2d(sliq,T,style=5,logflag="nn");
plot2d(sgaz,T,style=5);
xgrid(color("grey"));
xtitle("T (C)-S (kJ/K/kg)","S","T");
			
