一个高效并行全球海洋环流模式的设计与实现研究.pdf

一个高效并行全球海洋环流模式的设计与实现 宋君强 张理论 宇如聪 俞永强 (国防科大计算机学院高性能计算机研究所，长沙41Q173)(中科院大气物理研究所，北京110000) 一、引言 全球海洋环流模式是气候模拟的一个核心要素，它对于分析和预测气象状况、特别是 海洋气象状况具有特别重要的意义 海洋环流的数值计算往往被归结为一组时变偏微分 的数值求解.该类问题具有很强的物理背景，例如受流体的可压缩性、温度的垂直分布、 地球重力、地球 自转等因素的影响，海洋中存在多种性质不同的波动，包括重力内波、重力 外波、Ro}by波等，它们的时空尺度和传播速度等性质均有差别。另外，海洋环流模式和 大气模式、海冰模式间的祸合，特别是与大气模式间的藕合需要交换大量的数据资料。这 些特点决定海洋环流数值模拟具有海量存储、巨量计算的特点，对其采用并行计算是必T 的;我国在并行海洋模式的研究和实现方面尚是一项空白。本文实现了一个高效并行全 球海洋环流模式，其原始串行模式来源于中科院大气物理研究所第三代海洋环流模式 文中给出了并行计算的分析与实现方法。包括若干优化策略。优化后的并行海洋模式具 有较高的计算效率。 二、串行数值模式 中科院大气物理研究所的全球海洋环流模式要求解一组时变偏微分方程。包括动量 方程、连续方程、热力学方程、盐度方程和海水状态方程等。在该数值模式中包含了太阳 辐射、长波辐射等一系列参数化过程和基于热动力学模式的海冰计算(} icemodel)海 冰间的热交换在计算中并不是动态的。在对海洋原始方程的数值求解时，针对重力外波 传播速度快、运动尺度小，重力内波传播速度慢、运动尺度较大的不同特点，模式的动量方 程积分可分为两个相对独立的部分:二维外部模态即正压模态(externalorbaroturpic)和三 维内部模态即斜压模态(internalorbaloclinc)o 模式计算采用B一网格系统(B一grid)，与水平速度计算有关的格点称为V点，与温 度、盐度等模式变量有关的格点称为T点，虽然二者的经纬向网格下方法标所对应的地理 位置是不同的，但各自构成的网格子系统都可视为是规则的。在垂直方向上，垂直速度与 压力在相邻两层的交界面上计算，其他对变量的计算均视为在各层的中间层上进行。正 压模态、斜压模态以及热动力过程的数值积分采用了不同时间步长的蛙跳格式，为防止采 用蛙跳格式引起的数值解在奇偶时间步发散，时间积分在每月的第一步采用Euler前差格 式。正斜压模态的每个时间步都采用时间滤波以保证计算稳定。每天的积分过程大致如 下:首先是正斜压模态的祸合计算，正压模态积分NB步，时间步长为4ta,接着斜压模态 积分一步，其时间步长为4tc=NB乙tB;在正斜压模态计算N。步后，热盐过程积分一步， 其步长为4tts=Nc-,9tc.海洋模式的串行计算流程图1. 179 预处理过程 ‘ 告 亚 . 力 内 牛莫11初始化 波 积 分 I;llj史文们 对{】期循月 1 结束