Hi everyone:
With IDL 8.0-32bit, the following program does not show the filled contour over map.
But with IDL 7.1, this program runs successfully. I analized this code, and found that remove
the "LIMIT=[17.0000, 0.0000, 55.0000, 140.000]" of map_set, this program shows
the filled contour. (ps, this code also runs successfully with IDL 8.0-64bit version). So where
are the problems? Thanks.
pro test_map_contour
; Make a 10 degree latitude/longitude grid covering the Earth:
lat = REPLICATE(10., 37) # FINDGEN(19) - 90.
lon = FINDGEN(37) # REPLICATE(10, 19)
; Convert lat and lon to Cartesian coordinates:
X = COS(!DTOR * lon) * COS(!DTOR * lat)
Y = SIN(!DTOR * lon) * COS(!DTOR * lat)
Z = SIN(!DTOR * lat)
; Create the function to be plotted, set it equal
; to the distance squared from (1,1,1):
F = (X-1.)^2 + (Y-1.)^2 + (Z-1.)^2
; Create a plotting window
WINDOW, 0, TITLE='Mollweide Contour'
device,decomposed=0
loadct,13
;MAP_SET, /stereographic, 0, 0, /ISOTROPIC, $
; /HORIZON, /GRID, /CONTINENTS, $
; TITLE='Mollweide Contour'
MAP_SET, 90, 106., LIMIT=[17.0000, 0.0000, 55.0000, 140.000], $
/NOBORDER, /STEREOGRAPHIC
CONTOUR, F, lon, lat, NLEVELS=12, c_colors=indgen(12)*20,$
/OVERPLOT, /DOWNHILL, /FOLLOW, /cell_fill
map_continents
end
|