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.
Thanks for noticing. I'll submit a fix.