WRF-Chem wrfout_d01 输出文件太大有2.1G，如何控制输出变量

侠骨柔情

WRF-chem输出很大，有2.1G，每运行1一个时次需要3个小时，请问如何减少输出，提高运行效率？

astiny

wrf user guide第五章，run time io.
如果你wrfchem版本低于3.9那就只能改registry。
顺便强烈不建议删变量，不然你老板让你诊断什么东西你诊断不出来，还要重跑，追悔莫及。

侠骨柔情

astiny 发表于 2021-4-27 15:20
wrf user guide第五章，run time io.
如果你wrfchem版本低于3.9那就只能改registry。
顺便强烈不建议删变 ...

我用的是WRFV3.9.1，还想问一下，如何能跑快点儿，一个时次要将近3个小时，实在太慢了

astiny

侠骨柔情 发表于 2021-4-27 15:28
我用的是WRFV3.9.1，还想问一下，如何能跑快点儿，一个时次要将近3个小时，实在太慢了

没办法。唯一办法是提高单核性能。因为wrfchem不支持sm，所以网格有限、能用的cpu核数有限的情况下，只能提高单核性能，比如买主频快的cpu，或者如果开了超线程的话，把超线程关掉。

侠骨柔情

astiny 发表于 2021-4-27 15:43
没办法。唯一办法是提高单核性能。因为wrfchem不支持sm，所以网格有限、能用的cpu核数有限的情况下，只能 ...

好的，非常感谢！

侠骨柔情

楼主，你好，想问下，wrfchem输出的PM2.5和PM10均为零，其他变量比如no2、O3都有数值，这是怎么回事？我用的WRF3.9,namelist.input 设置如下:
&time_control
run_days                            = 3,
run_hours                           = 13,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2013,
start_month                         = 06,
start_day                           = 25,
start_hour                          = 12,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2013,
end_month                           = 06,
end_day                             = 29,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 21600,
input_from_file                     = .true.,
history_interval                    = 60,
frames_per_outfile                  = 1,
restart                             = .false.,
restart_interval                    = 43200,  
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
auxinput5_interval_m                = 60,
io_form_auxinput5                   = 2
auxinput6_interval_m                = 43200,
auxinput6_inname                    ='wrfbiochemi_d01',
io_form_auxinput6                   = 2,
debug_level                         = 0,
auxinput11_interval_s               = 120, 120, 120, 120,
auxinput11_end_h                    = 192, 30, 30, 30,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 301,
e_sn                                = 301,
e_vert                              = 35,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 6000,
dy                                  = 6000,
grid_id                             = 1,
parent_id                           = 0,
i_parent_start                      = 1,
j_parent_start                      = 1,
parent_grid_ratio                   = 1,
parent_time_step_ratio              = 1,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 10,
progn                               = 1,
ra_lw_physics                       = 1,
ra_sw_physics                       = 1,
radt                                = 10,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 1,
bl_pbl_physics                      = 1,
bldt                                = 1,
cu_physics                          = 1,
cudt                                = 0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,
num_urban_layers                    = 1200,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .false.
/

&fdda
obs_nudge_opt                       = 0,0,0,0,0
max_obs                             = 150000,
fdda_start                          =   0.,0.,     0.,     0.,     0.
fdda_end                            = 14400.
obs_nudge_wind                      = 1,1,1,0,0
obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp                      = 1,0,0,0,0
obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois                      = 1,1,1,1,1
obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy                           = 1080.,360.,90.,30.,180
obs_rinsig                          = 0.1,
obs_twindo                          = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi                            = 10,
obs_ionf                            = 2,2, 2, 2, 2,
obs_idynin                          = 0,
obs_dtramp                          = 40.,
obs_prt_freq                        = 10, 10, 10, 10, 10,
obs_prt_max                         = 10
obs_ipf_errob                       = .true.
obs_ipf_nudob                       = .true.
obs_ipf_in4dob                      = .true.
obs_ipf_init                        = .true.
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 2,
base_temp                           = 290.,
damp_opt                            = 0,
zdamp                               = 5000.,
dampcoef                            = 0.2,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 1,
chem_adv_opt                        = 1,
tracer_adv_opt                      = 1,
tke_adv_opt                         = 1,
/

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

&grib2
/

&chem
kemit = 1,
chem_opt = 5,
bioemdt = 30,
photdt = 30.,
chemdt = 10.0,
io_style_emissions = 2,
emiss_inpt_opt = 1,
emiss_opt = 13,
chem_in_opt = 1,
phot_opt = 0,
gas_drydep_opt = 1,
aer_drydep_opt = 1,
bio_emiss_opt = 0,
gas_bc_opt = 1,
gas_ic_opt = 1,
aer_bc_opt = 1,
aer_ic_opt = 1,
gaschem_onoff = 1,
aerchem_onoff = 1,
wetscav_onoff = 0,
cldchem_onoff = 0,
vertmix_onoff = 1,
chem_conv_tr = 0,
conv_tr_wetscav= 1,
conv_tr_aqchem = 1,
biomass_burn_opt = 0,
aer_ra_feedback = 0,
ne_area = 300,
opt_pars_out = 1,
have_bcs_chem = .true.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

侠骨柔情

astiny 发表于 2021-4-27 15:43
没办法。唯一办法是提高单核性能。因为wrfchem不支持sm，所以网格有限、能用的cpu核数有限的情况下，只能 ...

楼主，你好，想问下，wrfchem输出的PM2.5和PM10全部是零，其他变量比如no2、O3都有数值，这是怎么回事？我用的WRF3.9,namelist.input 设置如下:
&time_control
run_days                            = 3,
run_hours                           = 13,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2013,
start_month                         = 06,
start_day                           = 25,
start_hour                          = 12,
start_minute                        = 00,
start_second                        = 00,
end_year                            = 2013,
end_month                           = 06,
end_day                             = 29,
end_hour                            = 00,
end_minute                          = 00,
end_second                          = 00,
interval_seconds                    = 21600,
input_from_file                     = .true.,
history_interval                    = 60,
frames_per_outfile                  = 1,
restart                             = .false.,
restart_interval                    = 43200,  
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
auxinput5_interval_m                = 60,
io_form_auxinput5                   = 2
auxinput6_interval_m                = 43200,
auxinput6_inname                    ='wrfbiochemi_d01',
io_form_auxinput6                   = 2,
debug_level                         = 0,
auxinput11_interval_s               = 120, 120, 120, 120,
auxinput11_end_h                    = 192, 30, 30, 30,
/

&domains
time_step                           = 60,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 301,
e_sn                                = 301,
e_vert                              = 35,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
dx                                  = 6000,
dy                                  = 6000,
grid_id                             = 1,
parent_id                           = 0,
i_parent_start                      = 1,
j_parent_start                      = 1,
parent_grid_ratio                   = 1,
parent_time_step_ratio              = 1,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 10,
progn                               = 1,
ra_lw_physics                       = 1,
ra_sw_physics                       = 1,
radt                                = 10,
sf_sfclay_physics                   = 1,
sf_surface_physics                  = 1,
bl_pbl_physics                      = 1,
bldt                                = 1,
cu_physics                          = 1,
cudt                                = 0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
num_land_cat                        = 21,
sf_urban_physics                    = 0,
num_urban_layers                    = 1200,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .false.
/

&fdda
obs_nudge_opt                       = 0,0,0,0,0
max_obs                             = 150000,
fdda_start                          =   0.,0.,     0.,     0.,     0.
fdda_end                            = 14400.
obs_nudge_wind                      = 1,1,1,0,0
obs_coef_wind                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp                      = 1,0,0,0,0
obs_coef_temp                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois                      = 1,1,1,1,1
obs_coef_mois                       = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy                           = 1080.,360.,90.,30.,180
obs_rinsig                          = 0.1,
obs_twindo                          = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi                            = 10,
obs_ionf                            = 2,2, 2, 2, 2,
obs_idynin                          = 0,
obs_dtramp                          = 40.,
obs_prt_freq                        = 10, 10, 10, 10, 10,
obs_prt_max                         = 10
obs_ipf_errob                       = .true.
obs_ipf_nudob                       = .true.
obs_ipf_in4dob                      = .true.
obs_ipf_init                        = .true.
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 4,
diff_6th_opt                        = 2,
base_temp                           = 290.,
damp_opt                            = 0,
zdamp                               = 5000.,
dampcoef                            = 0.2,
khdif                               = 0,
kvdif                               = 0,
non_hydrostatic                     = .true.,
moist_adv_opt                       = 1,
scalar_adv_opt                      = 1,
chem_adv_opt                        = 1,
tracer_adv_opt                      = 1,
tke_adv_opt                         = 1,
/

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

&grib2
/

&chem
kemit = 1,
chem_opt = 5,
bioemdt = 30,
photdt = 30.,
chemdt = 10.0,
io_style_emissions = 2,
emiss_inpt_opt = 1,
emiss_opt = 13,
chem_in_opt = 1,
phot_opt = 0,
gas_drydep_opt = 1,
aer_drydep_opt = 1,
bio_emiss_opt = 0,
gas_bc_opt = 1,
gas_ic_opt = 1,
aer_bc_opt = 1,
aer_ic_opt = 1,
gaschem_onoff = 1,
aerchem_onoff = 1,
wetscav_onoff = 0,
cldchem_onoff = 0,
vertmix_onoff = 1,
chem_conv_tr = 0,
conv_tr_wetscav= 1,
conv_tr_aqchem = 1,
biomass_burn_opt = 0,
aer_ra_feedback = 0,
ne_area = 300,
opt_pars_out = 1,
have_bcs_chem = .true.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

哈库拉玛塔塔000

楼主，你的问题解决了吗？