We present the results of a 3D global magnetohydrodynamic simulation of the magnetosphere of Saturn for the period of Cassini's initial approach and entry into the magnetosphere. We compare calculated bow shock and magnetopause locations with the Cassini measurements. In order to match the measured locations we use a substantial mass source due to the icy satellites (~1 ¿ 1028 s-1 of water product ions). We find that the location of bow shock and magnetopause crossings are consistent with previous spacecraft measurements, although Cassini encountered the surfaces further from Saturn than the previously determined average location. In addition, we find that the shape of the model bow shock and magnetopause have smaller flaring angles than previous models and are asymmetric dawn-to-dusk. Finally, we find that tilt of Saturn's dipole and rotation axes results in asymmetries in the bow shock and magnetopause and in the magnetotail being hinged near Titan's orbit (~20 RS).