# M. Tluczykont, 2010
#
# demo script
# this script reads a fits file that contains data fields following the SLF definition.
import sys
import pyfits
import numpy as np

if len(sys.argv)>1:
  # open fits file for reading (read in hdus):
  hdulist = pyfits.open(sys.argv[1])
else:
  print 'usage: ' + sys.argv[0] + '<fits file>'

# print general info about file:
hdulist.info()

# print the primary header:
print hdulist[0].header
print hdulist[1].header

raw_input()

# read data
t_data = hdulist[1].data
for irow in range(0,len(t_data)):
  print t_data[irow]
  print "%s\n" % t_data[irow].field('experiment')

print 'read %s rows' % len(t_data)