Python script used to calculate phase saturation within a sample during injection of alternative phase (e.g. a dual phase system). Calculator relies on end phase saturation conditions of both fluids and intermediate saturation during injections determined from CT scanning.