;$Id: plater.pro,v 1.3 2021/12/11 21:56:39 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
;
;  to read data, run first ppower_all.
;
!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
;
;  For Run B1
;
xr=[10.,10240.]
it =[600,700,800,900,1100,1500,3000]
col=[  0,  0,  0, 55,  55, 122, 155]
yr1=[2d-10,7d-4] & yr2=[1d-23,7d-12]
;
;  For Run B10
;
xr=[100.,25600.]
it =[80,100,200,600,850,900]
col=[  0,  0,  0, 55,  55, 122, 155]
yr1=[2d-10,1d-2] & yr2=[1d-23,7d-13]
;
!p.multi=[0,2,1]
!x.title='!8k!6'
!y.title='!8E!6!dM!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='!8E!6!dGW!n(!8k!6)/!8k!6!u1.2!n'
!y.title='!8E!6!dGW!n(!8k!6)/!8k!6!u2!n'
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:*]^1.2
  oplot,col=col[i],k[1:*],grav1[1:*,it[i]]/k[1:*]^2
endfor
;
print,"$mv idl.ps ~/tex/GW/lowfreq/fig/plater.eps"
print,"$mv idl.ps ~/tex/GW/lowfreq/fig/plater_B10.eps"
END