# make postscript file with power-2 contours starting at $3
# in map.  Overlay the model on this plot.   
#

mapplot << eof
  style="g"
  mapfil = $2
  xrange 143, 589
  yrange 107, 556
  PLOTFIL "$1.ps/vcps"
  black = -0.015 white = 0.05
  palette=2
  ltype=1
  wedge
  title "",""
  annot = -1
  BEAM
  CHARSIZE=1.2,1.0
  colors = 1,0,0
  Z=0
  /
  overlay
  MAPFIL $1
  style ="c"
  xrange 143, 589
  yrange 107, 556
  black = 0 white = 1
  CHARSIZE=1.2,1.0
  colors = 1,0,0
  TITLE="",""
  topleft="Hydra A"
  topright="Radio + X-ray"
  ltype=4
  linewidth=3
  ANNOT=-1
  opaque
  BEAM
  levs=-1,1,2,4,8,16,32,64,128,256,512,1024,2048,4096
  plev=0
  clev=$3
  H0=75
  Z=0.05424
  BARLENGT=100000
  unit=arcmin
  /
eof