;$Id: pcomp_EEGW.pro,v 1.11 2019/11/11 13:03:27 roper Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=20,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; !p.charsize=1.7 !x.margin=[8.2,.5] !y.margin=[3.2,.5] bar='!20!s!A$!n!r!6' bar2='!20!s!u$!n!r!6' ; ; sample to compare time series form different directories ; dir1='M144_k10b_expA' dir2='M144_k10b_exp' ; dir1='M1152e3_exp6A' dir2='M1152e3_exp6' yrEEGW=[.0005,.30] yrhrms=[0,.16] ; dir1='M1152e_exp6k4_M4b' dir2='M1152e_exp6k4' dir3='M1152e_exp6k4_k60b' ;dir4='E288e_t11_M4e' ;(ac1) dir4='E1152e_t11_M4d_double_rep' ;(ac1) ;dir5='F288e_sig1_t11_M4' ;(hel1) dir5='F1152d2_sig1_t11_M2c_double' ;(hel1) ;dir6='F288a_sig1_t11b' ;(hel2) dir6='F1152a_sig1_t11d_double' ;(hel2) ; yrEEGW=[3e-12,1.2e-8] yrEEGW=[3e-12,1.2e-7] yrEEM=[3e-4,1.4e-1] yrhrms=[1e-23,1e-1] xr=[1.,1.27] xr=[1.,1.22] ; pc_read_ts,o=ts1,datadir='../'+dir1+'/data' pc_read_ts,o=ts2,datadir='../'+dir2+'/data' pc_read_ts,o=ts3,datadir='../'+dir3+'/data' pc_read_ts,o=ts4,datadir='../'+dir4+'/data' pc_read_ts,o=ts5,datadir='../'+dir5+'/data' pc_read_ts,o=ts6,datadir='../'+dir6+'/data' ; !p.multi=[0,1,2] !x.range=[1,10.4] ; !x.title='!6' ;!y.title='!13E!6!dM!n/!13E!6!drad!n' !y.title='!7X!6!dM,K!n' plot_io,ts1.t,ts1.EEM,xr=xr,yr=yrEEM oplot,ts2.t,ts2.EEM,col=122 oplot,ts3.t,ts3.EEM,col=55 oplot,ts4.t,ts4.EEK,li=2 oplot,ts5.t,ts5.EEM,li=2,col=122 oplot,ts6.t,ts6.EEM,li=2,col=55 ; xyouts,1.050,2.3e-2,'!6ini1',siz=siz xyouts,1.110,2.8e-3,'!6ini2',siz=siz,col=122 xyouts,1.005,8.0e-3,'!6ini3',siz=siz,col=55 ; xyouts,1.17,1.7e-3,'!6ac1',siz=siz xyouts,1.13,1.3e-2,'!6hel1',siz=siz,col=122 xyouts,1.02,1.0e-3,'!6hel2',siz=siz,col=55 ; ;!x.title='!8t!6' ;!y.title='!8h!6!drms!n' ;plot_io,ts1.t,ts1.hrms,xr=xr,yr=yrhrms ;oplot,ts2.t,ts2.hrms,col=122,li=2 ;oplot,ts3.t,ts3.hrms,col=55,li=2 ;oplot,ts4.t,ts4.hrms,col=155,li=2 ; ;!x.title=bar2+'!8t!6' !x.title='!8t!6' ;!y.title='!13E!6!dGW!n/!13E!6!drad!n' !y.title='!7X!6!dGW!n' plot_io,ts1.t,ts1.EEGW,xr=xr,yr=yrEEGW oplot,ts2.t,ts2.EEGW,col=122 oplot,ts3.t,ts3.EEGW,col=55 oplot,ts4.t,ts4.EEGW,li=2 oplot,ts5.t,ts5.EEGW,li=2,col=122 oplot,ts6.t,ts6.EEGW,li=2,col=55 ; xyouts,1.04,2.5e-9,'!6ini1',siz=siz xyouts,1.06,8.5e-12,'!6ini2',siz=siz,col=122 xyouts,1.06,8.0e-10,'!6ini3',siz=siz,col=55 ; xyouts,1.16,3.0e-8,'!6ac1',siz=siz xyouts,1.16,6.3e-9,'!6hel1',siz=siz,col=122 xyouts,1.08,1.2e-10,'!6hel2',siz=siz,col=55 ; print print,"$mv idl.ps ~/tex/roper/GW/fig/pcomp_EEGW.ps" print END