WRF数值气象预报模式系统在深腾6800上移植与测试.docVIP

WRF数值气象预报模式系统在深腾6800上移植与测试 　　摘 要：介绍了WRF中尺度数值气象预报模式与联想深腾6800高性能计算机系统，概述WRF模式在联想深腾6800上的移植与测试过程。采用3D数据集em_quarter_ss进行性能分析与正确性测试，并绘制气象图。分析串、并行状态的运算加速比与处理器个数间的关系，发现并行效率相当好，非常适用于机群系统。WRF在深腾6800上的成功移植与测试，也为进一步将WRF应用于中国国家网格奠定了研究基础。 　　关键词：天气研究和预报;深腾6800;中尺度;数值模拟;移植;测试;国家网格 　　中图分类号：TP311.5文献标志码：A 　　文章编号：1001-3695(2007)06-0245-03 　　 　　1 WRF模式特点与相关结构?? 　　WRF(Weather Research Forecast) 数值气象预报模式系统是美国气象界在继开发MM5数值模式之后，再一次联合开发的新一代中尺度气象研究与预报数值模式及同化系统［1］。对比2000年10月发布的第一版，2004年6月发布的WRF第二版具有更多特点，包括多重区域，从几公里到数千公里的灵活分辨率，多重嵌套网格，以及与之协调的三维变分同化系统3DVAR等。WRF数值模式采用高度模块化、并行化和分层设计技术，集成了迄今为止中尺度方面的研究成果［2］。?? 　　WRF为中尺度研究和业务数值气象预报提供共同框架，设计用于1－10km格距的模拟，与一般计算平台上典型的实时预报能力相适应；它也为理想化动力研究以及完整物理的数值气象预报与区域气候模拟提供共同框架。?? 　　与MM5数值模式相比，WRF的优点有：多维动态内核，更为合理的模式动力框架，先进的三维变分资料同化系统，可达几公里的水平分辨率及集合参数化物理过程方案等［3］。 　　图1为WRFv2模式与其他输入输出等模块的关系流图，描述了WRF模式的总体运行过程以及与其他模块的相互联系。 　　图1中，Part 1为外部数据集，即实际计算数据集，需经过WRF SI数据标准化预处理后，利用WRF进行实时预报与气象研究；Part 2为WRFv2模式提供的几个测试数据集；Part 3为标准输出与可视化模块，主要对数据进行后处理，并采用气象专用绘图软件GrADs对输出文件进行图形绘制。?? 　　对WRF内部模块按功能［6］进行归类，详见表1。 　　 　　2 深腾6800介绍?? 　　面向网格的联想深腾6800高性能计算机系统为COW结构，它包括265个四路结点机，1 060个主频为1.3GHz的安腾II处理芯片（其中1 024个处理机用于计算），分布式存储，内存总容量为2.6TB，磁盘存储总容量为80TB（其中光纤盘阵容量61TB）。主要节点采用LSF作业管理系统，高速连接网络为QsNet，点对点通信带宽大于us300Mbps，延迟时间小于7us。作为中国科学院e－science重要基础设施环境以及中国网格CNGrid的主要结点资源，深腾6800高性能计算机主要应用于大规模科学工程计算、数值模拟、商务计算和网格信息服务等多个领域。 　　 　　3 移植?? 　　移植过程是一个不断反复、螺旋上升的过程，主要对程序进行编译和调试。需要熟悉MPI、OpenMP并行编程以及c/fortran90混和编程技术，并需对WRF原理及模块间关系有一定了解［7,8］。下面总结移植过程，并辅以一些示例说明。 　　3.1 修改编译器?? 　　根据机器特点选择Intel Itanium2 ia64 Madison Linux DM－Parallel (e.g. mpp2 at PNNL), efc71 compiler (RSL, MPICH, nesting) 编译方案。由于使用深腾6800上的ecc/efc71编译器编译通不过，重新下载并安装和使用了Intel ecc/efc Compiler 7.0编译器，编译成功。?? 　　修改configure.wrf及各模块Makefile文件的CC/SC编译器为Intel ecc Compiler 7.0，FC/LD/SFC编译器为Intel efc Compiler 7.0。 　　3.2 修改编译规则?? 　　修改各Makefile文件及rules.make文件中的Fortran程序编译规则：Intel efc Compiler 编译器编译Fortran源程序为目标文件时，编译器只识别后缀名为f 90的文件，故需将*.F文件先编译为*.f 90，然后再生成*.o目标文件。涉及修改Makefile文件的模块有esmf_time_f90、io_esmf、io_grib1、io_int、io_mce1、io_netCDF、io_ph