if !d.name eq 'PS' then begin
  device,xsize=18,ysize=14,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
;  mv idl.ps ~/tex/gal/gaussianpot/fig/pspec512b3.ps
;
;  first run ppower_all (to read in the data), and then run pspec.
;
!p.multi=0
!p.charsize=1.8
siz=1.8
ii=[2,14,26,38]
i1=ii(0) & n=n_elements(ii) & i2=ii(n-1)
!x.title='!8k!6 / !8k!6!d1!n'
plot_oo,k,spec1(*,i1),xr=[1,256],yr=[1e-20,2e-4],thick=1,li=0
i2=ii(n_elements(ii)-1)
oplot,k,spec1(*,i2),thick=3,li=1
oplot,k,spec2(*,i2),thick=3,li=1
for j=0,n-1 do begin
  i=ii(j)
  oplot,k,spec1(*,i),thick=4,li=1
  oplot,k,spec2(*,i),thick=4
  print,t(i)
endfor
;
xyouts,70,1e-6,'!8E!6(!8k!6)',siz=siz
kk=[18,80] & oplot,kk,.0001/kk^2
xyouts,30,3e-9,'k!u-2!n',siz=siz
;
xyouts,3,1e-15,'!8M!6(!8k!6)',siz=siz
kk=[14,50] & oplot,kk,5e-19*kk^4.
xyouts,20,8e-13,'k!u4!n',siz=siz
;
;  arrow
;
yy1=1e-9 & yy2=1e-6 & yy3=1e-9 & siz=1.6 & fx=0.4
xx=2./.2 & arrow,/data,xx,yy1,xx,yy2,thick=3
xyouts,xx*fx,yy3,'!8R!6=0.2',siz=siz
;
END