Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged de...Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged density stratification of the whole computation region. But if there is great difference between the area-averaged and the local averaged density stratification, the error will be obvious. An example is given to show that the error from this method may be larger than that from no correction sometimes. The definition of local area is put forward. Then, four improved BPG difference schemes of subtracting the local averaged density stratification are designed to reduce the error. Two of them are for diagnostic calculation (density field is fixed), and the others are for prognostic calculation (density field is not fixed). The results show that the errors from these schemes all significantly decrease.展开更多
基金The Major State Basic Research Program of China under contract No. 2002412403the National Natural Science Foundation of China un-der contract No. 40306014.
文摘Much has been written of the error in computing the baroclinic pressure gradient (BPG) with sigma coordinates in ocean or atmospheric numerical models. The usual way to reduce the error is to subtract area-averaged density stratification of the whole computation region. But if there is great difference between the area-averaged and the local averaged density stratification, the error will be obvious. An example is given to show that the error from this method may be larger than that from no correction sometimes. The definition of local area is put forward. Then, four improved BPG difference schemes of subtracting the local averaged density stratification are designed to reduce the error. Two of them are for diagnostic calculation (density field is fixed), and the others are for prognostic calculation (density field is not fixed). The results show that the errors from these schemes all significantly decrease.
