function [p,T]=rhstat(r,h)
% compute reheat steam properties
rh=r*h;
p=27.06-1019.1*r-(1.735e-2)*h+1.228*rh;
T=-2013.4+189.65*r+1.963*h-(9.305e-2)*rh+459.67;
end

