;$Id: pearly_B7_comp.pro,v 1.2 2021/12/19 06:34:43 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
;
cwd,run
!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
;
dir1='512_1e2_1e4_4e8_1em4a'
dir2='1024_1e1_1e4_4e8_1em4a'
;
dir1='1024_1e1_1e4_4e8_1em4a'
dir2='2048_1e1_1e4_4e8_1em4a'
;
power,'_mag','_GWs',k=k1,spec1=spec1,spec2=grav1,i=n,tt=t1,/noplot,/lks,datatopdir='../'+dir1+'/data/'
power,'_mag','_GWs',k=k2,spec1=spec2,spec2=grav2,i=n,tt=t2,/noplot,/lks,datatopdir='../'+dir2+'/data/'
;
dk1=k1[1]
dk2=k2[1]
;
;  For Run B7
;
xr=[10.,25600.]
it1= [0  ,  1, 2, 4, 10, 16]
it2= [0  ,  1, 2, 4, 10, 16, 50,275]
it2= [0  ,  1, 2, 4, 10, 16]
col=[255,255,255, 55,122,155,188,188]
col=[  0,  0,  0, 55, 55,122,122,155]
;
xr=[10.,5120.]
it1= [ 2, 4, 8]
it2= [ 2, 4, 8]
col=[  0, 55,122]
yr1=[3d-16,2d-04] & yr2=[2d-14,2d-11]
;
xr=[5.,5120.]
it1= [ 0, 5, 20,100,162]
it2= [ 0, 5, 20,100,162]
col=[  0, 55,122,155,162]
yr1=[3d-16,2d-04] & yr2=[2d-14,2d-11]
;
th2=5
!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(it1)-1 do oplot,col=col[i],k1[1:*],spec1[1:*,it1[i]]/dk1
for i=0,n_elements(it2)-1 do oplot,col=col[i],k2[1:*],spec2[1:*,it2[i]]/dk2,li=1
;
!y.title='!8E!6!dGW!n(!8k!6)/!8k!6!u1.2!n'
!y.title='!8E!6!dGW!n(!8k!6) !8k!6!u1.8!n'
!y.title='!8E!6!dGW!n(!8k!6)'
plot_oo,xr,yr2,/nodata
for i=0,n_elements(it1)-1 do oplot,col=col[i],k1[1:*],grav1[1:*,it1[i]]/dk1
for i=0,n_elements(it2)-1 do oplot,col=col[i],k2[1:*],grav2[1:*,it2[i]]/dk2,li=1
;
print,"$mv idl.ps ~/tex/GW/lowfreq/fig/pearly_B7_comp.eps"
END