@parameters
default,isample,16
default,var,'/VAR3_'
spawn,'date >> rsample.log'
;
pc_read_kf,kf
pc_read_dim,obj=dim,/quiet
pc_read_param,/param2,obj=param,/quiet
nprocx=dim.nprocx
nprocy=dim.nprocy
nprocz=dim.nprocz
nxgrid=dim.nxgrid
nygrid=dim.nygrid
nzgrid=dim.nzgrid
;
nx=nxgrid/isample
ny=nygrid/isample
nz=nzgrid/isample
bb=fltarr(nx,ny,nz,3)
jj=fltarr(nx,ny,nz,3)
uu=fltarr(nx,ny,nz,3)
;
default,ipz1,0
default,ipz2,nprocz-1
for ipz=ipz1,ipz2 do begin
for ipy=0,nprocy-1 do begin
for ipx=0,nprocx-1 do begin
  l1=ipx*nx/nprocx & l2=(ipx+1)*nx/nprocx-1
  m1=ipy*ny/nprocy & m2=(ipy+1)*ny/nprocy-1
  n1=ipz*nz/nprocz & n2=(ipz+1)*nz/nprocz-1
  proc=ipx+ipy*nprocx+ipz*nprocx*nprocy
  rfile='data/proc'+str(proc)+var+str(isample)
  restore,rfile
  ;bb(l1:l2,m1:m2,n1:n2,*)=bb16
  ;jj(l1:l2,m1:m2,n1:n2,*)=jj16
  ;uu(l1:l2,m1:m2,n1:n2,*)=uu16
  ;bb(l1:l2,m1:m2,n1:n2,*)=bb2
  ;jj(l1:l2,m1:m2,n1:n2,*)=jj2
  ;uu(l1:l2,m1:m2,n1:n2,*)=uu2
  bb(l1:l2,m1:m2,n1:n2,*)=bb8
  jj(l1:l2,m1:m2,n1:n2,*)=jj8
  uu(l1:l2,m1:m2,n1:n2,*)=uu8
endfor
  print,rfile
endfor
endfor
;
save,file='data'+var+str(isample),bb,jj,uu
END