wrf运行问题（Water is zero in liquid phase）

affechen

photolysis_driver: called for domain            1
Water is zero in liquid phase
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:      76
Stopping in ASTEM_flux_wet
-------------------------------------------
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE:  <stdin>  LINE:      76
Stopping in ASTEM_flux_wet
-------------------------------------------
Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_OVERFLOW_FLAG IEEE_UNDERFLOW_FLAG IEEE_DENORMAL


Stopping in ASTEM_flux_wet这个问题把aerchem_onoff = 1关成0就能往下，但是似乎就不能模拟气溶胶了，该如何解决这个问题呢？

affechen

还是没能解决是不是排放清单的问题呢。。。。

affechen

我排放清单好像缺HONO数据。。。。

astiny

第六感告诉我应该不是排放源，而是参数设置的问题。你最好发一下你的namelist。

affechen

astiny 发表于 2018-4-12 13:18
第六感告诉我应该不是排放源，而是参数设置的问题。你最好发一下你的namelist。

麻烦您了~

&time_control            
run_days                 = 07,
run_hours                = 12,
run_minutes              = 0,
run_seconds              = 0,
start_year               = 2017,
start_month              = 11,
start_day                = 1,
start_hour               = 00,
start_minute             = 00,
start_second             = 00,
end_year                 = 2017,
end_month                = 11,
end_day                  = 08,
end_hour                 = 12,
end_minute               = 00,
end_second               = 00,
interval_seconds         = 7200,
input_from_file          = .false.,
history_interval         = 180,
frames_per_outfile       = 1000,
restart                  = .false.,
restart_interval         = 5000,
io_form_history          = 2,
io_form_restart          = 2,
io_form_input            = 2,
io_form_boundary         = 2,
debug_level              = 0,
/

&domains                 
eta_levels               = 1.000, 0.9937, 0.9874, 0.9811, 0.9748,
                0.9686, 0.9617, 0.9539, 0.9449, 0.9343,
                0.9217, 0.9063, 0.8875, 0.8647, 0.8369,
                0.8037, 0.7647, 0.7199, 0.6699, 0.6155,
                0.558, 0.4988, 0.439, 0.38, 0.3225,
                0.2674, 0.2149, 0.1656, 0.1194, 0.0764,
                0.0366, 0.000,
time_step                = 60,
time_step_fract_num      = 0,
time_step_fract_den      = 1,
max_dom                  = 1,
e_we                     = 160,
e_sn                     = 160,
e_vert                   = 32,
p_top_requested          = 21826.2,
num_metgrid_levels       = 32,
num_metgrid_soil_levels  = 4,
dx                       = 6000,
dy                       = 6000,
grid_id                  = 1,
parent_id                = 1,
i_parent_start           = 1,
j_parent_start           = 1,
parent_grid_ratio        = 1,
parent_time_step_ratio   = 1,
feedback                 = 1,
smooth_option            = 0,
/

&physics                 
mp_physics               = 10,
progn                    = 1,
ra_lw_physics            = 4,
ra_sw_physics            = 4,
radt                     = 30,
sf_sfclay_physics        = 5,
sf_surface_physics       = 2,
bl_pbl_physics           = 5,
bldt                     = 0,
cu_physics               = 5,
cudt                     = 0,
cugd_avedx               = 1,
isfflx                   = 1,
ifsnow                   = 1,
icloud                   = 1,
surface_input_source     = 1,
num_soil_layers          = 4,
sf_urban_physics         = 1,
mp_zero_out              = 2,
mp_zero_out_thresh       = 1.e-8,
cu_rad_feedback          = .true.,
cu_diag                  = 1,
slope_rad                = 0,
topo_shading             = 0,
maxiens                  = 1,
maxens                   = 3,
maxens2                  = 3,
maxens3                  = 16,
ensdim                   = 144,
/

&chem                    
kemit                    = 8,
chem_opt                 = 201,
bioemdt                  = 30,
photdt                   = 30.,
chemdt                   = 6.0,
io_style_emissions       = 1,
emiss_inpt_opt           = 102,
emiss_opt                = 10,
chem_in_opt              = 1,
phot_opt                 = 3,
gas_drydep_opt           = 1,
aer_drydep_opt           = 1,
depo_fact                = 0.25,
aer_op_opt               = 1,
bio_emiss_opt            = 3,
gas_bc_opt               = 1,
gas_ic_opt               = 1,
aer_bc_opt               = 1,
aer_ic_opt               = 1,
gaschem_onoff            = 1,
aerchem_onoff            = 1,
wetscav_onoff            = 1,
cldchem_onoff            = 0,
vertmix_onoff            = 1,
conv_tr_wetscav          = 1,
conv_tr_aqchem           = 0,
seas_opt                 = 0,
dust_opt                 = 0,
dmsemis_opt              = 1,
biomass_burn_opt         = 2,
plumerisefire_frq        = 30,
scale_fire_emiss         = .true.,
aer_ra_feedback          = 1,
n2o5_hetchem             = 1,
ne_area                  = 357,
opt_pars_out             = 1,
have_bcs_chem            = .true.,
chem_conv_tr             = 1,
have_bcs_upper           = .false.,
/

&fdda                    
/

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

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

&grib2                  
/

&namelist_quilt         
nio_tasks_per_group      = 0,
nio_groups               = 1,
/

affechen

astiny 发表于 2018-4-12 13:18
第六感告诉我应该不是排放源，而是参数设置的问题。你最好发一下你的namelist。

我想用MOZART-MOSAIC-4bin这个机制做chem_opt  只要关那个液相chem就好，神奇了

affechen

astiny 发表于 2018-4-12 13:18
第六感告诉我应该不是排放源，而是参数设置的问题。你最好发一下你的namelist。

新发现：我的wrfinput_d01倘若使用mozbc添加了mozart4化学初始边界条件，当aerchem_onoff = 1,就会出上面的错误，当aerchem_onoff = 0，就不出错，当我使用real.exe生成wrfinput_d01不用mozbc添加mozart4化学初始边界条件，当aerchem_onoff = 1就不会出错，但是似乎计算所得数据和我的监测数据一点都不相关。。。