;$Id: ppol_comp2.pro,v 1.5 2020/11/01 15:07:06 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=14,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
siz=1.0
cwd,run
@parameters
!p.charsize=1.0
!x.margin=[6.8,.5]
!y.margin=[3.2,.5]
!p.multi=[0,2,2]
;
dir1='F1152d2_sig1_t11_M2c_double'
dir2='F1152sig07c'
dir3='F1152sig05a'
dir4='F1152sig03a'
dir5='F1152sig01a'
dir6='F1152sigm01a'
dir7='F1152sigm1a'
dir8='F1152sigm0001a'
dir9='M1152e_exp6k4_M4b_rerun'
;dirB='M1152e_exp6k4_M4b_rerun_sig05'
;dirC='M1152e_exp6k4_M4b_rerun_sig01'
;dirD='M1152e_exp6k4_M4b_rerun_sig03'
;
dirA='M1152e_exp6k4'
dirB='M1152e_exp6k4_sig07'
dirC='M1152e_exp6k4_sig05'
dirD='M1152e_exp6k4_sig03'
dirE='M1152e_exp6k4_sig01'
;dirC='M1152e_exp6k4_sigm01b'
;dirC='M1152e_exp6k4_sig1d' ;(new)
;dirC='M1152e_exp6k4_sig1e' ;(new)
;dirC='M1152e_exp6k4_sig1f' ;(new)
;dirC='M1152e_exp6k4_sig1g' ;(new)
;dirC='F1024sig1t15a'
;
yr=[-1,1]
yr=[-0,1]
;------------------------------------------------------
!x.title='!6'
!y.title='!13P!6!dGW!n(!8k!6)'
restore,'../'+dirA+'/specm.sav'
plot_oi,k[1:*],grav2m[1:*]/grav1m[1:*],yr=yr
;
restore,'../'+dirB+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=55
;
loadct,6
restore,'../'+dirC+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=122
loadct,5
;
restore,'../'+dirD+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=155
;
restore,'../'+dirE+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=122
;
xyouts,siz=siz,2.1e2,0.85,'!7r!6=1'
xyouts,siz=siz,1.0e4,0.90,'!7r!6=0.7',col=55
loadct,6
xyouts,siz=siz,1.0e4,0.78,'!7r!6=0.5',col=122
loadct,5
xyouts,siz=siz,1.0e4,0.55,'!7r!6=0.3',col=155
xyouts,siz=siz,1.0e4,0.28,'!7r!6=0.1',col=122
;-----------------------------------------------------
!y.title='!6'
restore,'../'+dir1+'/specm.sav'
plot_oi,k[1:*],grav2m[1:*]/grav1m[1:*],yr=yr
;
restore,'../'+dir2+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=55
;
loadct,6
restore,'../'+dir3+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=122
loadct,5
;
restore,'../'+dir4+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=155
;
restore,'../'+dir5+'/specm.sav'
oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=122
;
;restore,'../'+dir6+'/specm.sav'
;oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=55
;
;restore,'../'+dir7+'/specm.sav'
;oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=55
;
;restore,'../'+dir8+'/specm.sav'
;oplot,k[1:*],grav2m[1:*]/grav1m[1:*],col=55
;
siz=0.8
xyouts,siz=siz,1.0e4,0.75,'!7r!6=1'
xyouts,siz=siz,7.0e3,0.52,'!7r!6=0.7',col=55
loadct,6
xyouts,siz=siz,1.1e4,0.22,'!7r!6=0.5',col=122
loadct,5
xyouts,siz=siz,6.0e2,0.65,'!7r!6=0.3',col=155
xyouts,siz=siz,8.0e2,0.48,'!7r!6=0.1',col=122
;xyouts,siz=siz,8.0e2,0.05,'!7r!6=0.01'
;xyouts,siz=siz,2.0e3,-.28,'!7r!6=-0.1'
;xyouts,siz=siz,7.0e3,-.65,'!7r!6=-1'
;
;==============================================================================
!x.title='!8k!6'
!y.title='!13P!6!dh!n(!8k!6)'
restore,'../'+dirA+'/specm.sav'
plot_oi,k[1:*],grah2m[1:*]/grah1m[1:*],yr=yr,/nodata
oplot,k[1:*],grah2m[1:*]/grah1m[1:*]
;
restore,'../'+dirB+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=55
;
loadct,6
restore,'../'+dirC+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=122
loadct,5
;
restore,'../'+dirD+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=155
;
restore,'../'+dirE+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=122
;
xyouts,siz=siz,2.3e2,0.85,'!7r!6=1'
xyouts,siz=siz,1.2e4,0.76,'!7r!6=0.7',col=55
loadct,6
xyouts,siz=siz,1.0e4,0.68,'!7r!6=0.5',col=122
loadct,5
xyouts,siz=siz,1.0e4,0.54,'!7r!6=0.3',col=155
xyouts,siz=siz,1.0e4,0.24,'!7r!6=0.1',col=122
;-----------------------------------------------------
!y.title='!6'
restore,'../'+dir1+'/specm.sav'
plot_oi,k[1:*],grah2m[1:*]/grah1m[1:*],yr=yr
;
restore,'../'+dir2+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=55
;
loadct,6
restore,'../'+dir3+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=122
loadct,5
;
restore,'../'+dir4+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=155
;
restore,'../'+dir5+'/specm.sav'
oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=122
;
restore,'../'+dir6+'/specm.sav'
;oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=55
;
restore,'../'+dir7+'/specm.sav'
;oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=55
;
restore,'../'+dir8+'/specm.sav'
;oplot,k[1:*],grah2m[1:*]/grah1m[1:*],col=55
;
siz=0.8
xyouts,siz=siz,6.0e3,0.52,'!7r!6=1'
;xyouts,siz=siz,7.0e3,0.52,'!7r!6=0.7',col=55
loadct,6
xyouts,siz=siz,1.9e2,0.42,'!60.5',col=122
loadct,5
xyouts,siz=siz,5.5e2,0.65,'!7r!6=0.3',col=155
xyouts,siz=siz,8.0e2,0.49,'!7r!6=0.1',col=122
;xyouts,siz=siz,8.0e2,0.05,'!7r!6=0.01'
;xyouts,siz=siz,2.0e3,-.28,'!7r!6=-0.1'
;xyouts,siz=siz,7.0e3,-.65,'!7r!6=-1'
;
print,'$mv idl.ps ~/tex/roper/helical/fig/ppol_comp2.ps'
END
;  sig   grav1mm     pol1     pol2      hel1     hel2   run
;-1.000  2.64e-06 -0.88109 -0.41315  -0.97212 -0.24926  F1152sigm1a
;-0.100  1.88e-06 -0.16662 -0.03742  -0.31665 -0.02251  F1152sigm01a
;-0.010  1.85e-06  0.01473 -0.00396  -0.03393 -0.00278  F1152sigm001a
;-0.001  1.85e-06  0.02997  0.00016  -0.00392  0.00074  F1152sigm0001a
; 0.001  1.86e-06  0.03117  0.00107   0.00275  0.00143  F1152sig0001a
; 0.010  1.86e-06  0.04504  0.00507   0.03255  0.00346  F1152sig001a
; 0.100  1.90e-06  0.22322  0.04499   0.31759  0.02252  F1152sig01a
; 0.300  2.03e-06  0.56261  0.12055   0.71544  0.05321  F1152sig03a
; 0.500  2.30e-06  0.75233  0.16868   0.87613  0.06953  F1152sig05a
; 0.700  2.53e-06  0.86410  0.26999   0.94389  0.16326  F1152sig07c
; 1.000  2.67e-06  0.89864  0.39000   0.96944  0.26901  F1152d2_sig1_t11_M2c_double