function [rho S Cp C ] = Constants()

rho = 1025;
S = 2.0;
C = 12.34*.0254;
Cp =(2.8750-0)*.0254;


end

