[ncl-talk] Unable to plot vector wind for a regional domain

[Geeta Geeta]

HiI have to plot 850hpa Vector wind for a regional domain (0-40N/50-100E). My dataset is FNL reanalysis. I tried following scripts. (png file is attached);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/wrf/WRFUserARW.ncl"; Open the file----------begina = addfile("~/Mesomodel/gfs/test.nc","r");printVarSummary(a);; -----Read data------Uwind = a->UGRD_P0_L100_GLL0(20,:,:); Read the variable to memoryVwind = a->VGRD_P0_L100_GLL0(20,:,:); Read the variable to memoryUwind at _FillValue = -999Vwind at _FillValue = -999printVarSummary(Uwind)printVarSummary(Vwind);; ----Create Plot-----------------wks   = gsn_open_wks("png", "wind850")res                             = Trueres at tiMainString                = "WIND 850hPa ";pltres                         = True;;               Map Resourcesmpres                           = Trueres at mpFillOn                    =Falseres at mpOutlineBoundarySets       ="National"res at mpGeophysicalLineColor      ="Black"res at mpMinLatF                   = 0.0res at mpMaxLatF                   = 40.0res at mpMinLonF                   = 50.0res at mpMaxLonF                   = 100.0res at mpDataBaseVersion           ="MediumRes";plot  = gsn_csm_map_ce(wks,res)        ;Draws a map;plot  = gsn_csm_vector_map(wks,Uwind,Vwind,True);plot  = gsn_csm_map_ce(wks,True);                                                                                                                             ;vecres                           =Trueres at vcMinFracLengthF             =1.0           ; Length of min vectorres at vcRefMagnitudeF              =3.0           ; Makes Vectors largerres at vcRefLengthF                 =0.045         ; ref vec Lengthlargerres at vcGlyphStyle                 ="CurlyVector" ; ref vec Lengthlargerres at vcMinDistanceF               =0.01          ; ref vec Lengthlargerres at vcMonoLineArrowColor         =False         ; ref vec Lengthlarger plot  = gsn_csm_vector_map(wks,Uwind,Vwind,True)end With this script, the problem is it is NOT Plotting the data for the regional domain.  

A little change in the above script (changes given below)plot  = gsn_csm_map_ce(wks,res)        ;Draws a map
draws a regional map BUT without data. 
Both figures attatched. 
Any suggestions............ Geeta.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150526/95be1d4b/attachment.html 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wind850.1.png
Type: image/png
Size: 4565 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150526/95be1d4b/attachment.png 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: wind850.png
Type: image/png
Size: 6311 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20150526/95be1d4b/attachment-0001.png