Global Fitting of the Response Surface via Estimating Multiple Contours of a Simulator