This paper is concerned with the mathematical model for the numerical simulation of supercritical surface flows. A boundary-fitted coordinate system was used to overcome the difficulties and inaccuracy associated with the determination of flow characteristics near the flow boundaries. The MacCormack scheme was applied for the solution of the transformed system of equations. Comparisions between computed results and experimental data show a satisfactory agreement.