<DIV>
<DIV><FONT size=4>Hi all:</FONT></DIV>
<DIV><FONT size=4> I am trying to make a plot with grib2 data(GFS),but it has warning and no picture appeared.</FONT></DIV>
<DIV><FONT size=4> 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</FONT></DIV>
<DIV><FONT size=4> I tried to use ncl_convert2nc,it also have the warning:</FONT></DIV>
<DIV>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<BR><FONT size=4>This is my data:</FONT></DIV>
<DIV><IMG src="cid:6A7E73CC@F326FE08.75484456.jpg" filesize="255656" naturalW="521" naturalH="432"></DIV>
<DIV><FONT size=4> How can I slove the problem?</FONT></DIV>
<DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"><FONT style="LINE-HEIGHT: 30px" size=2> <FONT size=4>Or do you guys have some advice?</FONT></FONT></SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"><FONT style="LINE-HEIGHT: 30px" size=4> Any information will be appreciated.</FONT></SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"><FONT size=4> This is my script:</FONT></SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px">load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/skewt_func.ncl"<BR>load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/wind_rose.ncl"<BR>;*********************************************************<BR>begin<BR> f = addfile("/public/home/huanglei/data/2015080200//gfs.t00z.pgrb2.0p50.f000.grb","r")<BR>;********************************************************************</SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"></SPAN> </DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"> prc = f->PWAT_P0_L200_GLL0<BR>;<BR>; The lat,lon coordinate points are only used in this<BR>; script to determine the area of the map to look at.<BR>; You could potentially use them to over this data<BR>; on a different map projection. See end of script.<BR>;<BR> lat2d = f->lat_0 ; (g5_x_0, g5_y_1)<BR> lon2d = f->lon_0 ; (g5_x_0, g5_y_1)</SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"> wks =gsn_open_wks("pdf","gfs_precip")<BR> gsn_define_colormap(wks,"BlAqGrYeOrReVi200")<BR> res = True<BR> <A href="mailto:res@gsnAddCyclic">res@gsnAddCyclic</A> = False<BR> <A href="mailto:res@gsnDraw">res@gsnDraw</A> = False ; don't draw yet<BR> <A href="mailto:res@gsnFrame">res@gsnFrame</A> = False ; don't advance frame yet<BR> <A href="mailto:res@mpOutlineOn">res@mpOutlineOn</A> = True ; Turn on map outlines<BR> <A href="mailto:res@mpGeophysicalLineColor">res@mpGeophysicalLineColor</A> = "Black"<BR> <A href="mailto:res@mpNationalLineColor">res@mpNationalLineColor</A> = "Black"<BR> <A href="mailto:res@mpDataSetName">res@mpDataSetName</A> = "Earth..4"<BR> <A href="mailto:res@mpDataBaseVersion">res@mpDataBaseVersion</A> = "MediumRes"<BR> <A href="mailto:res@mpOutlineSpecifiers">res@mpOutlineSpecifiers</A> = ("China:states")<BR> <A href="mailto:res@mpNationalLineThicknessF">res@mpNationalLineThicknessF</A>= 0.1<BR> <A href="mailto:res@mpOutlineDrawOrder">res@mpOutlineDrawOrder</A> = "PostDraw"<BR> <A href="mailto:res@mpLimitMode">res@mpLimitMode</A> = "LatLon" ; Limit the map view.<BR>; <A href="mailto:vcres@mpGeophysicalLineThicknessF">vcres@mpGeophysicalLineThicknessF</A> = 0.35 ; double the thickness of geophysical boundaries<BR>; <A href="mailto:vcres@mpNationalLineThicknessF">vcres@mpNationalLineThicknessF</A> = 0.5 ; double the thickness of national boundaries</SPAN></DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><SPAN style="LINE-HEIGHT: 24px"> <A href="mailto:res@pmTickMarkDisplayMode">res@pmTickMarkDisplayMode</A> = "Always"<BR> <A href="mailto:res@mpMinLatF">res@mpMinLatF</A> = 33 ;驽扮含搴?<BR> <A href="mailto:res@mpMaxLatF">res@mpMaxLatF</A> = 35 ;驽ぇ绾害<BR> <A href="mailto:res@mpMinLonF">res@mpMinLonF</A> = 107 ;驽扮诲害<BR> <A href="mailto:res@mpMaxLonF">res@mpMaxLonF</A> = 110 ;驽ぇ缁搴?</SPAN></DIV><SPAN style="LINE-HEIGHT: 24px">
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><BR>;************************************************************************************************</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"> <A href="mailto:res@cnFillOn">res@cnFillOn</A> = True ; turn on color<BR> <A href="mailto:res@gsnSpreadColors">res@gsnSpreadColors</A> = True ; use full colormap<BR> <A href="mailto:res@gsnSpreadColorStart">res@gsnSpreadColorStart</A> = 0<BR> <A href="mailto:res@gsnSpreadColorEnd">res@gsnSpreadColorEnd</A> = 199<BR> <A href="mailto:res@cnLinesOn">res@cnLinesOn</A> = False ; turn off contour lines<BR> <A href="mailto:res@cnLineLabelsOn">res@cnLineLabelsOn</A> = False ; tuen off line labels<BR> <A href="mailto:res@cnInfoLabelOn">res@cnInfoLabelOn</A> =False</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'">; <A href="mailto:res@cnLevelSelectionMode=%22ExplicitLevels">res@cnLevelSelectionMode="ExplicitLevels</A>"<BR>; <A href="mailto:res@cnLevels=(/14,18,22,26,30,34,38,42,46,50,52/">res@cnLevels=(/14,18,22,26,30,34,38,42,46,50,52/</A>)</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"> <A href="mailto:res@gsnLeftString">res@gsnLeftString</A> =""<BR> <A href="mailto:res@gsnRightString">res@gsnRightString</A> =""<BR> <A href="mailto:res@tiMainString">res@tiMainString</A> =""</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"> <A href="mailto:res@lbLabelBarOn">res@lbLabelBarOn</A> = False<BR> <A href="mailto:res@lbLabelStride">res@lbLabelStride</A> = 2 ; plot every other colar bar label<BR> <A href="mailto:res@lbOrientation">res@lbOrientation</A> = "vertical" ; vertical label bars</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"> </DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"> prc=smth9(prc,0.5,0.25,False)</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"><BR> plot = gsn_csm_contour_map(wks,prc,res)</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'">frame(wks)</DIV>
<DIV style="LINE-HEIGHT: 23px; FONT-FAMILY: 'lucida Grande', Verdana, 'Microsoft YaHei'"></DIV> end<BR> </SPAN></DIV></DIV>
<DIV></DIV>