WRF模式技术手册

1、附件（二）WRF模式技术手册目 录1. WRF模式概况32. WRF的预处理系统（WPS）42.1介绍42.2 各个WPS程序的功能52.2.1程序geogrid52.2.2程序ungrib52.2.3 metgrid程序62.3运行WPS72.4 WPS创建嵌套区域212.5 USGS和基于MODIS的地形数据之间的选择242.6 用于重力波拖曳方案的静态数据252.7 WPS中并行282.8检查WPS的输出292.9 WPS应用程序302.10 WRF区域向导352.11 把气象要素场写成过渡格式352.12 生成并编辑变量名列表（Vtables）382.13 Geogrid二进制格式写静

2、态数据气402.14 namellist变量的讲述432.15 GEOGRID.TBL选项的介绍502.16 指标选项的介绍532.17 METGRID.TBL选项介绍562.18 Geogrid和Mtgrid中可行的插值方案592.19 静态数据中的土地利用和土壤类别612.20 WPS输出场643.WRF-NMM初始化763.1介绍763.2 对实际数据进行初始化例子773.3运行real_nmm.exe784. WRF NMM 模式基本特性794.1引言794.2 WRF-NMM动力学框架804.3 物理参数化方案选项804.4业务配置864.5变量列表说明874.6如何运行WRF100

3、4.7多区域的设计1024.8实际数据测试案例：2005年1月23日00时到24日00时1054.9 WRF-NMM输出场的目录1064.10 WNM-NMM动力和物理过程的参考文献列表1161. WRF模式概况WRF(天气研究和预报)系统的核心NMM（非静力中尺度模式）由NOAA/NCEP发展而成。目前释放的版本为V3。WRF-NMM是灵活的，完美的大气模拟系统，具有易携带，高效，且可并行运算的特性。WRF-NMM广泛应用于从米到数千公里，包括：l 实时数值天气预报；l 预报研究；l 参数化研究；l 耦合模式应用；l 教学。WRF的前处理系统（WPS）用于实时的资料处理，功能包括：l 定义模

4、拟区域；l 插值地形资料（如地形，土表和土壤类型）到模拟区域；l 插值其他模式的资料到模拟区域和模式坐标。WRF-NMM的主要特征为：l 完全可压缩，具有静力选项的非静力模式；l Arakawa E-网格；l 水平传播的快波采用向前向后方案，垂直传播的声波采用隐式方法，水平平流采用Adams-Bashforth方案，垂直平流采用Crank-Nicholson方案；l 大量一二阶量的守恒，包括能量和位涡拟能；l 多物理选项；l 单方向和双向嵌套。WRF的后处理(WPP)：这个子程序可以用来处理WRF-ARW和WRF-NMM预报，设计如下：l 把模式垂直坐标差值到NWS标准输出层；l 把预报网格差

5、值到正常网格；l 计算诊断输出量；l 输出NWS和WMO标准的GRIB1。2. WRF的预处理系统（WPS）2.1介绍WRF前处理系统（WPS）是一个由三个程序组成的模块，这三个程序的作用是为真实数据模拟准备输入场。三个程序的各自用途为：geogrid确定模式区域并把静态地形数据插值到格点；ungrib从GRIB格式的数据中提取气象要素场；metgird则是把提取出的气象要素场水平插值到由geogrid确定的网格点上。把气象要素场垂直方向插值到WRF eta层则是WRF模块中的real程序的工作。上图给出了数据在WPS的三个程序之间的转换关系。正如图像所示，WPS里每个程序都会从一个共同的na

6、melist文件里读取参数。这个namelist文件按各个程序所需参数的不同分成了三个各自的记录部分及一个共享部分，它们分别定义了WPS系统所要用到的各种参数。被三个程序各自用到的表格文件没有在图中显示出来。尽管这些表格无需用户改动，但是这些个表格却提供了控制程序运行的额外信息。GEOGRID.TBL, METGRID.TBL,和Vtable文件将会在后文中被详细介绍。 C2 F4 w7 H4 9 R7 a1 d$ Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟安装WPS的步骤和安装WRF的步骤基本相同，都提供了

7、编译的选项，只是平台有所变化。当MPICH库及合适的库可以使用时，metgird和geogrid程序可以用分布式内存来编译，如果是这样操作，那当用户在设置大的模拟区域时就可以花更少的时间。但是ungrib程序却不能使用并行，因此只能用单CPU来操作。2.2 各个WPS程序的功能WPS包含三个独立的程序：geogrid, ungrib, 和 metgrid。下面给出三个主要程序的简要介绍。2.2.1程序geogridgeogrid的目的是确定模拟区域，及把各种地形数据集插值到模式格点上。模拟区域的确定是通过设置namelist.wps文件中的与“geogrid”有关的参数来实现的。除了计算经纬度

8、和地图每个格点的比例因子外，geogrid还会根据默认值来插值土壤类型、地表利用类型、地形高度、年平均深层土壤温度、月季植被覆盖、月季反照率、最大的积雪反照率及斜坡的类别。可以通过WRF的官方网站来下载这些场的全球数据集，而且这些资料可以被看成是不随时间改变的，因此只需下一次就可以了。很多数据集只有在某一特定精度上才能使用，但是其它的则可在30,2, 5, 和10中的任何一个精度使用；其中的意思是弧秒，而则是代表弧分（具体见下面的表格）。尽管一个与模拟区域分辨率接近的地形数据被应用后，这个插值后的地形场可能更具代表性，但是用户仍是无需下载所有可用的分辨率的数据（SJ认为还是全下载的好）。尽管如

9、此，那些希望应用拥有可以覆盖大范围的格点空间的区域的用户，还是愿意下载所有分辨率的静态地形数据。. TJ) F2 E1 R( f( n除了插值默认的静态数据，geogrid程序还可以插值进更多的连续的且不同种类的地形到模拟区域。可以通过应用表格文件GEOGRID.TBL来实现插值新的或额外的数据集到模拟区域。GEOGRID.TBL文件定义了所有可以被geogrid生成的地形场；它描述了插值一个地形场所需的方法，及所需数据所放的具体位置。由geogrid生成的文件的格式是WRF I/O API，因此可以通过选择NetCDF I/O格式使geogrid生成NetCDF格式的输出文件，以便更方便的用

10、一些外部软件ncview，NCL和最新版本的RIP4来实现可视化（画出地形图）。2.2.2程序ungribungrib程序读GRIB文件，“抽出”数据，然后把它们用一个简单的格式写出来，这种格式就是“过渡”格式（具体参考英文文献，及writing data to the intermediate format一节以获得关于这个格式的更多细节）。GRIB文件包含随时间变化的气象要素场，而且它是从其它区域或全球模式如NCEP的NAM或GFS模式而来的。ungrib程序可以读GRIB1的数据，当然，如果编译时选择了GRIB2的选项，则就可以读GRIB2的数据。* V. / : N5 ?6 H1 rM

11、eteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟GRIB文件包含的要素场多于启动WRF所需的数量。两种不同格式的GRIB用了不同的编码来确定变量和在GRIB文件中的层次。ungrib用这些编码表格Vtable（variable tables）来确定哪些场需要从GRIB文件里提取出来并写成过渡格式。关于编码的细节可以在WMO GRIB文档中找到。对于相同GRIB模式输出文件的各种Vtables可以在WPS主目录下的/ungrib/Variable_Tables/里。WPS为NAM104和212格点、NAM AWIP格式

12、、GFS、NCEP/NCAR再分析、RUC（气压坐标数据和混合坐标数据），AFWA的AGRMET地表模式输出，ECMWF及其它数据集。用户可以以其它Vtable做模板来为其它的模式输出创造自己的Vtable。Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟- Z% W$ u& g! ) z* PUngrib可以用三个用户可选格式中任何一个来写过渡数据，这三个格式是：WPS一个新的格式，它包含了对接下来的程序有用的额外信息；SIWRF之前使用过的过渡格式；MM5用来向MM5模式输入GRIB2数据。尽管WPS格式是被推

13、荐使用的，但是这三个中的任何一个格式都可以被用来启动WRF。Meteorologic2.2.3 metgrid程序metgrid程序的作用是把ungrib程序提取出的气象要素场水平插值到geogrid确定的模拟区域上。这个插值后的数据可以被WRF的real程序所识别并吸收。metgird插值的那些数据的时间段可以通过设置namelist.wps中share记录部分来调整，而且每个模拟区域（最外围区和嵌套区）的时间都要单独设置。与ungrib程序一样，metgird所处理的数据也是随时间改变的，因此每次做新的模拟时，都要运行metgird程序。! I% |/ Z* i/ zwww.mnmuc.o

14、rgMETGRID.TBL文件是用来控制如何把气象要素场进行插值的。METGRID.TBL文件为每个要素场都提供了一个区间，在这个区间里，可能会确定诸如要素场的插值方式、作为标记插值以及要素场所要插值的网格（如ARW的U，V；NMM的H，V）。由metgrid生成的文件的格式是WRF I/O API，因此可以通过选择NetCDF I/O格式使metgrid生成NetCDF格式的输出文件，以便更方便的用一些外部软件ncview，NCL和最新版本的RIP4来实现可视化（画出地形图）。2.3运行WPS运行WRF Preprocessing System（WPS）有如下三个步骤：+ i9 k4 x 0

15、 k% $ q! O1 利用geogrid模块确定一个模式的粗糙区域（最外围的范围），及其它嵌套区域；) y4 X) ?( _5 e6 o E2 利用ungrib把模拟期间所需的气象要素场从GRIB资料集中提取出来；# O+ N2 g9 L# 6 n5 F- P5 v mMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟3 利用metgird把上述的气象要素场（第二步所做的工作）水平插值到模式区域（第一步所做的工作）中。2 z. w( : $ A- L$ B) t当多个模拟在同一区域重复进行时，只需要做一次第一步的工

16、作即可（也就是说geogrid.exe所做出的地形资料geo_em.d0*.nc可以重复利用）；因此，只有随时间改变的数据才需要在每次模拟时用第二、三步来处理。类似的，如果在多次模拟中，气象数据是类似的，但是地形区域却不断改变的话，那第二步是可以省略的。下面是各个步骤的详细说明：5 q2 Z6 f: b# UStep 1: 确定模式区域- Y, _# C( n9 qMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟如果WPS安装成功，会在WPS根目录下出现三个可执行程序geogrid.exe, ungrib.exe和

17、metgird.exe（原程序在各自的同名子目录下）的链接。除了这三个可执行程序的链接，还有一个namelist.wps文件。下面是一个WPS主目录下各组件的列表$2 j$ r C4 s8 b?: v0 A/ s$ Hdrwxr-xr-x 2 4096 arch$ ( T. y( o/ CW U2 B* k( -rwxr-xr-x 1 1672 4 C+ w5 i% m. |7 Z5 X p-rwxr-xr-x 1 3510 compileI5 # A; S) w$ N4 -rw-r-r- 1 859

18、73 compile.output气象,数值,模式,气象数值模式,数值模式,气象模式, -rwxr-xr-x 1 4257 configure8 q! N2 M3 W2 y g气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛-rw-r-r- 1 2486 configure.wps% k& b8 a W# t1 r4 U7 W xdrwxr-xr-x 4 4096 geogridMeteorological Numerical Model Union of China (MNE6 V

19、- P8 lrwxrwxrwx 1 23 geogrid.exe - geogrid/src/geogrid.exe! d g0 ?( d# q- C) c9 w6 Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟-rwxr-xr-x 1 1328 link_grib.csh气象,数值,模式,气象数值模式,数值模式,气象模式,海drwxr-xr-x 3 4096 4 v( Q8 _w, b& W& 9 q# ! Ilrwxrwxrwx 1 23 metgrid.exe - m

20、etgrid/src/metgrid.exe# S B% U& a7 4 3 p气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛-rw-r-r- 1 1101 namelist.wps9 v P, l; c5 m5 F e-rw-r-r- 1 1987 namelist.wps.all_options$ U$ W- N- K: Z# + BMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟-rw-r-r-

21、1 1075 namelist.wps.global气象,数值,模式,气象数值模式,数值模式,气象模-rw-r-r- 1 652 namelist.wps.nmm w3 R( w: l) h5 h! oZ* Y/ P-rw-r-r- 1 4786 README8 B! m4 q% T2 |+ V2 c气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛drwxr-xr-x 4 4096 ungrib% t4 s1 o6 IA Y. G8 klrwxrwxrwx 1 21 ungrib.e

22、xe - ungrib/src/ungrib.exe+ L7 l p; C( V, |drwxr-xr-x 3 4096 util气象,模式的最外层区域和其它嵌套区域都是在namelist.list中geogrid记录里设置的，另外，在share记录里也有需要设置的参数。下面是关于这两个部分的一个模板，并且如果想获得更多关于每个参数的可能设置及其目的的信息，可以参考description of namelist variables&share2 t0 e+ ?4 O% c- A3 L- Z气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HY

23、COM,LICOM,同化,集合预报,模式联盟,模式论坛wrf_core = ARW,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模zmax_dom = 2,+ F3 U6 ND; u, N- M Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟start_date = 2008-03-24_12:00:00,2008-03-24_12:00:00,Meteorological end_date = 2008-03-24_18:00:00,2008-03-24_12:00:00,) i7 E%

24、x o5 L V6 FMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟interval_seconds = 21600,P0 X% i9 . W# f. Cio_form_geogrid = 2- |- Y! g5 Q) S) I* Z$ Y: t气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛/8 N6 |1 Z( G Q3 o# T% R+ |&geogrid- i. H-

25、e3 Y- H, V9 Dparent_id = 1, 1,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候parent_grid_ratio = 1, 3,; l9 K) / a6 g2 J( Ji_parent_start = 1, 31,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式, j_parent_start = 1, 17,* % l& eX5 . v5 s_we = 1, 1,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气e_we = 74, 112,. z8 / c& s) J( M. VMeteorological Nume

26、rical Model Union of China (MNMUC) 中国气象数值模式联盟s_sn = 1, 1,. n* j: L( 8 U6 e_sn = 61, 97, k( & e1 d! q! I! 1 X, t6 j geog_data_res = 10m,2m,Meteodx = 30000,Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模dy = 30000,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模map_proj = lam

27、bert,! w4 U+ I|% X+ ref_lat = 34.83,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模z2 Jref_lon = -81.03,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模truelat1 = 30.0,: t$ K) i; C/ b u气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛truelat2 = 60.0,% E/ g/ E8 v; N1 ( y; X气象,数值,模式,气象数值模式,数值模式,气象

28、模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛stand_lon = -98.,+ a1 u4 o- Q6 w6 fi8 t& 气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛geog_data_path = /mmm/users/wrfhelp/WPS_GEOG/9 7 j: n) O6 t/为了总结一些与geogrid有关的share记录中的典型改变，与WRF的动力核有关的wrf_core要首先被选择。如果WPS

29、要为ARW（Advanced Research WRF）模拟而运行，那wrf_core就设成ARW，如果要为NMM（Nonhydrostatic Mesoscale Model）模拟的话，则设成NMM。当选择好动力内核后，接下来选择max_dom，即区域（最外层的一个+嵌套数）的总数（当wrf_core=ARW）或者嵌套的层次（当wrf_core=NMM）。因为geogrid生成的仅仅是时间独立的数据，因此start_date,end_date, 和 interval_seconds这些参数将被其忽略。另外，还有一些可选的选项，如opt_output_from_geogrid_path，如果设

30、成默认值，则由geogrid.exe生成的地形文件将被放到当前工作目录（WPS的主目录），如果想放到别的目录下，则根据需要修改即可；io_form_geogrid则是设置地形数据输出格式的。在geogrid的记录部分，是关于模拟区域投影的设置，同时也设置了模式格点的大小和所在位置。模式所用的地图投影方式由map_proj来设置，其它用来设置投影的参数总结如下：* |0 T9 m! b9 T MMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟地图投影方式/与相应投影方式有关的参数变量气象,数值,模式,气象数值模式,数

31、值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛- z5 p1 m, u% Y8 D e3 + M2 a0 Ilambert truelat1F0 B$ j( G/ L/ U! Ctruelat2 (optional): fX- 8 m( h* v2 stand_lon1 d/ X+ i3 h9 j$ Lc/ mercator truelat1M1 a& x8 k) dMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模

32、式联盟polar truelat1# T7 _/ T X% / R& f气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛stand_lon+ H% d/ C+ w3 h: n. M% . xlat-lon pole_lat气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛/ S/ ) * ( d- r7 u9 r2 A, P9 ypole_)

33、m1 8 g5 F- W$ Ostand_lon% ? G$ k$ u; Ib/ U如果WRF是在一个局地区域里运行，那粗糙区域（最外层区域）位置则是通过ref_lat和ref_lon来定位，它们分别确定了粗糙区域的纬度和经度。如果也要处理嵌套区域，则它们的位置是通过i_parent_start 和j_parent_start来确定；更多关于设置嵌套区域的细节可以参考下一小节。接下来，粗糙区域的维数由dx和dy来确定，它们分别确定了x轴和y轴上标准格距的长度，而e_we和e_sn则分别给出了x轴（东西方向）和y轴（南北方向）上的格点数（具体设置方案可见后边description of name

34、list variables）；对应lambert, mercator, 和 polar投影方式，dx和dy的单位是米，对于lat-lon投影方式，dx和dy的单位则是度。对于嵌套区域，只有e_we和e_sn可以用来确定格点的维数，格点的dx和dy是不能被设定的，因为它们的值已经被parent_grid_ratio和parent_id所决定了，这两个参数分别确定了嵌套的上一级区域（父区域）的格点距离与嵌套的格点距离的比值和嵌套的上一级区域的ID。8 U+ |: ?% j; l7 v5 j) p0 r对于全球的模拟，粗糙区域（最外围区域）的覆盖范围就应该是全球，所以re

35、f_lat和ref_lon就不在被使用，而且dx和dy也不应该再被设置，这是因为这个格距将根据格点数被自动计算出来。同样需要注意的是，经纬度投影（map_proj = lat-lon）是WRF中唯一支持全球区域（全球模拟）的投影方式。: Z% v% 9 P* J# h4 y) d qn! m除了设置与模式区域的投影方式、位置和覆盖范围有关的参数外，静态数据集的路径也必须通过参数geog_data_path被正确的设置。同样，用户可能会通过设置geog_data_res参数来选择静态数据的分辨率以便于geogrid的差值，而分辨率的值要与GEOGRID.TBL中的数据分辨率吻合。如果静态数据集是

36、从WRF官网下载的，则可能的分辨率会包括30s, 2m,5m, 和 10m,分别对应的是30秒，2-，5-，10-分的数据。根据wrf_core所设值的不同，geogrid需要不同的GEOGRID.TBL与之对应，因为WPS会根据动力内核的不同插值出不同的网格。如果wrf_core=ARW，则应该使用GEOGRID.TBL.ARW，如果是wrf_core=NMM，则相应的要用GEOGRID.TBL.NMM。可以通过链接正确的GEOGRID.TBL文件来正确使用它，原始的文件放在geogrid目录里（另一种方法是通过opt_geogrid_tbl_path参数来设定文件的位置来使用它）。! F9

37、 Yz+ . o& 8 x气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛$ ls geogrid/GEOGRID.TBL9 l8 zW6 b* _) k6 x7 y lrwxrwxrwx 1 15 GEOGRID.TBL - GEOGRID.TBL.ARW) p1 a! N- |2 M$ R1 0 o) s; H- u如果想了解更多的关于各个参数的意义和可能的选项，可以参考description of the namelist variables。) $

38、c0 T& D0 Y$ d+ u8 T在namelist.wps中设置了合适的模拟和嵌套区域，geogrid.exe就可以处理地形文件了。在ARW的个例中，地形文件被命名为geo_em.d0N.nc，其中每个文件中的N代表嵌套的序号。当运行NMM个例时，粗糙区域（外围区）的地形数据文件的名字为geo_nmm.d01.nc，嵌套的文件则是geo_nmm_nest.l0N.nc，其中N代表嵌套的级数。还需要注意的是，文件的后缀会随参数io_form_geogrid设定的不同而不同。运行geogrid.exe的命令是：气象,数值,模式,气象数值模式,数值、$ ./geogrid.exewww.mnm

39、1 M A9 M1 W* n) & C: M Z$ F i当程序运行成功后，会出现如下信息 J2 . r! B/ O* q!& V9 d% o1 G & $ g) m, h- K! Successful completion of geogrid. !气象,数值,模式,气象数值模式,数值模式,气象模!2 g h4 6 i% Z: ?气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛并且WPS的主目录下会出现地形文件（或者出

40、现在opt_output_from_geogrid_path所设定的目录下）。如果没有出现，可以通过检查geogrid.log来找出可能的失败原因。如果想了解更多的关于检查geogrid输出结果的信息，可以参考2.20WPS 输出场一节。9 : M) g- b a: k5 J; L$ls8 B# e% z! A# J: G9 cdrwxr-xr-x 2 4096 arch4 J/ P. t) d5 b$ n& u4 7 r气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛-rwxr-x

41、r-x 1 1672 clean6 a1 M# t8 VZ0 b( O Z: -rwxr-xr-x 1 3510 5 S+ j1 k) U0 B6 B: x-rw-r-r- 1 85973 compile.output, v8 % V% O3 e2 z: UK7 W0 T-rwxr-xr-x 1 4257 configure气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛8 d: ! P! * r0 _: Y3 (

42、n: Z-rw-r-r- 1 2486 configure.wps5 # L6 P1 O # w4 k-rw-r-r- 1 geo_em.d01.ncMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟2 l+ m* B/ i 8 V- f-rw-r-r- 1 geo_em.d02.nc9 Q5 X: g# / W6 2 W气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛drwxr-xr-x 4 4096 g

43、eogrid1 g( L- P% J3 Ylrwxrwxrwx 1 23 geogrid.exe - geogrid/src/geogrid.exe# C$ R! ( k1 C( e3 b! K-rw-r-r- 1 11169 geogrid.log2 v( H p5 z% J, M8 $ w6 -rwxr-xr-x 1 1328 link_grib.csh! k9 ! ?8 d) Jdrwxr-xr-x 3 4096 metgridMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟- 5 D V+ |8 _lrwx

44、rwxrwx 1 23 metgrid.exe - metgrid/src/metgrid.exe* F0 T! 5 U9 L6 ( r4 z. -rw-r-r- 1 1094 namelist.wps9 l* e7 r2 I# S, m-rw-r-r- 1 1987 namelist.wps.all_options气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛! V: P, # A x* ?1 X-rw-r-r- 1 1075 namelist.wps

45、.global$ X# X; N/ . b( j% z) y% w/ T/ -rw-r-r- 1 652 6 S$ a6 U! T * z I: ( R5 -rw-r-r- 1 4786 README4 6 Y& p R! x+ & c# _# , f3 drwxr-xr-x 4 4096 ungrib: g& V) 7 y0 |* Rlrwxrwxrwx 1 21 ungrib.exe - ungrib/src/ungrib.exe* m6 T; I2 j% e5 M6 G% r v: drwxr-xr-x 3 4096 util2

46、A & e+ z0 g; A5 8 X4 I+ zMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟Step 2:利用ungrib从GRIB文件中提取气象要素场气象,数值当已经下载了GRIB格式的气象数据后，提取要素场以转成过渡格式的第一步所要做的参数设置涉及到“share”和“ungrib”两个部分。下面是关于这两部分的一个模板：气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛, x b2 |& f; C

47、/ w, H4 |. ?&V8 5 a0 % e6 w9 wrf_core = ARW,9 f% c4 : y3 R# F% % N气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛max_dom = 2,Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟7 K6 q8 q J) K0 s$ p9 sstart_date = 2008-03-24_12:00:00,2008

48、-03-24_12:00:00,- e: c% l* b- C0 K3 _气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛end_date = 2008-03-24_18:00:00,2008-03-24_12:00:00,气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛( S& P/ / % 2 ) f3 ointerval_seconds = 21600, L- L%

49、D9 C, |3 * N. X气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛io_form_geogrid = 2 z8 eZ8 S2 r$ ?2 LO! a5 Q7 l1 g8 N气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛/% Z l; L$ - L3 & |2 |&ungrib5 l5 K) Y. : _, EMeteorologica

50、l Numerical Model Union of China (MNMUC) 中国气象数值模式联盟out_format = WPS,! g& E+ w6 d, E) _prefix = FILE气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛& g; F6 Z# K* e% N& /7 Y! h: G! m- a在“share”部分，与ungrib有关的参数是粗糙区域（外围区）的开始和结束时间（start_date和end_date; 或者start_y

51、ear, start_month, start_day, start_hour, end_year, end_month, end_day和end_hour）以及气象数据文件的时间间隔（interval_seconds）。在“ungrib”部分，参数out_format是用来设置过渡数据文件的格式的；metgird程序可以读任何一个被ungrib支持的格式，因此虽然WPS是被推荐的格式,但是WPS，SI, 和MM5中的任何一个都可以作为选项被选择；另外，用户可以通过设置prefix参数来确定过渡数据文件的保存路径和前缀。例如，如果prefix被设成ARGRMET，那过渡数据的文件名就是AGRM

52、ET:YYYY-MM-DD_HH，其中YYYY-MM-DD_HH是文件中数据的真实的日期。# f6 F% 5 Q7 u- e3 k5 bdMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟在合理的修改了namelist.wps文件后，Vtable文件也是需要事先被提供的，并且GRIB文件也要被链接（拷贝）到ungrib.exe所期望的目录下（一般是WPS的主目录，如果WPS被安装成功的话）。如果WPS要处理气象数据，就要在Vtable的辅助下进行才行，即使Vtabel仅仅被象征性地以Vtable来命名并链接到ungr

53、ib.exe所期望的目录下（一般是WPS的主目录）。例如，如果是来自GFS模式的GRIB数据，可以如此处理：7 q/ w. f9 A5 x- Z1 e% T$ ln -s ungrib/Variable_Tables/Vtable.GFS Vtable气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛b4 I( W ) w& |3 lungrib程序将尝试以GRIBFILE.AAA,GRIBFILE.AAB, , GRIBFILE.ZZZ这样的文件名来读取GR

54、IB数据。为了能简化把GRIB数据链接到到合适的目录，并以这些名字来命名，程序提供了一个shell的脚本link_grib.csh。这个脚本通过读取一个以GRIB数据文件名为列表的命令行参数来完成上述工作。例如，如果GRIB数据被下载到一个/data/gfs的目录下，那这些文件可以用下面的方法来进行链接：8 9 e; V! 5 * J气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛$ ls /data/gfs+ R0 M4 z W+ Q; R- j! z-rw-r-r- 1 gfs_

55、12_00* q 9 ! o: m% U7 X% K6 e-rw-r-r- 1 gfs_12_06气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛 V i9 F( z: Z$ ./link_grib.csh /data/gfs/gfs*+ e$ j9 N/ . 上述操作完成后，WPS目录下的列表将是如下所示：6 U& T$ q9 a# |6 I% h气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合

56、预报,模式联盟,模式论坛drwxr-xr-x 2 4096 arch: Y. A* Q- b* M# l气象,数值,模式,气象数值模式,数值模式,气象模式,海洋模式,气候模式,WRF,MM5,Grapes,HYCOM,LICOM,同化,集合预报,模式联盟,模式论坛-rwxr-xr-x 1 1672 ! ?. b0 1 D1 p, t1 ?-rwxr-xr-x 1 3510 compile& s, A& I; _, NMeteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟-rw-r-r- 1 85973 compile.output# I( P) R6 x1 W& V- G, c+ P& Q. V, $ W-rwxr-xr-x 1 4257 configure, M p3 I6 c1 K) S-rw-r-r- 1 2486 configure.wps% G, S; u+ j2 f6 m-rw-r-r- 1 geo_em