add neumann test, add pressure modulation in darcy solver