from phaser import *
i = InputMR_DAT()
HKLIN = "beta_blip.mtz"
F = "Fobs"
SIGF = "Sigma"
i.setHKLI(HKLIN)
i.setLABI(F,SIGF)
r = runMR_DAT(i)
print r.logfile()
hkl = r.getMiller()
fobs = r.getF()
sigma = r.getSIGF()
nrefl = min(10,hkl.size())
print "Data read from: " , HKLIN
print "Spacegroup Name (Hall symbol) = %s (%s)" % \
  (r.getSpaceGroupName(), r.getSpaceGroupHall())
print "Unitcell = " , r.getUnitCell()
print "First ", nrefl , " reflections"
print "%4s %4s %4s %10s %10s" % ("H","K","L",F,SIGF)
for i in range(0,nrefl):
  print "%4d %4d %4d %10.4f %10.4f" % \
    (hkl[i][0],hkl[i][1],hkl[i][2],fobs[i],sigma[i])