沙瓦特指数的数值计算方案比较研究

比较了三种假相当位温(θse)计算式的精确度，分析了θse对沙瓦特指数(si)计算的重要影响，提出精度较高的si计算实用方案。该方案首先利用干绝热过程中位温和比湿守恒性质迭代求解抬升凝结高度处温压变量；计算起始气块的θse；然后基于湿绝热过程假相当位温的守恒性质，采用stackpole迭代方法计算湿绝热上升到500hpa的气块温度，进而计算si。研究发现，θse的误差是si计算误差的重要来源，采用准确的θse计算方案非常重要。经计算比较，发现θse的bolton公式和李任承公式计算精度较高，而经典的rossby公式计算得到的θse表现为系统性偏低，偏低量随气块温湿的增加而增大。相应地，采用bolton公式计算θse所得到的si之精度最高，与查表值703个算例之间的平均绝对误差仅0.45℃，李任承公式效果与之接近；而传统的rossby公式相应的平均绝对误差为0.66℃。与前人的si算例对比表明，本方案计算精度得到进一步提高，在预报业务中具有一定实用价值。