if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!x.margin=[6.5,.5]
!y.margin=[3.2,.5]
circ_sym,1.3,1
;file='lucky_idl.dat'
file='lucky.dat'
a=rtable(file,5)
t=a[0,*]
a1=a[3,*]
a2=a[4,*]
;;----erase bad ones from the table----------
;close,1
;openw,1,'lucky_new.dat'
;k=-1
;index_bad=make_array(1,n_elements(a1),/integer,value=0)
;for n=0,n_elements(a1)-2 do begin
;    bad_i=where(a1(n) lt a1(n+1:-1))
;    for m=0,n_elements(bad_i)-1 do begin
;	    if bad_i(m) gt -1 then begin
;	      k=k+1
;        index_bad(k)=n
;      endif
;    endfor
;endfor
;index_bad=index_bad(where(index_bad gt 0))
;index_bad=index_bad(uniq(index_bad))
;index_bad_size=size(index_bad)
;index_bad_size=index_bad_size[0]
;if index_bad_size ne 0 then begin
;  dims=size(a,/dimensions)
;  a_new=a[*, Where(~Histogram(index_bad, MIN=0, MAX=dims[1]-1), /NULL)]
;endif else begin
;  a_new=a
;endelse

;printf,1,a_new
;close,1
;
!x.title='!8t !6[s]'
!y.title='!8r !6[!7l!6m]'
;a1=a_new[3,*]
;a2=a_new[4,*]
plot,t,a1,/nodata,xr=[-20,max(t)+10],yr=[0,max(a1)+5]
oplot,t,a1,ps=4
oplot,t,a2,ps=5,col=55
;
cwd,run
default,texsubdir,''
;print,"mv idl.ps ~/tex/xiangyu/lucky/Figs/"
print,'mv idl.ps ~/tex/xiangyu/lucky/Figs/'+texsubdir+'pplucky_'+run+'.ps'
;
END