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 
72 namespace 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 ;
101  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 
119 namespace 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 }
Lorene prototypes.
Definition: app_hor.h:64
Standard units of space, time and mass.
Cmp pow(const Cmp &, int)
Power .
Definition: cmp_math.C:348
Standard electro-magnetic units.