from numarray import *
from numarray.random_array import *

size = 50000
bmean = 17.0
bsig = 2.0
bkmean = 2.0
bksig = 1.0
rmin = 0.1
rmax = 2.0


outhandle = open('simdata.d', 'w')
seed()
number = arange(size)
rah = randint(0,24,size)
ram = randint(0,60,size)
ras = uniform(0.0,60.0,size)
dd = 57.295*arccos(uniform(-1.0,1.0,size)) - 90.0
dd = dd.astype(Int8)
dm = randint(0,60,size)
ds = uniform(0.0,60.0,size)
bmag = normal(bmean,bsig,size)
bkcolour = normal(bkmean,bksig,size)
kmag = bmag - bkcolour
radio = uniform(rmin,rmax,size)

# Now loop through all the sources printing out their properties

for n in number:
    print >> outhandle, '%5d %5d %2d %5.2f %6d %2d %5.2f %8.2f %5.2f %5.3f' % (n, rah[n], ram[n], ras[n], dd[n], dm[n], ds[n], bmag[n], kmag[n], radio[n])
# Finally close your connection to the files.

outhandle.close()