I am having problems getting some filled contour charts of sea surface
temperature to look right on a new Linux machine. They look fine on
the old machine with the exact same program, and other types of contour
charts on the new machine are fine. The old machine had IDL Version
5.3 and the new machine has IDL Version 6.1.
The contouring and filling is correct, and the graphic size is correct
(the entire image size), but the graphic position and size of just the
contoured area (without the color bar etc) is incorrect. I suppressed
the lat/lon box using xstyle and ystyle=4 in the contour command since
it does not line up with the (correct) lat/lon box from PLOT.
Also, I can't get it to do the contouring, grid lines, and continents
all at once. I can either have it do just the continents and grid
lines, or just the contoured part without continents or grid lines.
When I take longitude and latitude off of the contour commands, the
size changes, but is still wrong. When I try to have contouring and
continents together, I get this error:
MAP_CONTINENTS: Map transform not established.
When I try to have contouring and grid lines, I get this error:
MAP_GRID: map_grid---Current ploting device must have mapping
coordinates
It's as if it is not recognizing MAP_SET?
I can get it to draw solid grid lines using PLOT since MAP_GRID does
not work (although I wanted dotted lines using MAP_GRID).
Any help in being able to solve these problems would be greatly
appreciated. Sometimes I also get an underflow error but not always.
Jill
Here is the code:
MAP_SET, 0, plon0, /NOBORDER, COLOR=drawColor, $
LIMIT=[lat_b,lon_b,lat_e,lon_e], $
POSITION=position
contour, arr, levels=c_levels, $
longitude, latitude, $
/cell_fill, /closed, c_colors=c_colors, $
position=position, $
xstyle=4, ystyle=4, $
min_value=ROUND(min_input), max_value=ROUND(max_input)
contour, arr, levels = c_levels, c_thick = .1,$
longitude, latitude, $
c_labels=c_labels, c_charsize=1.2, $
min_value=ROUND(min_input), max_value=ROUND(max_input), $
position=position, /overplot
; Draw the Continents.
if region EQ 'global50' then MAP_CONTINENTS, /COASTS, COLOR=drawColor
else $
if region EQ 'global100' then MAP_CONTINENTS, /COASTS, COLOR=drawColor
else $
if region EQ 'namerica' then MAP_CONTINENTS, /COASTS, COLOR=drawColor
else $
MAP_CONTINENTS, /CONT, /COUNTRIES, /USA, COLOR=drawColor
; Draw the grid lines.
MAP_GRID, color=drawColor, lons=lon_b, lats=lat_b, londel=londel,
latdel=latdel
PLOT, [lon_b, lon_e],[lat_b, lat_e], /NODATA, $
XStyle=1, YStyle=1, Xticks=num_lon_labels-1,
Yticks=num_lat_labels-1, $
POSITION=position, /NoErase, $
title=title, xtitle=xtitle, Color=drawColor, $
Charsize=char_size, $
YTICKV=[lat_labels], XTICKV=[lon_labels], $
; XTICKLEN=-0.01, YTICKLEN=-0.01
XTICKLEN=1, YTICKLEN=1
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
