Hi Vishali,
Here is a link where you can find solution for your problem
www.ncl.ucar.edu/Document/Graphics/hlures.shtml
You have to add .hluresfile in your home
Regards
--- En date de : Mer 9.6.10, Vishali P. <vishali.prat@gmail.com> a écrit :
De: Vishali P. <vishali.prat@gmail.com>
Objet: [ncl-talk] error in fuction code
À: ncl-talk@ucar.edu
Date: Mercredi 9 juin 2010, 11h54
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"
;************************************************
begin
;************************************************
; read in netCDF file
;************************************************
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@vcRefLengthF = 0.045 ; define length of vec ref
vcres@vcGlyphStyle = "CurlyVector" ; turn on curly vectors
;************************************************
; Choose a subregion ou fazer zoom
;************************************************
vcres@mpFillOn = False ; turn off gray fill
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
vector = gsn_csm_vector_map_ce(wks,u,v,vcres)
;************************************************
; add time stamp to plot
;************************************************
; size of text desired, and
; right string.
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 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 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 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 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 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 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 35 (T) IS NOT A LEGAL FUNCTION CODE
Please help
Vp
-----La pièce jointe associée suit-----
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Jun 9 07:48:18 2010
This archive was generated by hypermail 2.1.8 : Fri Jun 11 2010 - 09:33:24 MDT