@parameters
pc_read_dim,obj=dim,/quiet
nprocx=dim.nprocx
nprocy=dim.nprocy
nprocz=dim.nprocz
nxgrid=dim.nxgrid
nygrid=dim.nygrid
nzgrid=dim.nzgrid
;
print,nprocx,nprocy,nprocz
ipz=2
ipy=2
ipx=2
;
uu=fltarr(nxgrid,nygrid,3)
oo=fltarr(nxgrid,nygrid,3)
bb=fltarr(nxgrid,nygrid,3)
jj=fltarr(nxgrid,nygrid,3)
;
ipz=0
for ipy=0,nprocy-1 do begin
for ipx=0,nprocx-1 do begin
  l1=ipx*nxgrid/nprocx
  l2=(ipx+1)*nxgrid/nprocx-1
  m1=ipy*nygrid/nprocy
  m2=(ipy+1)*nygrid/nprocy-1
  proc=ipx+ipy*nprocx+ipz*nprocx*nprocy
  pc_read_var,/magic,/trimall,obj=var,proc=proc,/quiet,variables=['uu','oo','bb','jj']
  uu(l1:l2,m1:m2,*)=reform(var.uu(*,*,0,*))
  oo(l1:l2,m1:m2,*)=reform(var.oo(*,*,0,*))
  bb(l1:l2,m1:m2,*)=reform(var.bb(*,*,0,*))
  jj(l1:l2,m1:m2,*)=reform(var.jj(*,*,0,*))
  print,l1,l2,m1,m2
endfor
endfor
;
save,file='data/slice.sav',uu,oo,bb,jj
END