WRF

生生不息

我利用WRFCHEM模拟，第一次关闭chem运行纯气象模式，polar设置为true，运行real.exe,一直报错，请问这是为什么呢？该如何修改namelist.input? 以下是我的namelist.input和rsl.error.0000
&time_control
start_year                    = 2003,   2016,   2016,
start_month                   = 01,     08,     08,
start_day                     = 01,     26,     26,
! start_month                   = 09,     09,     09,
! start_day                     = 07,     07,     07,
start_hour                    = 00,     00,     00,
start_minute                  = 00,     00,     00,
start_second                  = 00,     00,     00,
end_year                      = 2003,   2016,   2016,
end_month                     = 01,     10,     10,
end_day                       = 03,     03,     03,
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,    60,
frames_per_outfile            = 9999,  9999,  9999,
frames_per_auxinput5          = 1, 1, 1,
! restart                       = .true.,
restart                       = .false.,
restart_interval              = 4320,
io_form_history               = 2,
io_form_restart               = 2,
io_form_input                 = 2,
io_form_boundary              = 2,
debug_level                   = 0,
auxinput6_inname              = 'wrfbiochemi_d<domain>',
! auxinput7_inname              = 'wrffirechemi_d<domain>_<date>',
auxinput5_interval_m          = 60, 60, 60,
! auxinput7_interval_m          = 60, 60, 60,
io_form_auxinput2             = 2,
io_form_auxinput4             = 2,
auxinput4_inname              = 'wrflowinp_d<domain>',
auxinput4_interval            = 360,360,360,
io_form_auxinput5             = 2,
io_form_auxinput6             = 2,
! io_form_auxinput7             = 2,
! io_form_auxinput8             = 0,
override_restart_timers       = .true.,
write_hist_at_0h_rst          = .true.,
frames_per_auxinput5          = 1, 1, 1,
frames_per_auxinput6          = 1, 1, 1,
! frames_per_auxinput7          = 1, 1, 1,
history_outname               =
auxinput1_inname              =
force_use_old_data            = T,
! iofields_filename             = 'my_file_d01.txt','my_file_d02.txt','my_file_d03.txt',
! ignore_iofields_warning       = .true.,
/


&domains
time_step                     = 120,
time_step_fract_num           = 0,
time_step_fract_den           = 1,
max_dom                       = 1,
e_we                          = 144, 181, 151,
e_sn                          = 90, 181, 151,
e_vert                        = 50,  50,  50,
num_metgrid_levels            = 27,
eta_levels   = 1.000, 0.990, 0.987, 0.979, 0.97,
                0.96, 0.949, 0.937, 0.924, 0.909,
                0.892, 0.873, 0.851, 0.826, 0.798,
                0.768, 0.736, 0.702, 0.666, 0.629,
                0.5915, 0.5536, 0.5153, 0.4773, 0.44,
                0.404, 0.3695, 0.3375, 0.3085, 0.2845,
                0.2645, 0.2465, 0.2305, 0.2165, 0.2035,
                0.1915, 0.1792, 0.1667, 0.1539, 0.1407,
                0.1272, 0.1134, 0.0995, 0.0855, 0.0713,
                0.0571, 0.0429, 0.0287, 0.0145, 0.000,
p_top_requested               = 5000,
num_metgrid_soil_levels       = 2,
dx                            = 277943.7, 9000, 3000,
dy                            = 222355.0, 9000, 3000,
grid_id                       = 1,   2,   3,
parent_id                     = 1,   1,   2,
i_parent_start                = 1,   71,  43,
j_parent_start                = 1,   35,  40,
parent_grid_ratio             = 1,   3,   3,
parent_time_step_ratio        = 1,   3,   3,
feedback                      = 1,
smooth_option                 = 2,
smooth_cg_topo                = .true.,
/


&physics
mp_physics                    = 2,   2,   2,
ra_lw_physics                 = 4,   4,   4,
ra_sw_physics                 = 4,   4,   4,
radt                          = 30,  30,  30,
sf_sfclay_physics             = 1,   1,   1,
sf_surface_physics            = 2,   2,   2,
bl_pbl_physics                = 1,   1,   1,
bldt                          = 0,   0,   0,
cu_physics                    = 5,   5,   5,
cu_rad_feedback               = .true., .true., .true.,
cudt                          = 0,   0,   0,
isfflx                        = 1,
ifsnow                        = 1,
icloud                        = 1,
surface_input_source          = 1,
num_soil_layers               = 4,
! mp_zero_out                   = 2,
! mp_zero_out_thresh            = 1.e-12,
sf_urban_physics              = 1,   1,   1,
tmn_update                    = 1,
maxiens                       = 1,
maxens                        = 3,
maxens2                       = 3,
maxens3                       = 16,
ensdim                        = 144,
sst_update                    = 1,
! usemonalb                     = .true.,
progn                         = 0,   0,   0,
! cu_diag                       = 1,   1,   1,
num_soil_layers               = 4,
num_land_cat                  = 21,
ghg_input                     = 0,
/

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


&dynamics
! rk_ord                        = 3,
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.01,   0.01,   0.01
khdif                         = 0,      0,      0,
kvdif                         = 0,      0,      0,
non_hydrostatic               = .true., .true., .true.,
moist_adv_opt                 = 2,      2,      2,
scalar_adv_opt                = 2,      2,      2,
chem_adv_opt                  = 2,      2,      2,
tke_adv_opt                   = 2,      2,      2,
! time_step_sound               = 4,      4,      4,
! 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,
! epssm                         = 0.5,    1,    0.1,    0.5,
/


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


&grib2
/

&namelist_quilt
nio_tasks_per_group           = 0,
nio_groups                    = 1,
/

&chem
kemit                         = 11,
ne_area                       = 100,
! chem_opt                      = 106,  106,  106,
chem_opt                      = 0,0,0
bioemdt                       = 30,   30,   30,
photdt                        = 30,   30,   30,
chemdt                        = 4,    4,    4,
io_style_emissions            = 1,
emiss_inpt_opt                = 1,    1,    1,
emiss_opt                     = 3,    3,    3,
chem_in_opt                   = 0,    0,    0,
phot_opt                      = 1,    1,    1,
gas_drydep_opt                = 1,    1,    1,
aer_drydep_opt                = 1,    1,    1,
bio_emiss_opt                 = 3,    3,    3,
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                 = 0,    0,    0,
cldchem_onoff                 = 0,    0,    0,
vertmix_onoff                 = 1,    1,    1,
chem_conv_tr                  = 1,    1,    1,
seas_opt                      = 2,   
dust_opt                      = 3,
biomass_burn_opt              = 2,    2,    2,
plumerisefire_frq             = 30,   30,   30,
have_bcs_chem                 = .false., .false., .false.,
aer_ra_feedback               = 0,    0,    0,
aer_op_opt                    = 3,3,3,1,    1,    1,
opt_pars_out                  = 0,
diagnostic_chem               = 0,    0,    0,
chemdiag                      = 0,    0,    1,
mozart_ph_diag                = 1,
/

这是rsl.error.0000文件
taskid: 0 hostname: login_a08
module_io_quilt_old.F        2931 T
Ntasks in X            1 , ntasks in Y            1
  Domain # 1: dx = ********* m
REAL_EM V4.4.2 PREPROCESSOR
git commit 6233639c599119e76fca17dba9ea211af53a0ba9
*************************************
Parent domain
ids,ide,jds,jde            1         144           1          90
ims,ime,jms,jme           -4         149          -4          95
ips,ipe,jps,jpe            1         144           1          90
*************************************
DYNAMICS OPTION: Eulerian Mass Coordinate
   alloc_space_field: domain            1 ,             2618776560  bytes allocated
d01 2003-01-01_00:00:00  Yes, this special data is acceptable to use: OUTPUT FROM METGRID V4.4
d01 2003-01-01_00:00:00  Input data is acceptable to use:
/met_em.d01.2003-01-01_00:00:00.nc
metgrid input_wrf.F first_date_input = 2003-01-01_00:00:00
metgrid input_wrf.F first_date_nml = 2003-01-01_00:00:00
d01 2003-01-01_00:00:00 Timing for input          0 s.
d01 2003-01-01_00:00:00          flag_soil_layers read from met_em file is  1
  *** Error in boundary condition specification
  boundary conditions at xs            1
  boundary conditions at xe            1
  boundary conditions at ys            2
  boundary conditions at ye            2
  boundary conditions logicals are
  periodic_x    F
  periodic_y    F
  symmetric_xs  F
  symmetric_xe  F
  symmetric_ys  F
  symmetric_ye  F
  open_xs       F
  open_xe       F
  open_ys       F
  open_ye       F
  polar         T
  nested        F
  specified     T
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:     165
*** Error in boundary condition specification
-------------------------------------------
application called MPI_Abort(MPI_COMM_WORLD, 1) - process 0
[unset]: write_line error; fd=-1 buf=:cmd=abort exitcode=1
:
system msg for write_line failure : Bad file descriptor

李知恩

家园好像以前有人问过这个问题：https://bbs.06climate.com/forum.php?mod=viewthread&tid=46030

还有我看楼主的num_metgrid_levels            = 27,
和eta_levels的层数似乎对不上？ 这么设置会报错吗？