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