How is the free surface calculated during seepage in OpenGeoSys? In which file can the relevant code be found?