[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