Abstract: New computational schemes of pressure gradient force are developed in this paper based on the coordinate transformation formula with difference-differential-difference consistence. In the mean time the so-called static subtraction method is suggested by using the initial field as the basic one. Several schemes are checked by some forecast experiments with real data. The results show that the new computational scheme is successful, but the initial static subtraction method is not as good as expected. The cause may be the errors in the vertical interpolations of the geostrophic wind.