;$Id: pom.pro,v 1.3 2020/11/21 15:50:50 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=10,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
siz=1.3
!p.charsize=1.7
!x.margin=[7.2,0.5]
!y.margin=[3.2,0.2]
;!x.title='!13E!6!d!8i!6!n/!13E!6!drad!n'
;!x.title='!7X!6!d!8i!6!n'
;
!p.multi=[0,2,1]
t1=1.4 & t2=2.
t1=1.102 & t2=2.
t1=1.002 & t2=1.1
t1=1.002 & t2=1.02
t1=1.0002 & t2=1.2
good=where(ts.t ge t1 and ts.t le t2)
nt=n_elements(good)
tmax=max(ts.t(good))
tmin=min(ts.t(good))
dtm=(tmax-tmin)/(nt-1)
tm=tmin+dtm*findgen(nt)
;
;plot,ts.t,ts.dt
;oplot,tm,ts.dt(good),col=55
;help,tm,ts.dt(good)
;oplot,ts.t(good),ts.t(good)*0+mean(ts.dt(good)),col=122
;
!x.title='!8t!6'
;fk=complex(ts.STrep2,ts.STimp2)
;fk=complex(ts.STrept,ts.STimpt)
;fk=complex(ts.hTrept,ts.hTimpt)
fk=complex(ts.gTrept,ts.gTimpt)
xr=[1,3] & yr=[0.,.0004]
xr=[1,1.1] & yr=[0.,.006]
xr=[1,1.3] & yr=[0.,.03]
xr=[1,1.3] & yr=[0.,3e-7] ;(h)
xr=[1,1.3] & yr=[0.,2e-5] ;(g)
plot,ts.t,abs(fk),xr=xr,yr=yr
oplot,ts.t(good),ts.t(good)*0+mean(abs(fk(good))),col=122
fko=fft(fk(good))
TT=nt*dtm
dom=2*!pi/TT
om=dom*findgen(nt)
!x.title='!7x!6'
;plot_oi,om,abs(fko),xr=[dom,1e4],yr=[3e-8,1.8e-5]
yr=[1e-7,6e-5]
yr=[3e-6,3e-3]
yr=[3e-11,6e-8] ;(h)
yr=[3e-9,6e-5] ;(g)
plot_oo,om,abs(fko),xr=[dom,6e3],yr=yr,/nodata
;plot_oo,om,abs(fko),xr=[dom,6e3],/nodata
oplot,om,abs(fko),col=55
;plot,om,abs(fko),xr=[dom,2e3]
help,om,abs(fko)
;
;-----------------------------------------------------------------------------
fk=complex(ts.STrep2,ts.STimp2)
fk=complex(ts.hTrep2,ts.hTimp2)
fk=complex(ts.gTrep2,ts.gTimp2)
;plot,ts.t,abs(fk)
fko=fft(fk(good))
TT=nt*dtm
dom=2*!pi/TT
om=dom*findgen(nt)
oplot,om,abs(fko),col=122
;plot,om,abs(fko),xr=[dom,2e3]
help,om,abs(fko)
xx=[50.,5000.] & oplot,xx,2e-3/xx^1
;
oplot,[1,1]*300*sqrt(3.),yr,li=1,col=55
oplot,[1,1]*700*sqrt(3.),yr,li=1,col=122
;
cwd,run
print,'$mv idl.ps ~/tex/rei/GW/fig/pom_'+run+'.ps'
END