On Apr 4, 4:26 pm, David Fanning <n...@idlcoyote.com> wrote:
> I typically use images instead of filled contours, too, in these
> situations, but I am very curious now to see the commands you
> used to create this PostScript file. I think it is a combination
> of a filled contour on a map projection that is causing the
> problem.
>
> I am curious to see if my solution to the gridding problem I
> am working on is also a solution to Matt's problem.
These are the relevant lines from my map plotting routine. The plot
above uses the version with the 'levels' set, but the scratches appear
whether the levels are specified or not:
ps_start, filename=filename
!p.multi=[0, 2, 1]
MAP_SET, limit=limit, position=[0.02, 0.05, 0.8, 0.95], $
title=title, charsize=charsize, /noborder, /isotropic
ctload, ct, brewer=brewer, reverse=reverse_ct
if keyword_set(levels) then $
fsc_contour, data, lon, lat, /overplot, /cell_fill, levels=levels,
missingvalue=missingvalue $
else $
fsc_contour, data, lon, lat, /overplot, /cell_fill,
nlevels=nlevels
ctload, 0
map_continents, thick=3
ctload, ct, brewer=brewer, reverse=reverse_ct
if keyword_set(levels) then $
fsc_colorbar, position=[0.93, 0.05, 0.98, 0.95], $
minrange=min(levels), maxrange=max(levels), /vertical,
$ ;ticknames=string(levels, format='(e8.1)')
title=cbtitle, divisions=n_elements(levels)-1, ticknames=labels
$
else $
fsc_colorbar, position=[0.93, 0.05, 0.98, 0.95], $
minrange=min(data, /nan), maxrange=max(data, /nan), /vertical, $
title=cbtitle
fsc_text, 0.05, 0.9, plot_label, /normal, charsize=1.5
ps_end
Cheers,
Matt
|