ln –sf /dgpfs/fs2/wrf_xp/grib/ grib
2) 根据数据来源模式具体情况，选择与此模式相应的变量列表文件。这里以GFS模式的
气象要素场为例：
ln –sf ungrib/Variable_Tables/Vtable.GFS Vtable
3) 连接GRIB码数据文件成GRIB.？？？格式 link_grib.csh grib/gfs 4) 编辑参数列表文件namelist.wps
&share
wrf_core = 'ARW', max_dom = 1,
start_date = '2006-08-16_12:00:00', end_date = '2006-08-16_18:00:00', interval_seconds = 21600 io_form_geogrid = 2,
opt_output_from_geogrid_path = './', debug_level = 0 /
&ungrib
out_format = 'WPS' /
提示：查看GRIB码数据的头信息可以用g1print.exe命令： util/g1print.exe GRIB.AAA
5) 创建运行ungrib.exe程序的作业卡
vi ungrib.cmd
创建后的ungrib.cmd的文件内容格式如下：
#################################################### #
# run the program ungrib.exe of WPS #
#################################################### #@ job_type = serial #@ class = interactive
#@ initialdir = /dgpfs/fs2/wrftc?/denglt/WPS #@ executable = ungrib.exe #@ input = /dev/null #@ output =ungrib.out #@ error =ungrib.err #@ notification = never
#@ wall_clock_limit = 00:10:00 #@ queue
6) 提交作业卡ungrib.cmd
6
llsubmit ungrib.cmd 提示：
a) ungrib的结果文件“FILE:…”可以用WPS的工具软件rd_intermediate.exe程序进行
查看，以验证数据的正确性。查看时，只需执行rd_intermediate.exe文件，并在这个可执行文件后面加上要查看的数据文件名。例如： util/rd_intermediate.exe FILE:2005-07-30_12
b) 用WPS的工具软件plotfmt.exe可以绘制中间文件的图形。具体操作是：
? 执行plotfmt.exe程序，并带上数据文件名作为输入参数。例如： util/plotfmt.exe FILE:2005-07-30_12 ? 运行图形查看软件idt
idt gmeta
? 弹出窗口里，点击“”键。
3． 运行metgrid.exe，对解码后的数据进行水平插值 1) 编辑参数列表文件namelist.wps
&share
wrf_core = 'ARW', max_dom = 1,
start_date = '2006-08-16_12:00:00', end_date = '2006-08-16_18:00:00', interval_seconds = 21600 io_form_geogrid = 2,
opt_output_from_geogrid_path = './', debug_level = 0 /
&metgrid
fg_name = './FILE' io_form_metgrid = 2,
opt_output_from_metgrid_path = './', opt_metgrid_ tbl_path = 'metgrid/', /
2) 创建运行metgrid.exe程序的作业卡 vi metgrid.cmd
创建后，metgrid.cmd文件的内容格式为：
#################################################### #
# run the program metgrid.exe of WPS #
#################################################### #@ job_type = parallel #@ class = normald
#@ initialdir = /dgpfs/fs2/wrftc?/denglt/WPS
7
#@ executable = metgrid.exe #@ input = /dev/null #@ output =metgrid.out #@ error =metgrid.err #@ node = 1
#@ tasks_per_node = 8 #@ node_usage = shared
#@ network.MPI = css0,shared,US #@ notification = never
#@ wall_clock_limit = 00:10:00 #@ queue
3） 提交作业卡metgrid.cmd llsubmit metgrid.cmd
运行成功后，在opt_output_from_metgrid_path设置的路径下会产生
met_em.d01.yyyy-mm-dd_hh:00:00.nc格式的文件。
三. 运行WRF 模式主体 1． 运行初始化程序real
1) 编辑参数列表文件namelist.input
WRF模式有专门运行real.exe和wrf.exe程序的路径（WRFV3/run或者WRFV3/test/em_real）。我们转到运行路径下
cd WRFV3/test/em_real/
然后就可以编辑参数列表文件
vi namelist.input
2) 连接real.exe的输入场数据
real.exe的输入场数据文件就是metgrid的输出文件，文件名格式为“met_em.d01.*”。在运行real.exe前，必须将这些文件拷贝或者连接到real.exe的运行路径 (test/em_real) 下面。
ln -sf ../../WPS/met_em.d01.* . vi real.exe
3) 创建运行real.exe程序的作业卡
提示：根据编译WRF模式时的选项情况，我们可知real.exe编译后是串行还是并行程序。根据real.exe的编译情况，写出相应的串行或者并行作业卡。 4) 提交作业
./real.exe
2． 运行WRF模式
1） 编辑参数列表文件namelist.input
进入到WRF的运行路径（WRFV2/run或者test/em_real）
cd WRFV3/run
8
编辑参数列表文件（一般是在real.exe运行完以后，wrf.exe的参数列表就不做任何变动）
vi namelist.input
2） 创建运行wrf.exe程序的作业卡job-submission.sh #!/bin/bash #PBS -q high
#PBS -l nodes=2:ppn=24 #PBS -j oe
cd $PBS_O_WORKDIR
cat \NP=`cat $PBS_NODEFILE | wc -l`
mpirun -np 48 -machinefile nodes wrf.exe
运行成功后，在real.exe的运行路径 (test/em_real) 下面会产生wrfinput_d01和wrfbdy_d01格式的文件。
提示：根据编译WRF模式时的选项情况，我们可知wrf.exe编译后是串行还是并行程序。根据wrf.exe的编译情况，写出相应的串行或者并行作业卡。 提交作业卡wrf.cmd
qsub job-submission.sh
运行成功后，在wrf.exe的运行路径 (test/em_real) 下面会产生wrfout_d01_*格式的文件。
四. 运行后处理程序 1． 编译安装程序 1） 获取程序
进入到自己的运行目录
cd /dgpfs/fs2/wrftc?/denglt
拷贝程序包
cp /u/wrf_xp/src/wrfpost20070829.tar
释放程序包
tar －xvf wrfpost20070829.tar
2） 配置编译环境
进入到后处理程序目录
cd wrfpost
编辑编译规则文件makefile，进行编译环境配置
vi makefile
3） 编译程序
9
在src目录下，编译MICAPS处理程序
make micaps
在src目录下，编译GrADS处理程序
make grads
2． 处理MICAPS数据
1) 编辑参数配置文件namelist.input
进入MICAPS处理路径
cd micaps
编辑参数配置文件
vi namelist.input
参数列表文件格式如下 !input_files
/dgpfs/fs2/wrftc?/bj15km/wvar_ubp_d01_2004-07-10_00:00:00 < --- 结果文件名列表 !end_input_files
!Global attributes
DATA_ROOT = \ ! 处理后的数据存放路径 interp_horizontal = 2, ! 水平插值方式
! 0: no interp, 1:convert to Lat/Lon grids
! 2: no interp, but output in lambert map project interp_vertical = 1, ! 垂直插值方式
! 0: no interp, 1: pressure levels, 2:height levels surface_levels = \ ! 垂直层次值
ntimes = \ ! 要从WRFOUT中读取的时次 grid_x = 228, ! 后处理的输出东西格点数 grid_y = 222, ! 后处理的输出东西格点数 debug = 0, ! 打印调试信息
lonlat_user = 0, ! 是否在模式区域中挖取一小块数据区进行输出 min_lat = 30, max_lat = 45, min_lon = 105, max_lon = 124.95, 2) 运行处理程序
proc.sh
运行成功后，在DATA_ROOT的路径下会存在一系列MICAPS格式的数据文件。
3． 处理GrADS数据（重点，可以讲一下GrADS操作方法）
10
1） 编辑参数列表文件namelist.input
进入GrADS处理路径
cd grads
编辑参数配置文件
vi namelist.input
2） 运行处理程序
proc.sh
运行成功后，在DATA_ROOT的路径下会存在一系列MICAPS格式的数据文件。
五. 应用T213场库数据运行WRF模式 1． 编译安装程序
1) 获取源程序包 (获取源程序代码可从WRF的官方网站下载)
cp /u/wrf_xp/src/ readbjgl20070522.tar.gz ./ 2) 解压源程序压缩包
gunzip readbjgl20070522.tar.gz 3) 释放源程序包
tar –xvf readbjgl20070522.tar 4) 进入释放后的源程序目录
cd readbjgl 5) 配置编译环境
vi makefile 6) 编译程序
make
2． 读取T213场库数据
3． 运行模式
11