$Id: pspecm_comp_lowfreq.pro,v 1.1 2022/01/22 06:06:35 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=6,yoffset=3
  !p.charthick=2 & !p.thick=2 & !x.thick=2 & !y.thick=2
end
;
;  mv idl.ps ~/tex/tina/GW/fig/pspecm_comp.ps
;  mv idl.ps ~/tex/tina/GW/fig/pspecm_comp2.ps
;  mv idl.ps ~/tex/prop/NASA/2018/GW/fig/pspecm_comp.ps
;
siz=1.0
!p.charsize=1.0
!p.multi=0
!x.margin=[6.8,.5]
!y.margin=[3.2,.5]
;
dir1='LowFreq1024sig1_k05_kf100a'
dir2='LowFreq1024sig1_k01_kf100a'
dir3='LowFreq1024sig0_k01_kf100a'
;
half='!s!u 1!n!r!s-!r!d 2!n'
circ_sym,0.7,1
file='specm.sav'
yr1=[1e-17,1e-6]
;
;-----------------------------------------------------------------------------
restore,'../'+dir1+'/'+file
;
!x.title='!8k!6/!8k!6!d1!n'
;!y.title='!8E!6(!8k!6)  and  !8f!6(!8k!6)'
!y.title='!6'
xr=[.09,266.]
plot_oo,k,spec1m,xr=xr,yr=yr1
;oplot,k,.5*k*abs(spec2m),li=1
;oplot,k,.5*k*spec2m,col=122,ps=8
;oplot,k,-.5*k*spec2m,col=55,ps=8
circ_sym,0.5,1
oplot,k,grav1m
;oplot,k,abs(grav2m),li=1
;oplot,k,grav2m,col=122,ps=8
;oplot,k,-grav2m,col=55,ps=8
;xyouts,280.,8e-6,siz=siz,'!6(a)'
;
xyouts,1.4,2e-7,siz=siz,'!8E!6!dM!n(!8k,!7g!6)'
xyouts,30.,8e-9,siz=siz,half+'!8kH!6!dM!n(!8k,!7g!6)',col=122
;
xyouts,160.,2e-12,siz=siz,'!8f!6(!8k,!7g!6)'
xyouts,3.0,1e-18,siz=siz,'!8g!6(!8k,!7g!6)!6'
;
xx=[3.,20] & oplot,xx,6e-13*xx^4.5
xx=[3.,20] & oplot,xx,4e-12/xx^0.0
xx=[.1,1.] & oplot,xx,2e-14*xx^1
;
;-----------------------------------------------------------------------------
restore,'../'+dir2+'/'+file
;
oplot,k,spec1m
oplot,k,grav1m
;
;-----------------------------------------------------------------------------
restore,'../'+dir3+'/'+file
;
oplot,k,spec1m,col=122
oplot,k,grav1m,col=122
;
END