change the code for calculating front tracking with surface tension and 2nd order terms