WRF-ChemV3.9.1 ，gasodesolver_rodas报错

闲云 · 发表于 2024-6-8 13:09:28

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由闲云于 2024-6-8 14:42 编辑

大家好,当我运行WRF-ChemV3.9.1时，发现报错gasodesolver_rodas failure no.1000，我尝试修改我的namelist.input，例如减小积分步长dt和修改max_dom=1，但又模拟一个小时后仍然报同样的错误。我怀疑我的排放源制作可能有问题，但是我用ncview查看，发现没有异常。如果我改变了该怎么办？

&time_control
run_days = 7,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2023,2023,
start_month = 06,06,
start_day = 12,12,
start_hour = 00, 00, 00,
start_minute = 00, 00, 00,
start_second = 00, 00, 00,
end_year = 2023,2023,
end_month = 06,06,
end_day = 19,19,
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, 30
frames_per_outfile = 1000, 1000, 1000
restart = .false.,
restart_interval = 1440,
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
auxinput5_interval_m = 60, 60, 60
auxinput6_inname = 'wrfbiochemi_d<domain>',
io_form_auxinput4 = 2,
io_form_auxinput5 = 2,
io_form_auxinput6 = 2,
frames_per_auxinput5 = 24, 24, 24
debug_level = 0
/
&domains
time_step = 36,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
parent_id = 1, 1,2
e_we = 74, 82
e_sn = 75, 73
i_parent_start = 1, 29
j_parent_start = 1, 26
parent_grid_ratio = 1, 3
grid_id = 1, 2, 3
parent_time_step_ratio = 1, 3, 3
p_top_requested = 5000,
num_metgrid_levels = 42,
num_metgrid_soil_levels = 4,
dx = 9000,3000
dy = 9000,3000
smooth_option = 0
num_metgrid_levels = 27
interp_type = 2
extrap_type = 2
t_extrap_type = 2
use_levels_below_ground = .true.
use_surface = .true.
lagrange_order = 1
zap_close_levels = 5000
lowest_lev_from_sfc = .false.
force_sfc_in_vinterp = 1
sfcp_to_sfcp = .false.
smooth_cg_topo = .false.
e_vert = 42,42,42,
eta_levels =1.0000, 0.9975, 0.9932, 0.9880, 0.9825, 0.9765, 0.9700, 0.9630,
0.9555, 0.9474, 0.9386, 0.9291, 0.9189, 0.9079, 0.8961, 0.8834,
0.8697, 0.8551, 0.8393, 0.8225, 0.8044, 0.7785, 0.7502, 0.7195, !24
0.6862, 0.6503, 0.6118, 0.5707, 0.5270, 0.4809, 0.4326, 0.3825,
0.3310, 0.2787, 0.2264, 0.1751, 0.1322, 0.0962, 0.0656, 0.0398,
0.0181, 0.0000
/
&physics
mp_physics = 2, 2, 2
mp_zero_out = 20
mp_zero_out_thresh = 1e-08
progn = 1, 1, 1
ra_lw_physics = 4, 4, 4
ra_sw_physics = 4, 4, 4
radt = 9, 3
sf_sfclay_physics = 1, 1, 1
sf_surface_physics = 2, 2, 2
sf_urban_physics = 1, 1, 1
bl_pbl_physics = 1, 1, 1
bldt = 0, 0, 0
cu_physics = 5, 0, 0
cudt = 0, 0, 0
!cu_diag = 1, 0, 0
sf_lake_physics = 1, 1, 1,
cu_rad_feedback = .false., .false., .false.
slope_rad = 0, 0, 0
topo_shading = 0, 0, 0
prec_acc_dt = 60, 60, 60
ishallow = 1,
ifsnow = 0,
! icloud = 1,
surface_input_source = 3,
isfflx = 1
surface_input_source = 1
num_soil_layers = 4
num_land_cat = 21
maxiens = 1
maxens = 3
maxens2 = 3
maxens3 = 16
ensdim = 144
num_urban_layers = 8000
/
&dynamics
rk_ord = 3
diff_opt = 1
km_opt = 4
damp_opt = 0
diff_6th_opt = 0, 0
diff_6th_factor = 0.12, 0.12
dampcoef = 0.01, 0.01
zdamp = 5000, 5000
w_damping = 0
base_temp = 290
khdif = 3000, 3000
kvdif = 0, 0
smdiv = 0.1, 0.1
emdiv = 0.01, 0.01
epssm = 0.1, 0.1
non_hydrostatic = .true.,.true.
pert_coriolis = .false.,.false.
top_lid = .false.,.false.
mix_isotropic = 0, 0
mix_upper_bound = 0.1, 0.1
tke_drag_coefficient = 0, 0
tke_heat_flux = 0, 0
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 = 1, 1
scalar_adv_opt = 1, 1
chem_adv_opt = 1, 1
tke_adv_opt = 1, 1
gwd_opt=1
/
&bdy_control
spec_bdy_width = 5
spec_zone = 1
relax_zone = 4
specified = .true.,.false.
spec_exp = 0
periodic_x = .false.
symmetric_xs = .false.
symmetric_xe = .false.
open_xs = .false.
open_xe = .false.
periodic_y = .false.
symmetric_ys = .false.
symmetric_ye = .false.
open_ys = .false.
open_ye = .false.
nested = .false.,.true.
polar = .false.
/
/
&grib2
/
&chem
kemit = 2
chem_opt = 10,10,
chem_in_opt = 1, 1, 1
chemdt = 0, 0, 0
emiss_opt = 3, 3, 3
io_style_emissions = 2
emiss_inpt_opt = 1,1,1, ! 101, 101
emiss_opt_vol = 0, 0, 0
emiss_ash_hgt = 20000.,
bio_emiss_opt = 3,3,
bioemdt = 30, 30, 30
ne_area = 666
biomass_burn_opt = 0, 0, 0
plumerisefire_frq = 60, 20, 20
!photol
phot_opt = 2, 2, 2
photdt = 30,30,30
dust_opt = 4,
dust_schme = 3,
dmsemis_opt = 0,
seas_opt = 0,
gas_drydep_opt = 1, 1, 1
aer_drydep_opt = 1, 1, 1
depo_fact = 0.25, 0.25, 0.25
have_bcs_chem = .true.,.true.,.true.
gas_bc_opt = 1, 1, 1
gas_ic_opt = 1, 1, 1
aer_bc_opt = 1, 1, 1
aer_ic_opt = 1, 1, 1
gaschem_onoff = 1, 1, 1
aerchem_onoff = 1, 1, 1
wetscav_onoff = 1, 1, 1
cldchem_onoff = 1, 1, 1
vertmix_onoff = 1, 1, 1
chem_conv_tr = 0, 0, 0
conv_tr_wetscav = 1, 1, 1
conv_tr_aqchem = 1, 1, 1
aer_ra_feedback = 1, 1, 1
chemdiag = 1, 1, 1
aer_op_opt = 1
opt_pars_out = 1
/
&namelist_quilt
!???
nio_tasks_per_group = 0
nio_groups = 1
/

复制代码

报错信息：

newnucbb dtnuc 3.60E+01
newnucbb ncnt 2050 0 0 0 0 0 0
coagbb ncntaa 2050 0 357 137 6 0 0 339 0 0
Timing for main: time 2023-06-12_01:00:00 on domain 1: 7.24536 elapsed seconds
Timing for Writing wrfout_d01_2023-06-12_01:00:00 for domain 1: 17.35099 elapsed seconds
mediation_integrate: med_read_wrf_chem_emissions: Read emissions for time 2023-06-12_01:00:00
d01 2023-06-12_01:00:00
d01 2023-06-12_01:00:00 *** gasodesolver_rodas failure no. 1
d01 2023-06-12_01:00:00 iregime, iok, i, j, k / t
d01 2023-06-12_01:00:00 1 -1001 6 2 1
d01 2023-06-12_01:00:00 3.6000000000E+03
d01 2023-06-12_01:00:00 inforodas_cur(1-6) =
d01 2023-06-12_01:00:00 1 16 8 0 7 0
d01 2023-06-12_01:00:00 tstart, tstop, abs_error, rel_error / temp, press, cair, cos_sza =
d01 2023-06-12_01:00:00 3.6000000000E+03 3.6360000000E+03 1.0000000000E+01 1.0000000475E-03
d01 2023-06-12_01:00:00 NaN NaN NaN -9.9000000000E+01
d01 2023-06-12_01:00:00 ngas_z, nrconst_nonzero =
d01 2023-06-12_01:00:00 71 74
d01 2023-06-12_01:00:00 l, name, cboxold, cboxnew for l=1,ngas_z
d01 2023-06-12_01:00:00 1 no NaN NaN
d01 2023-06-12_01:00:00 2 no2 NaN NaN
d01 2023-06-12_01:00:00 3 no3 NaN NaN
d01 2023-06-12_01:00:00 4 n2o5 NaN NaN
d01 2023-06-12_01:00:00 5 hono NaN NaN
d01 2023-06-12_01:00:00 6 hno3 NaN NaN
d01 2023-06-12_01:00:00 7 hno4 NaN NaN
d01 2023-06-12_01:00:00 8 o3 NaN NaN
d01 2023-06-12_01:00:00 9 o1d 0.0000000000E+00 0.0000000000E+00
d01 2023-06-12_01:00:00 10 o3p 0.0000000000E+00 0.0000000000E+00
d01 2023-06-12_01:00:00 11 oh NaN NaN
d01 2023-06-12_01:00:00 12 ho2 NaN NaN
d01 2023-06-12_01:00:00 13 h2o2 NaN NaN
d01 2023-06-12_01:00:00 14 co NaN NaN
d01 2023-06-12_01:00:00 15 so2 NaN NaN
d01 2023-06-12_01:00:00 16 h2so4 NaN NaN
d01 2023-06-12_01:00:00 17 nh3 NaN NaN
d01 2023-06-12_01:00:00 18 hcl NaN NaN
d01 2023-06-12_01:00:00 19 ch4 NaN NaN
d01 2023-06-12_01:00:00 20 c2h6 NaN NaN
d01 2023-06-12_01:00:00 21 ch3o2 NaN NaN
d01 2023-06-12_01:00:00 22 ethp NaN NaN
d01 2023-06-12_01:00:00 23 hcho NaN NaN
d01 2023-06-12_01:00:00 24 ch3oh NaN NaN
d01 2023-06-12_01:00:00 25 c2h5oh NaN NaN
d01 2023-06-12_01:00:00 26 ch3ooh NaN NaN
d01 2023-06-12_01:00:00 27 ethooh NaN NaN
d01 2023-06-12_01:00:00 28 ald2 NaN NaN
d01 2023-06-12_01:00:00 29 hcooh NaN NaN
d01 2023-06-12_01:00:00 30 par NaN NaN
d01 2023-06-12_01:00:00 31 aone NaN NaN
d01 2023-06-12_01:00:00 32 mgly NaN NaN
d01 2023-06-12_01:00:00 33 eth NaN NaN
d01 2023-06-12_01:00:00 34 olet NaN NaN
d01 2023-06-12_01:00:00 35 olei NaN NaN
d01 2023-06-12_01:00:00 36 tol NaN NaN
d01 2023-06-12_01:00:00 37 xyl NaN NaN
d01 2023-06-12_01:00:00 38 cres NaN NaN
d01 2023-06-12_01:00:00 39 to2 NaN NaN
d01 2023-06-12_01:00:00 40 cro NaN NaN
d01 2023-06-12_01:00:00 41 open NaN NaN
d01 2023-06-12_01:00:00 42 onit NaN NaN
d01 2023-06-12_01:00:00 43 pan NaN NaN
d01 2023-06-12_01:00:00 44 rcooh NaN NaN
d01 2023-06-12_01:00:00 45 rooh NaN NaN
d01 2023-06-12_01:00:00 46 c2o3 NaN NaN
d01 2023-06-12_01:00:00 47 ro2 NaN NaN
d01 2023-06-12_01:00:00 48 ano2 NaN NaN
d01 2023-06-12_01:00:00 49 nap NaN NaN
d01 2023-06-12_01:00:00 50 xo2 NaN NaN
d01 2023-06-12_01:00:00 51 xpar NaN NaN
d01 2023-06-12_01:00:00 52 isop NaN NaN
d01 2023-06-12_01:00:00 53 isoprd NaN NaN
d01 2023-06-12_01:00:00 54 isopp NaN NaN
d01 2023-06-12_01:00:00 55 isopn NaN NaN
d01 2023-06-12_01:00:00 56 isopo2 NaN NaN
d01 2023-06-12_01:00:00 57 dms NaN NaN
d01 2023-06-12_01:00:00 58 msa NaN NaN
d01 2023-06-12_01:00:00 59 dmso NaN NaN
d01 2023-06-12_01:00:00 60 dmso2 NaN NaN
d01 2023-06-12_01:00:00 61 ch3so2h NaN NaN
d01 2023-06-12_01:00:00 62 ch3sch2oo NaN NaN
d01 2023-06-12_01:00:00 63 ch3so2 NaN NaN
d01 2023-06-12_01:00:00 64 ch3so3 NaN NaN
d01 2023-06-12_01:00:00 65 ch3so2oo NaN NaN
d01 2023-06-12_01:00:00 66 ch3so2ch2oo NaN NaN
d01 2023-06-12_01:00:00 67 mtf NaN NaN
d01 2023-06-12_01:00:00 68 h2o NaN NaN
d01 2023-06-12_01:00:00 69 o2 NaN NaN
d01 2023-06-12_01:00:00 70 n2 NaN NaN
d01 2023-06-12_01:00:00 71 h2 NaN NaN
d01 2023-06-12_01:00:00 rconst for i=1,nrconst_nonzero
d01 2023-06-12_01:00:00 7.4897240847E-03 1.8635691702E-01 1.6598472139E-03 4.4300270474E-07
d01 2023-06-12_01:00:00 2.7440278245E-06 3.7815876567E-05 4.0247786092E-04 2.0091587430E-05
d01 2023-06-12_01:00:00 5.2978157328E-06 NaN NaN 2.2000000433E-10
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 2.2000000779E-11 NaN NaN NaN
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 NaN NaN 5.0000000181E-16 NaN
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 3.4999999860E-12 1.9999999365E-21 NaN NaN
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 2.2030435503E-05 3.6238096072E-05 9.9999999600E-12 NaN
d01 2023-06-12_01:00:00 4.0311911107E-06 3.7084710129E-06 NaN NaN
d01 2023-06-12_01:00:00 NaN NaN 1.0999999739E-12 2.4999999900E-12
d01 2023-06-12_01:00:00 NaN NaN NaN NaN
d01 2023-06-12_01:00:00 NaN 3.4452730233E-06 NaN NaN
d01 2023-06-12_01:00:00 NaN NaN NaN 3.9999999840E-12
d01 2023-06-12_01:00:00 NaN NaN
d01 2023-06-12_01:00:00
d01 2023-06-12_01:00:00 *** gasodesolver_lsodes failure no. 1
d01 2023-06-12_01:00:00 iregime, iok, i, j, k / t
d01 2023-06-12_01:00:00 1 -1 6 2 1
d01 2023-06-12_01:00:00 3.6000000000E+03

复制代码

闲云 · 发表于 2024-6-12 15:18:51

{:eb303:}{:eb303:}求大神帮助

JNU学术垃圾 · 发表于 2024-6-14 09:55:00

给楼主顶一下我也遇到这样的错误

		自动登录	找回密码
密码			立即注册

WRF-ChemV3.9.1 ，gasodesolver_rodas报错

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

相关帖子

回帖奖励 +5 金钱