Fixed OP to now return float values for 'occupancy' for BOTH fluids
Added final touches to Tortuosity
- If no fluid is sent, it get dry tortuosity between every pair of pores
- If a fluid IS sent, it will apply the 'throat.occupancy' to the connectivity matrix, to get the tortuosity in that phase.
- This method returns as a 2D Nt x Nt matrix that can be nicely visualized with plt.imshow