OP now takes an optional list of invasion pressures (inv_points=[list of presssures])
- If inv_points is given, the npts argument is ignored
- if inv_points is NOT given, then a list of N points is generated spanning the size range in the 'capillary_pressure' array.