Update hydrostatic pressure; calculate surface/domain top values from ghost cells, fill top ghost cell based on pressure domain top