[ncl-talk] warning: Entry (21) not found in code table file /public/software/mathlib/ncl_ncarg/6.1.0/gnu/lib/ncarg/grib2_codetables/ncep/4/4.2.0.0.table

[grace]

Hi all:
    I am trying to make a  plot with grib2 data(GFS),but it has warning and no picture appeared.
   warning: Entry (21) not found in code table file /public/software/mathlib/ncl_ncarg/6.1.0/gnu/lib/ncarg/grib2_codetables/ncep/4/4.2.0.0.table
     
   How can I slove the problem?
     Or do you guys have some advice?
    Any information will be appreciated.
    This is my 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/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"
;*********************************************************
begin
  f = addfile("/public/home/huanglei/data/2015080200//gfs.t00z.pgrb2.0p50.f000.grb","r")
;********************************************************************
  
   prc   = f->PWAT_P0_L200_GLL0
;
; The lat,lon coordinate points are only used in this
; script to determine the area of the map to look at.
; You could potentially use them to over this data
; on a different map projection. See end of script.
;
  lat2d = f->lat_0               ;  (g5_x_0, g5_y_1)
  lon2d = f->lon_0               ;  (g5_x_0, g5_y_1)
   wks =gsn_open_wks("pdf","gfs_precip")
  gsn_define_colormap(wks,"BlAqGrYeOrReVi200")
  res = True
  res at gsnAddCyclic  = False
  res at gsnDraw      = False                        ; don't draw yet
  res at gsnFrame     = False                        ; don't advance frame yet
  res at mpOutlineOn           = True         ; Turn on map outlines
  res at mpGeophysicalLineColor = "Black"
  res at mpNationalLineColor    = "Black"
  res at mpDataSetName         = "Earth..4"
  res at mpDataBaseVersion     = "MediumRes"
  res at mpOutlineSpecifiers = ("China:states")
  res at mpNationalLineThicknessF= 0.1
  res at mpOutlineDrawOrder    = "PostDraw"
  res at mpLimitMode = "LatLon"    ; Limit the map view.
;  vcres at mpGeophysicalLineThicknessF = 0.35      ; double the thickness of geophysical boundaries
;  vcres at mpNationalLineThicknessF    = 0.5      ; double the thickness of national boundaries
   res at pmTickMarkDisplayMode = "Always"
  res at mpMinLatF         = 33                         ;驽扮含搴?
  res at mpMaxLatF         = 35                        ;驽ぇ绾害
  res at mpMinLonF         = 107                        ;驽扮诲害
  res at mpMaxLonF         = 110                       ;驽ぇ缁搴?
 
;************************************************************************************************
  res at cnFillOn        = True                    ; turn on color
 res at gsnSpreadColors = True                    ; use full colormap
 res at gsnSpreadColorStart = 0
 res at gsnSpreadColorEnd   = 199
 res at cnLinesOn       = False                   ; turn off contour lines
 res at cnLineLabelsOn  = False                   ; tuen off line labels
 res at cnInfoLabelOn =False
 ; res at cnLevelSelectionMode="ExplicitLevels"
; res at cnLevels=(/14,18,22,26,30,34,38,42,46,50,52/)
   res at gsnLeftString =""
  res at gsnRightString =""
 res at tiMainString =""
  res at lbLabelBarOn = False
 res at lbLabelStride       = 2         ; plot every other colar bar label
 res at lbOrientation        = "vertical"         ; vertical label bars
  
    prc=smth9(prc,0.5,0.25,False)
 
 plot = gsn_csm_contour_map(wks,prc,res)
 frame(wks)
   end
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151112/e9a584e7/attachment.html