全球区域同化预报系统的使用及其数值模拟

《全球区域同化预报系统的使用及其数值模拟》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、全球区域同化预报系统的使用及其数值模拟*陈建萍（江西省环境预报中心，江西，南昌，330046）　　摘　要：介绍了中国气象科学研究院数值预报中心自主开发的GRAPESV2.1数值模式运行步骤，并利用该模式对2005年5月1日06～07时在江西南昌市区出现的30mm/h强降水过程,进行了数值模拟和诊断分析。结果发现，北方干冷空气南下，到达30°Ｎ附近后，与南方暖湿空气在南昌市区上空低层形成交汇，促使暖湿空气沿着锋面迅速抬升凝结，触发不稳定能量的释放，是造成南昌强降水产生的主要原因；GRAPES模式能很好地模拟出高低空环流形势场特征和主要天气系统，能提

2、供较准确的高分辨率诊断分析资料。　　关键词：GRAPES模式　强降水　数值模拟　诊断分析  中图分类号：P456.7；P458.1+21　文献标识码：A　文章编号：1007-9033（2005）02-00023-041　GRAPESV2.1简介及模式运行1.1　模式简介　　GRAPES（Global/RegionalAssimilationandPredictionSystem）—全球区域同化预报系统，是中国气象局负责研究开发的、研究与业务通用的数值天气预报系统。GRAPES系统是集常规与非常规变分同化、静力平衡与非静力平衡、全球与区域模式、科研

3、与业务应用、串行与并行计算、标准化与模块化程序、理想实验与实际预报等为一体，中小尺度与大尺度通用的先进数值预报系统。1.2　GRAPESV2.1模式的运行　　使用2005年4月29日20时的T213资料，以及2005年4月30日08时的探空资料，对2005年4月30日08时—5月1日08时的暴雨过程进行数值模拟。1.2.1　资料准备　　2005年4月29日20时的Ｔ２13资料（BJ*.*）放在grapes/data/9210_t213目录下，2005年4月30日08时的探空资料（*.ABJ）则放在grapes/data/9210_obs目录下。1

4、.2.2　探空资料预处理　　进入grapes/grapes_3dv/preproc/9210_obs_proc/datatran，使用命令pgf90-optrans0.exeptrans0.for和pgf90-optranstp.exeptranstp.for分别编译当前目录下的ptrans0.for和ptranstp.for；成功编译后，将在当前目录下分别生成ptrans0.exe和ptranstp.exe文件。然后使用命令catgrapes/data/9210_obs/*.ABJ>fn.dat把ABJ文件连接成１个文件，执行后当前路径下会生成

5、fn.dat。再在当前目录下把wsdat.dat的时间修改为05042300。最后分别执行ptrans0.exe和ptranstp.exe，将在grapes/grapes_3dv/data/input/GTS目录下生成TEMP文件。1.2.3　T213资料预处理　　首先进入grapes/9210_t213_proc目录，使用命令：9210_t213_proc.sh200504291248grapes/data/9210_t213/grapes/T213_DATA/，对T213原始报文资料进行解码。执行后将在grapes/T213_DATA目录下自

6、动生成2005042912目录，并在该目录下生成T2132005042912、T2132005042918……T2132005050112等9个文件。　　然后进入grapes/grapes_3dv/preproc/9210_t213_proc目录，修改namelist.t213中解码后T213资料的路径，即把record4中Ｔ213_file的路径修改为grapes/T213_DATA/2005042212/T2132005042300。　　最后运行make编译当前目录下的文件，成功编译后将生成t213proc.exe文件；执行t213proc.

7、exe后，将在grapes/grapes_3dv/data/input/backgrd目录下生成以xb.dat命名的背景场。1.2.4　分析系统编译连接及运行　　进入grapes/grapes_3dv/3dv/rundir，运行make_script脚本文件，编译整套分析系统，成功编译后有可执行文件Grapes3dvar.exe生成。运行时首先要编辑当前目录下的namelist.3dvar,修改record3中的ANALYSIS_DATE='2005043000'，以及record14中的file_dxa_binary='../../data/o

8、utput/dxa2005043000.dat'和file_xa_binary='../../data/output/xa2005043