GET_LUN, data

;
;--------Define a data structure and the parameters-----
ivar0=0
deltaivar=1
ivar1=14
;
restore,file='var_file/var0.sav'
power_snapshot, uu, eks=eksuu, ks=ksuu
nk = n_elements(x)/2-1
power_list=make_array(ivar1+1,n_elements(eksuu),/double,value=0)
power_mean=make_array(1,n_elements(eksuu),/double,value=0.)
save,file='ksuu.sav',ksuu
plot_oo, ksuu, eksuu, /xlog, /ylog, xrange=[1,nk]
;
for ivar=ivar0,ivar1 do begin
			restore,file='var_file/var'+str(ivar)+'.sav'
			; calculate power spectrum:
			power_snapshot, uu, eks=eksuu, ks=ksuu
			oplot, ksuu, eksuu
			print, 't=', t
      power_list[ivar,*]=eksuu
      power_mean=power_mean+eksuu
endfor
;
power_mean=power_mean/(ivar1+1)
save,file='power_list.sav',power_list
save,file='power_mean.sav',power_mean
;

end