p_nw_top Hp_nw_t(t)*rho_water*g_const*sigma p_nw_init p_nw_top+(8.34[cm]-y)*rho_nw*g_const p_w_init -rho_w*g_const*y p_w0 0.1[m]*sigma*rho_w*g_const*(t/170[h])