Various numerical methods for the two dimensional shallow water equations have been applied to the problems of flood routing, tidal circulation, storm surges, and atmospheric circulation. These methods are often based on the Alternating Direction Implicity(ADI) method. However, the ADI method results in inaccuracies for large time steps when dealing with a complex geometry or bathymetry. Since this method reduces the performance considerably, a fully implicit method developed by Wilders et al. (1998) is used to improve the accuracy for a large time step. Finite Difference Methods are defined on a rectangular grid. Two drawbacks of this type of grid are that grid refinement is not possibile locally and that the physical boundary is sometimes poorly represented by the numerical model boundary. Because of the second deficiency several purely numerical boundary effects can be involved. A boundary fitted curvilinear coordinate transformation is used to reduce these difficulties. It the curvilinear coordinate transformation is used to reduce these difficulties. If the coordinate transformation is orthogonal then the transformed shallow water equations are similar to the original equations. Therefore, an orthogonal coorinate transformation is used for defining coordinate system. A multigrid (MG) method is widely used to accelerate the convergence in the numerical methods. In this study, a technique using a MG method is proposed to reduce the computing time and to improve the accuracy for the orthogonal to reduce the computing time and to improve the accuracy for the orthogonal grid generation and the solutions of the shallow water equations.