Hi everybody, have you ever come across the following problem with IDL?
I trace contours lines (pressure fields) over the globe in polar
stereographic projection. Everything is fine except when the contour lines
happen to be exactly coincident with the latitude circles (i.e. circles
when seen in polar stereographic projection). In this case the contour
lines simply miss. This happens only when the keyword c_label is used (or
other keywords that impose the 'contour following' algorithm. When
no labels are used, everything is fine.
EXAMPLE:
map_set,90,limit=[20,-180,90,180],/stereo,/noerase,/isot ,/whole_map
map_continents,color=0
map_set,90,limit=[20,-180,80,180], /stereo,title=tit(j-1),/noerase $
,/ISOT,/whole_map
map_grid,latdel=20,londel=45,glinestyle=0,color=0
map_set,90,limit=[20,-180,90,180],/stereo,/noerase,/isot
contour,field,lon,lat,/overplot, levels=lev, $
c_linestyle=5*(LEV Lt 0.0),c_thick=thi,c_labels=lab
Now, I need the lines AND the labels. What can I do?
Thanks in advance
fabio
----------------------
Fabio D'Andrea
Laboratoire de Meteorologie Dynamique
Ecole Normale Superieure
24 Rue Lhomond 75231 PARIS CEDEX 05
Tel: +33 1 44 32 22 32
Fax: +33 1 43 36 83 92
|