f := gp.sat(gp)
gp.sat(gp) := f

Get/set the fluid saturation. This should be a value between 0.0 and 1.0. At 1.0, the gridpoint will be considered fully saturated and will begin to create pore pressure.

Returns:f - the current gridpoint saturation
Accepts:f - the new gridpoint saturation
Arguments:gp - a zone gridpoint pointer