; scatter plots
restore,'../polarization_comparisons_maps_10642.sav'
set_plot,'ps'
device,file='PAPER_polarization_model_comparisons_10642.eps',/encap,/color,/times,xsi=16,ysi=7

;!p.multi=[0,2,2]
!p.multi=[0,2,1]
!y.margin=[2,2]
!x.margin=[2,1]
!y.omargin=[2,0]
!x.omargin=[8,8]
nolab=[' ',' ',' ',' ',' ',' ',' ',' ',' ']
loadct,13
; maximum polarization, as f(mu) or something...too much to show!

; NOTE this is for wavelength-normalized polarization signal, 
; for final paper, remove titles. for now,helps recalling what's being plotted.

plot,MODELPcircMAX[*,*,5],CONVMODELPcircMAX[*,*,5],psym=4,xran=[0,0.85],yran=[0,0.45],/xst,/yst,$
   xtit='Max. Polariz. Fraction (Orig.)',ytit='Max. Polariz. Fraction (Model IVM)';,tit='Max Fract Circ'
   oplot,MODELPcircmax[*,*,1],CONVMODELPcircMAX[*,*,1],psym=4,col=0
for i=0,5 do begin  
   oplot,MODELPcircmax[*,*,5],IVMMODELPcircmax[*,*,5,i],psym=6,col=50+40*i,symsi=0.3
   oplot,MODELPcircmax[*,*,1],IVMMODELPcircmax[*,*,1,i],psym=6,col=50+40*i,symsi=0.3
endfor
oplot,5000*indgen(2),5000*indgen(2),lines=1

plot,MODELPlinMAX[*,*,1],CONVMODELPlinMAX[*,*,1],psym=4,xran=[0,0.4],yran=[0,0.15],/xst,/yst,$
	xtit='Max. Polariz Fraction (Orig.)',ytickn=nolab
	;,ytit='Convolved/IVM-sampled';,tit='Max Fract Linear'
oplot,MODELPlinMAX[*,*,5],CONVMODELPlinMAX[*,*,5],psym=4,col=0
for i=0,5 do begin 
   oplot,MODELPlinmax[*,*,5],IVMMODELPlinmax[*,*,5,i],psym=6,col=50+40*i,symsi=0.3
   oplot,MODELPlinmax[*,*,1],IVMMODELPlinmax[*,*,1,i],psym=6,col=50+40*i,symsi=0.3
endfor
oplot,5000*indgen(2),5000*indgen(2),lines=1
axis,yaxis=1,ytit='Max. Polariz. Fraction (Model IVM)'

; for clarity for the moment, only do max polarization...

;plot,MODELPcirctot[*,*,5],CONVMODELPcirctot[*,*,5],psym=4,xran=[0,0.35],yran=[0,0.35],/xst,/yst,$
;	xtit='Model',ytit='Convolved/IVM-sampled',tit='Integr Fract Circ'
;for i=0,5 do begin $ 
;oplot,MODELPcirctot[*,*,5],IVMMODELPcirctot[*,*,5,i],psym=6,col=50+40*i,symsi=0.3
;;oplot,MODELPcirctot[*,*,5],CONVMODELPcirctot[*,*,5],psym=4,col=0,symsi=0.1
;oplot,5000*indgen(2),5000*indgen(2)
;
;plot,MODELPlintot[*,*,5],CONVMODELPlintot[*,*,5],psym=4,xran=[0,0.015],yran=[0,0.015],/xst,/yst,$
;	xtit='Model',ytit='Convolved/IVM-sampled',tit='Integr Fract Linear'
;for i=0,5 do begin $ 
;oplot,MODELPlintot[*,*,5],IVMMODELPlintot[*,*,5,i],psym=6,col=50+40*i,symsi=0.3
;;oplot,MODELPlintot[*,*,5],CONVMODELPlintot[*,*,5],psym=4,col=0,symsi=0.1
;oplot,5000*indgen(2),5000*indgen(2)

device,/close
set_plot,'x'
reset

end