1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
|
# Copyright (c) 2017 Weitian LI <weitian@aaronly.me>
# MIT license
"""
Commonly used units and their conversions relations, as well as constants.
Astropy's units system is very powerful, but also very slow,
and may even be the speed bottleneck of the program.
This module provides commonly used units conversions by holding
them directly in a class, thus avoid repeated/unnecessary calculations.
"""
import astropy.units as au
import astropy.constants as ac
class Units:
"""
Commonly used units, especially in the CGS unit system.
"""
# Unit for electron momentum (p), thus its value is the Lorentz factor
# Unit: [g cm / s]
mec = ac.m_e.cgs.value*ac.c.cgs.value
# Energy of a still electron
# Unit: [erg]
mec2 = (ac.m_e * ac.c**2).to(au.erg).value
class UnitConversions:
"""
Commonly used units conversion relations.
Hold the conversion relations directly to avoid repeated/unnecessary
calculations.
"""
# Mass
Msun2g = au.solMass.to(au.g)
g2Msun = au.g.to(au.solMass)
# Time
Gyr2s = au.Gyr.to(au.s)
s2Gyr = au.s.to(au.Gyr)
# Length
kpc2m = au.kpc.to(au.m)
m2kpc = au.m.to(au.kpc)
Mpc2m = au.Mpc.to(au.m)
m2Mpc = au.m.to(au.Mpc)
kpc2cm = au.kpc.to(au.cm)
cm2kpc = au.cm.to(au.kpc)
Mpc2cm = au.Mpc.to(au.cm)
cm2Mpc = au.cm.to(au.Mpc)
Mpc2km = au.Mpc.to(au.km)
km2Mpc = au.km.to(au.Mpc)
kpc2km = au.kpc.to(au.km)
km2kpc = au.km.to(au.kpc)
km2cm = au.km.to(au.cm)
# Energy
keV2erg = au.keV.to(au.erg)
# Angle
rad2deg = au.rad.to(au.deg)
deg2rad = au.deg.to(au.rad)
rad2arcsec = au.rad.to(au.arcsec)
arcsec2rad = au.arcsec.to(au.rad)
arcmin2deg = au.arcmin.to(au.deg)
deg2arcmin = au.deg.to(au.arcmin)
arcsec2deg = au.arcsec.to(au.deg)
deg2arcsec = au.deg.to(au.arcsec)
arcmin2arcsec = au.arcmin.to(au.arcsec)
arcsec2arcmin = au.arcsec.to(au.arcmin)
# Temperature
eV2K = au.eV.to(ac.k_B*au.K)
keV2K = au.keV.to(ac.k_B*au.K)
class Constants:
"""
Commonly used constants, especially in the CGS unit system.
Astropy's constants are stored in SI units by default.
When request a constant in CGS unit system, additional (and slow)
conversions required.
"""
# Speed of light
c = ac.c.cgs.value # [cm/s]
# Atomic mass unit (i.e., a.m.u.)
u = ac.u.cgs.value # [g]
# Gravitational constant
G = ac.G.cgs.value # [cm^3/g/s^2]
# Electron charge
e = ac.e.gauss.value # [Fr] = [esu]
# Boltzmann constant
k_B = ac.k_B.cgs.value # [erg/K]
# Mean molecular weight
# Ref.: Ettori et al, 2013, Space Science Review, 177, 119-154, Eq.(6)
mu = 0.6
|
