Merge pull request #70 from andlaus/improve_doc

make the phase-presence docstring understandable for mere mortals
include capillary pressure in the PDEs

I'm neither sure that this is fully correct nor that I found all
occurences (so far, the output writing code is missing in this patch),
but it seems to work for SPE1...
use the reference pressure to calculate the transmiscibility factors