摘要
This study shows a new way to implement terrain-following s-coordinate in a numerical model,which does not lead to the well-known"pressure gradient force(PGF)"problem.First,the causes of the PGF problemare analyzedwith existing methods that are categorized into two different types based on the causes.Then,the new method that bypasses the PGF problem all together is proposed.By comparing these threemethods and analyzing the expression of the scalar gradient in a curvilinear coordinate system,this study finds out that only when using the covariant scalar equations of s-coordinate will the PGF computational form have one term in each momentum component equation,thereby avoiding the PGF problem completely.A convenient way of implementing the covariant scalar equations of s-coordinate in a numerical atmospheric model is illustrated,which is to set corresponding parameters in the scalar equations of the Cartesian coordinate.Finally,two idealized experimentsmanifest that the PGF calculated with the new method is more accurate than using the classic one.This method can be used for oceanic models as well,and needs to be tested in both the atmospheric and oceanic models.
基金
supported by the Knowledge Innovation Program of the Chinese Academy of Sciences(KZCX2-YW-Q11-04)
the National Basic Research Program of China(973 Program,Grant No.2011CB309704)
The second author was supported by the National Natural Science Foundation of China(NSFC)under Grant No.40875022,41175064 and 40633016.
