rho=fltarr(256,256,256)
vx=fltarr(256,256,256)
vy=fltarr(256,256,256)
vz=fltarr(256,256,256)
bx=fltarr(256,256,256)
by=fltarr(256,256,256)
bz=fltarr(256,256,256)
for i=0,7 do begin
i1=i*32
i2=i*32+31
aa=string('R000498QXP',i,format='(a10,i3.3)')
offset=0
offset=offset*256UL*256*32*4+4+offset*8
rho[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=1
offset=offset*256UL*256*32*4+4+offset*8
vx[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=2
offset=offset*256UL*256*32*4+4+offset*8
vy[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=3
offset=offset*256UL*256*32*4+4+offset*8
vz[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=4
offset=offset*256UL*256*32*4+4+offset*8
bx[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=5
offset=offset*256UL*256*32*4+4+offset*8
by[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
offset=6
offset=offset*256UL*256*32*4+4+offset*8
bz[*,*,i1:i2]=read_binary(aa,data_type=4,data_start=offset,data_dims=[256,256,32],endian='little')
endfor
vx/=rho
vy/=rho
vz/=rho
print, "min rho ", min(rho), "     max rho ", max(rho)
print, "min vx  ", min(vx),  "     max vx  ", max(vx)
print, "min vy  ", min(vy),  "     max vy  ", max(vy)
print, "min vz  ", min(vz),  "     max vz  ", max(vz)
print, "min bx  ", min(bx),  "     max bx  ", max(bx)
print, "min by  ", min(by),  "     max by  ", max(by)
print, "min bz  ", min(bz),  "     max bz  ", max(bz)
points=float(256)*256*256
rho_mean=total(rho)/points
rho_meansq=total(rho^2)/points
Ek=total((vx^2+vy^2+vz^2)*rho)/points
Bk=total(bx^2+by^2+bz^2)/points
print, "Kinetic Energy: ", Ek, "   Magnetic Energy: ", Bk
end