| Re: IDL 5.2.1 -- bug in filled contours? [message #17440] |
Wed, 27 October 1999 00:00 |
R.Bauer
Messages: 1424
Registered: November 1998
|
Senior Member |
|
|
Bob Yantosca wrote:
> Hi all,
>
> I've noticed a potential bug in IDL 5.2.1 when you try to overplot filled
> contours atop a map. Basically, what happens is that you get a
> small square of color that covers a tiny square centered on Lat=0,
> Lon=0, instead of covering the entire map area.
>
> See below for a test program I've devised to illustrate the error. The
> first plot (with regular color contours) comes out fine. When you use the
> /FILL keyword then you run into this bug. (I tried /CELL_FILL too, but
> that didn't make any difference).
>
> Per Martin Schultz's suggestion, I tried resetting the Y-axis indices so
> that the poles were -89.9 and +89.9 (instead of -90 and +90), but on my
> system (SGI-IRIX), that didn't help at all.
>
> Any suggestions?
>
> Thanks much,
>
> Bob Yantosca
> bmy@io.harvard.edu
>
> ------------------------------------------------------------ ----------
> pro Test_Contour
> loadct, 27 ; EOS-B Colortable
> !P.MULTI = [0, 1, 2, 0, 0] ; 2 plots per page
> TmpData = Dist(72, 46) ; Test array for plot
> XArr = findgen(72)*5.0 - 180.0 ; X-indices
> YArr = findgen(46)*4.0 - 90.0 ; Y-indices
> C_Levels = findgen(10)*4 ; Contour levels
> C_Colors = findgen(10)*10 ; Contour colors
>
> ; Call MAP_SET and CONTOUR with unfilled color contours
> Map_Set, /Grid, /Continents, Color=2, $
> Title='IDL 5.2 Color Contours (OK)'
>
> Contour, TmpData, XArr, YArr, $
> Levels=C_Levels, C_Colors=C_Colors, /OverPlot
>
> ; Call MAP_SET and CONTOUR with filled contours
> Map_Set, /Grid, /Continents, Color=2, /Advance, $
> Title='IDL 5.2 Filled Contours (/FILL)...screwed up'
>
> Contour, TmpData, XArr, YArr, $
> Levels=C_Levels, C_Colors=C_Colors, /OverPlot, /Fill
> end
>
/cell_fill works
I am using
IDL5.2.1 on WinNT
You should never use /fill with overplots on maps.
The failure on your platform is that's /cell_fill won't work.
R.Bauer
|
|
|
|
| Re: IDL 5.2.1 -- bug in filled contours? [message #17444 is a reply to message #17440] |
Tue, 26 October 1999 00:00  |
Charles Cavanaugh
Messages: 8
Registered: July 1997
|
Junior Member |
|
|
<!doctype html public "-//w3c//dtd html 4.0 transitional//en">
<html>
Bob,
<p>We ugraded, we saw, we reported. This bug is alive inside
the IDL demo as well.
<br>Fire up the demo, click on "Earth Sciences", click on "Mapping".
Look familiar?
<br>We grabbed a 5.3 beta version and the bug was still there.
<p>Charles
<br>
<br>
<p>Bob Yantosca wrote:
<blockquote TYPE=CITE>Hi all,
<p>I've noticed a potential bug in IDL 5.2.1 when you try to overplot filled
<br>contours atop a map. Basically, what happens is that you get
a
<br>small square of color that covers a tiny square centered on Lat=0,
<br>Lon=0, instead of covering the entire map area.
<p>See below for a test program I've devised to illustrate the error. The
<br>first plot (with regular color contours) comes out fine. When
you use the
<br>/FILL keyword then you run into this bug. (I tried /CELL_FILL
too, but
<br>that didn't make any difference).
<p>Per Martin Schultz's suggestion, I tried resetting the Y-axis indices
so
<br>that the poles were -89.9 and +89.9 (instead of -90 and +90), but on
my
<br>system (SGI-IRIX), that didn't help at all.
<p>Any suggestions?
<p>Thanks much,
<p>Bob Yantosca
<br>bmy@io.harvard.edu
<p> ------------------------------------------------------------ ----------
<br>pro Test_Contour
<br> loadct, 27 &nbs p; &nbs p; &nbs p; &nbs p;
; EOS-B Colortable
<br> !P.MULTI = [0, 1, 2, 0, 0] &nbs p;
; 2 plots per page
<br> TmpData = Dist(72, 46) &nb sp; &nb sp;
; Test array for plot
<br> XArr = findgen(72)*5.0 - 180.0
; X-indices
<br> YArr = findgen(46)*4.0 - 90.0
; Y-indices
<br> C_Levels = findgen(10)*4 &nb sp; &nb sp;
; Contour levels
<br> C_Colors = findgen(10)*10 &n bsp;
; Contour colors
<p> ; Call MAP_SET and CONTOUR with unfilled color contours
<br> Map_Set, /Grid, /Continents, Color=2, $
<br> Title='IDL 5.2 Color Contours (OK)'
<p> Contour, TmpData, XArr, YArr, $
<br> Levels=C_Levels, C_Colors=C_Colors,
/OverPlot
<p> ; Call MAP_SET and CONTOUR with filled contours
<br> Map_Set, /Grid, /Continents, Color=2, /Advance,
$
<br> Title='IDL 5.2 Filled Contours (/FILL)...screwed
up'
<p> Contour, TmpData, XArr, YArr, $
<br> Levels=C_Levels, C_Colors=C_Colors,
/OverPlot, /Fill
<br>end
<p>----
<br>Bob Yantosca &a mp;nbsp;
| Phone: (617) 495-4517 \ "Ooooh...floor pie!"
<br>Scientific Programmer | FAX: (617) 495-4551
\ -- Homer J. Simpson
<br>Atmospheric Chemistry | Email: bmy@io.harvard.edu +-------------------------
<br>Harvard University | WWW: <a href="http://www-as.harvard.edu/people/staff/bmy/">http://www-as.harvard.edu/people/staff/bmy/</a></blockquote>
<pre>--
Charles Cavanaugh
MOPITT & HIRDLS Programs
National Center for Atmospheric Research
Boulder, CO, USA
email: cavanaug@ucar.edu
phone: (303) 497-2919
fax : (303) 497-2920</pre>
</html>
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
