A grid combination method, based on a modified hyperbolic grid solver, is developed to generate tw+ dimensional grid systems in regions with prescribed grid distribution along all of the boundaries. The hyperbolic grid solver uses a predictor-corrector scheme and is a modified version of the hyperbolic grid generation method of Tai, Yih, and Soong, which applied the Roe first order upwind splitting to the Steger and Chaussee grid generation equations. Four grid sys-0 intersection angle between grid lines terns are first marched from each boundary toward its opposite boundary and then combined to give the a,o functions defined in As arc length, Eq.(12) A t , . . . = t i -ti-1,. . . user specified constant, Eq.(9) eigenvalue of C scaling factor, Eq.(ll) coordinates on computational domain t , 7