实用标准文档
x方向的格距，指不投影变形处的格距。
REAL*4 :: dy
y方向的格距，指不投影变形处的格距。
REAL*4 :: orient
与y轴平行的经度值，通常为投影不变形的经度值。
REAL*4 :: truelat1
投影不变形的第一纬度值。
REAL*4 :: truelat2
投影不变形的第二纬度值，只有兰勃特投影时才有第二不变形纬度值。
注意，第三版本的数据格式与第四版本的数据格式大部分相同，只是第四版本的数据比第三版本的数据多了两个字符型记录元：source和startloc。
(2) 运行MM5前处理程序，生成MM5的初始输入场和边界条件。 (3) 运行MM52WRF转换程序，将MM5的初始出入场和边界条件转换成WRF所需要
的输入数据文件real_input* 。
T213场库数据 readt213 MM5 grib_prep.exe gridgen_model.exhinterp.exe REGRID TERRAIN hinterp.exe INTERPF vinterp.exe vinterp.exe MM52WRF WRF 图2. WRF与T213嵌套实现流程示意图
文案大全
实用标准文档
3.4 编译嵌套程序
编译嵌套程序时，目前必须用64位形式编译，或者用编译选项“-qintsize=8”和“-qrealsize=8”来编译。另外，场库连接的库函数以及路径必须在编译连接时指定。IBM SP用户连接库为： ·场库在SP:/gpfs/fs2/nwp/gdbt213目录下 -L /u/nwp/gdbt213/lib -lpub
index_dic(1:28)='/u/nwp/gdbt213/dic/index_dic' ·场库在盘阵库/asm2/nwpgdb/t213gdb目录下 -L/u/nwp/gdbt213/lib -lpub_bx
index_dic(1:32)='/u/nwp/gdbt213/dic/index_nfs_dic’ 编译时，按照以下步骤进行：
(1) 进入到t213/mod目录下，输入make命令编译嵌套程序所需的模块程序； (2) 进入到t213/util目录下，输入make命令编译一些辅助子程序； (3) 回到t213/目录下，输入make命令编译主程序。
编译成功后，可产生一个可执行程序readt213.exe。
文案大全
实用标准文档
4. WRF模式系统的运行
4.1 理想大气方案
4.1.1 运行流程 运行WRF数据初始化程序ideal.exe
运行WRF模式主程序wrf.exe
wrfout
后处理
图3. WRF模式系统理想大气方案的运行流程
4.1.2 运行步骤
有多种理想大气试验方案可供试验运行。在编译好理想大气试验方案后，可以进入相应方案的目录下，如欧拉质量坐标x-z平面二维飚线方案对应的目录为./test/em_squall2d_x。
1. 设置该目录下的参数配置文件namelist.input；详细参数说明参见《WRF模式技术报告》中附录1。通过编辑namelist.input文件，可以修改积分的长度，输出的频次，计算范围的大小，时间步长，物理参数以前其他的一些参数。一些参数（比如时间步长，积分时间和输出时间，一些变化参数等）的改变，可以不需要重新执行“ideal.exe”程序。然而，当改变计算范围，边界条件，以及物理参数时，则必须重新执行“ideal.exe”程序。 2. 运行WRF模式的数据初始化程序ideal.exe。ideal.exe是串行运行的。如果要并行运行，则必须将并行节点和每个节点上的CPU数目都设置成1。运行正确后，可产生以下两个文件：
wrfbdy_d01 ： WRF的边界条件
wrfinput_d01 ： WRF的输入场
3. 运行WRF模式的主程序wrf.exe。wrf.exe可串行运行，也可并行运行。运行正确后，能产生以下形式的结果文件： wrfout_d01_000000
另外，还有一些程序运行时记录运行状态的中间记录文件。 4. WRF模式后处理
目前，已经初步地将WRF模式的输出结果转换成VIS5D、GrADS以及MICAPS的数据格式，显示WRF模式的输出结果。详细说明参见第四部分。
文案大全
实用标准文档
4.1.3 运行个例
通过编辑namelist.input文件，可以修改积分的长度，输出的频次，计算范围的大小，时间步长，物理参数以前其他的一些参数。一些参数（比如时间步长，积分时间和输出时间，一些变化参数等）的改变，可以不需要重新执行“ideal.exe”程序。然而，当改变计算范围，边界条件，以及物理参数时，则必须重新执行“ideal.exe”程序。
编辑原始程序代码，比如编辑dyn_eh/module_initialize_squall2d_x.F，已修改山体或圆形障碍物，或修改外围环境时，应该在WRF的顶级目录下重新运行“compile eh_squall2d_x”，重新编译程序代码。然后再象前面讲述的一样，分别执行“ideal.exe”和“wrf.exe”。
试验方案中，欧拉质量坐标x-z平面二维飚线方案采用Kessler微物理过程，粘性系数固定为300m2/s。并利用y方向的周期性边条件，从而使用三维的WRF模式模拟二维的飚线系统实际可行。在y方向，速度应该为零，并且没有其他物理量的计算输出。经过以上步骤地运行，得到结果文件有：
-rw-r--r-- 1 denglt numerica 7806948 Apr 14 06:33 wrfbdy -rw-r--r-- 1 denglt numerica 1630620 Apr 14 06:33 wrfinput -rw-r--r-- 1 denglt numerica 25498084 Apr 14 14:43 wrfout_01_000000
-rw-r--r-- 1 denglt numerica 7381448 Apr 14 14:05 wrfrst_01_000600
-rw-r--r-- 1 denglt numerica 7381448 Apr 14 14:43 wrfrst_01_001200
其中，wrfbdy和wrfinput两个文件是运行real.exe程序所产生的模式边界条件和初始场文件。wrfrst_01_000600和wrfrst_01_000600是模式“wrf.exe”运行中产生的中间文件，它的作用是使模式可以分各个时间段进行计算。wrfout_01_000000文件是模式的输出结果文件，从其中可以得到模式的各物理量的计算结果。图4是欧拉质量坐标二维(x-z)飚线方案模拟的水平风速的分布图。从模拟的图中，可以分析得出模式模拟的结果基本正确。
文案大全
实用标准文档
图4. 欧拉质量坐标二维(x-z)飚线方案模拟水平风速U分布
文案大全
WRF模式简易操作 - 中文指南设计 - 图文
