wrf-chem转出变量数值异常的问题

SG哭晕在厕所

最近在转一个雾霾的个例，本来PM2.5在近地面的正常浓度应该在200到300μg/m3间，但是转出来却只有0.002左右，PM10，SO2等几种污染物的情况类似，换了几种参数化方案都不行，很费解，有没有老司机帮忙解答一下这个问题。。。以下是input
&time_control            
run_days                 = 4,
run_hours                = 12,
run_minutes              = 0,
run_seconds              = 0,
start_year               = 2016,     2016,      
start_month              = 11,       11,         
start_day                = 02,       02,         
start_hour               = 00,       00,         
start_minute             = 00,       00,         
start_second             = 00,       00,         
end_year                 = 2016,     2016,      
end_month                = 11,       11,         
end_day                  = 06,       06,         
end_hour                 = 12,       12,         
end_minute               = 00,       00,         
end_second               = 00,       00,         
interval_seconds         = 21600,
input_from_file          = .true.,.true.,
history_interval         = 60,       60,         
frames_per_outfile       = 1,        1,         
restart                  = .false.,
restart_interval         = 5000,
auxinput5_inname         = 'wrfchemi_d<domain>_<date>'
auxinput5_interval_m     = 60,60,60,
frames_per_auxinput5     = 1,1,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,       
debug_level              = 0,
/

&domains                 
time_step                = 60,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 2,
s_we                     = 1,        1,      
e_we                     = 155,      130,      
s_sn                     = 1,        1,
e_sn                     = 152,      115,      
s_vert                   = 1,        1,
e_vert                   = 60,       60,
eta_levels               = 1.000, 0.996, 0.991, 0.985, 0.978,
                0.97, 0.961, 0.951, 0.94, 0.928,
                0.916, 0.904, 0.892, 0.88, 0.868,
                0.856, 0.844, 0.832, 0.82, 0.808,
                0.788, 0.768, 0.748, 0.723, 0.698,
                0.673, 0.643, 0.613, 0.583, 0.553,
                0.523, 0.493, 0.463, 0.433, 0.408,
                0.383, 0.363, 0.343, 0.323, 0.303,
                0.288, 0.273, 0.258, 0.243, 0.228,
                0.213, 0.198, 0.183, 0.168, 0.153,
                0.138, 0.123, 0.108, 0.093, 0.078,
                0.063, 0.048, 0.033, 0.018, 0.000,
p_top_requested          = 5000,      
num_metgrid_levels       = 32,       
num_metgrid_soil_levels  = 4,
dx                       = 27000,     9000,
dy                       = 27000,     9000,
grid_id                  = 1,        2,     
parent_id                = 0,        1,     
i_parent_start           = 1,       44,   
j_parent_start           = 1,       64,   
parent_grid_ratio        = 1,        3,     
parent_time_step_ratio   = 1,        3,
feedback                 = 1,
smooth_option            = 0,
use_levels_below_ground             = .true.
use_surface                         = .true.
lagrange_order                      = 1,
sfcp_to_sfcp                        = .true.
/

&physics                 
mp_physics               = 10,       10,   
ra_lw_physics            = 4,        4,     
ra_sw_physics            = 4,        4,
progn                    = 1,        1,   
radt                     = 27,       9,   
sf_sfclay_physics        = 4,        4,     
sf_surface_physics       = 2,        2,     
bl_pbl_physics           = 4,        4,     
bldt                     = 0,        0,     
cu_physics               = 3,        3,     
cu_diag                  = 1,        1,     
cudt                     = 0,        0,     
isfflx                   = 1,
ifsnow                   = 0,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 4,
sf_urban_physics                    = 0,     0,     
cu_rad_feedback                     = .true.,.true.,
/

&fdda                    
/

&dynamics               
rk_ord                   = 3,
w_damping                = 1,
diff_opt                 = 1,
km_opt                   = 4,
diff_6th_opt             = 0,        0,
diff_6th_factor          = 0.12      0.12,
damp_opt                 = 0,
base_temp                = 290.,
zdamp                    = 5000.,    5000.,   
dampcoef                 = 0.01,     0.01,     
khdif                    = 0,        0,            
kvdif                    = 0,        0,                  
non_hydrostatic          = .true.,   .true.,  
time_step_sound          = 4,        4,            
h_mom_adv_order          = 5,        5,            
v_mom_adv_order          = 3,        3,            
h_sca_adv_order          = 5,        5,            
v_sca_adv_order          = 3,        3,            
moist_adv_opt                 = 2,      2,         
scalar_adv_opt                = 2,      2,         
chem_adv_opt                  = 2,      2,         
tke_adv_opt                   = 2,      2,         
/

&bdy_control            
spec_bdy_width           = 5,
spec_zone                = 1,
relax_zone               = 4,
specified                = .true.,  .false.,
periodic_x               = .false.,  .false.,
symmetric_xs             = .false.,  .false.,
symmetric_xe             = .false.,  .false.,
open_xs                  = .false.,  .false.,
open_xe                  = .false.,  .false.,
periodic_y               = .false.,  .false.,
symmetric_ys             = .false.,  .false.,
symmetric_ye             = .false.,  .false.,
open_ys                  = .false.,  .false.,
open_ye                  = .false.,  .false.,
nested                   = .false.,   .true.,
/

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/
   

&chem
kemit                               = 1,
chem_opt                            = 11,          11,
bioemdt                             = 30,          30,
photdt                              = 30,          30,
chemdt                              = 0,          0,
io_style_emissions                  = 2,
emiss_opt                           = 0,          0,
emiss_opt_vol                       = 0,          0,
chem_in_opt                         = 0,          0,
emiss_inpt_opt                      = 1,          1,
phot_opt                            = 2,          2,
gas_drydep_opt                      = 1,          1,
aer_drydep_opt                      = 1,          1,
bio_emiss_opt                       = 1,          1,
dust_opt                            = 0,         
dmsemis_opt                         = 0,         
seas_opt                            = 2,         
gas_bc_opt                          = 1,          1,
gas_ic_opt                          = 1,          1,
aer_bc_opt                          = 1,          1,
aer_ic_opt                          = 1,          1,
gaschem_onoff                       = 1,          1,
aerchem_onoff                       = 1,          1,
wetscav_onoff                       = 1,          1,
cldchem_onoff                       = 1,          1,
vertmix_onoff                       = 1,          1,
chem_conv_tr                        = 1,          1,
conv_tr_wetscav                     = 1,          1,
conv_tr_aqchem                      = 1,          1,
biomass_burn_opt                    = 0          0,
plumerisefire_frq                   = 30,          30,
have_bcs_chem                       = .false., .false., .false.,
aer_ra_feedback                     = 1,         
aer_op_opt                          = 1,         
opt_pars_out                        = 1,         
diagnostic_chem                     = 0,         
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

astiny

源可能没给对。

ayzqs

不会是你单位没搞准吧？

SG哭晕在厕所

ayzqs 发表于 2016-12-5 17:23
不会是你单位没搞准吧？

我也挺怀疑的。。。但是ncdump里PM2.5的单位明明就是微克每立方米啊

astiny

可以把源plot一下，然后发一下。顺便ncdump一下单位。
chem_opt什么的应该不是问题。

astiny

ncdump一下源的单位。照理应该是mol/km2/h。

SG哭晕在厕所

astiny 发表于 2016-12-5 18:04
ncdump一下源的单位。照理应该是mol/km2/h。

源是指wrfchemi*吗？

SG哭晕在厕所

astiny 发表于 2016-12-5 18:04
ncdump一下源的单位。照理应该是mol/km2/h。

我看了一下，里面的单位是ug/m3 m/s，这什么意思。。。

astiny

SG哭晕在厕所 发表于 2016-12-5 18:54
我看了一下，里面的单位是ug/m3 m/s，这什么意思。。。

目测anthro_emis做的源吧233
这个只能用来给mozcart的那四个模式做源。或者用来做逐时次的，然后再在它基础上改。
因为单位不是mol那个，所以数值小。chem_opt读的时候只读数字不读单位。
解决方法。1.prep_chem_src.2找学长学姐要清华源。

SG哭晕在厕所

astiny 发表于 2016-12-5 19:06
目测anthro_emis做的源吧233
这个只能用来给mozcart的那四个模式做源。或者用来做逐时次的，然后再在它 ...

prep_chem_src是什么。。。。然后我输出wrfout是一个时次一个文件的，这算逐时次吗？