;$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