unites.h

/*
 *  Lorene's units
 *
 */

/*
 *   Copyright (c) 1999-2001 Eric Gourgoulhon
 *
 *   This file is part of LORENE.
 *
 *   LORENE is free software; you can redistribute it and/or modify
 *   it under the terms of the GNU General Public License as published by
 *   the Free Software Foundation; either version 2 of the License, or
 *   (at your option) any later version.
 *
 *   LORENE is distributed in the hope that it will be useful,
 *   but WITHOUT ANY WARRANTY; without even the implied warranty of
 *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 *   GNU General Public License for more details.
 *
 *   You should have received a copy of the GNU General Public License
 *   along with LORENE; if not, write to the Free Software
 *   Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */


/*
 * $Id: unites.h,v 1.4 2004/12/01 12:28:32 p_grandclement Exp $
 * $Log: unites.h,v $
 * Revision 1.4  2004/12/01 12:28:32  p_grandclement
 * Include math.h in unite.h
 *
 * Revision 1.3  2004/03/25 10:28:56  j_novak
 * All LORENE's units are now defined in the namespace Unites (in file unites.h).
 *
 * Revision 1.2  2004/03/22 13:12:43  j_novak
 * Modification of comments to use doxygen instead of doc++
 *
 * Revision 1.1.1.1  2001/11/20 15:19:27  e_gourgoulhon
 * LORENE
 *
 * Revision 1.2  2000/03/17  15:28:29  eric
 * Ajout de ggrav (G).
 *
 * Revision 1.1  1999/12/06  13:35:10  eric
 * Initial revision
 *
 *
 * $Header: /cvsroot/Lorene/C++/Include/unites.h,v 1.4 2004/12/01 12:28:32 p_grandclement Exp $
 *
 */

#include <math.h>

                  //--------------------------//
                  //  Standard LORENE units   //
                  //--------------------------//     

namespace Unites {
const double g_si = 6.6726E-11 ;     
const double c_si = 2.99792458E+8 ;  
const double rhonuc_si = 1.66E+17 ;  
const double km_si = 1.E+3 ;     
const double msol_si = 1.989E+30 ;   
const double mev_si = 1.6021892E-13 ;   

const double r_unit = 1.e4 ;  
const double v_unit = c_si ; 
const double rho_unit = rhonuc_si ; 
const double t_unit = r_unit/v_unit ; 
const double m_unit = rho_unit * pow(r_unit, 3.) ;  
const double g_unit = 1./(rho_unit*t_unit*t_unit) ; 
const double f_unit = 1./t_unit ;   
    
const double ggrav = g_si / g_unit ;  
const double qpig = 4 * M_PI * ggrav ; 
const double msol = msol_si/m_unit ; 
const double km = km_si/r_unit ;    

const double mevpfm3 = mev_si/( 1.66E-27 * v_unit *v_unit) *10 ;  

}
                  //-----------------------------//
                  //  Electro-mag LORENE units   //
                  //-----------------------------//     



namespace Unites_mag {
  using namespace Unites ;
  const double mu_si = 1.2566370614359173e-6 ;
  
  const double j_unit = 1e11 ; 

  const double mag_unit = mu_si * r_unit * j_unit / 1e9 ;
  const double elec_unit = mag_unit * c_si / 1e3 ;
  const double mu0 = mu_si * pow(j_unit ,2) * pow(r_unit,2) 
    / (rho_unit*pow(c_si,2)) ; 
}
  

---