WRF | WRF台风个例

一、台风介绍

2011年第17号强台风“纳沙”（英文名：NESAT，名字来源：柬埔寨，名字意义：渔夫）于2011年9月24日上午在西北太平洋洋面上生成，9月26日夜间和29日7时两度加强成为强台风，是2011年以来登陆中国强度最强的台风（以强台风级别登陆海南），也是6年来登陆海南的最强台风。

二、模拟范围

经度：105°E-125°E
纬度：12°N-24°N
时间：2011-09-28_06至2011-09-28_12
网格分辨率：9000 x 9000
本次模拟不涉及嵌套，作为示例时间只模拟6小时。

三、文件准备

1 文件夹准备

2 FNL文件

位置：/FNL

3 namelist.wps 文件

位置：/WPS

 &share
 wrf_core = 'ARW',
 max_dom = 1,
 start_date = 2011-09-28_06:00:00,'2006-08-16_12:00:00',
 end_date   = 2011-09-28_12:00:00,'2006-08-16_12:00:00',
 interval_seconds = 21600
 io_form_geogrid = 2,
/

&geogrid
 parent_id         =   1,   1,
 parent_grid_ratio =   1,   3,
 i_parent_start    =   1,  31,
 j_parent_start    =   1,  17,
 e_we              =  256, 112,
 e_sn              =  158,  97,
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! The default datasets used to produce the MAXSNOALB and ALBEDO12M
 ! fields have changed in WPS v4.0. These fields are now interpolated
 ! from MODIS-based datasets.
 !
 ! To match the output given by the default namelist.wps in WPS v3.9.1,
 ! the following setting for geog_data_res may be used:
 !
 ! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default', 
 !
 !!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
 !
 geog_data_res = '10m','2m',
 dx = 9000,
 dy = 9000,
 map_proj = 'lambert',
 ref_lat   = 18,
 ref_lon   = 115,
 truelat1  = 12,
 truelat2  = 24,
 stand_lon = 115,
 geog_data_path = '/home/daimu/Desktop/WRF/geog/geog'
/

&ungrib
 out_format = 'WPS',
 prefix = 'FILE',
/

&metgrid
 fg_name = 'FILE'
 io_form_metgrid = 2, 
/

4 namelist.input 文件

位置：/WRF/test/em_real

 &time_control
 run_days                            = 0,
 run_hours                           = 3,
 run_minutes                         = 0,
 run_seconds                         = 0,
 start_year                          = 2011, 2000, 2000,
 start_month                         = 09,   01,   01,
 start_day                           = 28,   24,   24,
 start_hour                          = 06,   12,   12,
 end_year                            = 2011, 2000, 2000,
 end_month                           = 09,   01,   01,
 end_day                             = 28,   25,   25,
 end_hour                            = 12,   12,   12,
 interval_seconds                    = 21600
 input_from_file                     = .true.,.true.,.true.,
 history_interval                    = 180,  60,   60,
 frames_per_outfile                  = 1, 1000, 1000,
 restart                             = .false.,
 restart_interval                    = 180,
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2
 /

 &domains
 time_step                           = 30,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 256,    112,   94,
 e_sn                                = 158,    97,    91,
 e_vert                              = 33,    33,    33,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 27,
 num_metgrid_soil_levels             = 4,
 dx                                  = 9000, 10000,  3333.33,
 dy                                  = 9000, 10000,  3333.33,
 grid_id                             = 1,     2,     3,
 parent_id                           = 0,     1,     2,
 i_parent_start                      = 1,     31,    30,
 j_parent_start                      = 1,     17,    30,
 parent_grid_ratio                   = 1,     3,     3,
 parent_time_step_ratio              = 1,     3,     3,
 feedback                            = 1,
 smooth_option                       = 0
 /

 &physics
 mp_physics                          = 3,     3,     3,
 ra_lw_physics                       = 1,     1,     1,
 ra_sw_physics                       = 1,     1,     1,
 radt                                = 30,    30,    30,
 sf_sfclay_physics                   = 1,     1,     1,
 sf_surface_physics                  = 2,     2,     2,
 bl_pbl_physics                      = 1,     1,     1,
 bldt                                = 0,     0,     0,
 cu_physics                          = 1,     1,     0,
 cudt                                = 5,     5,     5,
 isfflx                              = 1,
 ifsnow                              = 1,
 icloud                              = 1,
 surface_input_source                = 1,
 num_soil_layers                     = 4,
 sf_urban_physics                    = 0,     0,     0,
 /

 &fdda
 /

 &dynamics
 hybrid_opt                          = 2, 
 w_damping                           = 0,
 diff_opt                            = 1,      1,      1,
 km_opt                              = 4,      4,      4,
 diff_6th_opt                        = 0,      0,      0,
 diff_6th_factor                     = 0.12,   0.12,   0.12,
 base_temp                           = 290.
 damp_opt                            = 3,
 zdamp                               = 5000.,  5000.,  5000.,
 dampcoef                            = 0.2,    0.2,    0.2
 khdif                               = 0,      0,      0,
 kvdif                               = 0,      0,      0,
 non_hydrostatic                     = .true., .true., .true.,
 moist_adv_opt                       = 1,      1,      1,     
 scalar_adv_opt                      = 1,      1,      1,     
 gwd_opt                             = 1,
 /

 &bdy_control
 spec_bdy_width                      = 5,
 specified                           = .true.
 /

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,
 /

四 运行流程

1 WPS文件夹

step 1 : ./geogrid.exe

终端：

WPS文件夹-增加的文件：

step 2 : ./link_grib.csh …/FNL/fnl*

终端：

WPS文件夹-增加的文件：

step 3 : ./ungrib.exe

终端：

WPS文件夹-增加的文件：

step 4 : ./metgrid.exe

终端：

WPS文件夹-增加的文件：

2 WRF-test-em_real文件夹

step 1 : ln -sf …/…/…/WPS/met_em.d* .

终端：

WRF文件夹-增加的文件：

step 2 : ./real.exe

终端：

WRF文件夹-增加的文件：

step 3 : ./wrf.exe

终端：

WRF文件夹-增加的文件：

注意：如有问题，请查阅（用几个核就会出现几个文件）：