| Problems with /cell_fill [message #3096] |
Wed, 02 November 1994 06:31 |
polcher
Messages: 3
Registered: November 1994
|
Junior Member |
|
|
Hi
I am plotting a data set which covers the globe (+-177, +- 79) with IDL using the option
/cell_fill. It works very well for global projections (Mercator, Mollweide, Cylindrical)
when Greenwhich is in the middle of the picture.
When I rotate the globe in order to have another longitude in the middle the eastern part of
the map shows a problem. The filled contours are replaced by horizontal lines, this does not
happen for isolines !!! It seems to me that IDL has problems joining the data between +180
-180, the "periodicity" is not properly done. (A simple example follows)
A similar problem appears with projection which do not show the entire globe
(Stereographic, Orthographic, ...). The filled contours extend beyond the borders of the
globe.
Is there an easy patch to this problem ?
Sorry if this problem has already been posted in this news group.
Thanks for any help
Jan
EXAMPLE:
--------
pro test, lon=sellon
iim=64
jjm=50
lon=fltarr(iim,jjm)
lat=fltarr(iim,jjm)
FOR i=0,iim-1 DO lat(i,*)=-(findgen(jjm)/(jjm-1)*180-90)
FOR i=0,jjm-1 DO lon(*,i)=findgen(iim)/(iim-1)*360-180
car=fltarr(iim,jjm)
car=sin(findgen(iim,jjm))
map_set,0,sellon,0, /cyl, grid=1, continents=1
loadct, 12
contour, car, lon, lat , nlevels=10, /cell_fill, /overplot
map_continents
end
When we do "test, lon=0" it works and for "test, lon=120" it fails.
Or if /cyl is replaced by /stere it fails also.
--
Jan Polcher TEL: -33-1-44322243
Laboratoire de Meteorologie Dynamique / E.N.S. FAX: -33-1-43368392
24, rue Lhomond
75231 PARIS cedex 05, FRANCE
|
|
|
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
