import glob
import os
import numpy as np
import math

nx=512
Xray_r=np.zeros(nx**3)
filelist = glob.glob("jXRAY_*_dat")
for i,name in enumerate(filelist):
  f = open(name,'r')
  f.seek(12,os.SEEK_SET)
  xray=(np.fromfile(f,dtype='float64').reshape(256,256,256))*(4*math.pi)
  print "Reading %s"%name
  f.close()
  for ii in xrange(nx):
    for jj in xrange(nx):
      for kk in xrange(nx):
        index = ii*nx**2+jj*nx+kk
        ix = int(ii/2); iy = int(jj/2); iz = int(kk/2)
        Xray_r[index] = xray[ix][iy][iz]
  Xray_r.tofile(name)
  print "Writing %s"%name