mozart 更新化学初始场和侧边界场 问题

xyan88

各位，请问我用下载的mozart4geos5的数据来初始化初始场，PM2.5第0个时刻的结果为图所示，污染中心偏到甘肃那边去了 请问什么原因啊。多谢啦



namelist.input为

   &time_control
run_days                            = 0,
run_hours                           = 240,
run_minutes                         = 0,
run_seconds                         = 0,

start_year                          = 2016, 2016, 2999,
start_month                         = 08,   08,   06,
start_day                           = 14,   18,   11,
start_hour                          = 00,   12,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2016, 2016, 2999,
end_month                           = 08,   08,   06,
end_day                             = 24,   19,   12,
end_hour                            = 00,   12,   12,
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                       = 50000,
write_hist_at_0h_rst                = .true.

io_form_history                       = 2
io_form_restart                        = 2
io_form_input                          = 2
io_form_boundary                    = 2

debug_level                         =300 0
/



&domains
time_step                           = 180,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 160,    112,   94,
e_sn                                = 180,    97,    91,
e_vert                              = 27,    20,    20,
dx                                  = 20000, 20000, 6666.66,
dy                                  = 20000, 20000, 6666.66,
p_top_requested                     = 5000,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     30,    30,
j_parent_start                      = 1,     20,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 0,
smooth_option                       = 0
zap_close_levels                    = 50
interp_type                         = 1   
t_extrap_type                       = 3   
force_sfc_in_vinterp                = 0  
lagrange_order                      = 1   
lowest_lev_from_sfc                         = .true.      
/

&physics
mp_physics                          = 2,     2,     2,  
progn                               = 1,     1,         
ra_lw_physics                       = 4,     4,     1,
ra_sw_physics                       = 4,     4,     1,  
radt                                = 20,    20,    30,
sf_sfclay_physics                   = 1,     2,     1,
!naer                               = 1e9,            
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     2,     1,
bldt                                = 0,     0,     0,  
cudt                                = 0,     0,     5,  
cu_physics                          = 3,     0,     0,  
cu_diag                                               = 1,     0,        
shcu_physics                        =1,     
isfflx                              = 1,   
ifsnow                              = 1,   
icloud                              = 1,   
surface_input_source                = 1,   
num_soil_layers                     = 4,   
sf_urban_physics                    = 1,     1,     0,   
cugd_avedx                          = 1,  
mp_zero_out                         = 2,   
mp_zero_out_thresh                  = 1.e-8,
maxiens                             = 1,     
maxens                              = 3,     
maxens2                             = 3,     
maxens3                             = 16,   
ensdim                              = 144,   
cu_rad_feedback                     = .true.,
num_land_cat                        = 20,
/







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

&grib2
/

&chem
kemit                               = 20,
chem_opt                            = 9,        42,
!bioemdt                            = 30,       30,
photdt                              = 30,       30,
chemdt                              = 20,       2.,
io_style_emissions                  = 2,               
emiss_inpt_opt                      = 102,        1,   
emiss_opt                           = 4,        3,      
chem_in_opt                         = 1,        0,
phot_opt                            = 2,        1,
gas_drydep_opt                      = 1,        1,
aer_drydep_opt                      = 1,        1,
bio_emiss_opt                       = 0,        1,  !3
ne_area                             = 200,
dust_opt                            = 3,
dmsemis_opt                         = 0,
aer_op_opt                          = 1,       1,   1,
opt_pars_out                        = 1,
seas_opt                            = 0,
gas_bc_opt                          = 101,        1,
gas_ic_opt                          = 101,        1,
aer_bc_opt                          = 101,        1,
aer_ic_opt                          = 101,        1,
gaschem_onoff                       = 1,        1,
aerchem_onoff                       = 1,        1,
wetscav_onoff                       = 1,        0,
cldchem_onoff                       = 1,        0,
vertmix_onoff                       = 1,        1,
chem_conv_tr                        = 1,        1,
biomass_burn_opt                    = 0,        0,
plumerisefire_frq                   = 120,       30,
aer_ra_feedback                     = 0,        0,
have_bcs_chem                       = .true., .false.,
have_bcs_tracer                     = .true., .true.,.true.,
/

andrewsoong

0h的，那你有没有画一下原始数据看一下呢？最好把原始数据和wrf数据都画图看看。

astiny

可以把mozart4geos5那一时刻的场和PM2.5的排放源画一下看看。。。。

xyan88

andrewsoong 发表于 2016-12-11 22:34
0h的，那你有没有画一下原始数据看一下呢？最好把原始数据和wrf数据都画图看看。

你好，我发现是初始dust的原因，用的是CBMZ-MOSAIC_4bins.inp，如果不加oin_a01、oin_a02、oin_a03、oin_a04这四行就不会出现这种情况，但是不知道能不能直接把这四行直接删掉？或者有其他解决方案吗。

xyan88

astiny 发表于 2016-12-12 09:24
可以把mozart4geos5那一时刻的场和PM2.5的排放源画一下看看。。。。

你好，我发现是初始dust的原因，用的是CBMZ-MOSAIC_4bins.inp，如果不加oin_a01、oin_a02、oin_a03、oin_a04这四行就不会出现这种情况，但是不知道能不能直接把这四行直接删掉？或者有其他解决方案吗。

astiny

xyan88 发表于 2016-12-30 20:35
你好，我发现是初始dust的原因，用的是CBMZ-MOSAIC_4bins.inp，如果不加oin_a01、oin_a02、oin_a03、oin_ ...

肯定能啊。。。。

悠然一梦123

楼主，请问你这张图用的是全球排放源吗

xyan88

悠然一梦123 发表于 2017-3-16 18:36
楼主，请问你这张图用的是全球排放源吗

是的 数据可以在mozbc官网下载

helloClimate

xyan88 发表于 2017-3-18 09:04
是的 数据可以在mozbc官网下载

楼主你好~我最近正在学习WRF-Chem的使用，有一些东西还不太会，想向您请教一下，不知道方不方便留个联系方式呢？（QQ之类的最好了~）谢谢~

helloClimate

楼主我先问个比较傻的问题哈，你转换人为源污染数据的时候用的是prep_chem的程序还是anthro_emiss的程序呢？我感觉好像这两套做法都可以弄，但是不太懂区别。。。