Interpolation code. Plot contours