That’s not quite my intention, there is no specific global saturation for which I could create a liquid level (in which case your methods would work well).
Rather, for each cell I need to figure out where the liquid level would be, given an individual saturation (volume fraction) for that cell. Since my grids are unstructured and can contain polyhedron cells this needs to be done on a cell-by-cell analysis.