Based on the characteristics of the non-linear factor p in grid generation equation, a new procedure of control function which is used to control grid orthogonality is deduced theoretically and mathematically The numerical processes of present method in applications are discussed in detail. Numerical results show that the new procedure has strong influence on the grid orthogonality control.