沙瓦特指数一种迭代算法.PDFVIP

第 卷 期 应 用 气 象 学 报 ， ２０ ４ Ｖｏｌ．２０ Ｎｏ．４ 年 月　　 　 　 　 　 　 　 　　　　 　 ２００９ ８ ＪＯＵＲＮＡＬＯＦＡＰＰＬＩＥＤＭＥＴＥＯＲＯＬＯＧＩＣＡＬＳＣＩＥＮＣＥ Ａｕｕｓｔ２００９ 　 ｇ 沙瓦特指数的一种迭代算法 １） １） １） １） ２） 王学忠 胡邦辉 吕 梅 邹 力 倪东鸿 　 　 　 　 　 　 １）（解放军理工大学气象学院，南京 ２１１１０１） 　 ２）（南京信息工程大学江苏省气象灾害重点实验室，南京 ２１００４４） 　 摘 要 　　 为了提高沙瓦特指数的计算精度，在沙瓦特指数计算方案中引入二分法：先用于计算抬升凝结高度上的要素， 进而计算假相当位温，再依据假相当位温守恒性质用于计算 ５００ｈＰａ气块温度。与其他方案对比表明，该迭代算法 计算的沙瓦特指数精度较高；与查表法所得的气块温度对比表明，该迭代算法得到的气块温度的最大绝对误差为 ，平均误差（ ）与平均绝对误差（ ）数值接近；迭代算法与查表法之间以系统性偏差为主；此 １．３６℃ －０．６８℃ ０．６９℃ 外，该迭代算法计算量小，收敛速度快，具有推广应用价值。 关键词：沙瓦特指数；二分法；查表法；误差分析 其他诊断量的多级复合函数，在其计算过程中需要 引 言 考虑误差传递、放大效应，这就要求中间诊断量的计 　 算应达到较高精度。因此，提高计算精度的需求具 沙瓦特指数是 世纪 年代引入的一个稳定 有普遍性。在沙瓦特指数的计算方案中引入二分法 ２０ ５０ ［］ １ 度指数 ，定义为 ５００ｈＰａ等压面上环境温度与 可作为改进其他诊断量算法的参考。 ８５０ｈＰａ等压面上的湿空气块干绝热抬升，到达抬 升凝结高度后湿绝热上升至 ５００ｈＰａ时气块温度的 １ 二分法简介 　 差值。它反映了大气稳定度状况，当沙瓦特指数小 于 时，表示层结不稳定，且负值越大，不稳定程度 ０ 二分法的数学问题可描述为：当函数 （） ＝ 狓 狔 犳 越大，反之则是稳定的。该指数的提出在雷暴指数 在区间［，］上是单调函数，已知 且 ［ 犪犫 狔 狔