请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8567|回复: 8

WRF-Chem中止

[复制链接]

新浪微博达人勋

发表于 2018-7-6 10:51:29 | 显示全部楼层 |阅读模式

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

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

x
WRF-Chem3.6.1在跑2015年一个多月的模拟时突然中止,error文件显示如下:
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:      76
mieaer /refr/ outside range 1e-3 - 10 refr=     0.10016E+02
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0


似乎是数据溢出?偶尔会出现这种情况,有没有遇到的呢,或者知道问题所在的

namelist.input文件如下
&time_control
run_days                            = 37,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2015, 2999, 2999,
start_month                         = 08,   06,   06,
start_day                           = 14,   11,   11,
start_hour                          = 00,   12,   12,
start_minute                        = 00,   00,   00,
start_second                        = 00,   00,   00,
end_year                            = 2015, 2999, 2999,
end_month                           = 09,   06,   06,
end_day                             = 20,   12,   12,
end_hour                            = 00,   12,   12,
end_minute                          = 00,   00,   00,
end_second                          = 00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,
history_interval                    = 60,
frames_per_outfile                  = 19992, 1000, 1000,
restart                             = .false.,
restart_interval                    = 5000,
auxinput5_interval_m                = 86400,
auxinput7_interval_m                = 86400,  
auxinput8_interval_m                = 86400,
auxinput6_inname                    = 'wrfbiochemi_d<domain>',
auxinput7_inname                    = 'wrffirechemi_d<domain>',
auxinput8_inname                    = 'wrfchemi_gocart_bg_d<domain>',
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
io_form_auxinput4                   = 0
io_form_auxinput5                   = 2,
io_form_auxinput6                   = 2,
io_form_auxinput7                   = 2,
io_form_auxinput8                   = 2,
debug_level                         = 00
/

&domains
time_step                           = 30,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 1,
e_we                                = 60,
e_sn                                = 60,
e_vert                              = 31,
dx                                  = 12000, 20000, 6666.66,
dy                                  = 12000, 20000, 6666.66,
p_top_requested                     = 5000,
num_metgrid_levels                  = 27,
num_metgrid_soil_levels             = 4,
grid_id                             = 1,     2,     3,
parent_id                           = 0,     1,     2,
i_parent_start                      = 1,     30,    30,
j_parent_start                      = 1,     20,    30,
parent_grid_ratio                   = 1,     3,     3,
parent_time_step_ratio              = 1,     3,     3,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 2,     2,     2,
progn                               = 0,     0,     0,
naer                                = 1e9
ra_lw_physics                       = 1,     1,     1,
ra_sw_physics                       = 2,     2,     2,
radt                                = 30,    10,    10,
sf_sfclay_physics                   = 1,     1,     1,
sf_surface_physics                  = 2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,
bldt                                = 1,     0,     0,
cu_physics                          = 5,     5,     0,
cudt                                = 1,     1,     1,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
surface_input_source                = 1,
num_soil_layers                     = 4,
sf_urban_physics                    = 0,     0,     0,
maxiens                             = 1,
maxens                              = 3,
maxens2                             = 3,
maxens3                             = 16,
ensdim                              = 144,
cu_rad_feedback                     = .true.,
/

&fdda
/

&dynamics
w_damping                           = 1,
diff_opt                            = 1,      1,      1,
km_opt                              = 4,      4,      4,
diff_6th_opt                        = 0,      0,      0,
diff_6th_factor                     = 0.12,   0.12,   0.12,
base_temp                           = 290.
damp_opt                            = 0,
zdamp                               = 5000.,  5000.,  5000.,
dampcoef                            = 0.2,    0.2,    0.2
khdif                               = 0,      0,      0,
kvdif                               = 0,      0,      0,
non_hydrostatic                     = .true., .true., .true.,
moist_adv_opt                       = 1,      1,      1,     
scalar_adv_opt                      = 1,      1,      1,     
chem_adv_opt                        = 1,      1,      1,     
/

&bdy_control
spec_bdy_width                      = 5,
spec_zone                           = 1,
relax_zone                          = 4,
specified                           = .true., .false.,.false.,
nested                              = .false., .true., .true.,
/

&grib2
/

&chem
kemit                               = 1,
chem_opt                            = 303,        2,
bioemdt                             = 30,       30,
photdt                              = 30,       30,
chemdt                              = 2.,       2.,
io_style_emissions                  = 1,
emiss_inpt_opt                      = 1,        1,
emiss_opt                           = 5,        3,
chem_in_opt                         = 0,        0,
phot_opt                            = 1,        1,
gas_drydep_opt                      = 1,        1,
aer_drydep_opt                      = 1,        1,
bio_emiss_opt                       = 3,        1,
ne_area                              = 101,
dust_opt                            = 1,
dmsemis_opt                         = 1,
seas_opt                            = 1,
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,
biomass_burn_opt                    = 1,        0,
plumerisefire_frq                   = 30,       30,
aer_ra_feedback                     = 0,        0,
have_bcs_chem                       = .false., .false.,
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/



密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-6 14:33:27 | 显示全部楼层
你这是算复折射指数的时候溢出了吧。。。试试调一下光学相关的参数看看。或者也可能是chemdt的问题。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-7 15:45:59 | 显示全部楼层
astiny 发表于 2018-7-6 14:33
你这是算复折射指数的时候溢出了吧。。。试试调一下光学相关的参数看看。或者也可能是chemdt的问题。

请问负折射指数是哪个选项控制呢。      查了一下官方手册里给出的chemdt是1.5. 我试一下。
另外并不是每次模拟都会这样,是不是也和气象场和源清单有些关系。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-9 14:53:55 | 显示全部楼层
复折射指数应该只有在计算光学特性的时候才会调用到。所以可以尝试调一下phot_opt或者phot_dt.
chemdt的话可能还是和前面dt一样,0.5比较好一点。
和气象场应该没什么关系。和清单关系可能也不大(这一点我也不太清楚,是胡说的。)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-12 08:41:33 | 显示全部楼层
astiny 发表于 2018-7-9 14:53
复折射指数应该只有在计算光学特性的时候才会调用到。所以可以尝试调一下phot_opt或者phot_dt.
chemdt的话 ...

学习了,我试一下
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2022-3-20 23:55:14 | 显示全部楼层
您好,我也遇到了相同的问题,我想请问一下您解决了吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-3-22 14:41:50 | 显示全部楼层
啵啵不要香菜 发表于 2022-3-20 23:55
您好,我也遇到了相同的问题,我想请问一下您解决了吗?

Q: I am using chem_opt=300 with biomass burning emissions and GOCART background. I get this error when I turn on aerosol feedback:

-------------- FATAL CALLED ---------------

FATAL CALLED FROM FILE: LINE: 76

mieaer /refr/ outside range 1e-3 - 10 refr= 0.10003E+02

A: The error is due to the calculated refractive index (10.003) slightly exceeding the preset upper limit (10). The solution is to increase the upper limit a little bit. For example, to revise the upper limit into 15, make the following modification in WRFV3/chem/module_optical_averaging.F, and recompile the model.

Line 4678: if(abs(refr).gt.10.0.or.abs(refr).le.0.001)then

--> if(abs(refr).gt.15.0.or.abs(refr).le.0.001)then
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2022-3-22 22:20:15 | 显示全部楼层
bmbmnb 发表于 2022-3-22 14:41
Q: I am using chem_opt=300 with biomass burning emissions and GOCART background. I get this error  ...

好的谢谢您,我试一下!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2022-3-28 10:14:01 | 显示全部楼层
啵啵不要香菜 发表于 2022-3-22 22:20
好的谢谢您,我试一下!

我也是看到了这个网页,按照上面说的方法解决了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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