pro plot_xpol_amf,amfname,print=print

   if (keyword_set(print)) then begin
      dname = !d.name
      psav = !p
      set_plot,'printer'
      !p.charsize=2
   endif

   amf = readamf(amfname)

   plot,amf[0,*,0,0],amf[1,*,0,0],psym=-1,yran=[0.01,10],/ylog,/xlog,xran=[1,20],/xsty,$
         title='Parallel and Crossed 3C84 Amplitude',ytitle='Amplitude',xtitle='Frequency [GHz]'
   oplot,amf[0,*,0,0],amf[1,*,1,0],psym=-4
   oplot,amf[0,*,0,0],amf[1,*,0,1],psym=-1
   oplot,amf[0,*,0,0],amf[1,*,1,1],psym=-4

   plot,amf[0,*,0,0],amf[1,*,0,0]/amf[1,*,0,1],psym=-1,yran=[0.1,100],/ylog,/xlog,xran=[1,20],/xsty,$
         title='Parallel:Crossed 3C84 Ratio',ytitle='Ratio',xtitle='Frequency [GHz]'
   oplot,amf[0,*,0,0],amf[1,*,1,0]/amf[1,*,1,1],psym=-4

   if (keyword_set(print)) then begin
      device,/close
      set_plot,dname
      !p = psav
   endif
return
end