Aquila-Consortium/csiborgtools
0
0
Fork 0
mirror of https://github.com/Richard-Sti/csiborgtools.git synced 2024-12-22 17:18:02 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

change density to mass/mpc^3

Browse source
This commit is contained in:
rstiskalek 2022-11-06 10:06:40 +00:00
parent 5fa429ceb0
commit 859e3f4eb2
1 changed files with 7 additions and 7 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

14
csiborgtools/units/box_units.py
View file

@ -234,7 +234,7 @@ class BoxUnits:
def box2dens(self, density):
r"""
Convert density from box units to :math:`M_\odot / \mathrm{pc}^3`
Convert density from box units to :math:`M_\odot / \mathrm{Mpc}^3`
(with :math:`h=0.705`).
Parameters
@ -248,11 +248,11 @@ class BoxUnits:
Density in :math:`M_\odot / \mathrm{pc}^3`.
"""
return (density * self._unit_d / self._Msuncgs
* (units.pc.to(units.cm))**3)
* (units.Mpc.to(units.cm))**3)
def dens2box(self, density):
r"""
Convert density from :math:`M_\odot / \mathrm{pc}^3`
Convert density from :math:`M_\odot / \mathrm{Mpc}^3`
(with :math:`h=0.705`) to box units.
Parameters
@ -266,16 +266,16 @@ class BoxUnits:
Density in box units.
"""
return (density / self._unit_d * self._Msuncgs
/ (units.pc.to(units.cm))**3)
/ (units.Mpc.to(units.cm))**3)
def convert_from_boxunits(data, names, boxunits):
r"""
Convert columns named `names` in array `data` from box units to physical
units, such that
- length -> Mpc,
- mass -> solar mass,
- density -> solar mass per cubed kpc.
- length -> :math:`Mpc`,
- mass -> :math:`M_\odot`,
- density -> :math:`M_\odot / \mathrm{Mpc}^3`.
Any other conversions are currently not implemented. Note that the array
is passed by reference and directly modified, even though it is also
explicitly returned. Additionally centres the box coordinates on the