;$Id: ppspec_saveD.pro,v 1.4 2020/12/25 11:25:06 brandenb Exp $
if !d.name eq 'PS' then begin
  ;device,xsize=18,ysize=8,yoffset=3
  device,xsize=18,ysize=11,yoffset=3
  !p.charthick=1.6 & !p.thick=1.6 & !x.thick=1.6 & !y.thick=1.6
end
;
;  mv idl.ps ~/tex/roper/GW/fig/pspecm_scl2.ps
;  mv idl.ps ~/tex/prop/computing/pdc/18/fig/pspecm_scl3n.ps
;
@parameters
!p.charsize=1.4
!x.margin=[8.2,0.5]
!y.margin=[3.1,0.4]
!p.multi=0
si2=1.4
siz=1.6
;
file='spec_range.sav'
dir1='512_1e2_1e4_2p5e7_1em5a'  & i1=5 & col1=1   & fac1=0.8
dir2='512_1e2_1e4_2p5e7_2em5a'  & i2=1 & col2=55  & fac2=1.0
dir3='512_1e2_1e4_2p5e7_5em5b'  & i3=1 & col3=122 & fac3=0.3
dir4='512_1e2_1e4_2p5e7_1em4a2' & i4=23& col4=155 & fac4=1.
dir5='512_1e2_1e4_2p5e7_2em4a2' & i5=22& col5=122 & fac5=0.9
dir6='512_1e2_1e4_2p5e7_5em4a2' & i6=1 & col6=155 & fac6=1.
dir7='512_1e2_1e4_2p5e7_1em3a2' & i7=5 & col7=122 & fac7=1.0
dir8='512_1e2_1e4_2p5e7_2em5a2' & i8=1 & col8=122 & fac8=1.0
;
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k,t!6)'
yr=[1e-28,1e-22]*1.4
restore,'../'+dir1+'/'+file & i=i1
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
oplot,k[1:*],fac1*grav1_range(1:*,i)
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir1
;
restore,'../'+dir2+'/'+file & i=i2
oplot,k[1:*],fac2*grav1_range(1:*,i),col=col2
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir2
;
loadct,6
restore,'../'+dir3+'/'+file & i=i3
oplot,k[1:*],fac3*grav1_range(1:*,i),col=col3
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir3
loadct,5
;
restore,'../'+dir4+'/'+file & i=i4
oplot,k[1:*],fac4*grav1_range(1:*,i),col=col4
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir4
;
restore,'../'+dir5+'/'+file & i=i5
oplot,k[1:*],fac5*grav1_range(1:*,i),col=col5
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir5
;
xx=[220.,2400.] & oplot,xx,4.3e-26*(xx/100.)^2,li=3
xx=[1200.,4000.] & oplot,xx,1.2e-23/(xx/100.)^1.8,li=3
xyouts, 400.,1.6e-24,'!9A!8k!6!u2!n',siz=si2
xyouts,4500.,1.3e-26,'!9A!8k!6!u-1.8!n',siz=si2
;
xyouts,7000.,9.9e-28,'!6D1',siz=siz
xyouts,7000.,5.3e-26,'!6D2',siz=siz,col=55 & loadct,6
xyouts,6800.,5.0e-25,'!6D3',siz=siz,col=122 & loadct,5
xyouts,6600.,5.3e-24,'!6D4',siz=siz,col=155
xyouts,5000.,4.9e-23,'!6D5',siz=siz,col=122
;
cwd,run
print,'$mv idl.ps ~/tex/rei/GW/fig/ppspec_saveD.ps'
END