function pressure=lithostaticpressure(rho_ice,g,s,z);
%LITHOSTATICPRESSURE - compute the lithostatic pressure
%
%   Usage:
%      pressure=lithostacticpressure(rho_ice,g,surface,z)

pressure=rho_ice*g*(s-z);
