wrf跑2017全年，请问namelist.input这么设置合理吗？

physzw

3层嵌套，测试跑一天花费的时间约15分钟，每个wrfout文件的大小约为1.2G。
如果跑一整年，namelist.input这么设置合理吗？
尤其是frames_per_outfile 和 restart_interval应该如何设置？

&time_control
run_days                            = 365,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2017, 2017, 2017,
start_month                         = 01,   01,   01,
start_day                           = 01,   01,   01,
start_hour                          = 00,   00,   00,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2018, 2018, 2018,
end_month                           = 01,   01,   01,
end_day                             = 01,   01,   01,
end_hour                            = 00,   00,   00,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval      =  60,  60,   60,
frames_per_outfile   = 1000, 1000, 1000,
restart              = .false.,
restart_interval      = 45000,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                 = 216,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 3,
e_we                                = 121, 157, 166,
e_sn                                = 121, 163, 154,
e_vert                              = 30,    30,    30,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 36000, 12000,  4000,
dy                                  = 36000, 12000,  4000,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,  35, 51,
j_parent_start                      = 1,  35, 50,
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                = 3,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,     0,     0,
/

&fdda
/

&dynamics
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                            = 0,
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,     
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

packard

同意楼上，要思考思考，一次连续跑那么长时间的科学目的是什么。

目测楼主的模式并不是全球覆盖的，那么最终的结果就会很受输入文件影响。
而一般下载的输入资料，比如再分析资料，并不是一次独立连续跑下来的，而是一直在吃卫星的观测。
这样的话，卫星观测在影响输入资料，输入资料在影响你的WRF，所以你的WRF也就不是独立的预报，如果所谓提前一年报出什么事件，也是作弊（偷看卫星资料）。
所以呢，既然不独立了，也就不必一次连续跑下来。方法一是分开很多次来初始化，来跑，这样并不声称WRF能提前n个月的预报能力，做的科学问题可以是把再分析资料弄到更精细的模式里，等等。（比如：Du et al. 2014, https://journals.ametsoc.org/doi/pdf/10.1175/JCLI-D-13-00571.1  sect. 2a 他们每天的00UTC初始化模式，跑36小时，取后面的大概24小时）

方法二呢，就是做个全球模式（可以挺费CPU的），这样只要输入资料做初始条件，不用边界条件。但是呢，WRF长期预报还是一般需要海温的，否则会跑偏，要设置namelist里的sst_update。虽然也可以说作弊偷看海温，但除此以外没有偷看别的信息，还可以声称提前几个月的气候预报能力。（当然肯定不能做逐日天气预报了，E.N. Lorenz的混沌理论告诉我们的）

帅气小葛

我去，跑一年算个啥，我们院一个老师用WRF做东亚地区动力降尺度，都是直接跑30年，成果照样发在climate dynamics上。

小魔孤

跑全年时间太长了吧，后面的结果可能会很不合理，建议查一下文献，看看长时间的模拟怎样跑比较合理。

physzw

packard 发表于 2018-2-27 13:41
同意楼上，要思考思考，一次连续跑那么长时间的科学目的是什么。

目测楼主的模式并不是全球覆盖的，那么 ...

非常感谢您多次耐心指导。

physzw

小魔孤 发表于 2018-2-27 12:15
跑全年时间太长了吧，后面的结果可能会很不合理，建议查一下文献，看看长时间的模拟怎样跑比较合理。

非常感谢！我再仔细看看

qq531884069

刚好跟老板做过相关工作  说几点个人浅见
1.楼上的各位说的很正确   跑连续1年的目的是什么要明确
2.如果楼主是一定要跑一年  并且对数据精度要求比较高   建议进行数据同化(我们组就是做高分辨率资料同化的)     不同化数据直接跑一年   个人认为数据没有可参考性
个人浅见~

金陵蔡肉包

三层嵌套跑一年很占资源的，要考虑好。是否需要同化倒也不一定，主要在于模式结果用来干什么。

physzw

金陵蔡肉包 发表于 2018-3-2 08:37
三层嵌套跑一年很占资源的，要考虑好。是否需要同化倒也不一定，主要在于模式结果用来干什么。

谢谢！只要结果准确，占资源并不担心，因为新买的服务器整天也没事。
模式结果主要用来做大气预测。以前的高空数据都是买别人的。

physzw

qq531884069 发表于 2018-2-27 13:56
刚好跟老板做过相关工作  说几点个人浅见
1.楼上的各位说的很正确   跑连续1年的目的是什么要明确
2.如果 ...

谢谢指导！我现在就开始学习使用WRFDA同化

风之牧语

看楼主的目的是什么，如果是跑一年作为气候模拟的话，里面的一些设置需要简单修改更加合理（PS：这也不是绝对的）。如果是为了做气候模拟，建立搜一下WRF做气候模拟的相关设定，网上很多。现在用WRF气候模拟已经非常多了，几十年的都有。