Hi,
I'm using the contour function to draw values on a sub-region of a map.
The problem is that I obtain contour outside the limits I have defined.
I read the idl guide and went on this group but I didn't find the right
way to draw the contour. I have contour outside the limits only on one
side of my sub-region for the region between (-10.5E;-41.5E) and
latitude=57.5N. Hereafter my code:
nlevels=24
labels=fltarr(nlevels)
limits=[40,-60,80,0]
colors = fltarr(nlevels)
for i = 0, nlevels-1 do begin
colors(i) = 10 + i * (229-20) / (nlevels-1)
endfor
labels=[-1,0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19 ,20,21,22]
map_set, /cylindrical, color=white, /noerase, /continents,
Limit=limits, xmargin=5, ymargin=10
map_continents, color = black, mlinethick= 2.0
map_grid,/box_axes
x=sst_nadr[1,*] ; the longitude between -10.5E and -41.5E
y=sst_nadr[2,*] ; the latitude between 44.5N and 57.5N
sst=transpose(sst_nadr[3,*]) ; the values
contour,sst,x,y, levels = labels, c_colors = colors, /overplot,
/cell_fill, /irregular
Does anybody have an idea of the "bug"?
Thank you,
Cedric
|