Scheduled Downtime
On Friday 21 April 2023 @ 5pm MT, this website will be down for maintenance and expected to return online the morning of 24 April 2023 at the latest

wrfpress has null values at custom pressure levels

This post was from a previous version of the WRF&MPAS-A Support Forum. New replies have been disabled and if you have follow up questions related to this post, then please start a new thread from the forum home page.

vins_94

New member
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,
/
 
Dear Ming Chen,
I would like to make it clear that variables are not missing, but many of the values are null in the specified model levels(provided in &diags of namelist.input). wrfpress output is being created properly but some of the levels are having lot of missing values. whereas in the default levels of wrfout*, all the levels have proper values. Problem is that wrfout is extremely large to work with (~250mb/day of simulation, considering 15 years of files). I was interested only in some of the model levels by writing wrfpress* to history stream 23.
 
Hi,
When p_lev_diags is turned on, the following variables should be included in your wrfpress files:

p_pl,u_pl,v_pl,t_pl,rh_pl,ght_pl,s_pl,td_pl,q_pl

They are pressure level pressure, U, V, RH, GHT, Wind speed, Td, and water vapor mixing ratio.

Please check whether all the above variables exist in your wrfpress files.
 
Dear Chen,
All these variables are appearing in the output, but have missing values at specified pressure levels. While default wrfout file have all the values at all levels output properly.
Thank you
 
Hi,
I am a little perplexed about this issue. Please clarify that:
(1) you have all the variables of pressure, U, V, RH, GHT, Wind speed, Td, and water vapor mixing ratio in wrfpress files.
(2) The above variables have reasonable values in the pressure levels you specify in nameless.input
But
(3) some variables in wrfout (for example T2, Q2, HFX, QCLOUD, etc) are not available in wrfpress.
Please let me know whether this is the issue you have.
 
Top