The 2 dimond shapes on the contours in all the figures above are the locations, (14.5oS,104.5oW) and (42.5oS,166.5oW), did not converage to the solutions in the least-square procedure. Very high stopping tolerance, like 1E-10, is needed to obtain the solutions.