爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 6519|回复: 3

wrfchem是否只能每12小时restart一次?

[复制链接]

新浪微博达人勋

发表于 2018-4-18 18:40:09 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
跑wrfchem的时候,由于一些原因,需要每6小时restart一次,restart_interval设为360。排放源文件是用的00Z和12Z两个文件,每个文件12个时次,每天循环的。
然而,当我尝试在第6个小时restart的时候,rsl.error.0000中提示
HOURLY EMISSIONS UPDATE TIME    21600.0       0.0
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2016-
11-01_06:00:00
d01 2016-11-01_06:00:00  MMINLU error on input
            0 input_wrf: wrf_get_next_time current_date: 2016-11-01_06:00:00 Sta
tus =         -102
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     846
  ... Could not find matching time in input file
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
而若是从第12个小时restart则一切正常。
请问这可能是由于什么原因造成的?是因为排放源的原因,所以只能每12小时restart一次吗?
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 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
/
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-3-29 09:36:27 | 显示全部楼层
我也遇到了楼主的问题,有些时候也是6小时会报错,换成12小时就没问题
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-4-28 09:35:31 | 显示全部楼层
顶个先,谢谢楼主
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表