use phase pressure to compute FVF.

That is strange, since in function ``fluidDensity()``, the phase pressure are passed to compute ``fluidReciprocFVF()``, but in function ``computeAccum()``, reference pressure still be used to compute FVF.
After this change, a simple case is tested and the result improves a little.