check results to see if they look physical
The following changes where made in separate commits to separate the steps more easily:
- fix the 1p_2p test, as the boundary condition was wrong
- change the momentum coupling condition for multiphase flow to also reconstruct the coupling pressure with darcys law
- include the flux of the other phase for multiphase models in the coupling conditions
The results seem more dry with the changes, which seems logical as there should be more evaporation. The fluxes become 0 when we have a phase switch in the 2p2c model, which is also correct.