- 积分
- 5079
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2017-2-15
- 最后登录
- 1970-1-1
|
发表于 2018-5-20 15:20:58
|
显示全部楼层
本帖最后由 骨头333 于 2018-5-20 15:23 编辑
抱歉,我想请教您一下排放源文件是怎么做到包含了12个时次的呢
1.我根据https://ruc.noaa.gov/wrf/wrf-che ... alexercises003.html给出的练习进行训练
得到的wrfchemi_00z_d01和wrfchemi_12z_d01都只有一个时次
2.在进行exercise4使用global_emissions 数据作为排放清单的时候,按照官网给出的步骤,运行prepchem和convert_emiss之后已经能够成功得到三个nc文件wrfchemi_d01,wrffirechemi_d01和wrfchemi_gocart_bg_d01,可是我检查人为源 wrfchemi_d01发现每次也都是出现一个时次,这是不是意味着如果我的研究对象是一段时间,接着运行real和wrf的话就需要反复把这段时间中的每个小时都运行一次prepchem和convert_emiss呢
目前我的理解是io_style_emissions=1,frames_per_emissfile=12这两个会产生包含12时次的文件,io_style_emissions=2,frames_per_emissfile=1应该会产生一个时次的文件,然而我调试了多次,始终出现的文件都是一个时间点,而且文件大小是一样的。
exercise3中的namelist.input如下,还望您不吝赐教
run_days = 0,
run_hours = 12,
run_minutes = 0,
run_seconds = 0,
start_year = 2013, 2000, 2000,
start_month = 06, 01, 01,
start_day = 03, 24, 24,
start_hour = 12, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2013, 2000, 2000,
end_month = 06, 01, 01,
end_day = 03, 25, 25,
end_hour = 23, 12, 12,
end_minute = 00, 00, 00,
end_second = 00, 00, 00,
interval_seconds = 3600,
input_from_file = .true.,.true.,.true.,
history_interval = 60, 60, 60,
frames_per_outfile = 12, 12, 12,
restart = .false.,
restart_interval = 14400,
!auxinput4_inname = "wrflowinp_d<domain>"
!auxinput4_interval = 360., 360., 360.,
!io_form_auxinput4 = 2
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput2 = 2, ! to run ndown successfully
auxinput5_interval_m = 60, 60, 60, ! Anthropogenic emissions input
io_form_auxinput5 = 2,
debug_level = 0
/
&domains
time_step = 45,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
s_we = 1,
e_we = 81,
s_sn = 1,
e_sn = 51,
s_vert = 1,
e_vert = 31,
num_metgrid_levels = 27,
dx = 60000,
dy = 60000,
grid_id = 1,
parent_id = 0,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 0,
smooth_option = 0
numtiles = 1
p_top_requested = 5000,
/
&physics
mp_physics = 4,
ra_lw_physics = 4,
ra_sw_physics = 4,
radt = 5,
sf_sfclay_physics = 5,
sf_surface_physics = 1,
bl_pbl_physics = 5,
bldt = 0,
cu_physics = 5,
cudt = 0,
isfflx = 1,
ifsnow = 0,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat = 20,
sf_urban_physics = 1,
mp_zero_out = 2,
sst_update = 0
cu_rad_feedback = .true.
cu_diag = 1
ishallow = 1
/
&fdda
/
&dynamics
w_damping = 1,
diff_opt = 1,
km_opt = 4,
diff_6th_opt = 2,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000.,
dampcoef = 0.2,
khdif = 0,
kvdif = 0,
smdiv = 0.1,
emdiv = 0.01,
epssm = 0.1,
time_step_sound = 4,
h_mom_adv_order = 5,
v_mom_adv_order = 3,
h_sca_adv_order = 5,
v_sca_adv_order = 3,
non_hydrostatic = .true.,
moist_adv_opt = 2,
scalar_adv_opt = 2,
tke_adv_opt = 2,
do_avgflx_em = 1,
do_avgflx_cugd = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,
nested = .false.,
/
&grib2
/
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
&chem
chem_opt = 108, 108,
have_bcs_chem = .false., .true.,
chem_in_opt = 0, 0
gaschem_onoff = 1, 1,
aerchem_onoff = 1, 0,
vertmix_onoff = 1, 1,
chem_conv_tr = 1, 0,
gas_drydep_opt = 1, 1,
aer_drydep_opt = 1, 0,
chemdt = 0, 0,
bioemdt = 30, 30,
emiss_inpt_opt = 1, 1, !
emiss_opt = 3, 3, ! 16- for CO2/CO simulations
! 16- for CO2/CO simulations
kemit = 20,
io_style_emissions = 1, ! 2- Date/time specific emissions data files used
aircraft_emiss_opt = 0,
bio_emiss_opt = 0, 0, ! 2- includes biogenic emissions reference fields in wrfinput data file and modify values online based upon the weather
phot_opt = 1, 1, ! 1- uses Madronich photolysis (TUV)
photdt = 30, 30,
wetscav_onoff = 0, 0, ! Wet removal by resolved precipitation
cldchem_onoff = 0, 0, ! cloud chemistry turned off in the simulation, also see the chem_opt parameter
conv_tr_aqchem = 0, 0, ! CMAQ aqeuous chemistry for sub-grid clouds
conv_tr_wetscav = 0, 0, ! Wet removal by sub-grid precipitation
seas_opt = 0, ! no GOCART sea salt emissions, 2- MOSAIC or MADE/SORGAM sea salt emissions
dust_opt = 0, ! no GOCART dust emissions included
dmsemis_opt = 0, ! no GOCART dms emissions from sea surface
biomass_burn_opt = 0, 0, ! no biomass burning emissions
plumerisefire_frq = 30, 30, ! time interval for calling the biomass burning plume rise subroutine
gas_bc_opt = 1, 1, ! uses default boundary profile
gas_ic_opt = 1, 0, ! uses default initial condition profile
! 16 is for GHGs
aer_bc_opt = 1, 1,
aer_ic_opt = 1, 0,
aer_ra_feedback = 0, 0, ! feedback from the aerosols to the radiation schemes turned on, see also chem_opt parameter
opt_pars_out = 0, ! no optical properties output
/ |
|