小曾WRF自学日记(4)渐入佳境 ——WRF实例-运行WRF

WRF的运行

进入WRF/run文件夹，修改namelist.input

首先在WPS文件夹下，使用命令ncdump -h met_em.d01.2018-05-04_12:00:00.nc 查看最后的 global attributes部分。

// global attributes:
                :TITLE = "OUTPUT FROM METGRID V4.0" ;
                :SIMULATION_START_DATE = "2018-05-04_12:00:00" ;
                :WEST-EAST_GRID_DIMENSION = 150 ;
                :SOUTH-NORTH_GRID_DIMENSION = 150 ;
                :BOTTOM-TOP_GRID_DIMENSION = 38 ;
                :WEST-EAST_PATCH_START_UNSTAG = 1 ;
                :WEST-EAST_PATCH_END_UNSTAG = 149 ;
                :WEST-EAST_PATCH_START_STAG = 1 ;
                :WEST-EAST_PATCH_END_STAG = 150 ;
                :SOUTH-NORTH_PATCH_START_UNSTAG = 1 ;
                :SOUTH-NORTH_PATCH_END_UNSTAG = 149 ;
                :SOUTH-NORTH_PATCH_START_STAG = 1 ;
                :SOUTH-NORTH_PATCH_END_STAG = 150 ;
                :GRIDTYPE = "C" ;
                :DX = 45000.f ;
                :DY = 45000.f ;
                :DYN_OPT = 2 ;
                :CEN_LAT = 89.99999f ;
                :CEN_LON = 180.f ;
                :TRUELAT1 = 60.f ;
                :TRUELAT2 = 60.f ;
                :MOAD_CEN_LAT = 89.99999f ;
                :STAND_LON = 0.f ;
                :POLE_LAT = 90.f ;
                :POLE_LON = 0.f ;
                :corner_lats = 46.77416f, 46.77416f, 46.77416f, 46.77416f, 46.64143f, 46.64143f, 46.64143f, 46.64143f, 46.64143f, 46.64143f, 46.64143f, 46.64143f, 46.50927f, 46.50926f, 46.50926f, 46.50927f ;
                :corner_lons = -45.f, -135.f, 135.f, 45.f, -45.1929f, -134.8071f, 134.8071f, 45.19292f, -44.80708f, -135.1929f, 135.1929f, 44.80709f, -45.f, -135.f, 135.f, 45.f ;
                :MAP_PROJ = 2 ;
                :MMINLU = "MODIFIED_IGBP_MODIS_NOAH" ;
                :NUM_LAND_CAT = 21 ;
                :ISWATER = 17 ;
                :ISLAKE = 21 ;
                :ISICE = 15 ;
                :ISURBAN = 13 ;
                :ISOILWATER = 14 ;
                :grid_id = 1 ;
                :parent_id = 1 ;
                :i_parent_start = 1 ;
                :j_parent_start = 1 ;
                :i_parent_end = 150 ;
                :j_parent_end = 150 ;
                :parent_grid_ratio = 1 ;
                :sr_x = 1 ;
                :sr_y = 1 ;
                :NUM_METGRID_SOIL_LEVELS = 4 ;
                :FLAG_METGRID = 1 ;
                :FLAG_EXCLUDED_MIDDLE = 0 ;
                :FLAG_SOIL_LAYERS = 1 ;
                :FLAG_SNOW = 1 ;
                :FLAG_PSFC = 1 ;
                :FLAG_SM000007 = 1 ;
                :FLAG_SM007028 = 1 ;
                :FLAG_SM028100 = 1 ;
                :FLAG_SM100289 = 1 ;
                :FLAG_ST000007 = 1 ;
                :FLAG_ST007028 = 1 ;
                :FLAG_ST028100 = 1 ;
                :FLAG_ST100289 = 1 ;
                :FLAG_SLP = 1 ;
                :FLAG_SST = 1 ;
                :FLAG_MF_XY = 1 ;
                :FLAG_LAI12M = 1 ;
                :FLAG_VAR_SSO = 1 ;
}

根据内容修改namelist.input

time control部分

&time_control
 run_days                            = 5, 
 run_hours                           = 0,
 run_minutes                         = 0,
 run_seconds                         = 0,    %模型运行时间
 start_year                          = 2018,   
 start_month                         = 05,
 start_day                           = 01,
 start_hour                          = 12,
 end_year                            = 2018,
 end_month                           = 05,
 end_day                             = 06,
 end_hour                            = 12,   % 起始/结束时间
 interval_seconds                    = 21600  %输入文件的间隔时间，单位秒，此处为间隔6小时
 input_from_file                     = .true.,.true., %是否嵌套运行
 history_interval                    = 60,  60,    %输出文件间隔，单位为分，此处为每小时输出一个文件
 frames_per_outfile                  = 1, 1, %输出的out文件个数，如果out太大，可以切割输出
 restart                             = .false., %是否是重新运行
 restart_interval                    = 7200, %重新运行out文件的时间间隔，单位分钟
 % 以下全为输出文件、输入文件格式，全选2，表示格式为nc
 io_form_history                     = 2
 io_form_restart                     = 2
 io_form_input                       = 2
 io_form_boundary                    = 2

domains部分

与WPS的namelist和met_dem文件有关

&domain
 time_step                           = 90,
 time_step_fract_num                 = 0,
 time_step_fract_den                 = 1,
 max_dom                             = 1,
 e_we                                = 150,
 e_sn                                = 150,
 e_vert                              = 40,
 dzbot                               = 10,
 dzstretch_s                         = 1.3,
 dzstretch_u                         = 1.1,
 p_top_requested                     = 5000,
 num_metgrid_levels                  = 38,
 num_metgrid_soil_levels             = 4,
 dx                                  = 30000,
 dy                                  = 30000,
 grid_id                             = 1,
 parent_id                           = 1,
 i_parent_start                      = 1,
 j_parent_start                      = 1,
 parent_grid_ratio                   = 1,
 parent_time_step_ratio              = 1,
 feedback                            = 1,
 smooth_option                       = 0
 /
%全部与namelist中保持一致！

physics 部分

物理参数化方案的选择，注意查看WRF文档里，参数化方案的一些要求，

&physics
 physics_suite                       = 'CONUS'
 mp_physics                          = 2, %微物理方案，此处选择Lin
 cu_physics                          = 1, %积云参数化方案，此处选择Kain-Fritsch (new Eta) 方案
 ra_lw_physics                       = 1  %长波辐射方案，此处选择rrtn
 ra_sw_physics                       = 5, %短波辐射方案，此处选择Goddard
 bl_pbl_physics                      = 1, %边界层方案，此处选择YSU
 sf_sfclay_physics                   = 1, %地面层方案，应与边界层方案对应，此处为MM5
 sf_surface_physics                  = 2, %地面方案，这里选的Noah
 radt                                = 30, %与你的分辨率对应，dx=20km就写30
 bldt                                = 0, %推荐0
 cudt                                = 0, % 推荐0，除了  Kain-Fritsch参数化方案，这里我忘记改了
 icloud                              = 1,% 是否考虑云辐射效应，这里先不管，直接默认
 num_land_cat                        = 21,%与met_em文件对应
 sf_urban_physics                    = 2,%城市冠层方案，注意与边界层方案对应
 fractional_seaice                   = 1, %把海冰当作分数考虑，如果为0，则要么有冰有么没冰
 /

dynamic部分与其他

注意垂直层可能要加个etac上去

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

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

 &grib2
 /

 &namelist_quilt
 nio_tasks_per_group = 0,
 nio_groups = 1,