Code to fit surface ice properties to Viking Lander data