Hi! I am a student starting to learn wrf, and I'm running an example to learn the whole process of WRF.
However I met some problems that I do not know how to deal with. Could you please give me some help?
I have successfully preprocessed the data through WPS, and added chem、bio and anthro emission parts. When I tried to run wrf.exe, the error appeared as the following:
----------------- ERROR -------------------
namelist : NUM_LAND_CAT = 21
input files : NUM_LAND_CAT = 20 (from geogrid selections).
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 706
Mismatch between namelist and wrf input files for dimension NUM_LAND_CAT
-------------------------------------------
Warning: ieee_invalid is signaling
Warning: ieee_inexact is signaling
FORTRAN STOP
but when I change the namelist.input into: num_land_cat = 20 , it appears error like
----------------- ERROR -------------------
namelist : NUM_LAND_CAT = 20
input files : NUM_LAND_CAT = 21 (from geogrid selections).
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 706
Mismatch between namelist and wrf input files for dimension NUM_LAND_CAT
It seems that they always mismatch no matter I change , I don't know why
The followings are the namelist.wps and namelist.input I used.
namelist.wps
share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2015-10-01_12:00:00','2006-08-16_12:00:00',
end_date = '2015-10-08_18:00:00','2006-08-16_12:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 31,
j_parent_start = 1, 17,
e_we = 301, 112,
e_sn = 301, 97,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the MAXSNOALB and ALBEDO12M
! fields have changed in WPS v4.0. These fields are now interpolated
! from MODIS-based datasets.
!
! To match the output given by the default namelist.wps in WPS v3.9.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','default',
dx = 6000,
dy = 6000,
map_proj = 'lambert',
ref_lat = 38.0,
ref_lon = 116.0,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = 116.0,
geog_data_path = '/home/yxlv/geogv4.0/WPS_GEOG/'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
namelist.input
time_control
run_days= 8,
run_hours= 6,
run_minutes= 0,
run_seconds = 0,
start_year= 2015,
start_month= 10,
start_day = 01,
start_hour = 12,
start_minute = 00,
start_second= 00,
end_year = 2015,
end_month= 10,
end_day = 08,
end_hour = 18,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file= .true.,
history_interval = 60,
frames_per_outfile = 2,
restart= .false.,
restart_interval= 14400,
auxinput5_interval_m= 60,
auxinput6_interval_m = 86400,
auxinput6_inname = 'wrfbiochemi_d<domain>',
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput5 = 2
io_form_auxinput6 = 2
debug_level = 0
auxinput11_interval_s = 120, 120, 120, 120,
auxinput11_end_h = 888, 30, 30, 30,
/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 301,
e_sn = 301,
e_vert = 30,
dx = 6000,
dy = 6000,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
grid_id = 1,
parent_id = 0,
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 = 6,
progn = 0,
naer = 1e9
ra_lw_physics = 5,
ra_sw_physics = 5,
radt = 6,
sf_sfclay_physics = 2,
sf_surface_physics = 2,
bl_pbl_physics = 2,
bldt = 1,
cu_physics = 0,
cudt = 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat =21
sf_urban_physics = 0,
num_urban_layers = 1200,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false.,
/
&fdda
obs_nudge_opt= 1,1,0,0,0
max_obs= 150000,
fdda_start = 0., 0., 0., 0., 0.
fdda_end = 99999., 99999., 99999., 99999., 99999.
obs_nudge_wind = 1,1,1,0,0
obs_coef_wind = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp = 0,0,0,0,0
obs_coef_temp = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois = 1,1,1,1,1
obs_coef_mois = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy = 1080.,360.,90.,30.,180
obs_rinsig = 0.1,
obs_twindo = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi = 10,
obs_ionf = 2, 2, 2, 2, 2,
obs_idynin = 0,
obs_dtramp = 40.,
obs_prt_freq = 10, 10, 10, 10, 10,
obs_prt_max = 10
obs_ipf_errob = .true.
obs_ipf_nudob = .true.
obs_ipf_in4dob = .true.
obs_ipf_init = .true.
/
&dynamics
w_damping = 1,
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,
tracer_adv_opt = 1,
tke_adv_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,
nested = .false.,
&grib2
/
&chem
kemit = 2,
chem_opt = 1,
bioemdt = 30,
photdt = 15,
chemdt = 2.0,
frames_per_emissfile = 840,
emiss_inpt_opt = 2,
emiss_opt = 1
bb_emiss_opt = 1,
io_style_emissions = 2
chem_in_opt = 0,
phot_opt = 1,
drydep_opt = 1,
bio_emiss_opt = 3,
aer_ra_feedback = 1,
so2_haze_react = 10,
fe3_ion_factor = 1.0,
mn2_ion_factor = 1.0,
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I would appreciate your help deeply.
Regards
However I met some problems that I do not know how to deal with. Could you please give me some help?
I have successfully preprocessed the data through WPS, and added chem、bio and anthro emission parts. When I tried to run wrf.exe, the error appeared as the following:
----------------- ERROR -------------------
namelist : NUM_LAND_CAT = 21
input files : NUM_LAND_CAT = 20 (from geogrid selections).
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 706
Mismatch between namelist and wrf input files for dimension NUM_LAND_CAT
-------------------------------------------
Warning: ieee_invalid is signaling
Warning: ieee_inexact is signaling
FORTRAN STOP
but when I change the namelist.input into: num_land_cat = 20 , it appears error like
----------------- ERROR -------------------
namelist : NUM_LAND_CAT = 20
input files : NUM_LAND_CAT = 21 (from geogrid selections).
-------------- FATAL CALLED ---------------
FATAL CALLED FROM FILE: <stdin> LINE: 706
Mismatch between namelist and wrf input files for dimension NUM_LAND_CAT
It seems that they always mismatch no matter I change , I don't know why
The followings are the namelist.wps and namelist.input I used.
namelist.wps
share
wrf_core = 'ARW',
max_dom = 2,
start_date = '2015-10-01_12:00:00','2006-08-16_12:00:00',
end_date = '2015-10-08_18:00:00','2006-08-16_12:00:00',
interval_seconds = 21600
io_form_geogrid = 2,
/
&geogrid
parent_id = 1, 1,
parent_grid_ratio = 1, 3,
i_parent_start = 1, 31,
j_parent_start = 1, 17,
e_we = 301, 112,
e_sn = 301, 97,
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
! The default datasets used to produce the MAXSNOALB and ALBEDO12M
! fields have changed in WPS v4.0. These fields are now interpolated
! from MODIS-based datasets.
!
! To match the output given by the default namelist.wps in WPS v3.9.1,
! the following setting for geog_data_res may be used:
!
! geog_data_res = 'maxsnowalb_ncep+albedo_ncep+default', 'maxsnowalb_ncep+albedo_ncep+default',
!
!!!!!!!!!!!!!!!!!!!!!!!!!!!! IMPORTANT NOTE !!!!!!!!!!!!!!!!!!!!!!!!!!!!
!
geog_data_res = 'default','default',
dx = 6000,
dy = 6000,
map_proj = 'lambert',
ref_lat = 38.0,
ref_lon = 116.0,
truelat1 = 30.0,
truelat2 = 60.0,
stand_lon = 116.0,
geog_data_path = '/home/yxlv/geogv4.0/WPS_GEOG/'
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/
namelist.input
time_control
run_days= 8,
run_hours= 6,
run_minutes= 0,
run_seconds = 0,
start_year= 2015,
start_month= 10,
start_day = 01,
start_hour = 12,
start_minute = 00,
start_second= 00,
end_year = 2015,
end_month= 10,
end_day = 08,
end_hour = 18,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file= .true.,
history_interval = 60,
frames_per_outfile = 2,
restart= .false.,
restart_interval= 14400,
auxinput5_interval_m= 60,
auxinput6_interval_m = 86400,
auxinput6_inname = 'wrfbiochemi_d<domain>',
io_form_history = 2
io_form_restart = 2
io_form_input = 2
io_form_boundary = 2
io_form_auxinput5 = 2
io_form_auxinput6 = 2
debug_level = 0
auxinput11_interval_s = 120, 120, 120, 120,
auxinput11_end_h = 888, 30, 30, 30,
/
&domains
time_step = 30,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 301,
e_sn = 301,
e_vert = 30,
dx = 6000,
dy = 6000,
p_top_requested = 5000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
grid_id = 1,
parent_id = 0,
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 = 6,
progn = 0,
naer = 1e9
ra_lw_physics = 5,
ra_sw_physics = 5,
radt = 6,
sf_sfclay_physics = 2,
sf_surface_physics = 2,
bl_pbl_physics = 2,
bldt = 1,
cu_physics = 0,
cudt = 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
num_land_cat =21
sf_urban_physics = 0,
num_urban_layers = 1200,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
cu_rad_feedback = .false.,
/
&fdda
obs_nudge_opt= 1,1,0,0,0
max_obs= 150000,
fdda_start = 0., 0., 0., 0., 0.
fdda_end = 99999., 99999., 99999., 99999., 99999.
obs_nudge_wind = 1,1,1,0,0
obs_coef_wind = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_temp = 0,0,0,0,0
obs_coef_temp = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_nudge_mois = 1,1,1,1,1
obs_coef_mois = 6.E-4,6.E-4,6.E-4,6.E-4,6.E-4
obs_rinxy = 1080.,360.,90.,30.,180
obs_rinsig = 0.1,
obs_twindo = 0.6666667,0.6666667,0.6666667,0.6666667,0.6666667,
obs_npfi = 10,
obs_ionf = 2, 2, 2, 2, 2,
obs_idynin = 0,
obs_dtramp = 40.,
obs_prt_freq = 10, 10, 10, 10, 10,
obs_prt_max = 10
obs_ipf_errob = .true.
obs_ipf_nudob = .true.
obs_ipf_in4dob = .true.
obs_ipf_init = .true.
/
&dynamics
w_damping = 1,
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,
tracer_adv_opt = 1,
tke_adv_opt = 1,
/
&bdy_control
spec_bdy_width = 5,
spec_zone = 1,
relax_zone = 4,
specified = .true.,
nested = .false.,
&grib2
/
&chem
kemit = 2,
chem_opt = 1,
bioemdt = 30,
photdt = 15,
chemdt = 2.0,
frames_per_emissfile = 840,
emiss_inpt_opt = 2,
emiss_opt = 1
bb_emiss_opt = 1,
io_style_emissions = 2
chem_in_opt = 0,
phot_opt = 1,
drydep_opt = 1,
bio_emiss_opt = 3,
aer_ra_feedback = 1,
so2_haze_react = 10,
fe3_ion_factor = 1.0,
mn2_ion_factor = 1.0,
&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/
I would appreciate your help deeply.
Regards