# magnitude limit
vmaglim=16.5
# assumed rotation velocity at 8 kpc
vrotsun=220.0
# position of the sun
xsun=-8.0
# input file name
inputfilei='../ubgaiaerr/gaiaei-out.dat'
inputfiled='../ubgaiaerr/gaiaed-out.dat'
# reading binary data
f=FortranFile(inputfilei)
nset=f.readInts()
print ' input Number of stars =',nset
f.close()
f=FortranFile(inputfiled)
i=0
rdata=np.zeros((nset,27))
while i < nset:
rdata[i,:]=f.readReals('d')
i+=1
f.close()
# parallax (mas) -> arcsec
paradet=rdata[:,8]/1000.0
# only chose positive parallax and < 16 kpc away
sindx=np.where(paradet > 1.0/(16.0*1.0e3))
# counting using parallax
vmaglim=18.5
# assumed rotation velocity at 8 kpc
vrotsun=220.0
# assumed M_v of RC, V-I colour
mvrc=1.2714
virc=1.0
# solar metallicity
zsun=np.float(0.019)
f=FortranFile(inputfile)
i=0
rdata=np.zeros((nset,10))
while i < nset:
rdata[i,:]=f.readReals('d')
i+=1
f.close()
print ' Number of stars=',nset
# set date
# position (kpc)
xs=rdata[:,0]
ys=rdata[:,1]
zs=rdata[:,2]
# velocity (km/s)
vxs=rdata[:,3]
vys=rdata[:,4]
vzs=rdata[:,5]
# particle mass (Msun)
ms=rdata[:,6]
# age (Gyr)
vmaglim=16.5
# assumed rotation velocity at 8 kpc
vrotsun=220.0
# assumed M_v of RC, V-I colour
mvrc=1.2714
virc=1.0
# solar metallicity
zsun=np.float(0.019)
f=FortranFile(inputfile)
i=0
rdata=np.zeros((nset,10))
while i < nset:
rdata[i,:]=f.readReals('d')
i+=1
f.close()
print ' Number of stars=',nset
# set date
# position (kpc)
xs=rdata[:,0]
ys=rdata[:,1]
zs=rdata[:,2]
# velocity (km/s)
vxs=rdata[:,3]
vys=rdata[:,4]
vzs=rdata[:,5]
# particle mass (Msun)
ms=rdata[:,6]
# age (Gyr)
vmaglim = 18.5
# assumed rotation velocity at 8 kpc
vrotsun = 220.0
# assumed M_v of RC, V-I colour
mvrc = 1.2714
virc = 1.0
# solar metallicity
zsun = np.float(0.019)
f = FortranFile(inputfile)
i = 0
rdata = np.zeros((nset, 10))
while i < nset:
rdata[i, :] = f.readReals('d')
i += 1
f.close()
print ' Number of stars=', nset
# set date
# position (kpc)
xs = rdata[:, 0]
ys = rdata[:, 1]
zs = rdata[:, 2]
# velocity (km/s)
vxs = rdata[:, 3]
vys = rdata[:, 4]
vzs = rdata[:, 5]
# particle mass (Msun)
ms = rdata[:, 6]
# age (Gyr)
