;$Id: ppspec3_save.pro,v 1.1 2020/12/22 15:58:09 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.0
!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_5em6a' & i1=2 & fac1=1.
dir2='512_1e2_1e4_4e8_1em5a' & i2=7 & col2=55 & fac2=1.
dir3='512_1e2_1e4_4e8_2em5a' & i3=1 & col3=122 & fac3=3.2
dir4='512_1e2_1e4_4e8_5em5a' & i4=1 & col4=155 & fac4=1.
dir5='512_1e2_1e4_4e8_1em4a2' & i5=5 & col5=122 & fac5=1.4
;
!x.title='!8k!6'
!y.title='!8E!6!dGW!n(!8k,t!6)'
yr=[1e-28,1e-22]
yr=[1d-27,1e-7]
restore,'../'+dir1+'/'+file & i=i1
plot_oo,k[1:*],k(1:*),yr=yr,/nodata
oplot,k[1:*],spec1_range(1:*,i)
;
restore,'../'+dir2+'/'+file & i=i2
oplot,k[1:*],spec1_range(1:*,i),col=col2
;
loadct,6
restore,'../'+dir3+'/'+file & i=i3
oplot,k[1:*],fac3*spec1_range(1:*,i),col=col3
loadct,5
;
restore,'../'+dir4+'/'+file & i=i4
oplot,k[1:*],spec1_range(1:*,i),col=col4
;
restore,'../'+dir5+'/'+file & i=i5
oplot,k[1:*],fac5*spec1_range(1:*,i),col=col5
;
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.,4e-27,'!9A!8k!6!u-2!n',siz=si2
;
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.9e-23,'!6B7',siz=siz,col=122
;
cwd,run
print,'$mv idl.ps ~/tex/rei/GW/fig/ppspec_save.ps'
END