error in fuction code

*dear ncl users
I am a beginner in ncl.*

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"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
  a = addfile("/home/Vp/Documents/uwnd.1998.nc","r")
  print(a)
  b = addfile("/home/Vp/Documents/vwnd.1998.nc","r")
;************************************************
; read in zonal [u] and meridional [v] winds (July)
;************************************************
  u = a->uwnd(1,0,:,:)
  v = b->vwnd(1,0,:,:) ; Get u, v, time (1),level (1000hpa),latitude(-90:90)
and longitude(0:360) data.

  printVarSummary(u)
  printVarSummary(v)

;************************************************
; create plot
;************************************************
  wks = gsn_open_wks("ps","vector1") ; open a ps file

  vcres = True ; plot mods desired
  vcres@gsnFrame = False ; so we can draw time
stamp
  vcres@vcRefAnnoOrthogonalPosF = -1.0 ; move ref vector up
  vcres@vcRefMagnitudeF = 10.0 ; define vector ref mag
  vcres@vcRefLengthF = 0.045 ; define length of vec
ref
  vcres@vcGlyphStyle = "CurlyVector" ; turn on curly vectors
  vcres@vcMinDistanceF = 0.017

;************************************************
; Choose a subregion ou fazer zoom
;************************************************
  vcres@mpFillOn = False ; turn off gray fill
  vcres@mpOutlineBoundarySets = "National" ; turn on country
boundaries
  vcres@mpGeophysicalLineColor = "Navy" ; color of cont. outlines
  vcres@mpGeophysicalLineThicknessF = 1.5 ; thickness of outlines

  vcres@mpMaxLatF = 0 ;maximum latitude
  vcres@mpMinLatF = -45 ;minimum latitude
  vcres@mpMaxLonF = 55 ;maximum longitude
  vcres@mpMinLonF = 5 ;minimum longitude

  vector = gsn_csm_vector_map_ce(wks,u,v,vcres)
;************************************************
; add time stamp to plot
;************************************************
  infoTimeStamp(wks,0.025,"uwnd.1998.nc") ; arguments are workstation
                                           ; size of text desired, and
                                           ; right string.
  frame (wks)
end

*
I am using this script to plot the wind vectors using ncep/reanalysis data
1998 and I have this strange message*

 PLCHHQ - CHARACTER NUMBER 13 (W) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 14 (e) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 15 (d) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 17 (J) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 18 (u) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 19 (n) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 22 (9) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 34 (M) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 35 (T) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 13 (W) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 14 (e) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 15 (d) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 17 (J) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 18 (u) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 19 (n) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 22 (9) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 34 (M) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 35 (T) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 13 (W) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 14 (e) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 15 (d) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 17 (J) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 18 (u) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 19 (n) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 22 (9) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 34 (M) IS NOT A LEGAL FUNCTION CODE
 PLCHHQ - CHARACTER NUMBER 35 (T) IS NOT A LEGAL FUNCTION CODE

*
Please help

Vp*

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk