WRF官网案例运行01

本案例为WRF官网2001年1月的默认案例（官网参考链接）

安装好WRF之后开始进行WRF模式运行

WRF的预处理系统（WPS）

熟悉一下wps三个程序组成模块的作用（为模拟准备输入场）

（1）geogrid:确定模式区域并把静态地理数据插值到格点。

（2）ungrib:从GRIB格式的数据中提取气象要素场。

（3）metgrid:把提取出的气象要素场水平插值到由geogrid确定的网格点上。

0 运行geogrid.exe

首先我们要从WRF官网上下载静态地理数据
Static geography data
WRF建模系统能够创建理想的模拟，尽管大多数用户对实际数据案例感兴趣。要启动实际数据案例，必须创建域在地球上的物理位置以及该位置的静态信息。
这需要一个数据集，其中包括地形和土地用途分类等字段。移至Build_WRF目录，下载文件并解压缩。解压缩文件为“ geog”文件夹，将其重命名为“ WPS_GEOG”。
The tar file can be found at:
geog_complete.tar.bz2
geog_10m.tar.gz

不wget了，直接下载放到Build_WRF文件夹下面解压

tar -zxvf geog_10m.tar.gz

把解压出来的geog文件夹改名字成WPS_GEOG

mv geog WPS_GEOG

The directory infomation is given to the geogrid program in the namelist.wps file in the &geogrid section. The complete data expands to approximately 10 GB. This data allows a user to run the geogrid.exe program.

修改namelist.wps静态地理数据的打开路径

xiaomo@ubuntu:~/Build_WRF$ cd WPS
xiaomo@ubuntu:~/Build_WRF/WPS$ nano namelist.wps

终端在WPS目录下输入nano namelist.wps (非root用户需要输入密码）或者vi namelist.wps 进入编辑模式。

将max_dom后的嵌套修改为1。

将geog_data_path后面的路径修改为静态地理数据的储存路径
我的是把geog_data_path = ‘/glade/p/work/wrfhelp/WPS_GEOG/’
改成geog_data_path = ‘/home/xiaomo/Build_WRF/WPS_GEOG/’

保存后退出（ctrl+s,ctrl+x)

接着输入

./geogrid.exe

完成后会生成一个静态文件geo_em.d01.nc （如果有错误可以查看geogrid.log寻找error）

这时候可能遇到地理数据不全的情况：

xiaomo@ubuntu:~/Build_WRF/WPS$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1
ERROR: Could not open /home/xiaomo/Build_WRF/WPS_GEOG/orogwd_10m/con/index
application called MPI_Abort(MPI_COMM_WORLD, 22077) - process 0

这时候就缺啥补啥，哪里不会点哪里哈哈哈，如上面的例子就是缺少orogwd_10m文件

去https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.html这个网址下载相应数据即可

我少了不少，网又不好，下载了好久┭┮﹏┭┮，看下面这一个个的，还只是一部分！！！

xiaomo@ubuntu:~/Build_WRF/WPS$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1
ERROR: Could not open /home/xiaomo/Build_WRF/WPS_GEOG/lai_modis_10m/index
application called MPI_Abort(MPI_COMM_WORLD, 21953) - process 0
xiaomo@ubuntu:~/Build_WRF/WPS-4.1$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1
ERROR: Could not open /home/xiaomo/Build_WRF/WPS_GEOG/maxsnowalb_modis/index
application called MPI_Abort(MPI_COMM_WORLD, 22056) - process 0
xiaomo@ubuntu:~/Build_WRF/WPS-4.1$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1
ERROR: Could not open /home/xiaomo/Build_WRF/WPS_GEOG/orogwd_10m/con/index
application called MPI_Abort(MPI_COMM_WORLD, 22077) - process 0
xiaomo@ubuntu:~/Build_WRF/WPS$ ./geogrid.exe

好在最后成功了：

xiaomo@ubuntu:~/Build_WRF/WPS$ ./geogrid.exe
Parsed 28 entries in GEOGRID.TBL
Processing domain 1 of 1Processing XLAT and XLONGProcessing MAPFACProcessing F and EProcessing ROTANGProcessing LANDUSEFCalculating landmask from LANDUSEF ( WATER =  17  21 )Processing HGT_MProcessing SOILTEMPProcessing SOILCTOPProcessing SCT_DOMProcessing SOILCBOTProcessing SCB_DOMProcessing ALBEDO12MProcessing GREENFRACProcessing LAI12MProcessing SNOALBProcessing CONProcessing VARProcessing OA1Processing OA2Processing OA3Processing OA4Processing OL1Processing OL2Processing OL3Processing OL4Optional fields not processed by geogrid:VAR_SSO (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/varsso_10m/')LAKE_DEPTH (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/lake_depth/')URB_PARAM (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/NUDAPT44_1km/')FRC_URB2D (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/urbfrac_nlcd2011/')IMPERV (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/nlcd2011_imp_ll_9s/')CANFRA (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/nlcd2011_can_ll_9s/')EROD (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/erod/')CLAYFRAC (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/clayfrac_5m/')SANDFRAC (priority=1, resolution='default', path='/home/xiaomo/Build_WRF/WPS_GEOG/sandfrac_5m/')!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of geogrid.        !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

完成后会生成一个静态文件geo_em.d01.：

或者用ncview 软件查看.

#安装ncview
xiaomo@ubuntu:~/Build_WRF/WPS$ sudo apt install ncview
[sudo] password for xiaomo:
Reading package lists... Done
Building dependency tree
Reading state information... Done
The following packages were automatically installed and are no longer required:linux-headers-5.4.0-42 linux-headers-5.4.0-42-genericlinux-image-5.4.0-42-generic linux-modules-5.4.0-42-genericlinux-modules-extra-5.4.0-42-generic
Use 'sudo apt autoremove' to remove them.
The following additional packages will be installed:libaec0 libhdf5-103 libnetcdf15 libsz2 libudunits2-0 libudunits2-data
The following NEW packages will be installed:libaec0 libhdf5-103 libnetcdf15 libsz2 libudunits2-0 libudunits2-data ncview
0 upgraded, 7 newly installed, 0 to remove and 209 not upgraded.
Need to get 434 kB/2,179 kB of archives.
After this operation, 8,539 kB of additional disk space will be used.
Do you want to continue? [Y/n] y
Get:1 http://us.archive.ubuntu.com/ubuntu focal/universe amd64 ncview amd64 2.1.8+ds-3build2 [434 kB]
Get:1 http://us.archive.ubuntu.com/ubuntu focal/universe amd64 ncview amd64 2.1.8+ds-3build2 [434 kB]
Get:1 http://us.archive.ubuntu.com/ubuntu focal/universe amd64 ncview amd64 2.1.8+ds-3build2 [434 kB]
Get:1 http://us.archive.ubuntu.com/ubuntu focal/universe amd64 ncview amd64 2.1.8+ds-3build2 [434 kB]
Get:1 http://us.archive.ubuntu.com/ubuntu focal/universe amd64 ncview amd64 2.1.8+ds-3build2 [434 kB]
Fetched 15.6 kB in 5min 11s (50 B/s)
Selecting previously unselected package libaec0:amd64.
(Reading database ... 230851 files and directories currently installed.)
Preparing to unpack .../0-libaec0_1.0.4-1_amd64.deb ...
Unpacking libaec0:amd64 (1.0.4-1) ...
Selecting previously unselected package libsz2:amd64.
Preparing to unpack .../1-libsz2_1.0.4-1_amd64.deb ...
Unpacking libsz2:amd64 (1.0.4-1) ...
Selecting previously unselected package libhdf5-103:amd64.
Preparing to unpack .../2-libhdf5-103_1.10.4+repack-11ubuntu1_amd64.deb ...
Unpacking libhdf5-103:amd64 (1.10.4+repack-11ubuntu1) ...
Selecting previously unselected package libnetcdf15:amd64.
Preparing to unpack .../3-libnetcdf15_1%3a4.7.3-1_amd64.deb ...
Unpacking libnetcdf15:amd64 (1:4.7.3-1) ...
Selecting previously unselected package libudunits2-data.
Preparing to unpack .../4-libudunits2-data_2.2.26-5_all.deb ...
Unpacking libudunits2-data (2.2.26-5) ...
Selecting previously unselected package libudunits2-0:amd64.
Preparing to unpack .../5-libudunits2-0_2.2.26-5_amd64.deb ...
Unpacking libudunits2-0:amd64 (2.2.26-5) ...
Selecting previously unselected package ncview.
Preparing to unpack .../6-ncview_2.1.8+ds-3build2_amd64.deb ...
Unpacking ncview (2.1.8+ds-3build2) ...
Setting up libaec0:amd64 (1.0.4-1) ...
Setting up libudunits2-data (2.2.26-5) ...
Setting up libsz2:amd64 (1.0.4-1) ...
Setting up libhdf5-103:amd64 (1.10.4+repack-11ubuntu1) ...
Setting up libnetcdf15:amd64 (1:4.7.3-1) ...
Processing triggers for libc-bin (2.31-0ubuntu9) ...
Processing triggers for man-db (2.9.1-1) ...
Processing triggers for sgml-base (1.29.1) ...
Setting up libudunits2-0:amd64 (2.2.26-5) ...
Setting up ncview (2.1.8+ds-3build2) ...
Processing triggers for libc-bin (2.31-0ubuntu9) ...#查看
xiaomo@ubuntu:~/Build_WRF/WPS$ ncview geo_em.d01.nc
Ncview 2.1.8 David W. Pierce  8 March 2017
http://meteora.ucsd.edu:80/~pierce/ncview_home_page.html
Copyright (C) 1993 through 2015, David W. Pierce
Ncview comes with ABSOLUTELY NO WARRANTY; for details type `ncview -w'.
This is free software licensed under the Gnu General Public License version 3; type `ncview -c' for redistribution details.Note: could not open file /home/xiaomo/.ncviewrc for reading
Warning: Cannot convert string "-*-lucida-bold-r-*-*-14-*-*-*-*-*-*-*" to type FontStruct
X connection to :0 broken (explicit kill or server shutdown).

查看结果：

1 运行./ungrib.exe

在与WPS同级别目录下创建DATA文件夹

将从官网下载的JAN00_GFS_FNL.tar.gz文件放在其目录下.

解压文件

tar -zxvf JAN00_GFS_FNL.tar.gz

链接Variable_Table文件，运行ungrib.exe,解码GRIB数据，提取气象场。

首先通过一个数据文件运g1print

xiaomo@ubuntu:~/Build_WRF/WPS$ ./util/g1print.exe ../DATA/JAN00/fnl_20000124_12_00.grib1

#运行结果：
Copen: File = ../DATA/JAN00/fnl_20000124_12_00.grib1
Fortran Unit = 0
UNIX File descriptor: 3----------------------------------------------------rec GRIB GRIB  Lvl  Lvl  Lvl         Time      FcstNum Code name  Code one  two                   hour
----------------------------------------------------1   7 HGT      100 1000    0  2000-01-24_12:00 + 002   7 HGT      100  975    0  2000-01-24_12:00 + 003   7 HGT      100  950    0  2000-01-24_12:00 + 004   7 HGT      100  925    0  2000-01-24_12:00 + 005   7 HGT      100  900    0  2000-01-24_12:00 + 006   7 HGT      100  850    0  2000-01-24_12:00 + 007   7 HGT      100  800    0  2000-01-24_12:00 + 008   7 HGT      100  750    0  2000-01-24_12:00 + 009   7 HGT      100  700    0  2000-01-24_12:00 + 0010   7 HGT      100  650    0  2000-01-24_12:00 + 0011   7 HGT      100  600    0  2000-01-24_12:00 + 0012   7 HGT      100  550    0  2000-01-24_12:00 + 0013   7 HGT      100  500    0  2000-01-24_12:00 + 0014   7 HGT      100  450    0  2000-01-24_12:00 + 0015   7 HGT      100  400    0  2000-01-24_12:00 + 0016   7 HGT      100  350    0  2000-01-24_12:00 + 0017   7 HGT      100  300    0  2000-01-24_12:00 + 0018   7 HGT      100  250    0  2000-01-24_12:00 + 0019   7 HGT      100  200    0  2000-01-24_12:00 + 0020   7 HGT      100  150    0  2000-01-24_12:00 + 0021   7 HGT      100  100    0  2000-01-24_12:00 + 0022   7 HGT      100   70    0  2000-01-24_12:00 + 0023   7 HGT      100   50    0  2000-01-24_12:00 + 0024   7 HGT      100   30    0  2000-01-24_12:00 + 0025   7 HGT      100   20    0  2000-01-24_12:00 + 0026   7 HGT      100   10    0  2000-01-24_12:00 + 0027  11 TMP      100 1000    0  2000-01-24_12:00 + 0028  11 TMP      100  975    0  2000-01-24_12:00 + 0029  11 TMP      100  950    0  2000-01-24_12:00 + 0030  11 TMP      100  925    0  2000-01-24_12:00 + 0031  11 TMP      100  900    0  2000-01-24_12:00 + 0032  11 TMP      100  850    0  2000-01-24_12:00 + 0033  11 TMP      100  800    0  2000-01-24_12:00 + 0034  11 TMP      100  750    0  2000-01-24_12:00 + 0035  11 TMP      100  700    0  2000-01-24_12:00 + 0036  11 TMP      100  650    0  2000-01-24_12:00 + 0037  11 TMP      100  600    0  2000-01-24_12:00 + 0038  11 TMP      100  550    0  2000-01-24_12:00 + 0039  11 TMP      100  500    0  2000-01-24_12:00 + 0040  11 TMP      100  450    0  2000-01-24_12:00 + 0041  11 TMP      100  400    0  2000-01-24_12:00 + 0042  11 TMP      100  350    0  2000-01-24_12:00 + 0043  11 TMP      100  300    0  2000-01-24_12:00 + 0044  11 TMP      100  250    0  2000-01-24_12:00 + 0045  11 TMP      100  200    0  2000-01-24_12:00 + 0046  11 TMP      100  150    0  2000-01-24_12:00 + 0047  11 TMP      100  100    0  2000-01-24_12:00 + 0048  11 TMP      100   70    0  2000-01-24_12:00 + 0049  11 TMP      100   50    0  2000-01-24_12:00 + 0050  11 TMP      100   30    0  2000-01-24_12:00 + 0051  11 TMP      100   20    0  2000-01-24_12:00 + 0052  11 TMP      100   10    0  2000-01-24_12:00 + 0053  39 V VEL    100 1000    0  2000-01-24_12:00 + 0054  39 V VEL    100  975    0  2000-01-24_12:00 + 0055  39 V VEL    100  950    0  2000-01-24_12:00 + 0056  39 V VEL    100  925    0  2000-01-24_12:00 + 0057  39 V VEL    100  900    0  2000-01-24_12:00 + 0058  39 V VEL    100  850    0  2000-01-24_12:00 + 0059  39 V VEL    100  800    0  2000-01-24_12:00 + 0060  39 V VEL    100  750    0  2000-01-24_12:00 + 0061  39 V VEL    100  700    0  2000-01-24_12:00 + 0062  39 V VEL    100  650    0  2000-01-24_12:00 + 0063  39 V VEL    100  600    0  2000-01-24_12:00 + 0064  39 V VEL    100  550    0  2000-01-24_12:00 + 0065  39 V VEL    100  500    0  2000-01-24_12:00 + 0066  39 V VEL    100  450    0  2000-01-24_12:00 + 0067  39 V VEL    100  400    0  2000-01-24_12:00 + 0068  39 V VEL    100  350    0  2000-01-24_12:00 + 0069  39 V VEL    100  300    0  2000-01-24_12:00 + 0070  39 V VEL    100  250    0  2000-01-24_12:00 + 0071  39 V VEL    100  200    0  2000-01-24_12:00 + 0072  39 V VEL    100  150    0  2000-01-24_12:00 + 0073  39 V VEL    100  100    0  2000-01-24_12:00 + 0074  52 R H      100 1000    0  2000-01-24_12:00 + 0075  52 R H      100  975    0  2000-01-24_12:00 + 0076  52 R H      100  950    0  2000-01-24_12:00 + 0077  52 R H      100  925    0  2000-01-24_12:00 + 0078  52 R H      100  900    0  2000-01-24_12:00 + 0079  52 R H      100  850    0  2000-01-24_12:00 + 0080  52 R H      100  800    0  2000-01-24_12:00 + 0081  52 R H      100  750    0  2000-01-24_12:00 + 0082  52 R H      100  700    0  2000-01-24_12:00 + 0083  52 R H      100  650    0  2000-01-24_12:00 + 0084  52 R H      100  600    0  2000-01-24_12:00 + 0085  52 R H      100  550    0  2000-01-24_12:00 + 0086  52 R H      100  500    0  2000-01-24_12:00 + 0087  52 R H      100  450    0  2000-01-24_12:00 + 0088  52 R H      100  400    0  2000-01-24_12:00 + 0089  52 R H      100  350    0  2000-01-24_12:00 + 0090  52 R H      100  300    0  2000-01-24_12:00 + 0091  52 R H      100  250    0  2000-01-24_12:00 + 0092  52 R H      100  200    0  2000-01-24_12:00 + 0093  52 R H      100  150    0  2000-01-24_12:00 + 0094  52 R H      100  100    0  2000-01-24_12:00 + 0095  41 ABS V    100 1000    0  2000-01-24_12:00 + 0096  41 ABS V    100  975    0  2000-01-24_12:00 + 0097  41 ABS V    100  950    0  2000-01-24_12:00 + 0098  41 ABS V    100  925    0  2000-01-24_12:00 + 0099  41 ABS V    100  900    0  2000-01-24_12:00 + 00100  41 ABS V    100  850    0  2000-01-24_12:00 + 00101  41 ABS V    100  800    0  2000-01-24_12:00 + 00102  41 ABS V    100  750    0  2000-01-24_12:00 + 00103  41 ABS V    100  700    0  2000-01-24_12:00 + 00104  41 ABS V    100  650    0  2000-01-24_12:00 + 00105  41 ABS V    100  600    0  2000-01-24_12:00 + 00106  41 ABS V    100  550    0  2000-01-24_12:00 + 00107  41 ABS V    100  500    0  2000-01-24_12:00 + 00108  41 ABS V    100  450    0  2000-01-24_12:00 + 00109  41 ABS V    100  400    0  2000-01-24_12:00 + 00110  41 ABS V    100  350    0  2000-01-24_12:00 + 00111  41 ABS V    100  300    0  2000-01-24_12:00 + 00112  41 ABS V    100  250    0  2000-01-24_12:00 + 00113  41 ABS V    100  200    0  2000-01-24_12:00 + 00114  41 ABS V    100  150    0  2000-01-24_12:00 + 00115  41 ABS V    100  100    0  2000-01-24_12:00 + 00116  41 ABS V    100   70    0  2000-01-24_12:00 + 00117  41 ABS V    100   50    0  2000-01-24_12:00 + 00118  41 ABS V    100   30    0  2000-01-24_12:00 + 00119  41 ABS V    100   20    0  2000-01-24_12:00 + 00120  41 ABS V    100   10    0  2000-01-24_12:00 + 00121 154 O3MR     100  100    0  2000-01-24_12:00 + 00122 154 O3MR     100   70    0  2000-01-24_12:00 + 00123 154 O3MR     100   50    0  2000-01-24_12:00 + 00124 154 O3MR     100   30    0  2000-01-24_12:00 + 00125 154 O3MR     100   20    0  2000-01-24_12:00 + 00126 154 O3MR     100   10    0  2000-01-24_12:00 + 00127 222 5WAVH    100  500    0  2000-01-24_12:00 + 00128  33 U GRD    100 1000    0  2000-01-24_12:00 + 00129  33 U GRD    100  975    0  2000-01-24_12:00 + 00130  33 U GRD    100  950    0  2000-01-24_12:00 + 00131  33 U GRD    100  925    0  2000-01-24_12:00 + 00132  33 U GRD    100  900    0  2000-01-24_12:00 + 00133  33 U GRD    100  850    0  2000-01-24_12:00 + 00134  33 U GRD    100  800    0  2000-01-24_12:00 + 00135  33 U GRD    100  750    0  2000-01-24_12:00 + 00136  33 U GRD    100  700    0  2000-01-24_12:00 + 00137  33 U GRD    100  650    0  2000-01-24_12:00 + 00138  33 U GRD    100  600    0  2000-01-24_12:00 + 00139  33 U GRD    100  550    0  2000-01-24_12:00 + 00140  33 U GRD    100  500    0  2000-01-24_12:00 + 00141  33 U GRD    100  450    0  2000-01-24_12:00 + 00142  33 U GRD    100  400    0  2000-01-24_12:00 + 00143  33 U GRD    100  350    0  2000-01-24_12:00 + 00144  33 U GRD    100  300    0  2000-01-24_12:00 + 00145  33 U GRD    100  250    0  2000-01-24_12:00 + 00146  33 U GRD    100  200    0  2000-01-24_12:00 + 00147  33 U GRD    100  150    0  2000-01-24_12:00 + 00148  33 U GRD    100  100    0  2000-01-24_12:00 + 00149  33 U GRD    100   70    0  2000-01-24_12:00 + 00150  33 U GRD    100   50    0  2000-01-24_12:00 + 00151  33 U GRD    100   30    0  2000-01-24_12:00 + 00152  33 U GRD    100   20    0  2000-01-24_12:00 + 00153  33 U GRD    100   10    0  2000-01-24_12:00 + 00154  34 V GRD    100 1000    0  2000-01-24_12:00 + 00155  34 V GRD    100  975    0  2000-01-24_12:00 + 00156  34 V GRD    100  950    0  2000-01-24_12:00 + 00157  34 V GRD    100  925    0  2000-01-24_12:00 + 00158  34 V GRD    100  900    0  2000-01-24_12:00 + 00159  34 V GRD    100  850    0  2000-01-24_12:00 + 00160  34 V GRD    100  800    0  2000-01-24_12:00 + 00161  34 V GRD    100  750    0  2000-01-24_12:00 + 00162  34 V GRD    100  700    0  2000-01-24_12:00 + 00163  34 V GRD    100  650    0  2000-01-24_12:00 + 00164  34 V GRD    100  600    0  2000-01-24_12:00 + 00165  34 V GRD    100  550    0  2000-01-24_12:00 + 00166  34 V GRD    100  500    0  2000-01-24_12:00 + 00167  34 V GRD    100  450    0  2000-01-24_12:00 + 00168  34 V GRD    100  400    0  2000-01-24_12:00 + 00169  34 V GRD    100  350    0  2000-01-24_12:00 + 00170  34 V GRD    100  300    0  2000-01-24_12:00 + 00171  34 V GRD    100  250    0  2000-01-24_12:00 + 00172  34 V GRD    100  200    0  2000-01-24_12:00 + 00173  34 V GRD    100  150    0  2000-01-24_12:00 + 00174  34 V GRD    100  100    0  2000-01-24_12:00 + 00175  34 V GRD    100   70    0  2000-01-24_12:00 + 00176  34 V GRD    100   50    0  2000-01-24_12:00 + 00177  34 V GRD    100   30    0  2000-01-24_12:00 + 00178  34 V GRD    100   20    0  2000-01-24_12:00 + 00179  34 V GRD    100   10    0  2000-01-24_12:00 + 00180  11 TMP      116   30    0  2000-01-24_12:00 + 00181  52 R H      116   30    0  2000-01-24_12:00 + 00182  51 SPF H    116   30    0  2000-01-24_12:00 + 00183  33 U GRD    116   30    0  2000-01-24_12:00 + 00184  34 V GRD    116   30    0  2000-01-24_12:00 + 00185  11 TMP      103 1829    0  2000-01-24_12:00 + 00186  11 TMP      103 2743    0  2000-01-24_12:00 + 00187  11 TMP      103 3658    0  2000-01-24_12:00 + 00188  33 U GRD    103 1829    0  2000-01-24_12:00 + 00189  33 U GRD    103 2743    0  2000-01-24_12:00 + 00190  33 U GRD    103 3658    0  2000-01-24_12:00 + 00191  34 V GRD    103 1829    0  2000-01-24_12:00 + 00192  34 V GRD    103 2743    0  2000-01-24_12:00 + 00193  34 V GRD    103 3658    0  2000-01-24_12:00 + 00194   1 PRES       1    0    0  2000-01-24_12:00 + 00195  54 P WAT    200    0    0  2000-01-24_12:00 + 00196  52 R H      200    0    0  2000-01-24_12:00 + 00197   7 HGT        7    0    0  2000-01-24_12:00 + 00198  11 TMP        7    0    0  2000-01-24_12:00 + 00199   1 PRES       7    0    0  2000-01-24_12:00 + 00200 136 VW SH      7    0    0  2000-01-24_12:00 + 00201 131 LFT X      1    0    0  2000-01-24_12:00 + 00202 157 CAPE       1    0    0  2000-01-24_12:00 + 00203 156 CIN        1    0    0  2000-01-24_12:00 + 00204 132 4LFTX      1    0    0  2000-01-24_12:00 + 00205 157 CAPE     116  180    0  2000-01-24_12:00 + 00206 156 CIN      116  180    0  2000-01-24_12:00 + 00207   7 HGT        6    0    0  2000-01-24_12:00 + 00208  11 TMP        6    0    0  2000-01-24_12:00 + 00209   1 PRES       6    0    0  2000-01-24_12:00 + 00210   7 HGT        1    0    0  2000-01-24_12:00 + 00211   2 PRMSL    102    0    0  2000-01-24_12:00 + 00212  52 R H      108   44  100  2000-01-24_12:00 + 00213  52 R H      108   72   94  2000-01-24_12:00 + 00214  52 R H      108   44   72  2000-01-24_12:00 + 00215  52 R H      108   33  100  2000-01-24_12:00 + 00216  13 POT      107 9950    0  2000-01-24_12:00 + 00217  11 TMP      107 9950    0  2000-01-24_12:00 + 00218  39 V VEL    107 9950    0  2000-01-24_12:00 + 00219  52 R H      107 9950    0  2000-01-24_12:00 + 00220  10 TOZNE    200    0    0  2000-01-24_12:00 + 00221  33 U GRD      7    0    0  2000-01-24_12:00 + 00222  33 U GRD      6    0    0  2000-01-24_12:00 + 00223  33 U GRD    107 9950    0  2000-01-24_12:00 + 00224  34 V GRD      7    0    0  2000-01-24_12:00 + 00225  34 V GRD      6    0    0  2000-01-24_12:00 + 00226  34 V GRD    107 9950    0  2000-01-24_12:00 + 00227  11 TMP        1    0    0  2000-01-24_12:00 + 00228 144 SOILW    112    0   10  2000-01-24_12:00 + 00229 144 SOILW    112   10  200  2000-01-24_12:00 + 00230  11 TMP      112    0   10  2000-01-24_12:00 + 00231  11 TMP      112   10  200  2000-01-24_12:00 + 00232  65 WEASD      1    0    0  2000-01-24_12:00 + 00233  81 LAND       1    0    0  2000-01-24_12:00 + 00234  91 ICE C      1    0    0  2000-01-24_12:00 + 00235  11 TMP      105    2    0  2000-01-24_12:00 + 00236  51 SPF H    105    2    0  2000-01-24_12:00 + 00237 221 HPBL       1    0    0  2000-01-24_12:00 + 00238  71 T CDC    244    0    0  2000-01-24_12:00 + 00239   1 PRES     243    0    0  2000-01-24_12:00 + 00240   1 PRES     242    0    0  2000-01-24_12:00 + 00241  52 R H      105    2    0  2000-01-24_12:00 + 00242  33 U GRD    105   10    0  2000-01-24_12:00 + 00243  34 V GRD    105   10    0  2000-01-24_12:00 + 00244  27 GP A     100 1000    0  2000-01-24_12:00 + 00245  27 GP A     100  500    0  2000-01-24_12:00 + 00246 230 5WAVA    100  500    0  2000-01-24_12:00 + 00***** End-Of-File on C unit   3

链接Variable_Table文件

ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable

通过脚本./link_grib.csh链接GRIB数据

./link_grib.csh ../DATA/JAN00/fnl_2000012

Note: For explanation of these variables see Chapter 3 of the WRF User’s Guide, or our Best Practice WPS Namelist web page for detailed descriptions of the namelist variables, and suggestions for best practices.

修改namelist.wps静态地理数据的开始和结束时间

终端在WPS目录下输入：

nano namelist.wps

将namelist.wps中的时间：

start_date = '2006-08-16_12:00:00','2006-08-16_12:00:00',end_date   = '2006-08-16_18:00:00','2006-08-16_12:00:00',

更改为以下：

**start_date = '2000-01-24_12:00:00',
end_date = '2000-01-25_12:00:00',**

start_date = '2000-01-24_12:00:00','2000-01-24_12:00:00',
end_date   = '2000-01-25_12:00:00','2000-01-25_12:00:00',

另外注意这两项对不对

**interval_seconds = 21600,prefix = 'FILE',**

运行ungrib.exe

./ungrib.exe

(这里有可能报错./ungrib.exe: error while loading shared libraries: libpng12.so.0: cannot open shared object file: No such file or directory 其实就是你的libpng没有安装到它默认的路径下。解决办法如下：找到你的libpng12.so.0的安装路径，一般实在一个/lib文件夹下面。终端在WPS目录下输入以下命令：

vim /etc/ld.so.conf
#按i进入编辑模式将libpng12.so.0的安装路径输入进去.
这是我自己的安装路径 eg:/home/xiaomo/Build_WRF/LIBRAREIES/grib2/lib 然后按住Esc再按shift+zz保存退出

接着输入命令

ldconfig

这样这个error就解决了。另一个error好象是DATA NOT FOUND，我没有遇到，可以去气象家园看看相关解决办法的帖子。

再运行

./ungrib.exe

结果如下

./ungrib.exe运行结果截取部分如下：

**********
Deleting temporary files created by ungrib...
**********Deleting file:  ./PFILE:2000-01-24_12
Deleting file:  ./PFILE:2000-01-24_18
Deleting file:  ./PFILE:2000-01-25_00
Deleting file:  ./PFILE:2000-01-25_06
Deleting file:  ./PFILE:2000-01-25_12**********
Done deleting temporary files.
**********!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of ungrib.   !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

这将创建以下文件：

FILE:2000-01-24_12

FILE:2000-01-24_18

FILE:2000-01-25_00

FILE:2000-01-25_06

FILE:2000-01-25_12

同时解码成功后会在WPS下面生成以GRIBFILE.为开头的文件。eg:GRIBFILE.AAA GRIBFILE.AAB GRIBFILE.AAC

运行此操作时，将自动创建一个ungrib.log文件。它包含有关输入文件中的字段以及这些字段在哪些级别可用的信息。如果此步骤中出现问题，它也有助于进行故障排除。细读此日志以确保找到 Vtable 中的所有字段。

2 运行./metgrid.exe

将输入数据插值到我们的模型域(metgrid.exe)

./metgrid.exe

The text below should appear on your screen during run time. Look for "Successful completion of metgrid " at the end of the run, indicating that all went well during the execution. A metgrid.log file containing more information will also be created, and can be useful for troubleshooting purposes if the run is not successful.

下面的文本应该在运行时出现在屏幕上。在运行结束时查看“成功完成metgrid”，表明在执行期间一切都进展顺利。还将创建一个包含更多信息的metgrid.log文件，如果运行不成功，这个文件对于排除故障非常有用。

xiaomo@ubuntu:~/Build_WRF/WPS$ ./metgrid.exe
Processing domain 1 of 1Processing 2000-01-24_12FILEProcessing 2000-01-24_18FILEProcessing 2000-01-25_00FILEProcessing 2000-01-25_06FILEProcessing 2000-01-25_12FILE
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
!  Successful completion of metgrid.  !
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
Note: The following floating-point exceptions are signalling: IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL

此步骤将创建以下文件

met_em.d01.2000-01-24_12:00:00.nc

met_em.d01.2000-01-24_18:00:00.nc

met_em.d01.2000-01-25_00:00:00.nc

met_em.d01.2000-01-25_06:00:00.nc

met_em.d01.2000-01-25_12:00:00.nc

（至此，前处理完毕）

3 运行WRF

进入WRF模式的主体，终端进入WRF目录下

cd test/em_real

将上一步生成的met_em.文件链接到WRF运行的路径下

ln -sf ../../../WPS/met_em.d01.2000-01* .

修改namelist.input（参照官网参数修改wrf.php），可以在文件夹下直接打开namelist.input修改，也可以用vi命令来编辑（namelist.input的存放路径在/WRF/test/em_real下）

Note 1: You should not need to make any changes, as these are the default settings, but check to ensure they are set as below

Note 2: As we are only interested in domain 1, you can ignore all values in other columns, and the model will do the same.

Note3: For explanation of these variables, in addition to some suggestions for best practices, see our Best Practice WRF Namelist page or Chapter 5 of the WRF User’s Guide.
run_days = 0,

run_hours = 12,

run_minutes = 0,

run_seconds = 0,

start_year = 2000,

start_month = 01,

start_day = 24,

start_hour = 12,

end_year = 2000,

end_month = 01,

end_day = 25,

end_hour = 12,

interval_seconds = 21600

input_from_file = .true.,

history_interval = 180,

frames_per_outfile = 1000,

time_step = 180,

max_dom = 1,

s_we = 1,

e_we = 74,

s_sn = 1,

e_sn = 61,

s_vert = 1,

e_vert = 33,

num_metgrid_levels = 27

dx = 30000,

dy = 30000,

#执行real.exe,初始化模式

./real.exe

#检查有没有出错

tail rsl.error.0000

(根据具体error修改namelist.input参数,然后再执行./real.exe。再检查有没有出错）

出错了：

xiaomo@ubuntu:~/Build_WRF/WRF/test/em_real$ tail rsl.error.0000metgrid input_wrf.F first_date_nml = 2000-01-24_12:00:00
d01 2000-01-24_12:00:00 Timing for input          0 s.
d01 2000-01-24_12:00:00          flag_soil_layers read from met_em file is  1
Max map factor in domain 1 =  1.03. Scale the dt in the model accordingly.
ERROR in psfc: flag_psfc =  1, flag_soilhgt =  1, flag_slp =  0, sfcp_to_sfcp = F
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:    1396
not enough info for a p sfc computation
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0

可见ERROR in psfc: flag_psfc = 1, flag_soilhgt = 1, flag_slp = 0, sfcp_to_sfcp = F

解决方法：
在namelist.input文件的“domains”部分添加下面这一条即可
sfcp_to_sfcp = .true.
再次运行，成功。

（运行成功后会生成wrfbdy_d01 ， wrfinput_d01两个文件）

运行wrf.exe

./wrf.exe

(完了后检查rsl.out.0000文件）

运行成功后会生成wrfout_d01_2000-01-24_12:00:00文件

Quick check of the contents of this file:

a. Use the ncdump utility:

ncdump -h wrfout_d01_2000-01-24_12:00:00
ncdump -v Times wrfout_d01_2000-01-24_12:00:00
(to see which forecast times are in the file)

b. Use the ncview utility.

下次写用WRF Domain Wizard来运行官网案例O(∩_∩)O