computer program for the calculation of steady three dimensional adiabatic flows in supersonic rocket engine nozzles. The fluid is assumed to be an inviscid perfect gas and the flow is assumed to be shock free. The governing equations for this flow model consist of the continuity equation, the momentum equations, and the speed of sound equation. The steady three dimensional supersonic flows is computed using a pentahedral bicharacteristic method. An inverse marching method based on the network of wavelines and the pathline is employed. The spline surface fit is used in the iteration procedure to evaluate the location of the unknown point in the new reference plane and its flow properties-computer program for the calculation of steady three dimensional adiabatic flows in supersonic rocket engine nozzles. The fluid is　 assumed to be an inviscid perfect gas and the flow is assumed to be shock free. The governing equations for this flow model consist of the continuity equation, the momentum equations, and the speed of sound equation. The steady three dimensional supersonic flows is computed using a pentahedral bicharacteristic method. An inverse marching method based on the network of wavelines and the pathline is employed. The spline surface fit is used in the iteration procedure to evaluate the location of the unknown point in the new reference plane and its flow properties