;$Id: plater2.pro,v 1.1 2021/12/10 01:50:47 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=7,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
!p.charsize=1.1
!x.margin=[7.8,1.5]
!y.margin=[3.2,0.3]
!x.title='!6'
!p.multi=[0,1,3]
si2=1.0
siz=1.3
;
yr=[1d-56,1d-10]
yr=[1d-56,1d-40]
;
xr=[10.,10240.]
;
it= [0  ,  1, 10, 30, 60,100,200,300]
col=[255,255,255, 55,122,155,188,188]
col=[  0,  0,  0, 55, 55,122,122,155]
yr1=[3d-36,2d-28] & yr2=[1d-48,4d-35]
;
it =[600,700,800,900,1100,1500,3000]
col=[  0,  0,  0, 55,  55, 122, 155]
yr1=[2d-10,7d-4] & yr2=[2d-10,7d-8]
;
!p.multi=[0,2,1]
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k!6)'
plot_oo,xr,yr1,/nodata
for i=0,n_elements(it)-1 do begin
  print,t[it[i]] 
  oplot,col=col[i],k[1:*],spec1[1:*,it[i]]
endfor
;
!y.title='!8k!6!u2!n!8E!6!dGW!n(!8k!6)'
plot_oo,xr,yr2,/nodata
for i=0,n_elements(it)-1 do begin
  print,t[it[i]] 
  oplot,col=col[i],k[1:*],grav1[1:*,it[i]]*k[1:*]^2
endfor
;
print,"$mv idl.ps ~/tex/GW/lowfreq/fig/plater2.eps"
END