change waterdensity to 1040kg/m^3, in addition, init water pressure in hydrostaticIC should not be added to rogh