# M. Tluczykont 2010
#
# 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>'

columns = hdulist[1].get_coldefs()
#print type(obj)
#print obj.__class__
#print obj
#print obj.name
print '# column definitions (name/units):',
verbose = 0
if verbose:
  # print general info about file:
  hdulist.info()
  # print the primary header:
  print hdulist[0].header
  print hdulist[1].header
else:
  for col in columns:
    print col.name,"/",col.unit,' | ',
print '#'

#print "Starting file contents dump (press enter):"
#raw_input()

# read data
t_data = hdulist[1].data
for irow in range(0,len(t_data)):
  for field in t_data[irow]:
    print "\t%10s" % str(field),# comma prevents \n being printed
  print "\n",			# end of line, explicitly print a \n