天河二号试算快报——全球正压大气浅水波模式-上海超级计算中心.PDF

高性能计算应用 67 天河二号试算快报 ——全球正压大气浅水波模式  杨超 中国科学院软件研究所 北京100190  薛巍 付吴桓 尤洋 王欣亮 徐阳彤 甘霖 清华大学 北京100084 前言 关注。大气模式作为上述研究的核心量化工具，其 特殊性和重要性决定了它是目前地球各圈层模拟系 2013年6月13日至14日，在天河二号正式发布前 夕，受国防科学技术大学邀请，中国科学院软件研 统中最复杂、最关键也是发展最完善的模块。观测 究所和清华大学联合对天河二号进行了一次应用试 和数学物理模型的不断完善，数值模拟的不确定性 算。 使得发展超高分辨率大气模式成为当前该领域最重 要的研究方向。 本次试算所使用的测试程序是Scala-Core- gSWE-E1.3模式，其通用处理器版本由中国科学院软 然而，分辨率的提高毫无疑问地导致计算量大 件研究所研制；混合架构优化版本由清华大学和中 幅攀升。以天气/气候模式为例，如果水平分辨率 国科学院软件研究所联合完成。该模式用于求解全 提高10倍，计算量将增加数百乃至上千倍。这使得 球正压大气浅水波方程，其主要特点如下： 天气、气候模式成为当前高性能计算领域公认最重 • 基于立方球网格，避免极区所带来的计算精 要的大型应用之一，也是未来E级计算系统的一大 度和负载不均等问题； 驱动，因此得到世界上主要超级计算机中心、计算 • 二阶有限体积离散格式，严格保持质量守 机提供商和研究机构的高度重视。针对这样的计算 恒，并能够有效处理任意复杂地形； 挑战，大气模式可扩展并行算法和性能优化的研究 得到了广泛的重视。比如美国国家大气研究中心支 • 支持主流的显式Runge-Kutta时间积分方 法； 持了动力框架研究计划HOMME，其发展的谱元方 • 构建于PETSc数值计算库之上，充分利用其 法已经在CAM中得到应用，并行可扩展性迈向10万 高级数据结构和MPI通信库； 处理器量级，而针对CPU-GPU混合并行架构的新 • 采用新型水管流方法隐藏立方球网格上的邻 一代谱元动力框架也在积极研制当中。欧洲中尺度 居通信； 天气预报中心长期不懈地改进其发展谱方法，已经 • 支持CPU、CPU-GPU、CPU-MIC等多种架构 可以支持数千甚至上万处理器的有效并行，正在积 的大型并行计算系统。 极研制基于Intel MIC的新一代动力框架。在日本， 与Scala-Core-gSWE-E1.3模式有关的参考文献包 K Computer的拥有者日本理化研究所成立专门的研 括[1]-[6]。特别指出，基于该模式，中国科学院软件 究小组支持超高分辨率大气模式的研制，以及在K 研究所和清华大学的团队曾联合在天河一号上开展 computer上的开发和优化。当前，高分辨率大气模式 了CPU-GPU异构算法的研究工作，该工作已发表于正在从万核有效模拟向百万核模拟攀升。 PPoPP’13上（参见