 To obtain the coordinates of CV vertices at free surface, we have more unknowns and less equations due to single volumetric flow rate for each cell.Hence the CVs are defined by the cell face centers rather than vertices and vertices are obtained by interpolation.This gives a tridiagonal system for 2D and can be solved using TDMA method.