;PRO concatenate

journ
f=findfile('vv_hd96*.save')
s=size(f)
nf=s[1]
s=size(f)
nf=s[1]
print,'nf=',nf
print,f

nt=0
for i=0,nf-1 do begin $
  restore,f[i] &$
  s=size(timet) &$
  nt=nt+s[1] &$
endfor
print,'nt=',nt
s=size(vy1)
nx=s[1]
nz=s[2]
print,nx,nz
vy01t=fltarr(nx,nz,nt)
vx01t=vy01t
vz01t=vy01t
times=fltarr(nt)
nxt=times

tt=0
;tt=nt0
for i=0,nf-1 do begin $
  print,f[i] &$
  restore,/ver,f[i] &$
  s=size(vy01) &$
  nx=s[1] &$
  nz=s[2] &$
  nt1=s[3] &$
  print,'nt1=',nt1 &$
  timet[0]=timet[1]-0.6 &$
  vx01t[0:nx-1,0:nz-1,tt:tt+nt1-1]=vx01 &$
  vy01t[0:nx-1,0:nz-1,tt:tt+nt1-1]=vy01 &$
  vz01t[0:nx-1,0:nz-1,tt:tt+nt1-1]=vz01 &$
  times[tt:tt+nt1-1]=timet &$
  nxt[tt:tt+nt1-1]=nx &$
  tt=tt+nt1 &$
  print,'tt=',tt &$
  plot,times,yst=3 &$
endfor

 save,/verb,file='vyt_hd96_747-1355.save',vy01t,times,nxt
 save,/verb,file='vxt_hd96_747-1355.save',vx01t,times,nxt
 save,/verb,file='vzt_hd96_747-1355.save',vz01t,times,nxt
print,'times=',min(times,max=max),max

flush,!journal
;exit
;END