;$Id: pearly2.pro,v 1.1 2021/12/10 17:30:04 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=[8.8,0.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 =[600,700,800,900,1000,1200,1400]
col=[255,255,255, 55,122 ,155 ,188]
yr1=[1d-12,2d-2]
yr2=[1d-23,2d-10]
;
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-42,4d-29]
;
!p.multi=[0,2,1]
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k!6)/!8k!6!u3.9!n'
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]]/k[1:*]^3.9
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/pearly2.eps"
END