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
license: mit |
# 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] |
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