;$Id: pspecm.pro,v 1.23 2021/01/21 01:33:31 roper Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/tina/GW/fig/pspecm.ps
;  mv idl.ps ~/tex/prop/NSF/2017/fig/pspecm.ps
;  mv idl.ps ~/tex/prop/NASA/2018/GW/fig/pspecm.ps
;
fact1=1e-2 & fact2=1e-4
fact1=1e-2 & fact2=1e-2
siz=1.8
cwd,run
@parameters
!p.charsize=1.3
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
half='!s!u 1!n!r!s-!r!d 2!n'
circ_sym,0.7,1
file='specm.sav'
restore,file
;
!x.title='!8ck!6/!8H!6!d*!n'
!x.title='!8k!6'
!y.title='!8k!dH!nE!6!dGW!n/!13E!6!drad!n  and  !8k!dH!nE!6!dM!n/!13E!6!drad!n'
!y.title='!8E!6!dGW!n !8k!dH!n!6/!13E!6!drad!n  and  !8E!6!dM!n !8k!dH!n!6/!13E!6!drad!n'
!y.title='!7X!6!dGW!n(!8k!6)/!8k!6  and  !7X!6!dM!n(!8k!6)/!8k!6'
!y.title='!8E!6!dGW!n(!8k!6)  and  !8E!6!dM!n(!8k!6)'
xr=[k[1],max(k)]
plot_oo,k,fact1*spec1m,xr=xr,yr=yr
; oplot,k,spec1m+kine1m,col=55
;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,fact2*grav1m
;oplot,k,abs(grav2m),li=1
; oplot,k,grav2m,col=122,ps=8
; oplot,k,-grav2m,col=55,ps=8
;
;oplot,k,k^2*grah1m,li=2
;AB: removed because of Arthur
;
@postproc_specm
;
print,'$mv idl.ps ~/tex/roper/GW/fig/pspecm_'+run+'.ps'
print,'$mv idl.ps ~/tex/tina/Proposals/NASA-LPS-2020/trunk/pspecm_'+run+'.ps'
END