[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

[David Brown]

The warning concerns another parameter in the GRIB file for which no
name or description was available at the time of the release of 6.1.0.
This information is available for more recent versions of NCL and if
you upgraded NCL to 6.3.0 you would no longer get the warning.
However, since you are reading a different parameter this warning has
nothing to do with your problem. It appears as if your issue may
simply be that you have set gsnDraw to False. This prevents the the
gsn_csm_contour_map routine from actually drawing the plot, requiring
you to draw it explicitly. You have two choices:

1) comment out or remove the line:
  res at gsnDraw = False
or 2) add an explicit draw command:
 plot = gsn_csm_contour_map(wks,prc,res)
 draw(plot)
 frame(wks)

If this does not fix your problem, then you need to examine your data
variable, prc. Use printVarSummary, printMinMax, etc. Let us know if
you still have problems.
 -dave




On Thu, Nov 12, 2015 at 12:25 AM, grace <313695096 at qq.com> wrote:
> 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
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>