; $Id: rslice_para.pro,v 1.2 2022/01/22 20:20:43 brandenb Exp $
;
;  reads xy slices
;  this routine is not very general and has been adjusted to this
;  particular purpose. Once data are processed, continue with
;  pcomp_last_QCD.pro.
;
@parameters
dir='data/proc'
file_slice='/slice_ggT.xy' & prefix='ggT_'
file_slice='/slice_hhT.xy' & prefix='hhT_'
;
precision=1d0
precision=1.
t=0.*precision
slice_ypos=0.*precision
pc_read_dim,obj=dim
nx=dim.nx/dim.nprocx
ny=dim.ny/dim.nprocy
nz=dim.nz/dim.nprocz
xy_slice=fltarr(nx,ny)*precision
slice_ypos=0.
;
;  open all files for x-direction
;
for iprocx=1,dim.nprocx do begin
  close,iprocx
  openr,iprocx,dir+str(iprocx-1)+file_slice,/f77
  print,'opened ',iprocx
endfor
;
icount=0
while not eof(1) do begin
  for iprocx=1,dim.nprocx do begin
    readu,iprocx,xy_slice,t,slice_ypos
    if iprocx eq 1 then xy_slice_=xy_slice else xy_slice_=[xy_slice_,xy_slice]
  endfor
  if icount mod 1000 eq 0 then print,icount,t
  f=reform(xy_slice_[*,0])
  if icount eq 0 then tt=t else tt=[tt,t]
  if icount eq 0 then ff=f else ff=[ff,f]
  ;plot,xy_slice(11,*),yr=[-1,1]*1e-3
  ;tvscl,xy_slice_
  icount=icount+1
  ;wait,.1
end
;
nt=icount
fff=reform(ff,nx*dim.nprocx,nt)
close,1
;
;print,824-512
;     312
;IDL> help,fff[*,312:*]
;print,"tvscl,fff[*,312:*]"
;i1_fff=nt-nomspec
i1_fff=0
fff_last=fff[*,i1_fff:*]
tt_last=tt[i1_fff:*]
nx=n_elements(f)
dx=2*!pi/nx
x=dx*(findgen(nx)-.5*nx+.5)
save,file=prefix+'last.sav',fff_last,tt_last,x
;<Expression>    DOUBLE    = Array[1024, 512]
END