;$Id: calc_ko.pro,v 1.11 2022/10/19 02:47:02 brandenb Exp $
if !d.name eq 'PS' then begin
  device,xsize=18,ysize=19,yoffset=3
  !p.charthick=3 & !p.thick=3 & !x.thick=3 & !y.thick=3
  col0=255
endif else begin
  col0=0
endelse
;
!p.charsize=1.7
!x.margin=[7.8,.5]
!y.margin=[3.2,.5]
;
;  number of panels
;
npanels=2
npanels=1
;
@parameters
restore,pfile+'.sav'
;
fff_om=fff
s=size(fff)
nk=s[1]
nt=s[2]
;
;  time range
;
it1=500 & it2=nt-1 & nom=451  ;(K512irro_k6_ramp1c_rep)
it1=0 & it2=nt-1 & ikmax=10 & fmin=1e-8 & nom=287  ;(K512irrok6c_rep1_cont)
it1=500 & it2=nt-1 & nom=623  ;(K512sig0_k6_ramp1a_cont)
it1=500 & it2=nt-1 & nom=550  ;(K512irro_k6_dtf0p05h2)
;
for ik=0,nk-1 do begin
  fff_om[ik,*]=fft(fff[ik,*])
endfor
;
!x.title='!8k!6'
!y.title='!7x!6'
kk=dk*findgen(nk)
;
if npanels eq 2 then begin
  !p.multi=[0,2,1]
  contour,alog(abs(fff[0:ikmax,it1:it2])>fmin),/fil,nlev=30,kk[0:ikmax],tt[it1:it2]
endif else begin
  !p.multi=0
endelse
;
ikmax=10
dt=tt[1]-tt[0]
print,'dt=',dt
t_interval=dt*nt
dom=2.*!pi/t_interval
om=dom*findgen(nt)
contour,alog(abs(fff_om[0:ikmax,0:nom])>fmin),/fil,nlev=30,kk[0:ikmax],om[0:nom],/iso
oplot,kk,kk,col=255
oplot,kk,kk/sqrt(3.),col=255,li=3
;
cwd,run
print,'$mv idl.ps ~/GitHub/Axel/VortexGen/fig/calc_ko_'+run+'.eps'
END