;$Id: pol_analytical_nS_-3.67_nA_-3.67.pro,v 1.1 2020/08/14 14:46:56 roper Exp $ if !d.name eq 'PS' then begin ;device,xsize=18,ysize=8,yoffset=3 device,xsize=18,ysize=16,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; dir_ext='nS_-3.67_nA_-3.67/' yr=[-1,1] dir='~/tex/roper/helical/data/' dir=dir+dir_ext k_ast=2*!DPI*100 file=dir+'h01.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] plot_oi,k2,P,yr=yr oplot,k2,P,thick=4,color=60 ;nel = n_elements(k2) ;xyouts,k2[119],P[119],'!6o',col=122 file=dir+'h02.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h03.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h04.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h05.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h06.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h07.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h08.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h09.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P file=dir+'h1.txt' c1=rtable(file,2,head=1) k2=c1[0,*]*k_ast P=c1[1,*] oplot,k2,P,thick=4,col=122 !x.title='!8k!6' !y.title='!13P!6!dGW!n(!8k!6)' xyouts,1e4,-0.4,'!8n!dA!n!6=-11/3!n',siz=siz xyouts,1e4,-0.2,'!8n!dS!n!6=-11/3!n',siz=siz END