Should now be in sync with cfs_tpfa_residual C interface. Simple well gravity model implemented. More flexibility in well gravity models would be a natural future extension.