fix bug on initial surface caculation and add Tneighbor. A much more stable simulation than before, the avgTemp is almost the same between each iteration
merge conflict resolve