


LITHOSTATICPRESSURE - compute the lithostatic pressure
Usage:
pressure=lithostacticpressure(rho_ice,g,surface,z)

0001 function pressure=lithostacticpressure(rho_ice,g,s,z); 0002 %LITHOSTATICPRESSURE - compute the lithostatic pressure 0003 % 0004 % Usage: 0005 % pressure=lithostacticpressure(rho_ice,g,surface,z) 0006 0007 pressure=rho_ice*g*(s-z);