关于WRF-Chem中的上边界条件问题

shiyu

大家好，不知道是否有老师和同学在WRF-Chem中添加过上边界条件（upper boundary condition），我在模拟过程中仅能模拟10天就会出现tropopause_climate: Warning: some tropopause levels still NOTFOUND的错误然后任务中断。想请教一下做过的前辈，非常感谢。

这是我的报错信息：
2018-07-10_16:00:00
**WARNING** Time in input file not being checked **WARNING**
photolysis_driver: called for domain            1
photolysis_driver: called for domain            2
tropopause_climate: Warning: some tropopause levels still NOTFOUND
tropopause_climate: Warning: some tropopause levels still NOTFOUND
tropopause_climate: Warning: some tropopause levels still NOTFOUND

Program received signal SIGSEGV: Segmentation fault - invalid memory reference.

Backtrace for this error:
#0  0x2B18AC1966F7
#1  0x2B18AC196D3E
#2  0x2B18ACC2826F
#3  0x2022B72 in __module_upper_bc_driver_MOD_upper_bc_driver
#4  0x19FDBC3 in chem_driver_
#5  0x1905B69 in solve_interface_
#6  0x47915F in __module_integrate_MOD_integrate
#7  0x4052A3 in __module_wrf_top_MOD_wrf_run

shiyu

这是我的namelist.input，非常感谢大家

&time_control
run_days                            = 30,
run_hours                           = 0,
run_minutes                         = 0,
run_seconds                         = 0,
start_year                          = 2018, 2018, 2018, 2018,
start_month                         = 07,   07,   07,   07,
start_day                           = 01,   01,   01,   01,
start_hour                          = 00,   00,   00,   00,
start_minute                        = 00,   00,   00,   00,
start_second                        = 00,   00,   00,   00,
end_year                            = 2018, 2018, 2018, 2018,
end_month                           = 07,   07,   07,   07,
end_day                             = 31,   31,   31,   31,
end_hour                            = 00,   00,   00,   00,
end_minute                          = 00,   00,   00,   00,
end_second                          = 00,   00,   00,   00,
interval_seconds                    = 21600
input_from_file                     = .true.,.true.,.true.,.true.,
history_interval                    = 60,   60,   60,    60,
frames_per_outfile                  = 1000, 1000, 1000,  1000,     
restart                             = .false.,
restart_interval                    = 3600,
io_form_history                     = 2
io_form_restart                     = 2
io_form_input                       = 2
io_form_boundary                    = 2
debug_level                         = 0
/

&domains
time_step                           = 150,
time_step_fract_num                 = 0,
time_step_fract_den                 = 1,
max_dom                             = 2,
e_we                                = 140,    103
e_sn                                = 120,    100,  
e_vert                              = 59,     59,
p_top_requested                     = 100,
num_metgrid_levels                  = 32,
num_metgrid_soil_levels             = 4,
dx                                  = 27000, 9000,  3000,  1000,
dy                                  = 27000, 9000,  3000,  1000,
feedback                            = 1,
smooth_option                       = 0
/

&physics
mp_physics                          = 2,     2,     2,     2,
mp_zero_out                         = 2,
mp_zero_out_thresh                  = 1.e-8
ra_lw_physics                       = 4,     4,     4,     4,
ra_sw_physics                       = 4,     4,     4,     4,
radt                                = 27,    9,     3,     1,   
sf_sfclay_physics                   = 1,     1,     1,     1,
sf_surface_physics                  = 2,     2,     2,     2,
bl_pbl_physics                      = 1,     1,     1,     1,
bldt                                = 0,     0,     0,     0,
cu_physics                          = 5,     5,     0,     0,
cudt                                = 0,     0,     0,     0,
isfflx                              = 1,
ifsnow                              = 1,
icloud                              = 1,
num_soil_layers                     = 4,   
sf_urban_physics                    = 1,     1,     1,     1,
num_land_cat                        = 24,
/

&fdda
/

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

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

&chem
kemit                               = 1,
chem_opt                            = 112,        112,   112,   112,
photdt                              = 30,         30,    30,    30,
chemdt                              = 1.5,        1.5,   1.5,   1.5,
io_style_emissions                  = 1,
emiss_opt                           = 8,          8,     8,     8,
emiss_inpt_opt                      = 111,        111,   111,   111,
chem_in_opt                         = 1,          1,     1,     1,
phot_opt                            = 1,          1,     1,     1,
gas_drydep_opt                      = 1,          1,     1,     1,
aer_drydep_opt                      = 1,          1,     1,     1,
dust_opt                            = 3,            
depo_fact                           = 0.25,
bio_emiss_opt                       = 0,          0,    0,    0,
gas_bc_opt                          = 1,          1,    1,    1,
gas_ic_opt                          = 1,          1,    1,    1,
aer_bc_opt                          = 1,          1,    1,    1,
aer_ic_opt                          = 1,          1,    1,    1,
gaschem_onoff                       = 1,          1,    1,    1,
aerchem_onoff                       = 1,          1,    1,    1,
wetscav_onoff                       = 0,          0,    0,    0,
cldchem_onoff                       = 0,          0,    0,    0,
vertmix_onoff                       = 1,          1,    1,    1,
chem_conv_tr                        = 0,          0,    0,    0,
conv_tr_wetscav                     = 1,          1,    1,    1,
conv_tr_aqchem                      = 0,          0,    0,    0,
biomass_burn_opt                    = 0,          0,    0,    0,
have_bcs_chem                       = .true., .true., .true., .true.,
have_bcs_upper                      = .true., .true., .true., .true.,
aer_op_opt                          = 1,
opt_pars_out                        = 0,
aer_ra_feedback                     = 1,
chemdiag                            = 0,          0,      0,     0,
fixed_ubc_press                     = 1.,         1.,     1.,    1.,
fixed_ubc_inname                    = 'ubvals_b40.20th.track1_1996-2005.nc'
/



&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

Lucia007

请问博主解决了吗

刘书华

请问博主解决了吗

zy18981925350

请问博主解决了吗

flower.

shiyu 发表于 2020-1-11 15:51
这是我的namelist.input，非常感谢大家

&time_control

p_top_requested=1000试试,另外，time_step 改小一点