| map_continents,/hires problem in IDL 5.0.3 [message #10497] |
Wed, 10 December 1997 00:00  |
f055
Messages: 29
Registered: April 1995
|
Junior Member |
|
|
I have just installed IDL 5.0.3 on our Digital Alpha (Dig. UNIX 4.0) having
previously had IDL 5.0.2. I got it via ftp, and since the high resolution
coastlines only come on CD I had to copy my 5.0.2 hires coastline files to
my 5.0.3 directory.
When using the hires coastlines, however, 5.0.3 doesn't work properly.
Depending upon what region is being used, various bits are missing. Try
the following:
window,ysize=850
!p.multi=[0,1,2,0,0]
;
; Global works fine
;
map_set,/advance
map_continents
map_set,/advance
map_continents,/hires
;
; Regional works fine
;
map_set,/advance,limit=[0,-90,90,90]
map_continents
map_set,/advance,limit=[0,-90,90,90]
map_continents,/hires
;
; Where is Greenland, Iceland and Svalbard in the hi-res version?
;
map_set,/advance,limit=[30,-50,80,50]
map_continents
map_set,/advance,limit=[30,-50,80,50]
map_continents,/hires
;
; Where is Britain?
;
map_set,/advance,limit=[45,-15,65,15]
map_continents
map_set,/advance,limit=[45,-15,65,15]
map_continents,/hires
If I try map_continents,/hires,/coasts then I get Britain, but also many
lakes that I do not want.
I don't think that this problem is related to my using 5.0.2 coastlines with
5.0.3 since the release notes do not indicate a change in the files. But,
map_continents.pro has been changed slightly. Try
UNIX> diff $IDL_DIR/lib/map_continents.pro $IDL5.0.2_DIR/lib
and you get
1c1
< ; $Id: map_continents.pro,v 1.20 1997/10/01 16:38:40 dave Exp $
---
> ; $Id: map_continents.pro,v 1.19 1997/03/26 00:15:35 dave Exp $
177c177
< ((ndx.lonmin lt lonmax) and (ndx.lonmax gt lonmin)), count)
---
> ((ndx.lonmin lt lonmax) OR (ndx.lonmax gt lonmin)), count)
180c180
< subs = where((ndx.lonmin lt lonmax) and (ndx.lonmax gt lonmin), count)
---
> subs = where((ndx.lonmin lt lonmax) OR (ndx.lonmax gt lonmin), count)
186c186
<
---
>
Two ORs have been changed to ANDs. If I try IDL 5.0.3 but compile the old
map_continents
IDL> .run <path of IDL5.0.2>/lib/map_continents.pro
then the above examples work fine.
Has anyone else encountered problems like this. Can I just stick with the
old map_continents.pro, or have the ORs been changed to ANDs for a reason?
Tim
......................... Dr Tim Osborn . t.osborn@uea.ac.uk
.... ___/.. __ /.. /.. /. Senior Research Associate . phone:01603 592089
... /..... /. /.. /.. /.. Climatic Research Unit . fax: 01603 507784
.. /..... __/.. /.. /... School of Environmental Sciences.
. /..... /\ ... /.. /.... University of East Anglia .
____/.._/..\_..____/..... Norwich NR4 7TJ .
......................... UK .
|
|
|
|
| Re: map_continents [message #11432 is a reply to message #10497] |
Thu, 16 April 1998 00:00  |
Harald Frey
Messages: 41
Registered: March 1997
|
Member |
|
|
Harald Frey wrote:
>>> Please try
>>>
>>> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
>>> map_continents
>>>
>>> and you see a nice map of the northern hemisphere continents. But now
>>> limit the region
>>>
>>> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
>>> map_continents
>>>
>>> and not any of the existing coastlines will appear. Any reason?
>>>
>>> Harald Frey
>>> SSL-UC Berkeley
>>> hfrey@ssl.berkeley.edu
>>
I asked RSI about this problem and it turned out that it is really a bug in
IDL 5.0.3. They sent me a new version of the map_continents procedure and
after replacing it for the old version it now works fine.
Thank you, RSI support for the help.
Any interested user can get the corrected map_continents.pro upon request.
At the moment it should still be on their ftp site, but it may not stay
there for a long time.
Harald Frey
|
|
|
|
| Re: map_continents [message #11452 is a reply to message #10497] |
Tue, 14 April 1998 00:00 |
dmarshall
Messages: 23
Registered: December 1997
|
Junior Member |
|
|
Looks OK in student version 5.0.1 , in Win95.
Dave.
In article <MPG.f985fe2273a38ab98977e@news.frii.com>, davidf@dfanning.com (David Fanning) writes:
> Harald Frey (hfrey@ssl.berkeley.edu) writes:
>
>> Please try
>>
>> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
>> map_continents
>>
>> and you see a nice map of the northern hemisphere continents. But now
>> limit the region
>>
>> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
>> map_continents
>>
>> and not any of the existing coastlines will appear. Any reason?
>
> Old version of IDL?
>
> This code looks fine on my Windows NT 4 machine with IDL
> 5.0.3. The continents look correct to me.
>
> Cheers,
>
> David
> -----------------------------------------------------------
> David Fanning, Ph.D.
> Fanning Software Consulting
> E-Mail: davidf@dfanning.com
> Phone: 970-221-0438
> Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
| Re: map_continents [message #11459 is a reply to message #10497] |
Mon, 13 April 1998 00:00 |
David Foster
Messages: 341
Registered: January 1996
|
Senior Member |
|
|
Harald Frey wrote:
>
> Please try
>
> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
> map_continents
>
> and you see a nice map of the northern hemisphere continents. But now
> limit the region
>
> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
> map_continents
>
> and not any of the existing coastlines will appear. Any reason?
I get the same problem as you. Sorry, no solutions.
IDL> help, !version, /struct
** Structure !VERSION, 5 tags, length=40:
ARCH STRING 'sparc'
OS STRING 'sunos'
OS_FAMILY STRING 'unix'
RELEASE STRING '5.0.3'
BUILD_DATE STRING 'Nov 3 1997'
Dave
--
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~
David S. Foster Univ. of California, San Diego
Programmer/Analyst Brain Image Analysis Laboratory
foster@bial1.ucsd.edu Department of Psychiatry
(619) 622-5892 8950 Via La Jolla Drive, Suite 2240
La Jolla, CA 92037
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~
|
|
|
|
| Re: map_continents [message #11462 is a reply to message #10497] |
Mon, 13 April 1998 00:00 |
Harald Frey
Messages: 41
Registered: March 1997
|
Member |
|
|
Harald Frey wrote:
> Martin Schultz wrote:
>
>> Harald Frey wrote:
>>>
>>> Please try
>>>
>>> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
>>> map_continents
>>>
>>> and you see a nice map of the northern hemisphere continents. But now
>>> limit the region
>>>
>>> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
>>> map_continents
>>>
>>> and not any of the existing coastlines will appear. Any reason?
>>>
>>> Harald Frey
>>> SSL-UC Berkeley
>>> hfrey@ssl.berkeley.edu
>>
>> just tested your example: works for me ! I see North America and
>> Greenland as well as a tiny littl ebit of Russia and the UK.
>>
>> What platform and version are you using ?
>
> I use
>
>> IDL> print,!version
>> { sparc sunos unix 5.0.3 Nov 3 1997}
>>
> which is the latest version I know. Is it possible this is some problem
> with the
> monitor? I first thought the problem might be caused by using the low
> resolution continent data only. However, it did not disappear after
> installing the high-res data set.
>
> Harald
Just let me add some more observations about this problem.
I get the grid lines, but I don't get the coast lines. I tries the same on
another computer
IDL> print,!version
{ alpha vms 3.5.1}
And it works just fine there with the full outline of the coastlines. So. I
don't really know if it is a problem of IDL or of SPARC UNIX.
Harald
|
|
|
|
| Re: map_continents [message #11463 is a reply to message #10497] |
Mon, 13 April 1998 00:00 |
Harald Frey
Messages: 41
Registered: March 1997
|
Member |
|
|
Martin Schultz wrote:
> Harald Frey wrote:
>>
>> Please try
>>
>> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
>> map_continents
>>
>> and you see a nice map of the northern hemisphere continents. But now
>> limit the region
>>
>> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
>> map_continents
>>
>> and not any of the existing coastlines will appear. Any reason?
>>
>> Harald Frey
>> SSL-UC Berkeley
>> hfrey@ssl.berkeley.edu
>
> just tested your example: works for me ! I see North America and
> Greenland as well as a tiny littl ebit of Russia and the UK.
>
> What platform and version are you using ?
I use
> IDL> print,!version
> { sparc sunos unix 5.0.3 Nov 3 1997}
>
which is the latest version I know. Is it possible this is some problem
with the
monitor? I first thought the problem might be caused by using the low
resolution continent data only. However, it did not disappear after
installing the high-res data set.
Harald
|
|
|
|
| Re: map_continents [message #11468 is a reply to message #10497] |
Sat, 11 April 1998 00:00 |
Martin Schultz
Messages: 515
Registered: August 1997
|
Senior Member |
|
|
Harald Frey wrote:
>
> Please try
>
> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
> map_continents
>
> and you see a nice map of the northern hemisphere continents. But now
> limit the region
>
> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
> map_continents
>
> and not any of the existing coastlines will appear. Any reason?
>
> Harald Frey
> SSL-UC Berkeley
> hfrey@ssl.berkeley.edu
just tested your example: works for me ! I see North America and
Greenland as well as a tiny littl ebit of Russia and the UK.
What platform and version are you using ?
Martin.
--
------------------------------------------------------------ -------
Dr. Martin Schultz
Department for Earth&Planetary Sciences, Harvard University
186 Pierce Hall, 29 Oxford St., Cambridge, MA-02138, USA
phone: (617)-496-8318
fax : (617)-495-4551
e-mail: mgs@io.harvard.edu
IDL-homepage: http://www-as.harvard.edu/people/staff/mgs/idl/
------------------------------------------------------------ -------
|
|
|
|
| Re: map_continents [message #11471 is a reply to message #10497] |
Fri, 10 April 1998 00:00 |
davidf
Messages: 2866
Registered: September 1996
|
Senior Member |
|
|
Harald Frey (hfrey@ssl.berkeley.edu) writes:
> Please try
>
> map_set,70.,270.,limit=[50.,0.,90.,360.],/grid,/isotropic
> map_continents
>
> and you see a nice map of the northern hemisphere continents. But now
> limit the region
>
> map_set,70.,270.,limit=[50.,180.,90.,360.],/grid,/isotropic
> map_continents
>
> and not any of the existing coastlines will appear. Any reason?
Old version of IDL?
This code looks fine on my Windows NT 4 machine with IDL
5.0.3. The continents look correct to me.
Cheers,
David
-----------------------------------------------------------
David Fanning, Ph.D.
Fanning Software Consulting
E-Mail: davidf@dfanning.com
Phone: 970-221-0438
Coyote's Guide to IDL Programming: http://www.dfanning.com/
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
