LORENE
unites.h
1/*
2 * Lorene's units
3 *
4 */
5
6/*
7 * Copyright (c) 1999-2001 Eric Gourgoulhon
8 *
9 * This file is part of LORENE.
10 *
11 * LORENE is free software; you can redistribute it and/or modify
12 * it under the terms of the GNU General Public License as published by
13 * the Free Software Foundation; either version 2 of the License, or
14 * (at your option) any later version.
15 *
16 * LORENE is distributed in the hope that it will be useful,
17 * but WITHOUT ANY WARRANTY; without even the implied warranty of
18 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
19 * GNU General Public License for more details.
20 *
21 * You should have received a copy of the GNU General Public License
22 * along with LORENE; if not, write to the Free Software
23 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
24 *
25 */
26
27
28/*
29 * $Id: unites.h,v 1.8 2015/03/17 14:20:00 j_novak Exp $
30 * $Log: unites.h,v $
31 * Revision 1.8 2015/03/17 14:20:00 j_novak
32 * New class Hot_eos to deal with temperature-dependent EOSs.
33 *
34 * Revision 1.7 2014/10/13 08:52:37 j_novak
35 * Lorene classes and functions now belong to the namespace Lorene.
36 *
37 * Revision 1.6 2014/06/30 14:34:57 j_novak
38 * Update of the values of some constants (G, M_sol and MeV).
39 *
40 * Revision 1.5 2014/05/13 10:06:12 j_novak
41 * Change of magnetic units, to make the Lorene unit system coherent. Magnetic field is now expressed in Lorene units. Improvement on the comments on units.
42 *
43 * Revision 1.4 2004/12/01 12:28:32 p_grandclement
44 * Include math.h in unite.h
45 *
46 * Revision 1.3 2004/03/25 10:28:56 j_novak
47 * All LORENE's units are now defined in the namespace Unites (in file unites.h).
48 *
49 * Revision 1.2 2004/03/22 13:12:43 j_novak
50 * Modification of comments to use doxygen instead of doc++
51 *
52 * Revision 1.1.1.1 2001/11/20 15:19:27 e_gourgoulhon
53 * LORENE
54 *
55 * Revision 1.2 2000/03/17 15:28:29 eric
56 * Ajout de ggrav (G).
57 *
58 * Revision 1.1 1999/12/06 13:35:10 eric
59 * Initial revision
60 *
61 *
62 * $Header: /cvsroot/Lorene/C++/Include/unites.h,v 1.8 2015/03/17 14:20:00 j_novak Exp $
63 *
64 */
65
66#include <cmath>
67
68 //--------------------------//
69 // Standard LORENE units //
70 //--------------------------//
71
72namespace Lorene {
80 namespace Unites {
81 const double g_si = 6.6738E-11 ;
82 const double c_si = 2.99792458E+8 ;
83 const double kB_si = 1.3806488E-23 ;
84 const double rhonuc_si = 1.66E+17 ;
85 const double km_si = 1.E+3 ;
86 const double msol_si = 1.9885E+30 ;
87 const double mev_si = 1.602176565E-13 ;
88
89 const double r_unit = 1.e4 ;
90 const double v_unit = c_si ;
91 const double rho_unit = rhonuc_si ;
92 const double t_unit = r_unit/v_unit ;
93 const double m_unit = rho_unit * pow(r_unit, 3.) ;
94 const double g_unit = 1./(rho_unit*t_unit*t_unit) ;
95 const double f_unit = 1./t_unit ;
96
97 const double ggrav = g_si / g_unit ;
98 const double qpig = 4 * M_PI * ggrav ;
99 const double msol = msol_si/m_unit ;
100 const double km = km_si/r_unit ;
102 const double mevpfm3 = mev_si/( rho_unit * v_unit *v_unit) *1.e45 ;
104 const double m_u_mev = rho_unit / 1.e44 *c_si*c_si / mev_si ;
105 }
106
107
108 //----------------------------------//
109 // Electro-magnetic LORENE units //
110 //----------------------------------//
111
112
113
119namespace Unites_mag {
120 using namespace Unites ;
121 const double mu_si = 1.2566370614359173e-6 ;
122
123 const double j_unit = 1e11 ;
124
126 const double mag_unit = rho_unit*v_unit*v_unit/ (r_unit * j_unit) ;
128 const double elec_unit = mag_unit * v_unit ;
129
131 const double mu0_unit = rho_unit*v_unit*v_unit / (pow(j_unit ,2)*pow(r_unit,2));
133 const double mu0 = mu_si / mu0_unit ;
134}
135
136}
Cmp pow(const Cmp &, int)
Power .
Definition cmp_math.C:348
Lorene prototypes.
Definition app_hor.h:64
Standard electro-magnetic units.
Standard units of space, time and mass.