ERROR: Mismatch between namelist and global attribute NUM_METGRID_SOIL_LEVELS

[Mengyao]

Hi! I am a student using WRF(v4.1). 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(v4.1). When I tried to run real.exe, the error appeared as the following:


----------------- ERROR -------------------
namelist : num_metgrid_soil_levels = 4
input files : NUM_METGRID_SOIL_LEVELS = 2 (from met_em files).
d01 2000-08-01_00:00:00 ---- ERROR: Mismatch between namelist and global attribute NUM_METGRID_SOIL_LEVELS
NOTE: 1 namelist vs input data inconsistencies found.


The followings are the namelist.wps and namelist.input I used.

----------------------------------------------------------------------------------------------------------------------------------------------------------------

&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2000-08-01_00:00:00',
end_date = '2000-08-01_18:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
/

&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 701,
e_sn = 401,
dx = 0.1,
dy = 0.1,
map_proj = 'lat-lon',
pole_lat = 90
pole_lon = 0
stand_lon = 105
ref_lat = 35,
ref_lon = 105,
geog_data_res = 'modis_5m_lake+soilgrids+5m',
geog_data_path = '/work/06956/mguo/WPS_static_data/WPS_GEOG'
/

&ungrib
out_format = 'WPS',
prefix = 'FILE',
/

&metgrid
fg_name = 'FILE'
io_form_metgrid = 2,
/

--------------------------------------------------------------------------------------------------------------------------------------------------------------

&time_control
run_days = 0,
run_hours = 18,
run_minutes = 0,
run_seconds = 0,
start_year = 2000,
start_month = 08,
start_day = 01,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2000,
end_month = 08,
end_day = 01,
end_hour = 18,
end_minute = 00,
end_second = 00,
interval_seconds = 21600
input_from_file = .true.,.true.,.true.,
history_interval = 360, 60, 60,
frames_per_outfile = 1000, 1000, 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
time_step = 66,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 701,
e_sn = 401,
e_vert = 40
p_top_requested = 6000,
num_metgrid_levels = 27,
num_metgrid_soil_levels = 4,
dx = 11117.75,
dy = 11117.75,
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,
sfcp_to_sfcp =.true.,
/

&physics
mp_physics = 3,
ra_lw_physics = 1,
ra_sw_physics = 1,
radt = 11,
sf_sfclay_physics = 1,
sf_surface_physics = 1,
bl_pbl_physics = 1,
bldt = 0,
cu_physics = 1,
cudt = 5,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
sf_urban_physics = 0,
num_land_cat = 21,
/

&fdda
/

&dynamics
w_damping = 0,
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 = 6000., 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,
/

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

&grib2
/

&namelist_quilt
nio_tasks_per_group = 0,
nio_groups = 1,
/

---------------------------------------------------------------------------------------------------------------------------------------------------------------


I would appreciate your help deeply.


Regards,

Mengyao Guo

 

[Ming Chen]

Mengyao,
Please see the line below to &domains of your namelist.input)
NUM_METGRID_SOIL_LEVELS = 2

The rerun real.exe.

 

[Mengyao]

Thanks for your help Ming!
I actully want to get the soil data with 4 layers. So I was wondering if there is a way to fix it.
I have used ncdump met_em.d01.2000-08-01_00:00:00.nc |grep num_st and ncdump met_em.d01.2000-08-01_00:00:00.nc |grep num_sm for troubleshooting. But results showed num_st_layers = 2 and num_sm_layers = 2 .

By the way, I use FNL data as input data. And I have linked the Vtable.GFS.new to Vtable.

Are there any other possible methods to fix it?

Thank you!

 

[Ming Chen]

Mengyao,
FNL data only has two layers of soil moisture and temperature. I don't think you can change that.
If possible, can you switch to GFS? GFS provides 4-layer of soil information.
Vtable.GFS works for both FNL and GFS data.

 

[Mengyao]

Thank you, Ming! I forgot to check the message these days and see your suggestion just now.

I tried to change the GRIB1 FNL data to GRIB2 yesterday, and it finally worked! So maybe the GRIB2 data from 2007.12.06 could provide the 4-layer soil information as well?

 

[Ming Chen]

That is possible, although I am not 100% sure. Thanks for the update.