* add a checktwophasesolution to check also the discontious pressure, the codes is ugly now