LegendContourPlot[exp_, rangex_, rangey_, opts___]::usage =
"LegendContourPlot[exp,rangex,rangey,opts] Does a contour Plot with \
a colour bar a la IDL";
LegendContourPlot[exp_, rangex_, rangey_, opts___] :=
With[{func = extractPureFunction[exp], size = 300,
cf = ColorFunction /. FilterRules[{opts}, Options[ContourPlot]] /.
ColorFunction -> Automatic},
Row[Show[#, ImageSize -> {Automatic, 300},
ImagePadding -> 20] & /@ {#1,
colorbar[Through[{Min, Max}[#2]], cf]} & @@
Reap[ContourPlot[
func[x, y], {x, rangex[[2]], rangex[[3]]}, {y, rangey[[2]],