Added output at edges in profiles_edges.c
Scalar is now CWC
Fixed surface discontinuity via NewCells before ComputeVelocityVector
Added interpolation at top boundary
Working with strong tides and wind.
Working SFBay 3D and all other test cases
Works when 2d and 1d cases are the same
Working Newton iter. with huge dt