;$Id: ppol_comp3.pro,v 1.3 2021/01/21 01:33:31 roper 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' ; dir_ext1='nS_-3.67_nA_-3.67/' dir_ext2='nS_-3.67_nA_-4.67/' dir_ext3='nS_-4.33_nA_-4.33/' dir='~/tex/roper/helical/data/' dir_ext1=dir+dir_ext1 dir_ext2=dir+dir_ext2 dir_ext3=dir+dir_ext3 yr=[-1,1] yr=[-0,1.2] ; -- -- -- ; -- -- -- panel for initial field case (plotting P=\XiGW/\OmGW) ; -- -- -- ;------------------------------------------------------ int_k0 = 1000 int_k1 = 3000 !x.title='!6' !y.title='!13P!6!dGW!n(!8k!6)' restore,'../'+dirA+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] polA=int_tabulated(kaux,polaux) plot_oi,k[1:*],pol2m,yr=yr ; restore,'../'+dirB+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] polB=int_tabulated(kaux,polaux) oplot,k[1:*],pol2m,col=55 ; loadct,6 restore,'../'+dirC+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] polC=int_tabulated(kaux,polaux) oplot,k[1:*],pol2m,col=122 loadct,5 ; restore,'../'+dirD+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] polD=int_tabulated(kaux,polaux) oplot,k[1:*],pol2m,col=155 ; restore,'../'+dirE+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] polE=int_tabulated(kaux,polaux) oplot,k[1:*],pol2m,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 dir_ext = dir_ext3 ; Analytical (compensated) for initial field case (plotting P=\XiGW/\OmGW) thick_an=.5 k_ast=2*!DPI*100 file=dir_ext+'h01.txt' c1=rtable(file,2,head=1) fack01=1. k2=c1[0,*]*k_ast/fack01 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol01=int_tabulated(k2aux,P2aux) facP01=polE/pol01 oplot,k2,P*facP01,thick=thick_an,color=60 file=dir_ext+'h03.txt' c1=rtable(file,2,head=1) fack03=1.15 k2=c1[0,*]*k_ast/fack03 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol03=int_tabulated(k2aux,P2aux) facP03=polD/pol03 oplot,k2,P*facP03,thick=thick_an,color=60 file=dir_ext+'h05.txt' c1=rtable(file,2,head=1) fack05=1.25 k2=c1[0,*]*k_ast/fack05 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol05=int_tabulated(k2aux,P2aux) facP05=polC/pol05 oplot,k2,P*facP05,thick=thick_an,color=60 file=dir_ext+'h07.txt' c1=rtable(file,2,head=1) fack07=1.35 k2=c1[0,*]*k_ast/fack07 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol07=int_tabulated(k2aux,P2aux) facP07=polB/pol07 oplot,k2,P*facP07,thick=thick_an,color=60 file=dir_ext+'h1.txt' c1=rtable(file,2,head=1) fack1=1.5 k2=c1[0,*]*k_ast/fack1 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol1=int_tabulated(k2aux,P2aux) facP1=1.04*polA/pol1 oplot,k2,P*facP1,thick=thick_an,color=60 print,'values for initial field:' print,'factors k: ',fack01,fack03,fack05,fack07,fack1 print,'factors P: ',facP01,facP03,facP05,facP07,facP1 ; -- -- -- ; -- -- -- panel for driven field case (plotting P=\XiGW/\OmGW) ; -- -- -- ;----------------------------------------------------- int_k0 = 500 int_k1 = 2000 !y.title='!6' restore,'../'+dir1+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] plot_oi,k[1:*],pol2m,yr=yr good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] poldir1=int_tabulated(kaux,polaux) ; restore,'../'+dir2+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] oplot,k[1:*],pol2m,col=55 good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] poldir2=int_tabulated(kaux,polaux) ; loadct,6 restore,'../'+dir3+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] oplot,k[1:*],pol2m,col=122 good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] poldir3=int_tabulated(kaux,polaux) loadct,5 ; restore,'../'+dir4+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] oplot,k[1:*],pol2m,col=155 good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] poldir4=int_tabulated(kaux,polaux) ; restore,'../'+dir5+'/specm.sav' pol2m=grav2m[1:*]/grav1m[1:*] oplot,k[1:*],pol2m,col=122 good=where(k[1:*] GT int_k0) kaux=k[good] polaux=pol2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polaux=polaux[good] poldir5=int_tabulated(kaux,polaux) ; ;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' ; panel for analytical (compensated) driven field case (plotting P=\XiGW/\OmGW) dir_ext = dir_ext2 k_ast=2*!DPI*100 file=dir_ext+'h01.txt' c1=rtable(file,2,head=1) fack01_dir=1.2 k2=c1[0,*]*k_ast/fack01_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol01=int_tabulated(k2aux,P2aux) facP01_dir=1.03*poldir5/pol01 oplot,k2,P*facP01_dir,thick=thick_an,color=60 file=dir_ext+'h03.txt' c1=rtable(file,2,head=1) fack03_dir=1.2 k2=c1[0,*]*k_ast/fack03_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol03=int_tabulated(k2aux,P2aux) facP03_dir=poldir4/pol03 oplot,k2,P*facP03_dir,thick=thick_an,color=60 file=dir_ext+'h05.txt' c1=rtable(file,2,head=1) fack05_dir=1.2 k2=c1[0,*]*k_ast/fack05_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol05=int_tabulated(k2aux,P2aux) facP05_dir=.95*poldir3/pol05 oplot,k2,P*facP05_dir,thick=thick_an,color=60 file=dir_ext+'h07.txt' c1=rtable(file,2,head=1) fack07_dir=1.2 k2=c1[0,*]*k_ast/fack07_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol07=int_tabulated(k2aux,P2aux) facP07_dir=.9*poldir2/pol07 oplot,k2,P*facP07_dir,thick=thick_an,color=60 file=dir_ext+'h1.txt' c1=rtable(file,2,head=1) fack1_dir=1.2 k2=c1[0,*]*k_ast/fack1_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] pol1=int_tabulated(k2aux,P2aux) facP1_dir=.9*poldir1/pol1 oplot,k2,P*facP1_dir,thick=thick_an,color=60 print,'values for driven field:' print,'factors k: ',fack01_dir,fack03_dir,fack05_dir,fack07_dir,fack1_dir print,'factors P: ',facP01_dir,facP03_dir,facP05_dir,facP07_dir,facP1_dir ; ;============================================================================== ; -- -- -- ; -- -- -- panel for initial field case (plotting P_h=Ah/Sh) ; -- -- -- int_k0 = 1000 int_k1 = 3000 !x.title='!8k!6' !y.title='!13P!6!dh!n(!8k!6)' restore,'../'+dirA+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhA=int_tabulated(kaux,polhaux) plot_oi,k[1:*],polh2m,yr=yr,/nodata oplot,k[1:*],polh2m ; restore,'../'+dirB+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhB=int_tabulated(kaux,polhaux) oplot,k[1:*],polh2m,col=55 ; loadct,6 restore,'../'+dirC+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhC=int_tabulated(kaux,polhaux) oplot,k[1:*],polh2m,col=122 loadct,5 ; restore,'../'+dirD+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhD=int_tabulated(kaux,polhaux) oplot,k[1:*],polh2m,col=155 ; restore,'../'+dirE+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhE=int_tabulated(kaux,polhaux) oplot,k[1:*],polh2m,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 ; -- -- -- analytical (compensated) initial field case (plotting P_h=Ah/Sh) dir_ext = dir_ext3 k_ast=2*!DPI*100 file=dir_ext+'h01.txt' c1=rtable(file,2,head=1) fackh01=1.1 k2=c1[0,*]*k_ast/fackh01 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh01=int_tabulated(k2aux,P2aux) facPh01=polhE/polh01 oplot,k2,P*facPh01,thick=thick_an,color=60 file=dir_ext+'h03.txt' c1=rtable(file,2,head=1) fackh03=1.5 k2=c1[0,*]*k_ast/fackh03 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh03=int_tabulated(k2aux,P2aux) facPh03=polhD/polh03*1.04 oplot,k2,P*facPh03,thick=thick_an,color=60 file=dir_ext+'h05.txt' c1=rtable(file,2,head=1) fackh05=1.5 k2=c1[0,*]*k_ast/fackh05 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh05=int_tabulated(k2aux,P2aux) facPh05=1.03*polhC/polh05 oplot,k2,P*facPh05,thick=thick_an,color=60 file=dir_ext+'h07.txt' c1=rtable(file,2,head=1) fackh07=1.5 k2=c1[0,*]*k_ast/fackh07 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh07=int_tabulated(k2aux,P2aux) facPh07=1.03*polhB/polh07 oplot,k2,P*facPh07,thick=thick_an,color=60 file=dir_ext+'h1.txt' c1=rtable(file,2,head=1) fackh1=1.5 k2=c1[0,*]*k_ast/fackh1 P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh1=int_tabulated(k2aux,P2aux) facPh1=1.06*polhA/polh1 oplot,k2,P*facPh1,thick=thick_an,color=60 print,'values for initial field:' print,'factors kh: ',fackh01,fackh03,fackh05,fackh07,fackh1 print,'factors Ph: ',facPh01,facPh03,facPh05,facPh07,facPh1 ;----------------------------------------------------- ; -- -- -- ; -- -- -- panel for driven field case (plotting P_h=Ah/Sh) ; -- -- -- int_k0 = 500 int_k1 = 2000 !y.title='!6' restore,'../'+dir1+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] plot_oi,k[1:*],polh2m,yr=yr good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhdir1=int_tabulated(kaux,polhaux) ; restore,'../'+dir2+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] oplot,k[1:*],polh2m,col=55 good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhdir2=int_tabulated(kaux,polhaux) ; loadct,6 restore,'../'+dir3+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] oplot,k[1:*],polh2m,col=122 good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhdir3=int_tabulated(kaux,polhaux) loadct,5 ; restore,'../'+dir4+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] oplot,k[1:*],polh2m,col=155 good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhdir4=int_tabulated(kaux,polhaux) ; restore,'../'+dir5+'/specm.sav' polh2m=grah2m[1:*]/grah1m[1:*] oplot,k[1:*],polh2m,col=122 good=where(k[1:*] GT int_k0) kaux=k[good] polhaux=polh2m[good] good=where(kaux[1:*] LT int_k1) kaux=kaux[good] polhaux=polhaux[good] polhdir5=int_tabulated(kaux,polhaux) ; 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' ; -- -- -- analytical (compensated) driven field case (plotting P_h=Ah/Sh) dir_ext = dir_ext2 k_ast=2*!DPI*100 file=dir_ext+'h01.txt' c1=rtable(file,2,head=1) fackh01_dir=1.15 k2=c1[0,*]*k_ast/fackh01_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh01=int_tabulated(k2aux,P2aux) facPh01_dir=polhdir5/polh01 oplot,k2,.96*P*facPh01_dir,thick=thick_an,color=60 file=dir_ext+'h03.txt' c1=rtable(file,2,head=1) fackh03_dir=1.2 k2=c1[0,*]*k_ast/fackh03_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh03=int_tabulated(k2aux,P2aux) facPh03_dir=.93*polhdir4/polh03 oplot,k2,P*facPh03_dir,thick=thick_an,color=60 file=dir_ext+'h05.txt' c1=rtable(file,2,head=1) fackh05_dir=1.2 k2=c1[0,*]*k_ast/fackh05_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh05=int_tabulated(k2aux,P2aux) facPh05_dir=.91*polhdir3/polh05 oplot,k2,P*facPh05_dir,thick=thick_an,color=60 file=dir_ext+'h07.txt' c1=rtable(file,2,head=1) fackh07_dir=1.2 k2=c1[0,*]*k_ast/fackh07_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh07=int_tabulated(k2aux,P2aux) facPh07_dir=.91*polhdir2/polh07 oplot,k2,P*facPh07_dir,thick=thick_an,color=60 file=dir_ext+'h1.txt' c1=rtable(file,2,head=1) fackh1_dir=1.3 k2=c1[0,*]*k_ast/fackh1_dir P=c1[1,*] good=where(k2 GT int_k0) k2aux=k2[good] P2aux=P[good] good=where(k2aux LT int_k1) k2aux=k2aux[good] k2aux=k2aux[rem_dup(k2aux)] P2aux=P2aux[good] P2aux=P2aux[rem_dup(k2aux)] polh1=int_tabulated(k2aux,P2aux) facPh1_dir=.95*polhdir1/polh1 oplot,k2,P*facPh1_dir,thick=thick_an,color=60 print,'values for driven field:' print,'factors kh: ',fackh01_dir,fackh03_dir,fackh05_dir,fackh07_dir,fackh1_dir print,'factors Ph: ',facPh01_dir,facPh03_dir,facPh05_dir,facPh07_dir,facPh1_dir ; 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