;$Id: pe_comp_above.pro,v 1.2 2020/11/11 03:30:58 brandenb Exp $ if !d.name eq 'PS' then begin device,xsize=18,ysize=16,yoffset=3 !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3 end ; siz=1.7 !p.charsize=1.7 !x.margin=[8.3,.5] !y.margin=[3.2,.5] ; dir1='F1152a_sig0_t11_M4_ramp01b' & li1=1 & col1=122 dir2='F1152a_sig0_t11_M4_ramp02a' & li2=2 & col2=122 dir3='F1152a_sig0_t11_M4_ramp05a' & li3=3 & col3=55 dir4='F1152a_sig0_t11_M4_ramp1a' & li4=4 & col4=122 dir5='F1152a_sig0_t11_M4_ramp2a' & li5=0 ; taver5=2.2 ; default,iread,0 if iread eq 0 then begin 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' iread=1 endif ; default,kin,0 default, r_pe,[1.0,1.64] default,taver1,1.15 default,kf,600. default,xr_pe,1. circ_sym,1.3,1 fo="(a,e8.2,e10.2,i3,f6.2,i5,2x,a)" !p.multi=[0,1,2] ; t1a=1.14 & t1b=1.27 & t1c=1.2 t2a=1.14 & t2b=1.32 & t2c=1.3 t3a=1.14 & t3b=1.4 & t3c=1.5 t4a=1.14 & t4b=1.5 & t4c=1.7 t5a=1.14 & t5b=1.7 & t5c=2.2 ; ;----------------------------------------------------------------------------- ; panel 1 ; !x.title='!6' !y.title='!7X!6!dM!n(!8t!6)' xr_pe=[1.,3.24] yr=[0.,.043] yr=[1e-3,.1] th=6 ; run=dir1 if kin then EEKM=ts1.EEK else EEKM=ts1.EEM ;plot,ts1.t,EEKM,yr=yr,xr=xr_pe,li=li1 plot_io,ts1.t,EEKM,yr=yr,xr=xr_pe,li=li1 imax=(where(EEKM eq max(EEKM)))[0] n2=n_elements(EEKM)-1 EEKMmax=mean(EEKM[(imax-1)>0:(imax+1)0:(imax+1)0:(imax+1)0:(imax+1)0:(imax+1) 1) ',EEKMmax,EEGWm,kin,tmax,kf,run,fo=fo good=where(ts1.t ge t1c) & tgood=ts1.t(good) pp1=linfit(tgood,ts1.EEGW(good)) & mm1=mean(ts1.EEGW(good)) ;oplot,tgood,s*(tgood*pp1[1]+pp1[0]),li=li1,th=th ; loadct,6 oplot,ts2.t,s*ts2.EEGW,li=li2,col=col2 good=where(ts2.t ge taver1) EEGWm=mean(ts2.EEGW(good)) ;oplot,ts2.t(good),ts2.t(good)*0.+EEGWm,col=122 print,'--> 2) ',EEKMmax,EEGWm,kin,tmax,kf,run,fo=fo good=where(ts2.t ge t2c) & tgood=ts2.t(good) pp2=linfit(tgood,ts2.EEGW(good)) & mm2=mean(ts2.EEGW(good)) ;oplot,tgood,s*(tgood*pp2[1]+pp2[0]),li=li2,col=col2,th=th loadct,5 ; oplot,ts3.t,s*ts3.EEGW,li=li3,col=col3 good=where(ts3.t ge taver1) EEGWm=mean(ts3.EEGW(good)) ;oplot,ts3.t(good),ts3.t(good)*0.+EEGWm,col=122 print,'--> 3) ',EEKMmax,EEGWm,kin,tmax,kf,run,fo=fo good=where(ts3.t ge t3c) & tgood=ts3.t(good) pp3=linfit(tgood,ts3.EEGW(good)) & mm3=mean(ts3.EEGW(good)) ;oplot,tgood,s*(tgood*pp3[1]+pp3[0]),li=li3,col=col3,th=th ; oplot,ts4.t,s*ts4.EEGW,li=li4,col=col4 good=where(ts4.t ge taver1) EEGWm=mean(ts4.EEGW(good)) ;oplot,ts4.t(good),ts4.t(good)*0.+EEGWm,col=122 print,'--> 4) ',EEKMmax,EEGWm,kin,tmax,kf,run,fo=fo good=where(ts4.t ge t4c) & tgood=ts4.t(good) pp4=linfit(tgood,ts4.EEGW(good)) & mm4=mean(ts4.EEGW(good)) ;oplot,tgood,s*(tgood*pp4[1]+pp4[0]),li=li4,col=col4,th=th ; oplot,ts5.t,s*ts5.EEGW,li=li5 good=where(ts5.t ge taver5) EEGWm=mean(ts5.EEGW(good)) oplot,ts5.t(good),ts5.t(good)*0.+s*EEGWm,col=122 print,'--> 5) ',EEKMmax,EEGWm,kin,tmax,kf,run,fo=fo good=where(ts5.t ge t5c) & tgood=ts5.t(good) pp5=linfit(tgood,ts5.EEGW(good)) & mm5=mean(ts5.EEGW(good)) ;oplot,tgood,s*(tgood*pp5[1]+pp5[0]),li=li5,col=col5,th=th ; xyouts,1.07,.225,siz=siz,'!6(b)' ; !p.multi=0 print,'$convert idl.ps /D/Print/Tuusula_GW/pe_comp.png' print,'$mv idl.ps ~/tex/tina/duration/fig/pe_comp.ps' print,'$mv idl.ps ~/tex/tina/CircPol/fig/pe_comp.ps' ; print,p1 print,p2 print,p3 print,p4 print,p5 print,mm1,pp1 print,mm2,pp2 print,mm3,pp3 print,mm4,pp4 print,mm5,pp5 END