comp.lang.idl-pvwave archive
Messages from Usenet group comp.lang.idl-pvwave, compiled by Paulo Penteado

Home » Public Forums » archive » (cG) windbarb problem
Show: Today's Messages :: Show Polls :: Message Navigator
E-mail to friend 
Return to the default flat view Create a new topic Submit Reply
Re: (cG) windbarb problem [message #88162 is a reply to message #88160] Tue, 25 March 2014 06:41 Go to previous messageGo to previous message
natha is currently offline  natha
Messages: 482
Registered: October 2007
Senior Member
I fixed the routine by using the CLIP keyword and by modifying the lines 245 and 246
x1 = clip[0] > (xx[j] + sindr * sr) < clip[2]
y1 = clip[1] > (yy[j] + cosdr * sr * aspect) < clip[3]:
by
x1 = clip[0] > (xx[j]) < clip[2]
y1 = clip[1] > (yy[j]) < clip[3]

Now everything is working properly

device, dec=0

!p.background=255
!p.color=0

limit=[20.,-128.,53.,-65.]

dim=[800,800]

window, 0, xsize=dim[0], ysize=dim[1]

map_set, 0, 0, limit=limit, position=[0,0,1,1], /noborder, /isotropic, /mercator

map_continents, /countries, /coasts

uv_box=!map.uv_box
aspect_ratio=(uv_box[2]-uv_box[0])/(uv_box[3]-uv_box[1])

data=intarr(100,100)
temp=map_image(data,xstart,ystart,xsize,ysize,latmin=limit[0 ],latmax=limit[2],$
lonmin=limit[1],lonmax=limit[3],/compress)

clip=[xstart,ystart,xstart+xsize,ystart+ysize]/float([dim,di m])

windbarb, [-110], [35], [50.], [45.], length=0.3, color=[0,0,0], clip=clip, aspect=1.

plots, [-110], [35], psym=1
[Message index]
 
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Read Message
Previous Topic: Trouble plotting GSHHS shapefile
Next Topic: a GDL Workshop in June 2014

-=] Back to Top [=-
[ Syndicate this forum (XML) ] [ RSS ] [ PDF ]

Current Time: Thu Oct 09 19:53:27 PDT 2025

Total time taken to generate the page: 0.00186 seconds