Many changes done until here. Code is running for annular, but not correct yet. Problems with motion of points neighbour to annular interface. Advances in computation of curvature for periodic points, but requires validation.