wrfpress has null values at custom pressure levels

Topics specifically related to any of the WRF runtime options, such as SST update, DFI, nudging, restarts, global runs, adaptive time step, ndown.exe, tc.exe, moving nest, etc.
Post Reply
vins_94
Posts: 8
Joined: Fri May 21, 2021 2:55 pm

wrfpress has null values at custom pressure levels

Post by vins_94 » Fri May 21, 2021 3:13 pm

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,
/

Ming Chen
Posts: 1597
Joined: Mon Apr 23, 2018 9:42 pm

Re: wrfpress has null values at custom pressure levels

Post by Ming Chen » Mon May 24, 2021 9:29 pm

Can you specify what variables are missing in your wrfpress file?
WRF Help Desk

vins_94
Posts: 8
Joined: Fri May 21, 2021 2:55 pm

Re: wrfpress has null values at custom pressure levels

Post by vins_94 » Wed Jun 09, 2021 7:57 am

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.

Ming Chen
Posts: 1597
Joined: Mon Apr 23, 2018 9:42 pm

Re: wrfpress has null values at custom pressure levels

Post by Ming Chen » Mon Jun 14, 2021 11:14 pm

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.
WRF Help Desk

vins_94
Posts: 8
Joined: Fri May 21, 2021 2:55 pm

Re: wrfpress has null values at custom pressure levels

Post by vins_94 » Wed Jun 16, 2021 2:59 pm

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

Ming Chen
Posts: 1597
Joined: Mon Apr 23, 2018 9:42 pm

Re: wrfpress has null values at custom pressure levels

Post by Ming Chen » Tue Jun 22, 2021 10:16 pm

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.
WRF Help Desk

Post Reply

Return to “Special Running Options”