为了体现次网格尺度能量升尺度转换过程中存在的不确定性,文中将随机动能补偿(Stochastic Kinetic Energy Backscatter,SKEB)方案应用于GRAPES(Global/Regional Assimilation and Prediction System)全球集合预报系统(GRAPES-GEPS),以...为了体现次网格尺度能量升尺度转换过程中存在的不确定性,文中将随机动能补偿(Stochastic Kinetic Energy Backscatter,SKEB)方案应用于GRAPES(Global/Regional Assimilation and Prediction System)全球集合预报系统(GRAPES-GEPS),以更好地表征模式误差并且增大集合离散度。使用的SKEB方案基于具有一定时、空相关特征的随机型以及由数值扩散导致的局地动能耗散率来构造随机流函数强迫。并根据流函数与水平风速旋转分量的关系,将SKEB方案中的流函数强迫转化为适用于GRAPES全球模式的水平风速扰动。结果表明,SKEB方案的使用一方面能够提高GRAPES对大气动能谱的模拟能力;另一方面能够改善GRAPES-GEPS的集合离散度与集合平均误差的关系,增加了集合离散度,并在一定程度上减小了集合平均误差,尤其是在热带地区这种改进更为显著。而且该方案使得热带地区连续分级概率评分(CRPS评分)显著减小。就降水预报而言,从Brier评分与相对作用特征面积(AROC,Area under the Relative Operating Characteristics)的结果来看,SKEB方案有助于改善中国地区小雨[0.1 mm,10 mm)、中雨[10 mm,25 mm)与大雨[25 mm,50 mm)量级降水的概率预报技巧,而对暴雨[50 mm,∞)量级降水预报技巧影响很小(24 h降水量)。总体上,模式扰动随机动能补偿方案提高了GRAPES-GEPS的概率预报技巧。展开更多
GRAPES(Global/Regional Assimilation and PrEdiction System)模式动力框架中垂直方向变量的跳层设置采用Charney-Phillips分布,在整层上进行位温、水物质的计算,物理过程中在半层上对其进行处理。这样在GRAPES模式中,进入物理过程之...GRAPES(Global/Regional Assimilation and PrEdiction System)模式动力框架中垂直方向变量的跳层设置采用Charney-Phillips分布,在整层上进行位温、水物质的计算,物理过程中在半层上对其进行处理。这样在GRAPES模式中,进入物理过程之前和物理过程计算完毕之后,都要采用线性插值进行整层和半层之间物理量的转换。由于线性插值精度欠佳,为提高上述反馈过程的精度,并保证水物质的正定性。该研究引入样条插值,并在水物质的插值过程中进行保单调处理,有效减小了位温场、水物质场的预报偏差,并提升了模式的综合预报性能。展开更多
随着高性能计算机技术的发展和应用,并行计算已成为保证数值天气预报模式业务运行时效的关键技术之一。目前高性能计算机计算能力已达到每秒千万亿次浮点计算,系统中处理器数量也早已达十万甚至更多,如此巨大的计算资源对应用软件系统...随着高性能计算机技术的发展和应用,并行计算已成为保证数值天气预报模式业务运行时效的关键技术之一。目前高性能计算机计算能力已达到每秒千万亿次浮点计算,系统中处理器数量也早已达十万甚至更多,如此巨大的计算资源对应用软件系统的设计也提出了挑战。数值天气预报软件系统要充分利用高性能计算机提供的计算资源,必须依靠并行计算方法,这包括适合计算问题的可扩展并行算法的设计、合适的数据分配方案以及良好的任务负载平衡方案。作为中国新一代数值天气预报格点模式,GRAPES(Global and Regional Assimilation and PrEdiction System)设计的最终目标是一个科研/业务通用,区域/全球通用模式。作为一个格点模式,GRAPES的并行计算具有与欧洲中期数值预报研究中心谱模式并行计算不同的特点,GRAPES的并行计算采用了经典的水平网格数据划分。但对于全球的GRAPES模式,由于采用拉格朗日差分方案,模式极地及附近区域格点与格点之间距离的减小,使得模式并行计算在采用简单的经纬网格划分方式实现时,必须考虑极地区域并行计算跨越多个处理器时导致的频繁通讯解决途径。本研究提出了利用消息传递组通讯实现全球格点模式并行计算的一种方法,其核心思想是将极点附近一定区域内的处理器按纬向划归不同的处理器组。文中还给出了该实现方法的任务分配算法,提出了改进的任务分配负载平衡方案。在中国气象局高性能计算机IBM-cluster1600上的测试表明,算法具有较好的可扩展性,其负载平衡方案改善了计算的绝对墙钟时间,使并行计算效率提高10%以上。模式的准业务运行结果表明计算墙钟时间基本可以满足数值预报业务的实时性要求。展开更多
文摘GRAPES(Global/Regional Assimilation and PrEdiction System)模式动力框架中垂直方向变量的跳层设置采用Charney-Phillips分布,在整层上进行位温、水物质的计算,物理过程中在半层上对其进行处理。这样在GRAPES模式中,进入物理过程之前和物理过程计算完毕之后,都要采用线性插值进行整层和半层之间物理量的转换。由于线性插值精度欠佳,为提高上述反馈过程的精度,并保证水物质的正定性。该研究引入样条插值,并在水物质的插值过程中进行保单调处理,有效减小了位温场、水物质场的预报偏差,并提升了模式的综合预报性能。
文摘随着高性能计算机技术的发展和应用,并行计算已成为保证数值天气预报模式业务运行时效的关键技术之一。目前高性能计算机计算能力已达到每秒千万亿次浮点计算,系统中处理器数量也早已达十万甚至更多,如此巨大的计算资源对应用软件系统的设计也提出了挑战。数值天气预报软件系统要充分利用高性能计算机提供的计算资源,必须依靠并行计算方法,这包括适合计算问题的可扩展并行算法的设计、合适的数据分配方案以及良好的任务负载平衡方案。作为中国新一代数值天气预报格点模式,GRAPES(Global and Regional Assimilation and PrEdiction System)设计的最终目标是一个科研/业务通用,区域/全球通用模式。作为一个格点模式,GRAPES的并行计算具有与欧洲中期数值预报研究中心谱模式并行计算不同的特点,GRAPES的并行计算采用了经典的水平网格数据划分。但对于全球的GRAPES模式,由于采用拉格朗日差分方案,模式极地及附近区域格点与格点之间距离的减小,使得模式并行计算在采用简单的经纬网格划分方式实现时,必须考虑极地区域并行计算跨越多个处理器时导致的频繁通讯解决途径。本研究提出了利用消息传递组通讯实现全球格点模式并行计算的一种方法,其核心思想是将极点附近一定区域内的处理器按纬向划归不同的处理器组。文中还给出了该实现方法的任务分配算法,提出了改进的任务分配负载平衡方案。在中国气象局高性能计算机IBM-cluster1600上的测试表明,算法具有较好的可扩展性,其负载平衡方案改善了计算的绝对墙钟时间,使并行计算效率提高10%以上。模式的准业务运行结果表明计算墙钟时间基本可以满足数值预报业务的实时性要求。