Re: can't plot vectors

Vishali,

I don't think this is an NCL problem, but it would help if you could provide your script and data.

--Mary

On Jun 28, 2010, at 2:32 AM, Vishali P. wrote:

> Adam,
>
> Thanks for help, I've tried to change those values , and I even succeeded to cut, but the problem is that the vectors are going all to the same direction. I've also tried to not make zoom of the map, but there's no change on vectors direction (see "vector_2"). Is there any other solution?
>
> vcres@vcRefMagnitudeF = 10.0 ; define vector ref mag
> vcres@vcRefLengthF = 0.00004 ; define length of vec ref
>
>
> Please help
>
> Vishali P
>
>
>
>
> On Wed, Jun 23, 2010 at 2:48 PM, Adam Phillips <asphilli@ucar.edu> wrote:
> VP,
> vcres@vcRefMagnitudeF = 10.0 ; define vector ref mag
> vcres@vcRefLengthF = 0.045 ; define length of vec
>
> I believe the curly vectors are being drawn, but they are too long. Try
> adjusting the resources I pasted above downward.. At first guess I would
> cut them in half... Let us know if you still can't figure it out... Adam
>
> > *Dear NCL users.
> >
> > I am using this script to plot the wind vectors using ncep/reanalysis data
> > 1998. I've tried to plot wind components, but in the plot they aren't
> > going
> > the way that it should be. Something strange is happening.
> >
> >
> > Any help is strongly appreciated
> >
> > Vp*
> >
> > ___________________________________________________________________________________
> >
> > 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("/uwnd.1998.nc","r")
> > print(a)
> > b = addfile("/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.
> >
> > speed = sqrt(u^2+v^2)
> >
> > printVarSummary(u)
> > printVarSummary(v)
> >
> > ;************************************************
> > ; create plot
> > ;************************************************
> > wks = gsn_open_wks("ps","vector_") ; 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
> > _______________________________________________
> > ncl-talk mailing list
> > List instructions, subscriber options, unsubscribe:
> > http://mailman.ucar.edu/mailman/listinfo/ncl-talk
> >
>
>
>
> <vector_2.jpeg>_______________________________________________
> 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 Mon Jun 28 17:19:56 2010