基于新型LE水平跳点网格的阴阳网格平流方案设计

为了模拟球面平流传输过程,本文基于球面阴阳重叠网格设计了一种两时间层半拉格朗日平流方案.该方案在球面坐标下采用新型的LE水平跳点网格,同时针对阴阳网格重叠区,采用了不同插值方法进行比较分析,且进行了相关的理想数值试验对方案设计效果进行评估.数值试验表明方案设计是成功的,阴阳网格重叠区平流对插值方案比较敏感;半拉格朗日方案能较好地模拟球面刚体平流和变形涡旋的结构、位置及演变过程,并具有较好的数值稳定性和较高的数值精度.     

基于国家气象观测站逐日降水格点数据的交叉检验误差分析

利用2006-2008年2 403个国家气象观测站地面雨量计的逐日降水量资料,采用与网格点最近的观测站有、无降水确定该网格点有、无降水和Barnes插值方案确定网格点降水大小的混合插值方案,得到全国空间分辨率为0.1°×0.1°（约10km×10km）的逐日降水量格点数据,在此基础上通过交叉检验方法统计格点数据的误差,...     

守恒型有理函数插值半拉格朗日平流方案及性能分析

通过多种理想试验对正定、保形守恒型有理函数插值半拉格朗日平流方案分别在平面直角坐标以及阴阳网格球面坐标中进行了计算性能分析,并采用多种误差模对守恒型有理函数插值半拉格朗日平流方案的网格收敛性进行评估。结果表明,采用分段有理函数插值的守恒型半拉格朗日平流方案可以有效消除不连续分布处的数值振荡、保证正定性,物理场平滑分布时维持1-2阶收敛速度;而在不连续点或大梯度区域以及应用分维技术的多维算法都会通过有理函数的降阶特性,影响平流计算的收敛阶数,并且,在球面坐标中受球面曲率的影响,守恒型有理函数插值半拉格朗日平流算法的网格收敛速度有所降低。     

广西冬季最低气温的小网格分析

研究各重现期广西冬季最低气温的小网格气候分析方法。采用网格距为 0 .1× 0 .1经纬度的小网格 ,从地理信息系统 GIS中读取海拔高度资料 ,用回归分析方法求气候方程 ,利用 Gr ADS绘图系统提供的 oacres内部插值函数求取格点余差插值 ,利用改进的 Gr ADS绘图系统绘制小网格气候分析图     

雷达径向速度资料同化中不同坐标转换方案的对比试验

利用江苏省气象局与美国强风暴实验室联合开发的高精度数值分析及预报系统(Precision Weather Analysis and Forecast System,PWAFS)对雷达资料同化中径向速度资料的两种坐标转换方案进行对比分析。Grid方案将雷达径向速度资料通过最小二乘法从极坐标映射到模式三维网格;Tilt方案将雷达径向速度资料通过双线性插值在水平方向插值至标量水平网格,但在垂直方向不进行插值,保留在雷达仰角对应的高度上。两种方案对反射率资料的处理均是插值到模式三维网格点。Grid方案在近雷达处进行平滑,在远雷达处进行插值,会导致低层数据平滑,Tilt方案减少了雷达径向风观测垂直插值引发的误差,更多的保留了雷达观测的特性。本研究分别通过龙卷、大风及梅雨锋暴雨个例对这两种方案的同化结果进行对比分析。龙卷个例中Grid方案得到了部分虚假的较大的同化风场,Tilt方案结果清楚展示了龙卷发生位置的回波及流场的精细结构。大风个例中两种方案得到的最大风速值差3 m·s^-1,Tilt方案的结果更接近观测最大风速值,且得到的大风速区分布更符合观测。梅雨锋暴雨个例中Grid方案对东北及西南两个区域的大风速区均未能很好的反映,Tilt方案得到的水平风速大值区范围明显优于Grid方案。在靠近雷达中心的低层,观测资料密集,Tilt方案能够更好的反应实际大气状态。但是因为缺乏其他观测资料进行验证,两种方案的效果还需要利用数值预报或其他方法进行对比。     

基于Surfer Automation接口的气象等值线图的绘制

气象等值线图是由许多离散的站点数据经过网格点插值,再通过逐个网格判断来绘制等值线,加上颜色渲染而生成的图形。通过研究Surfer软件提供的Automation接口,用VB编程实现快速绘制出各种气象等值线图,可自定义绘制并可实时修改图层属性,定制出形式多样的等值线图,实现了Surfer功能的完美封装。     

阴阳网格上质量守恒计算性能分析

质量守恒数值计算是球面准均匀阴阳网格构造全球大气环流模式的重要条件,也是提高阴阳网格应用质量的重要技术手段。本文针对通量形式平流方程,在球面坐标上采用多种理想数值试验对阴阳网格上的三种守恒计算方案和边界插值非守恒计算方案进行了比较检验。发现,质量守恒方案不仅对全球数值积分重要,还影响数值计算精度,满足局地守恒条件的全球强迫守恒方法可以获得较高的精度;网格内质量均匀分布的阴阳网格边界通量一致性守恒强迫计算方案,实现了在不增加计算误差条件下保证局地和全球守恒的目的,且具有很小的计算负担,可以作为阴阳网格上全球质量强迫守恒的有效计算方案;而网格质量的线性分布可以有效提高阴阳网格的数值积分计算精度,但在一定程度上会增加计算负担。     

EC细网格对黑龙江省日极端气温预报的误差分析

本文采用2013、2014(1-10月)年的EC细网格2 m气温预报资料,通过双线性插值方法对黑龙江省83个站点进行插值,得出EC细网格模式对于日极端气温的预报误差并进行分析。结果表明,EC细网格模式本身对于最高气温的预报准确率高于最低气温,预报误差具有显著的季节性变化特征,但不同的起报时间对误差的影响并不明显。EC细网格模式的预报误差符合正态分布,具有可订正性。     

“准拉格朗日法”和“欧拉法”的数学一致性与三次插值函数算法时间积分方案

从预报方程组通式和欧拉算符出发,用泰勒级数展开,推导出时空间微商余项为二阶、四阶完全预报方程组。可以证明,同阶时空间微商余项的准拉格朗日法和向前差分欧拉法时间积分方案,具有相同的物理意义和数学一致性,相比之下,传统中央差欧拉法时间积分方案是"简单格式"与"增大不可预测计算误差"并存。进而讨论用"三次插值函数"实现二阶时空间微商余项准拉格朗日法、或同阶向前差分欧拉法(可"二选一"),它们应该分别替代"双线性插值"准拉格朗日法和传统时空间中央差欧拉法,因前二者时空间微商余项及计算精度高于后二者。所以,三次插值函数算法可将准拉格朗日法和欧拉法时间积分方案、以及CFL判据统一起来。由于三次插值函数具有对原函数"变量场"的数学定律"收敛性"和二阶可导"最优性",且一次"三次插值函数"运算,即具对网格变量场二阶可导拟合"等价性":不仅拟合变量场斜率、还拟合其曲率和挠率。并因周期"三次插值函数",可作全球变量场"三次插值函数"二阶可导拟合,实现全球"三次"数值模式,并且可按变量场曲率判断,作变量场局域或单点平滑,保持"三次"模式时间积分的稳定性。     

可变网格的样条客观分析方法

本文根据我国气象工作的实际需要,采用样条插值方法,建立了一种网格和区域都可以变化的客观分析方法。曾对该方案中有关参数进行了大量的试验,指出了在客观分析中采用样条插值一些值得注意的问题。该方案不但能适应多方面工作的需要,而且能显著地提高工作效率。所做的实例分析,效果也是较好的。