;$Id: pxyaver.pro,v 1.2 2008/06/23 16:05:56 brandenb Exp $
;
;  reads xyaverage file and down-samples the data amount
;
pc_read_xyaver,o=xya
pc_read_grid,o=grid
pc_read_dim,o=dim
pc_read_param,o=param,/param2
zzz=grid.z(dim.n1:dim.n2)
Omega=param.Omega
cs0=param.cs0
;
nevery=40
nevery=1
s=size(xya.bxmz)
nz=s[1]
nt=s[2]
;
ntout=nt/nevery
print,'ntout=',ntout
print,'nz=',nz
it2=nt-1
it1=nt-nevery*ntout
;
tt=reform((reform(xya.t(it1:it2),nevery,ntout))(nevery-1,*))
bxmz=reform((reform(xya.bxmz(*,it1:it2),nz,nevery,ntout))(*,nevery-1,*))
bymz=reform((reform(xya.bymz(*,it1:it2),nz,nevery,ntout))(*,nevery-1,*))
;
contour,transpose(bymz),tt,zzz,nlev=30,/fil
;
save,file='xyaver.sav',bxmz,bymz,tt,zzz,Omega,cs0
END