;$Id: ppspec_save.pro,v 1.7 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=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
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.6
!x.margin=[8.2,0.5]
!y.margin=[3.1,0.4]
!p.multi=0
si2=1.4
siz=1.
;
file='spec_range.sav'
dir1='512_1e2_1e4_4e8_1em6a' & i1=21 & fac1=1.
dir2='512_1e2_1e4_4e8_2em6a' & i2=5 & fac2=1.
dir3='512_1e2_1e4_4e8_5em6a' & i3=2 & fac3=1.
dir4='512_1e2_1e4_4e8_1em5a' & i4=7 & col4=55 & fac4=1.
dir5='512_1e2_1e4_4e8_2em5a' & i5=1 & col5=122 & fac5=3.2
dir6='512_1e2_1e4_4e8_5em5a' & i6=1 & col6=155 & fac6=1.
dir7='512_1e2_1e4_4e8_1em4a2'& i7=5 & col7=122 & fac7=1.4
;dir8='512_1e2_1e4_4e8_2em4a' & i8=1 & col8=122 & fac8=1.
;dir9='512_1e2_1e4_4e8_5em4a' & i9=1 & col9=122 & fac9=13000.
dir0='512_1e2_1e4_4e8_1em3a' & i0=1 & col0=122 & fac0=1.0
;
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k,t!6)'
yr=[1e-28,1e-22]
yr=[1e-29,1e-21]
restore,'../'+dir3+'/'+file & i=i3
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
oplot,k[1:*],grav1_range(1:*,i)
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir3
;
restore,'../'+dir1+'/'+file & i=i1
oplot,k[1:*],grav1_range(1:*,i),li=1
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:*],grav1_range(1:*,i),li=2
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir2
;
restore,'../'+dir4+'/'+file & i=i4
oplot,k[1:*],grav1_range(1:*,i),col=col4
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir4
;
loadct,6
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
loadct,5
;
restore,'../'+dir6+'/'+file & i=i6
oplot,k[1:*],grav1_range(1:*,i),col=col6
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir6
;
restore,'../'+dir7+'/'+file & i=i7
oplot,k[1:*],fac7*grav1_range(1:*,i),col=col7
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir7
;
;restore,'../'+dir8+'/'+file & i=i8
;oplot,k[1:*],fac8*grav1_range(1:*,i),col=col8
;kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir8
;
;restore,'../'+dir9+'/'+file & i=i9
;oplot,k[1:*],fac9*grav1_range(1:*,i),col=col9
;kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir9
;
restore,'../'+dir0+'/'+file & i=i0
oplot,k[1:*],fac0*grav1_range(1:*,i),li=1
kGW=total(grav1_range(1:*,i))/total(grav1_range(1:*,i)/k[1:*]) & print,'kGW=',kGW,'  ',dir0
;
xx=[220.,1400.] & oplot,xx,1.3e-25*(xx/100.)^2,li=3
xx=[820.,7400.] & oplot,xx,7e-24/(xx/100.)^2,li=3
xyouts, 400.,5e-24,'!9A!8k!6!u2!n',siz=si2
xyouts,4500.,3e-27,'!9A!8k!6!u-2!n',siz=si2
;
xyouts,1000.,3.4e-28,'!6B1',siz=siz
xyouts,7000.,3.1e-28,'!6B3',siz=siz
xyouts,7000.,1.4e-26,'!6B4',siz=siz,col=55 & loadct,6
xyouts,6800.,1.9e-25,'!6B5',siz=siz,col=122 & loadct,5
xyouts,6600.,2.0e-24,'!6B6',siz=siz,col=155
xyouts,5000.,1.4e-23,'!6B7',siz=siz,col=122
xyouts,3600.,1.9e-22,'!6B10',siz=siz
;
cwd,run
print,'$mv idl.ps ~/tex/rei/GW/fig/ppspec_save.ps'
END