; read in the highest resolution image first
img_vla=readfits('../J1428p3526.fits',hdr_vla)
CX=281 ; center of source in X dimension
CY=282 ; center of source in Y dimension
HW=100 ; half width of box to extract

HEXTRACT,img_vla,hdr_vla,CX-HW,CX+HW,CY-HW,CY+HW

img_350=readfits('../../059/set_16301_smooth.fits',hdr_350) ; read in sharc image
hastrom,img_350,hdr_350,img_350p,hdr_350p,hdr_vla,MISSING=0 ; make image match size/shape of VLA

loadct,0                                                    ; black and white color table
gamma_ct,1.0,/CURRENT                                       ; normal gamma stretch
img_350p=bytscl(img_350p,min=-0.05,max=0.07)                ; choose plotting range
set_plot, 'ps'
device , filename='sharc.eps', /encap, xsize=10., ysize = 10., $
         yoffset = 1., BITS_PER_PIXEL=8, COLOR=1
imcontour, img_350p, hdr_350p, levels=0, xtitle=' ', ytitle=' ',$ ; plot the axis
           charsize=1.2,charthick=3,/nodata,subtitle=' ',COLOR=0,$
           XSTYLE=4,YSTYLE=4
tvimage,img_350p,/overplot,/keep_aspect; display the image

  ; overlay radio contours
imcontour, img_vla, hdr_vla, levels=[2e-4,3e-4,4e-4],xtitle=' ', ytitle=' ',$
  charsize=1.2,charthick=3,/noerase,subtitle=' ',C_THICK=3,C_COLOR=0,COLOR=0, $
  XTHICK=2,YTHICK=2

XYOUTS,10,10,'350 + VLA contours',charsize=1.5,charthick=3
DEVICE, /close