;$Id: prelhel_comp.pro,v 1.4 2010/06/20 15:24:08 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=12,ysize=20,yoffset=3
  device,xsize=12,ysize=17,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/arakel/kinhel/fig/prelhel_comp_standard.ps
;
!p.charsize=2.2
s=.14
;
;  The safest way to detime t1 (called t0 in the paper)
;  is to cut and paste this value from time_series.dat after the #.
;
dir1='64b_pert1' & t1=5755.57  & tit1='!6Re=40'
dir2='64a_pert1' & t2=6342.69  & tit2='!6Re=75'
;dir3='128b_pert1'& t3=516.28   & tit3='!6Re=135'
dir3='128b_pert2'& t3=1159.57  & tit3='!6Re=135'
;dir4='256b_pert2' & t4=1308.41  & tit4='!6Re=400'
;
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'
;
!x.title='!6'
!x.range=[-50,100]
!y.range=[-.3,.7]
!p.multi=[0,1,4]
!p.multi=[0,1,3]
plot,s*(ts1.t-t1),ts1.oum/ts1.urms^2,tit=tit1
oplot,!x.range,[0,0],li=3
oplot,[0,0],[-1,1],li=1

plot,s*(ts2.t-t2),ts2.oum/ts2.urms^2,tit=tit2
oplot,!x.range,[0,0],li=3
oplot,[0,0],[-1,1],li=1
;
;!y.title='                       '+ $
!y.title='                                                       '+ $
  '!6<!7x!8.u!6>/!8k!6!df!n<!8u!6!u2!n>'
!x.title='(!8t-t!6!d*!n) !8u!6!drms!n!8k!6!df!n'
plot,s*(ts3.t-t3),ts3.oum/ts3.urms^2,tit=tit3
oplot,!x.range,[0,0],li=3
oplot,[0,0],[-1,1],li=1
;
!y.title='!6'
!x.title='(!8t-t!6!d*!n) !8u!6!drms!n!8k!6!df!n'
;plot,s*(ts4.t-t4),ts4.oum/ts4.urms^2,tit=tit4
;oplot,!x.range,[0,0],li=3
;oplot,[0,0],[-1,1],li=1
;
END