小曾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,