;PRO bflux

fb=findfile('bb*.save')
fi=findfile('int_h5*.save')
fv=findfile('vv*.save')
fb=fb[1:*]
fi=fi[1:*]
fv=fv[1:*]

n=10
f2=fb[0:n]
f1=fb[n+1:*]
fb=[f1,f2]
f2=fi[0:n]
f1=fi[n+1:*]
fi=[f1,f2]
f2=fv[0:n]
f1=fv[n+1:*]
fv=[f1,f2]
print,fb
print,fi
print,fv

s=size(fb)
nf=s(1)
print,'nf=',nf
nt=0
for i=0,nf-1 do begin $
  restore,fi[i] &$
  s=size(timet) &$
  nt=nt+s(1) &$
endfor
print,'nt=',nt
s=size(int)
nx=s[1]
nz=s[2]
intt=fltarr(nx,nz,nt)
;bxt=intt
;byt=intt
;bzt=intt
;vxt=intt
;vyt=intt
;vzt=intt
times=fltarr(nt)
tmp=intt

tt=0
for i=0,nf-1 do begin $
  print,fi[i] &$
; restore,/ver,fv[i] &$
; restore,/ver,fb[i] &$
  restore,/ver,fi[i] &$
  help,int,bx1,vx1,timet &$
  s=size(int) &$
  nt1=s[3] &$
  print,'nt1=',nt1 &$
  if (s[1] ne nx) then begin &$
    tmp=int & int=fltarr(nx,nz,nt1) &$
    for t=0,nt1-1 do int[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
  endif &$
  intt[*,*,tt:tt+nt1-1]=int &$
; print,'Bx' &$
; if (s[1] ne nx) then begin &$
;   tmp=bx1 & bx1=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do bx1[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; bxt [*,*,tt:tt+nt1-1]=bx1 &$
; print,'By' &$
; if (s[1] ne nx) then begin &$
;   tmp=by1 & by1=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do by1[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; byt [*,*,tt:tt+nt1-1]=by1 &$
; print,'Bz' &$
; if (s[1] ne nx) then begin &$
;   tmp=bz1 & bz1=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do bz1[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; bzt [*,*,tt:tt+nt1-1]=bz1 &$
; print,'Vx' &$
; if (s[1] ne nx) then begin &$
;   tmp=vx01 & vx01=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do vx01[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; vxt [*,*,tt:tt+nt1-1]=vx01 &$
; print,'Vy' &$
; if (s[1] ne nx) then begin &$
;   tmp=vy01 & vy01=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do vy01[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; vyt [*,*,tt:tt+nt1-1]=vy01 &$
; print,'Vz' &$
; if (s[1] ne nx) then begin &$
;   tmp=vz01 & vz01=fltarr(nx,nz,nt1) &$
;   for t=0,nt1-1 do vz01[*,*,t]=spline2d(reform(tmp[*,*,t]),nx,nz) &$
; endif &$
; vzt [*,*,tt:tt+nt1-1]=vz01 &$
  times[tt:tt+nt1-1]=timet &$
  tt=tt+nt1 &$
  print,'tt=',tt &$
  plot,times &$
endfor

; fix glitches in intt

;save,/verb,file='intt_bbt_vvt_h5b.save',intt,bxt,byt,bzt,vxt,vyt,vzt,times 
save,/verb,file='intt_h5.save',intt,times 
print,'times=',min(times,max=max),max

exit
;END