用WRF-chem做出来的2.5

Gumho

这是我的namelist.input
&time_control
run_days                            = 31,
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                            = 2017, 2017, 2017,
end_month                           = 02,   02,  02,
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                    = 0,
io_form_history                     = 2,
io_form_restart                     = 2,
io_form_input                       = 2,
io_form_boundary                    = 2,
debug_level                         = 0,
auxinput5_inname                    = 'wrfchemi_d01',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
auxinput5_interval_m                = 21600,
auxinput7_interval_m                = 21600,
auxinput8_interval_m                = 21600,
io_form_auxinput5                   = 2,
/

&dfi_control
/

&domains
time_step                           = 120,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
s_we                                =   1,  
e_we                                = 74,   88,   91,
s_sn                                =   1,
e_sn                                = 61,   70,   82,
e_vert                              = 31,   31,   31,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 27000,  9000, 1000,
dy                                  = 27000,  9000, 1000,
grid_id                             = 1,   2,   3,
parent_id                           = 1,   1,   2,
i_parent_start                      = 1,   31,  20,
j_parent_start                      = 1,   23,  13,
parent_grid_ratio                   = 1,   3,   9,
parent_time_step_ratio              = 1,   3,   9,
p_top_requested                     = 5000,
feedback                            = 1,
smooth_option                       = 0,
zap_close_levels                    = 50,
interp_type                         = 1,
t_extrap_type                       = 2,
force_sfc_in_vinterp                = 0,
use_levels_below_ground             = .true.,
use_surface                         = .true.,
lagrange_order                      = 1,
/
sfcp_to_sfcp                        = .true.,

&physics
mp_physics                          = 2,
progn                               = 0,
ra_lw_physics                       = 1,   1,   1,
ra_sw_physics                       = 2,   2,   2,
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                          = 3,   3,   3,
cu_diag                             = 1,
cudt                                = 5,   5,   5,
ishallow                            = 0,
isfflx                              = 1,
ifsnow                              = 0,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,   0,   0,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-12,         
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,   
cu_rad_feedback                     = .true.,
num_land_cat                        = 24,
/


&fdda
/

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

&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,
/

&chem
kemit                               = 1,
chem_opt                            = 301,
bioemdt                             = 0,
photdt                              = 30,
chemdt                              = 30,
io_style_emissions                  = 1,
emiss_opt                           = 5,      0,
emiss_inpt_opt                      = 1,
emiss_opt_vol                       = 0,      0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 0,      0,
phot_opt                            = 2,      2,
gas_drydep_opt                      = 1,      1,                     
aer_drydep_opt                      = 1,      1,
bio_emiss_opt                       = 1,      1,
dust_opt                            = 1,
dmsemis_opt                         = 1,
seas_opt                            = 1,
depo_fact                           = 0.25,
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                       = 0,      0,
cldchem_onoff                       = 0,      0,
vertmix_onoff                       = 1,      1,
chem_conv_tr                        = 1,      1,
conv_tr_wetscav                     = 0,      0,
conv_tr_aqchem                      = 0,      0,
biomass_burn_opt                    = 1,      1,
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,
/

这是我用grads查看的生成的2017-01-01.ctl的数据
ga-> open 2017-01-01.ctl
Scanning description file:  2017-01-01.ctl
Data file 2017-01-01.dat is open as file 1
LON set to 102.162 124.662
LAT set to 15.09 31.6306
LEV set to 1 1
Time values set: 2017:1:1:0 2017:1:1:0
E set to 1 1
Notice: Implied interpolation for file 2017-01-01.ctl
Interpolation will be performed on any data displayed from this file
ga-> q dim
Default file number is: 1
X is varying   Lon = 102.162 to 124.662   X = 1 to 186
Y is varying   Lat = 15.09 to 31.6306   Y = 1 to 137
Z is fixed     Lev = 1  Z = 1
T is fixed     Time = 00Z01JAN2017  T = 1
E is fixed     Ens = 1  E = 1
ga->

这是我提取数据写的脚本
'reinit'

'open /home/hp/li/2018-10-01.ctl'

'set timelab off'
'set grads off'

latmin='22.567'
latmax='23.256'
lonmin='113.466'
lonmax='114.298'

'set t 1 744'

'define var0=tloop(aave(pm2_5_dry,lon= 'lonmin',lon= 'lonmax',lat= 'latmin',lat= 'latmax'))'
'd var0'
file="/home/hp/li/pm2_5_dry.xls"
rc=write(file,result)
;

提取出来了数据是只有一个小时的

希望各位大神可以帮帮忙。

Gumho

已解决！在第一步修改prep.....inp那里的倒数几行的project-RAMS-grib....那里改为YES！