In this paper, Laplace-Beltrami equations are used to generate orthogonal curvilinear grids on the sphere for ocean models. In addition to overcoming the pole-problem, the grid configuration has quasi-uniform cell-size on the whole sphere. Some quantities such as the grid length along two directions, the angle deviation from orthogonality, the area of the cell to evaluate the quality of the grid, which demonstrate the grid produced is fit to be a model grid on which the finite difference method or finite volume method can be implemented for numerical simulating of global atmosphere and ocean dynamics on large scale.