Re: vectors on map-- partvelvec? [message #80719] |
Wed, 04 July 2012 04:11 |
anil
Messages: 34 Registered: August 2009
|
Member |
|
|
On Jul 3, 4:10 pm, David Fanning <n...@idlcoyote.com> wrote:
> anil writes:
>> I am trying to plot a map and on top of it, the geostrophic velocities
>> from a data source.
>> I am using partvelvec.pro . It works fine when i do not do any
>> positioning, but this time the plot is awkward.
>> When i set a position like;
>> Map_Set,/grid, Limit=[40 ,27,48, 42],position=[0.1,0.15,0.9,0.9] ;
>> for example, and i do ,
>> partvelvec,u,v,lon,lat,/over
>> then i end up with some of the data being plotted on the map but some
>> are missing !!!??!!.
>> The missing part corresponds ~ to the second half of the data file( an
>> ascii file containing longitude, latitude,u,v and dates).
>> Why is this happening? I know the partvelvec function uses cg*
>> routines from Coyote library, but I have it installed in my libray.
>> Why is it plotting everything, when position is not set but it is not
>> when the position is set?
>> I checked the ascii files a thousand times but nothing seems to be
>> wrong, just an ordinary ascii file.
>
> The good news is, this is not possible. The bad news is,
> the hallucinations are getting worse. :-(
>
> What happens if you try this without a map projection?
>
> Cheers,
>
> David
>
> --
> David Fanning, Ph.D.
> Fanning Software Consulting, Inc.
> Coyote's Guide to IDL Programming:http://www.idlcoyote.com/
> Sepore ma de ni thui. ("Perhaps thou speakest truth.")
It seems that i have forgot to close the device and set_plot at the
end of the code. That was the medication to the hallucinations :).
Thanks David.
|
|
|
Re: vectors on map-- partvelvec? [message #80726 is a reply to message #80719] |
Tue, 03 July 2012 06:10  |
David Fanning
Messages: 11724 Registered: August 2001
|
Senior Member |
|
|
anil writes:
> I am trying to plot a map and on top of it, the geostrophic velocities
> from a data source.
> I am using partvelvec.pro . It works fine when i do not do any
> positioning, but this time the plot is awkward.
> When i set a position like;
> Map_Set,/grid, Limit=[40 ,27,48, 42],position=[0.1,0.15,0.9,0.9] ;
> for example, and i do ,
> partvelvec,u,v,lon,lat,/over
> then i end up with some of the data being plotted on the map but some
> are missing !!!??!!.
> The missing part corresponds ~ to the second half of the data file( an
> ascii file containing longitude, latitude,u,v and dates).
> Why is this happening? I know the partvelvec function uses cg*
> routines from Coyote library, but I have it installed in my libray.
> Why is it plotting everything, when position is not set but it is not
> when the position is set?
> I checked the ascii files a thousand times but nothing seems to be
> wrong, just an ordinary ascii file.
The good news is, this is not possible. The bad news is,
the hallucinations are getting worse. :-(
What happens if you try this without a map projection?
Cheers,
David
--
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.idlcoyote.com/
Sepore ma de ni thui. ("Perhaps thou speakest truth.")
|
|
|