[ncl-talk] NCEP local GRIB2 table
Paula Freeman - NOAA Affiliate
paula.freeman at noaa.gov
Mon Jul 2 10:23:10 MDT 2018
Hello,
I am a complete newbie to NCL, I am grateful for any help you can provide.
I am trying to plot the accumulated 6-hr precip from a 1-degee GFS GRIB2
file.
ncl_filedump lists the grid I want as APCP_P8_L1_GLL0_acc6h:
$ ncl_filedump gfs.t12z.pgrb2.1p00.f012.grib2
<...>
float APCP_P8_L1_GLL0_acc6h ( lat_0, lon_0 )
center : US National Weather Service - NCEP (WMC)
production_status : Operational products
long_name : Total precipitation
units : kg m-2
_FillValue : 1e+20
grid_type : Latitude/longitude
parameter_discipline_and_category : Meteorological products,
Moisture
parameter_template_discipline_category_number : ( 8, 0, 1,
8 )
level_type : Ground or water surface
level : 0
type_of_statistical_processing : Accumulation
statistical_process_duration : 6 hours (ending at forecast time)
forecast_time : 12
forecast_time_units : hours
initial_time : 06/25/2018 (12:00)
Here's my simple script:
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
;---Open file and read data
filename = "/u/Paula.Freeman/test/gfs.t12z.pgrb2.1p00.f012.grib2"
a = addfile(filename,"r")
p = a->APCP_P8_L1_GLL0_acc6h
;---Start the graphics
wks = gsn_open_wks("png","/u/Paula.Freeman/test/gfs_f012_precip")
;---Set some resources
res = True
res at cnFillOn = True ; turn on color fill
res at gsnMaximize = True ; maximize plot in frame
plot = gsn_csm_contour_map(wks,p,res)
end
I get the following error:
$ ncl gfs_precip_p06.ncl
Copyright (C) 1995-2012 - All Rights Reserved
University Corporation for Atmospheric Research
NCAR Command Language Version 6.1.0
The use of this software is governed by a License Agreement.
See http://www.ncl.ucar.edu/ for more details.
warning: Entry (21) not found in code table file
/usrx/local/NCL_NCARG-6.1.0_gcc446/OpenNDAPenabled/lib/
ncarg/grib2_codetables/ncep/4/4.2.0.0.table
When I look at the table listed above, there is no entry for 21, just this
placeholder:
# 21-191:21-191:Reserved:-:
# 192-254:192-254:Reserved for Local Use:-:
So it appears that these tables are not updated for NCEP's local variables.
How can I direct NCL to an updated version of the tables? My local GRIB2
expert has pointed me to a local copy.
I am running this script on the WCOSS at NCEP:
$ uname -a
Linux g14a1 2.6.32-573.48.1.el6.x86_64 #1 SMP Thu Sep 21 11:19:32 EDT
2017 x86_64 x86_64 x86_64 GNU/Linux
Thanks,
Paula Freeman
--
Paula Freeman
Software Developer
Cyberdata Technologies, Inc.
NOAA/NWS/NCEP/NCO
NCWCP, Cube 1034
Onsite: (301) 683-3837
Teleworking: (301)283-3596
Cell: (202)709-1537
Paula.Freeman at noaa.gov
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20180702/a276d6e8/attachment.html>
More information about the ncl-talk
mailing list