Hi,
I am trying to run WRF with following namelist options. But when I used aux_hist23 option with custom pressure levels for wrfpress_d01_<domain><date>. The output has lot of missing values where as, wrfout* file has proper values for u and v variables. when I try to extract them for particular pressure levels, some of the values go missing in wrfpress. Is this behaviour normal??. I have compared the pressure levels between metemfiles and my custom levels in &diags section in namelist. everything seems correct. Please help.
NAMELIST WPS
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2009-12-31_12:00:00',
end_date = '2009-12-31_18:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = './',
debug_level = 0,
/
&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 127,
e_sn = 141,
geog_data_res = '30s',
dx = 30000,
dy = 30000,
map_proj = 'mercator',
ref_lat = 21.474,
ref_lon = 79.796,
truelat1 = 21.474,
truelat2 = 0,
stand_lon = 79.796,
geog_data_path = '../WPS_GEOG',
opt_geogrid_tbl_path = './',
ref_x = 63.5,
ref_y = 70.5,
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/DATA/user4/LA_project/',
opt_metgrid_tbl_path = './',
/
&mod_levs
press_pa = 201300 , 200100 , 100000 ,
95000 , 90000 ,
85000 , 80000 ,
75000 , 70000 ,
65000 , 60000 ,
55000 , 50000 ,
45000 , 40000 ,
35000 , 30000 ,
25000 , 20000 ,
15000 , 10000 ,
5000 , 1000
/
&domain_wizard
grib_data_path = '/media/user/Research/ERAint',
grib_vtable = 'Vtable.ERA-interim.pl',
dwiz_name =eraintfine
dwiz_desc =
dwiz_user_rect_x1 =2710
dwiz_user_rect_y1 =545
dwiz_user_rect_x2 =3137
dwiz_user_rect_y2 =997
dwiz_show_political =true
dwiz_center_over_gmt =true
dwiz_latlon_space_in_deg =10
dwiz_latlon_linecolor =-8355712
dwiz_map_scale_pct =25.0
dwiz_map_vert_scrollbar_pos =463
dwiz_map_horiz_scrollbar_pos =2613
dwiz_gridpt_dist_km =30.0
dwiz_mpi_command =null
dwiz_tcvitals =null
dwiz_bigmap =Y
/
namelist.input
&time_control
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2008,
start_month = 01,
start_day = 01,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2008,
end_month = 01,
end_day = 01,
end_hour = 18,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file = .true.,
history_interval = 360,
frames_per_outfile = 1000,
restart = .false.,
restart_interval = 43200,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
nocolons = .false.,
io_form_auxinput4 = 2 !? 11 on cray...parallel netcdf
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 360
output_diagnostics = 1,
auxhist3_outname = 'wrfxtrm_d<domain>_<date>',
io_form_auxhist3 = 2,
auxhist3_interval = 1440,
frames_per_auxhist3 =1,
auxhist23_outname ='wrfpress_d<domain>_<date>'
io_form_auxhist23 = 2, !? 11 on cray...parallel netcdf
auxhist23_interval = 360,
frames_per_auxhist23 = 4,
auxhist7_outname ='CTPframework_d<domain>_<date>'
io_form_auxhist7 = 2, !? 11 on cray...parallel netcdf
auxhist7_interval = 180,
frames_per_auxhist7 = 8,
iofields_filename = "additional_fields.txt"
diag_print = 1,
/
&diags
p_lev_diags = 1,
num_press_levels = 19,
press_levels = 95000, 92500, 90000, 87500, 85000, 82500, 80000, 77500, 75000, 72500, 70000, 67500, 65000, 50000, 40000, 30000, 20000, 10000, 5000
use_tot_or_hyd_p = 2,
/
&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 127,
e_sn = 141,
e_vert = 40,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 30000,
dy = 30000,
grid_id = 1,
parent_id = 1,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 0,
smooth_option = 0,
sfcp_to_sfcp = .true.,
/
&physics
sst_update = 1,
mp_physics = 6, ! WSM5 scheme
ra_lw_physics = 1, ! RRTM scheme
ra_sw_physics = 1, ! Dudhia scheme
radt = 30, !minutes=dx [km]
sf_sfclay_physics = 1, ! Monin-Obukhov
sf_surface_physics = 4, ! NOAHMP
bl_pbl_physics = 1, ! Yonsei University
bldt = 0,
topo_wind = 0, ! only needed for YSU 1, 1, ! Jimenez method
cu_physics = 1, ! Kain-Fritsch
cudt = 5,
cu_rad_feedback = .true. ! add sub-grid cloud effect to the optical depth in radiation
kfeta_trigger = 2,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
mp_zero_out = 0,
sf_urban_physics = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
slope_rad = 0,
topo_shading = 0,
num_land_cat = 21,
iz0tlnd = 1,
shcu_physics = 0 ! needs to be 0 with KF schemes
sf_ocean_physics = 0
hail_opt = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
/
&noah_mp ! josipas settings
dveg = 3 ! dynamic vegetation option
opt_crs = 1 ! stomata resistance
opt_sfc = 1 ! surface layer drag coefficient
opt_btr = 2 ! soil moisture factor for stomatal resistance
opt_run = 3 ! TOPMODEL with groundwater
opt_rad = 3 ! radiative transfer scheme
/
&fdda
/
&dynamics
w_damping = 1, !1 for operational
diff_opt = 1,
km_opt = 4,
gwd_opt = 0,
diff_6th_opt = 2,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000.,
dampcoef = 0.2
khdif = 0
kvdif = 0
non_hydrostatic = .true.
moist_adv_opt = 1,
scalar_adv_opt = 1,
epssm = 0.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,
/
I am trying to run WRF with following namelist options. But when I used aux_hist23 option with custom pressure levels for wrfpress_d01_<domain><date>. The output has lot of missing values where as, wrfout* file has proper values for u and v variables. when I try to extract them for particular pressure levels, some of the values go missing in wrfpress. Is this behaviour normal??. I have compared the pressure levels between metemfiles and my custom levels in &diags section in namelist. everything seems correct. Please help.
NAMELIST WPS
&share
wrf_core = 'ARW',
max_dom = 1,
start_date = '2009-12-31_12:00:00',
end_date = '2009-12-31_18:00:00',
interval_seconds = 21600,
io_form_geogrid = 2,
opt_output_from_geogrid_path = './',
debug_level = 0,
/
&geogrid
parent_id = 1,
parent_grid_ratio = 1,
i_parent_start = 1,
j_parent_start = 1,
e_we = 127,
e_sn = 141,
geog_data_res = '30s',
dx = 30000,
dy = 30000,
map_proj = 'mercator',
ref_lat = 21.474,
ref_lon = 79.796,
truelat1 = 21.474,
truelat2 = 0,
stand_lon = 79.796,
geog_data_path = '../WPS_GEOG',
opt_geogrid_tbl_path = './',
ref_x = 63.5,
ref_y = 70.5,
/
&ungrib
out_format = 'WPS',
prefix = 'FILE',
/
&metgrid
fg_name = 'FILE',
io_form_metgrid = 2,
opt_output_from_metgrid_path = '/DATA/user4/LA_project/',
opt_metgrid_tbl_path = './',
/
&mod_levs
press_pa = 201300 , 200100 , 100000 ,
95000 , 90000 ,
85000 , 80000 ,
75000 , 70000 ,
65000 , 60000 ,
55000 , 50000 ,
45000 , 40000 ,
35000 , 30000 ,
25000 , 20000 ,
15000 , 10000 ,
5000 , 1000
/
&domain_wizard
grib_data_path = '/media/user/Research/ERAint',
grib_vtable = 'Vtable.ERA-interim.pl',
dwiz_name =eraintfine
dwiz_desc =
dwiz_user_rect_x1 =2710
dwiz_user_rect_y1 =545
dwiz_user_rect_x2 =3137
dwiz_user_rect_y2 =997
dwiz_show_political =true
dwiz_center_over_gmt =true
dwiz_latlon_space_in_deg =10
dwiz_latlon_linecolor =-8355712
dwiz_map_scale_pct =25.0
dwiz_map_vert_scrollbar_pos =463
dwiz_map_horiz_scrollbar_pos =2613
dwiz_gridpt_dist_km =30.0
dwiz_mpi_command =null
dwiz_tcvitals =null
dwiz_bigmap =Y
/
namelist.input
&time_control
run_days = 0,
run_hours = 0,
run_minutes = 0,
run_seconds = 0,
start_year = 2008,
start_month = 01,
start_day = 01,
start_hour = 00,
start_minute = 00,
start_second = 00,
end_year = 2008,
end_month = 01,
end_day = 01,
end_hour = 18,
end_minute = 00,
end_second = 00,
interval_seconds = 21600,
input_from_file = .true.,
history_interval = 360,
frames_per_outfile = 1000,
restart = .false.,
restart_interval = 43200,
io_form_history = 2,
io_form_restart = 2,
io_form_input = 2,
io_form_boundary = 2,
debug_level = 0,
nocolons = .false.,
io_form_auxinput4 = 2 !? 11 on cray...parallel netcdf
auxinput4_inname = "wrflowinp_d<domain>"
auxinput4_interval = 360
output_diagnostics = 1,
auxhist3_outname = 'wrfxtrm_d<domain>_<date>',
io_form_auxhist3 = 2,
auxhist3_interval = 1440,
frames_per_auxhist3 =1,
auxhist23_outname ='wrfpress_d<domain>_<date>'
io_form_auxhist23 = 2, !? 11 on cray...parallel netcdf
auxhist23_interval = 360,
frames_per_auxhist23 = 4,
auxhist7_outname ='CTPframework_d<domain>_<date>'
io_form_auxhist7 = 2, !? 11 on cray...parallel netcdf
auxhist7_interval = 180,
frames_per_auxhist7 = 8,
iofields_filename = "additional_fields.txt"
diag_print = 1,
/
&diags
p_lev_diags = 1,
num_press_levels = 19,
press_levels = 95000, 92500, 90000, 87500, 85000, 82500, 80000, 77500, 75000, 72500, 70000, 67500, 65000, 50000, 40000, 30000, 20000, 10000, 5000
use_tot_or_hyd_p = 2,
/
&domains
time_step = 120,
time_step_fract_num = 0,
time_step_fract_den = 1,
max_dom = 1,
e_we = 127,
e_sn = 141,
e_vert = 40,
p_top_requested = 5000,
num_metgrid_levels = 38,
num_metgrid_soil_levels = 4,
dx = 30000,
dy = 30000,
grid_id = 1,
parent_id = 1,
i_parent_start = 1,
j_parent_start = 1,
parent_grid_ratio = 1,
parent_time_step_ratio = 1,
feedback = 0,
smooth_option = 0,
sfcp_to_sfcp = .true.,
/
&physics
sst_update = 1,
mp_physics = 6, ! WSM5 scheme
ra_lw_physics = 1, ! RRTM scheme
ra_sw_physics = 1, ! Dudhia scheme
radt = 30, !minutes=dx [km]
sf_sfclay_physics = 1, ! Monin-Obukhov
sf_surface_physics = 4, ! NOAHMP
bl_pbl_physics = 1, ! Yonsei University
bldt = 0,
topo_wind = 0, ! only needed for YSU 1, 1, ! Jimenez method
cu_physics = 1, ! Kain-Fritsch
cudt = 5,
cu_rad_feedback = .true. ! add sub-grid cloud effect to the optical depth in radiation
kfeta_trigger = 2,
isfflx = 1,
ifsnow = 1,
icloud = 1,
surface_input_source = 1,
num_soil_layers = 4,
mp_zero_out = 0,
sf_urban_physics = 0,
maxiens = 1,
maxens = 3,
maxens2 = 3,
maxens3 = 16,
ensdim = 144,
slope_rad = 0,
topo_shading = 0,
num_land_cat = 21,
iz0tlnd = 1,
shcu_physics = 0 ! needs to be 0 with KF schemes
sf_ocean_physics = 0
hail_opt = 1,
tmn_update = 1,
lagday = 150,
sst_skin = 1,
/
&noah_mp ! josipas settings
dveg = 3 ! dynamic vegetation option
opt_crs = 1 ! stomata resistance
opt_sfc = 1 ! surface layer drag coefficient
opt_btr = 2 ! soil moisture factor for stomatal resistance
opt_run = 3 ! TOPMODEL with groundwater
opt_rad = 3 ! radiative transfer scheme
/
&fdda
/
&dynamics
w_damping = 1, !1 for operational
diff_opt = 1,
km_opt = 4,
gwd_opt = 0,
diff_6th_opt = 2,
diff_6th_factor = 0.12,
base_temp = 290.
damp_opt = 3,
zdamp = 5000.,
dampcoef = 0.2
khdif = 0
kvdif = 0
non_hydrostatic = .true.
moist_adv_opt = 1,
scalar_adv_opt = 1,
epssm = 0.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,
/