Jo Bovy jobovy

These Python scripts can be used to reproduce the figures in Bovy & Tremaine (2012): On the local dark matter density.

Figure 1 is produced by:

python plot_dvcapprox.py

while Figure 2 is produced by:

python plot_sigmaz.py

	# The DF of a tidal stream peppered with impacts
	import copy
	import hashlib
	import numpy
	from scipy import integrate, special, stats, optimize, interpolate, signal
	from galpy.df import streamdf, streamgapdf
	from galpy.df.streamgapdf import _rotation_vy
	from galpy.util import bovy_conversion, bovy_coords
	class streampepperdf(streamdf):
	"""The DF of a tidal stream peppered with impacts"""

	#Module to read and process the Koposov et al. (2010) data
	import os
	import numpy
	from galpy.util import bovy_coords
	_DATADIR='../data/'
	# Koposov et al. (2010) transformation matrix
	_TKOP= numpy.zeros((3,3))
	_TKOP[0,:]= [-0.4776303088,-0.1738432154,0.8611897727]
	_TKOP[1,:]= [0.510844589,-0.8524449229,0.111245042]
	_TKOP[2,:]= [0.7147776536,0.4930681392,0.4959603976]