WRF-Chem 求助： exceeded lsodes failure limit = 1000 error

小米_millet

大家好， 我在 run wrf.exe 出现以下问题

当 wrffirechem_d01*  文件中所有变量的值都是 0  的时候可以正常运行；当 wrffirechemi_d01* 文件中变量有非0值的时候，就会报错

-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:      76
*** exceeded lsodes failure limit =        1000


这个错误信息来自于module_cbmz.F。

有rsl.error.0013 文件可以看到， 是从 2006-10-01_04:48:00 开始出错。积分开始不稳定， iok 一直小于 0， gasodesolver_rodas failure no. 开始递增直至超过 1000 然后报错。

d01 2006-10-01_04:48:00  *** gasodesolver_rodas failure no.           1
d01 2006-10-01_04:48:00 iregime, iok, i, j, k / t
d01 2006-10-01_04:48:00            1       -1001          41          71           1
d01 2006-10-01_04:48:00   1.7280000000E+04
d01 2006-10-01_04:48:00 inforodas_cur(1-6) =
d01 2006-10-01_04:48:00            1          16           8           0           7           0
d01 2006-10-01_04:48:00 tstart, tstop, abs_error, rel_error / temp, press, cair, cos_sza =
d01 2006-10-01_04:48:00   1.7280000000E+04  1.7928000000E+04  1.0000000000E+01  1.0000000475E-03
d01 2006-10-01_04:48:00                NaN               NaN  3.9062062569E-05 -9.9000000000E+01
d01 2006-10-01_04:48:00 ngas_z, nrconst_nonzero =
d01 2006-10-01_04:48:00           71          74
d01 2006-10-01_04:48:00 l, name, cboxold, cboxnew for l=1,ngas_z
d01 2006-10-01_04:48:00   1 no                         NaN               NaN
d01 2006-10-01_04:48:00   2 no2                        NaN               NaN
d01 2006-10-01_04:48:00   3 no3                        NaN               NaN
d01 2006-10-01_04:48:00   4 n2o5                       NaN               NaN
d01 2006-10-01_04:48:00   5 hono                       NaN               NaN
d01 2006-10-01_04:48:00   6 hno3                       NaN               NaN
d01 2006-10-01_04:48:00   7 hno4                       NaN               NaN
d01 2006-10-01_04:48:00   8 o3                         NaN               NaN
d01 2006-10-01_04:48:00   9 o1d           0.0000000000E+00  0.0000000000E+00
d01 2006-10-01_04:48:00  10 o3p           0.0000000000E+00  0.0000000000E+00
d01 2006-10-01_04:48:00  11 oh                         NaN               NaN
d01 2006-10-01_04:48:00  12 ho2                        NaN               NaN



求助！ 请问是因为WRF-chem配置的不恰当造成的模式计算不稳定吗？有人用过 chem_opt = 10 这个选项，知道如何配置其它参数使其运行吗？

namelist.input 中的配置如下：


   &physics
mp_physics                          = 10,     2,     2,
progn                               = 1,
ra_lw_physics                       = 4,     1,     1,
ra_sw_physics                       = 4,     1,     1,
radt                                = 27,    30,    30,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     1,
bl_pbl_physics                      = 9,     1,     1,
bldt                                = 0,     0,     0,
cu_physics                          = 7,     5,     0,
cudt                                = 0,     0,     0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
sst_update                          = 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
cu_rad_feedback                     = .false.,
num_land_cat                        = 21,
/

&fdda
grid_fdda                           = 1,  0,  0,  0,
gfdda_inname                        = "wrffdda_d<domain>",
gfdda_end_h                         = 9999, 9999, 96, 96,
gfdda_interval_m                    = 360, 360, 60, 60,
fgdt                                = 0,  0,  0,  0,
if_no_pbl_nudging_uv                = 1,  1,  0,  0,
if_no_pbl_nudging_t                 = 1,  1,  1,  1,
if_no_pbl_nudging_q                 = 1,  1,  0,  0,
if_zfac_uv                          = 0,  0,  0,  0,
  k_zfac_uv                          = 10, 10, 10, 10,
if_zfac_t                           = 0,  0,  0,  0,
  k_zfac_t                           = 10, 10, 10, 10,
if_zfac_q                           = 1,  0,  0,  0,
  k_zfac_q                           = 10, 10, 10, 10,
guv                                 = 0.0003, 0.0000, 0.0001, 0.0001,
gt                                  = 0.0003, 0.0000, 0.0001, 0.0001,
gq                                  = 0.0001, 0.0000, 0.0001, 0.0001,
if_ramping                          = 0,
dtramp_min                          = 60.0,
/


&dynamics
rk_ord                              = 3,
w_damping                           = 1,
diff_opt                            = 1,
km_opt                              = 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,      0,      0,
scalar_adv_opt                      = 2,      0,      0,
chem_adv_opt                        = 2,      0,      0,
tke_adv_opt                         = 2,      0,      0,
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,
/


&chem
kemit                               = 1,
chem_opt                            = 10,
bioemdt                             = 30,
photdt                              = 30,
chemdt                              = 10,
io_style_emissions                  = 2,
emiss_opt                           = 4,          5,
emiss_inpt_opt                      = 101,
emiss_opt_vol                       = 0,          0,
emiss_ash_hgt                       = 20000.,
chem_in_opt                         = 0,          0,
phot_opt                            = 0,          0,
gas_drydep_opt                      = 1,          1,
aer_drydep_opt                      = 1,          1,
bio_emiss_opt                       = 0,          0,
dust_opt                            = 0,
dmsemis_opt                         = 0,
seas_opt                            = 0,
gas_bc_opt                          = 1,          0,
gas_ic_opt                          = 1,          0,
aer_bc_opt                          = 1,          0,
aer_ic_opt                          = 1,          0,
gaschem_onoff                       = 1,          0,
aerchem_onoff                       = 1,          0,
wetscav_onoff                       = 1,          0,
cldchem_onoff                       = 1,          0,
vertmix_onoff                       = 1,          0,
biomass_burn_opt                    = 1,          0,
plumerisefire_frq                   = 60,         0,
have_bcs_chem                       = .false., .false., .false.,
aer_ra_feedback                     = 1,
aer_op_opt                          = 5,
opt_pars_out                        = 1,
diagnostic_chem                     = 0,
chem_conv_tr                        = 0,        1,
/

小米_millet

修改初始时间（start_day），从第二天开始运行就不会出现这个问题。有没有遇到类似问题的小伙伴，知道应该如何debug? 谢谢

小米_millet

修改了区域的范围（e_we, e_sn, CEN_LAT, CEN_LON), 模式积分就变稳定了。

小米_millet

修改了区域的范围（e_we, e_sn, CEN_LAT, CEN_LON), 模式积分就变稳定了。

772412162

小米_millet 发表于 2017-4-19 13:23
修改了区域的范围（e_we, e_sn, CEN_LAT, CEN_LON), 模式积分就变稳定了。

请问楼主当时从什么范围改成了什么范围后积分开始稳定的呢？万分感谢

郭小小

楼主您好，请问你这个是啥原因造成的，和你一摸一样的错误，改了区域还是不行。

xibaokaixin

一样的问题，怎么改都不对，楼主解决了吗

Sun丶Yang

郭小小 发表于 2017-12-2 16:54
楼主您好，请问你这个是啥原因造成的，和你一摸一样的错误，改了区域还是不行。

一样的问题 请问你解决了吗？这个error原因是什么 或者怎么解决呢

芦花浅水

我也有这样的问题，是namelist.input的问题，我的是因为设置没有边界初始条件，然后list里面把chem_in_opts设置成1，出现这样的问题，将这个变量改成0，即可解决问题。

yangziningQX

请问能请教一下关于这个error怎么解决的问题嘛，我的namelist.input里面chem_in_opts设置的就是0