Include option for liquid water in mean profiles in hydrostatic pressure calculation (2nd order only for now)
Don't use simplified exner() function
Unified version of calchydropress for 2nd and 4th order scheme
Change in variable name from dn to rho
Allow for saturation at surface in pressure calculation