WRF实例运行系列（1）

## WRF实例运行1

****单层嵌套模型****

step1：定义模型模拟域

• 下载静态地理数据并解压（以本人计算机实践为例）

下载地址WRFV4：https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog.html#specifi

下载地址WRFV3：

https://www2.mmm.ucar.edu/wrf/users/download/get_sources_wps_geog_V3.html

解压后存储地址：$HOME/Public/Build_WRF/WPS_GEOG

可以不用全部下载，运行geogrid时，程序告诉你缺什么就可以下载什么。 

• 修改namelist.wps内容

cd ~/Public/Build_WRF/WRFV4/WPS
vi namelist.wps

修改namelist.wps中的&geogrid和&share中与时间无关的字段，均可按照自己想设置的域范围进行修改。

&share

 wrf_core = 'ARW',
 max_dom = 1, #单层粗域处理

&geogrid  #定义模拟域投影及所有模型网格尺寸和位置。具体每项的含义可以参考前文“the WRF preprocessing system(WPS)介绍”

 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  31,  
 j_parent_start    =   1,  17,  
 e_we              =  74, 112,  
 e_sn              =  61,  97,  

geog_data_res = 'default','default',
 dx = 30000,
 dy = 30000,
 map_proj = 'lambert',
 ref_lat   =  34.83,
 ref_lon   =  -81.03,
 truelat1  =  30.0,
 truelat2  =  60.0,
 stand_lon = -98.0,
 geog_data_path = '/home/hsy/Public/Build_WRF/WPS_GEOG/'
最主要的是一定要输入WPS_GEOG的绝对路径，注意每行后面都有','作为结尾。

• 运行geogrid

./geogrid.exe

运行成功后。 

nc文件可以使用ncview查看。

 ncview geogrid_output/20000124/geo_em.d01.nc 

 

step2：用ungrib从网格数据中提取气象字段

• 连接气象数据，使用./link_grib.csh脚本连接到自己解压数据后的存储地址。（本文使用的是20000124-20000125的gfs数据，可自行下载，也可以去博主资源主页https://download.csdn.net/download/A18040554844/18937434下载。）

./link_grib.csh ../../../grid_data/JAN00_GFS_FNL/JAN00/

• 连接Vtable。（本文使用的是Vtable.GFS，可根据下载的数据选择不同的Vtable，更可以根据自己下载的数据源创建新的Vtable。）

每个Vtable包含7个或11个字段，具体取决于Vtable是用于GRIB Edition 1数据源还是用于GRIB Edition 2数据源。Vtable的字段可以分为三类:描述如何在GRIB文件中标识数据的字段，描述ungrib和metgrid程序如何标识数据的字段，以及特定于GRIB Edition 2的字段。具体如何创建Vtable可以自行根据官网教程设置。

ln -sf ungrib/Variable_Tables/Vtable.GFS Vtable 

• 修改namelist.wps中与时间相关的share字段以及ungrib字段。

&share#修改与模拟时间相关的部分

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

&ungrib
 out_format = 'WPS',
 prefix = 'ungrib_output/20000124/FILE',#输出的临时文件位置为ungrib_output/20000124/，FILE为文件前缀。

• 运行ungrib。

./ungrib.exe

运行成功后。 

step3：用metgrid水平插值气象数据

• 修改namelist.wps中&metgrid部分，这里可以再次检查一下&share中时间设置是否恰当。

&metgrid
 fg_name = 'ungrib_output/20000124/FILE',
 io_form_metgrid = 2,
 opt_output_from_metgrid_path = 'metgrid_output/'        

• 运行metgrid。

./metgrid.exe

 运行成功后。

 附单层粗域namelist.wps设置以供参考(WRF4.0)。

step4：运行real.exe。

• 连接metgrid产生的nc数据。

cd ../WRF/test/em_real/
ls -l ../../../WPS/metgrid_output/met_em.d01.2000-01-2* 
ln -sf ../../../WPS/metgrid_output/met_em.d01.2000-01-2* .

连接成功后。

• 修改namelist.input内容

主要修改&time_control和&domains部分。

 &time_control
 run_days                            = 0,
 run_hours                           = 24,
 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.,.true.,.true.,
 history_interval                    = 180,  60,  60,
 frames_per_outfile                  = 1000, 1000, 1000,
 restart                             = .false.,
 restart_interval                    = 5000,
 io_form_auxinput2                   = 2,
 io_form_history                     = 2,
 io_form_restart                     = 2,
 io_form_input                       = 2,
 io_form_boundary                    = 2,
 /

 &domains
 time_step                           = 180,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 74,   ,   ,
 e_sn                                = 61,   ,   ,
 e_vert                              = 30,    30,    30,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 30000, 10000,  3333.33,
 dy                                  = 30000, 10000,  3333.33,
 grid_id                             = 1,
 parent_id                           = 0,
 i_parent_start                      = 1,    ,    ,
 j_parent_start                      = 1,    ,    ,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0

• 运行real

./real.exe
#运行完以后，立马检查real输出。
tail rsl.error.0000
#or
tail rsl.out.0000
#or 
vi rsl.error.0000

real运行报错，证明namelist.input设置与气象数据输入不同，请耐心根据报错内容修改namelist.input中&domain字段，直至出现SUCCESS COMPLETE REAL_EM INIT以及
 wrfbdy_d01和 wrfinput_d01两个文件。

step5：运行wrf模型。

• 运行wrf

./wrf.exe

运行成功后，会出现SUCCESS COMPLETE WRF提示，并且出现wrfout_d0*文件。