大气数值模式设计中的一些问题.pdfVIP

大气数值模式设计中的一些问题 序洞贤 (国家气象中,u，北京100081) 下面讨论3个问题，即模式的构造、模式各部份之间的协调和不连续问题。 一、模式的构造 通常，用“尺度分析”等方法构造好模式并进行离散后，虽其所引进的误差早已为人能 熟知，但它对方程中其他项的作用却很少有人讨论。分析表明:当方程中各项量级相差 悬殊时，如分辩率取法不当，某项的截断误差可能大于或相当于某些重要的项，从而，影响 预报效果。在取近似时也有类似情况。另一方面，尺度分析往往只用来分析方程的动力 部份，对非绝热物理过程，还没有见到像尺度分析那样的方法。有人建议用统计方法对其 相对重要性进行分析，但具体细节还需要研究。下面我们将对误差的影响和物理过程相 对重要性分析分别进行讨论。 1.误差的影响 (1)取近似的影响 以取重力加速度为常数为例。按公式[[I] g=g,(1+0.0053026W，一0.0000079.229)(1一a,z) (1) e.04 1一海平面 0.03 2-5000. 0.02 000Gm 0.01 -0.01 刃.02 一 03 -0.04 一 05 一 06 图1各等高面上，E随甲的交化 如定义E二g-9.8(以幽一2为单位)根据计算，可以得到图1,这里go=9.7031ma-2 a,=3.14x10-7.-1从图中可以看出:在赤道上空的10000.高度,E可以达到 15 · 5x10-2M-2，一般也有lo-I.-I,这和在中一Y尺度运动中aw/at的量级相当fzl也和一 些强对流中asv/at的量级相近。显然.不考虑9随甲和:的变化，要描写好非静力作用是 困难的。 (2)截断误差的影响 用中心差商代替微商，垂直运动方程可以写成 SPH_工m (2) aatr，其他项和截断误差=一告 乱 p8z一8一告、一告Ed rz, 0 《z簇介， (3) 兀一rzr, rZ ‘z . 并定义 已殊二EH/(aPH/az)， 则 一E1a二 1apH relEb (4) N p az 至于一玛巾，如设P，在:=o,zT处为零，则P可‘表示为 P=三Pmsintmtz 定义 7e1(Ed)=(Ed).,/(即/az)m, P1Ea，一E“d)二，·告az (5) 其中()m是()中第m个成份;玩和F,d各是p的静压力的垂直差商和气压偏差的垂 直差商的截断误差。于是，从静力方程的解，如一(aPH/az)净为lo.-I,则可以得到在各 不同高度曰与和}EH/p1随AZ的变化(表1).;还可以得出，在不同I(ap/az)/pI的情况 下，在M二to时，