A bunch of wrf.exe mistakes in two WRF-Chem Versions

[Ren JY]

Dear all,
I'm a freshman of WRF model and last month I just solved a problem about WPS running(link:http://forum.mmm.ucar.edu/phpBB3/viewtopic.php?f=32&t=5365). Now I've got new problems.
At the very beginning, I ran real.exe successfully with my met* files metgridded from ERA-I data under WRF-Chem Ver3.5, but wrf.exe stopped in 30 seconds with one mistake called "Segmentation fault-Invalid memory reference". According to the bug report, I found after setting ra_lw_physics&ra_sw_physics = 0, wrf.exe could run with no pause.
However I need radiation, so I downloaded and compiled WRF-Chem Ver3.9.1 to see if it works. Then it gave me new error messages about auxinput5. I added "io_form_auxinput5 = 2" in &time control, ran and still fall. The third error is:

d01 2018-07-01_00:00:00 Input data processed for aux input 4 for domain 1
3 input_wrf: wrf_get_next_time current_date: 2018-07-01_00:00:00 Status = -11
d01 2018-07-01_00:00:00 ---- ERROR: Could not find matching time in input file
NOTE: 1 namelist vs input data inconsistencies found.
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 1148
NOTE: Please check and reset these options

Every time I searched the keywords about my mistake at this forum and read all the related topics, but it seems this time searching doesn't help. All my attempts running wrf.exe will be stopped in about 30 sec, leaving a small-sized wrfout file. During steps before executing wrf.exe, I didn't receive error reports. Actually I turned off all the other options in my namelist (no chemistry either). Am I wrong in namelist settings? Do I have to redownload the reanalysis data from rda.ucar? Or maybe I need to complie a pure WRF at the latest version? Any suggestion is urgely welcomed. Thanks!
Grateful,
Ren

 

[douglowe]

You shouldn't need to use WRF-Chem to get the long and short wave radiation options to work, so probably you have a WRF, not WRF-Chem, issue.

I've had a quick look through your namelist.input file - the one thing that jumps out at me is this:
P_TOP_REQUESTED = 1000000.000 ,

This is the pressure altitude at which you want the top of the domain to be. The default value is 5000 Pa. Your value of 1,000,000 Pa would be more than atmospheric pressure at sea level, so it's likely that your model levels are almost flat on the ground surface. Try a value of 5000 again, and see if that solves your run-time problems.

 

[Ren JY]

Hi Douglowe,
Thanks for your suggestion! I asked for another namelist.input from my classmate to conduct the wrf running and it worked. Since there is quite much difference including the p_top_requested = 5,000, I'm not really sure if that is the key to my problem. Next I'll try to do control experiments to find the exact reason. Here is my new namelist.input.

Ren

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

&domains
time_step = 72,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
s_we = 1, 1, 1,
e_we = 432, 193, 133,
s_sn = 1, 1, 1,
e_sn = 372, 106, 97,
s_vert = 1, 1, 1,
e_vert = 38, 33, 33,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 12000, 10000, 3333.33,
dy = 12000, 10000, 3333.33,
grid_id = 1, 2, 3,
parent_id = 0, 1, 2,
i_parent_start = 1, 31, 30,
j_parent_start = 1, 17, 30,
parent_grid_ratio = 1, 3, 3,
parent_time_step_ratio = 1, 3, 3,
feedback = 1,
smooth_option = 2

/

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

&fdda
/

&dynamics
w_damping = 0,
diff_6th_opt = 0, 0, 0,
diff_6th_factor = 0.12, 0.12, 0.12,
diff_opt = 1,
km_opt = 4,
base_temp = 290.
damp_opt = 3,
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,
gwd_opt = 0,
/

&bdy_control
spec_bdy_width = 12,
specified = .true.
/

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/