;$Id: ppower_all_save_H1152.pro,v 1.1 2018/03/23 20:22:49 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/roper/GW/fig/ppower_all_save_H1152.ps
;
!p.charsize=1.7
!x.margin=[5.2,.5]
!y.margin=[3.1,.5]
siz=1.7
;
@parameters
!p.title='!6'
!x.title='!8k!6'
!y.title='!6'
;
;  determine normalization
;
xr=[1,max(k)]
yr=[1e-5,3e3]
;
restore,'../H1152a/spec_save.sav'
k1=k & k1[0]=1. & k1=1./k1 & k1[0]=0.
EEM=total(spec1_save)
kM=EEM/total(k1*spec1_save)
EEGW0=EEM^2/(8.*kM^2)
spec1_save0=EEM/kM
grav1_save0=EEGW0/kM
plot_oo,k,spec1_save/spec1_save0,xr=xr,yr=yr
oplot,k,grav1_save/grav1_save0,col=122
oplot,k,grah1_save/grav1_save0*k^2,col=55
print,spec1_save0
;
fact=1.7
restore,'../H1152b/spec_save.sav'
;k1=k & k1[0]=1. & k1=1./k1 & k1[0]=0.
;EEM=total(spec1_save)
;kM=EEM/total(k1*spec1_save)
;EEGW0=EEM^2/(8.*kM^2)
;spec1_save0=EEM/kM
;grav1_save0=EEGW0/kM
oplot,k,fact*spec1_save/spec1_save0,li=1
oplot,k,fact*grav1_save/grav1_save0,col=122,li=1
oplot,k,fact*grah1_save/grav1_save0*k^2,col=55,li=1
print,spec1_save0
;
END