;$Id: rslice.pro,v 1.3 2023/08/11 08:11:08 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=12,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
end
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
;  show two different slices
file_slice1='data/slice_bb3.xy3'
;
t=0.
nx=1024L
ny=1024L
nz=256L
slice1=fltarr(nx,ny)
;
close,1 & openr,1,file_slice1,/f77
;
jump1=6 & jump2=2
jump1=0 & jump2=0
;
icount=0
!p.multi=0
for it=0,jump1 do readu,1,slice1,t
while not eof(1) do begin
    readu,1,slice1,t
    print,icount,t
    ;contour,slice1,/fil,nlev=30,tit=str(t)
    ;tvscl,slice1
    power2d_s,slice1,wav=k,spe=e,/noplo
    plot_oo,k[1:*],e[1:*]
    oplot,k[1:*],1e-6*k[1:*]^3
    if icount eq 0 then begin
      tt=t
      ee=e
    endif else begin
      tt=[tt,t]
      ee=[ee,e]
    endelse
    for it=0,jump2 do readu,1,slice1,t
    icount=icount+1
    wait,.6
endwhile
;
close,1
;
nk=n_elements(k)
nt=n_elements(tt)
ee=reform(ee,nk,nt)
save,file='2D_spectra.sav',tt,ee,k,nk,nt
;
print,"$mv idl.ps ~/tex/axel/vertical/fig/pbzmmax.eps"
END