30char et_rot_mag_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Etoile/et_rot_mag.C,v 1.24 2014/10/13 08:52:58 j_novak Exp $" ;
117#include "et_rot_mag.h"
118#include "utilitaires.h"
208 if ( withbphi == 0 ) {
343 ost <<
"Electromagnetic quantities" << endl ;
344 ost <<
"----------------------" << endl ;
347 ost <<
"***************************" << endl ;
348 ost <<
"**** perfect conductor ****" << endl ;
349 ost <<
"***************************" << endl ;
350 ost <<
"Prescribed charge : " <<
Q*j_unit*
pow(r_unit,3)/v_unit
354 ost <<
"***************************" << endl ;
355 ost <<
"**** isolator ****" << endl ;
356 ost <<
"***************************" << endl ;
358 ost <<
"Prescribed current amplitude : " <<
a_j*j_unit
359 <<
" [A/m2]" << endl ;
360 ost <<
"Magnetic Momentum : " <<
MagMom()/1.e32
361 <<
" [10^32 Am^2]" << endl ;
362 ost <<
"Radial magnetic field polar value : " <<
366 ost <<
"Tangent magnetic field equatorial value : " <<
371 ost <<
"Central magnetic field values : " <<
372 Magn()(0)(0,0,0,0)*mag_unit/1.e9
375 ost <<
"Radial electric field polar value : " <<
380 ost <<
"Radial electric field equatorial value : " <<
385 ost <<
"Magnetic/fluid pressure : "
386 << 1/(2*mu_si)*(
pow(
Magn()(0)(0,0,0,0),2) +
388 pow(
Magn()(2)(0,0,0,0),2))*mag_unit*mag_unit
389 / (
press()(0,0,0,0)*rho_unit*
pow(v_unit,2)) << endl ;
390 ost <<
"Computed charge : " <<
Q_comput() <<
" [C]" << endl ;
391 ost <<
"Interior charge : " <<
Q_int() <<
" [C]" << endl ;
394 ost <<
"Gyromagnetic ratio : " <<
GyroMag() << endl ;
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
void sauve(FILE *) const
Save in a file.
Equation of state base class.
Class for magnetized (isolator or perfect conductor), rigidly rotating stars.
virtual ~Et_rot_mag()
Destructor.
double Q_int() const
Computed charge from the integration of charge density over the star (i.e.
Tenseur Srr_em
rr component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame....
int conduc
Flag: conduc=0->isolator, 1->perfect conductor.
Cmp j_phi
-component of the current 4-vector
double Q_comput() const
Computed charge deduced from the asymptotic behaviour of At [SI units].
void operator=(const Et_rot_mag &)
Assignment to another Et_rot_mag.
Et_rot_mag(Map &mp_i, int nzet_i, bool relat, const Eos &eos_i, const int cond)
Standard constructor.
Cmp A_phi
-component of the electromagnetic potential 1-form divided by .
double a_j
Amplitude of the curent/charge function.
Tenseur Spp_em
component of the electromagnetic stress 3-tensor, as measured in the Eulerian frame.
double MagMom() const
Magnetic Momentum in SI units.
Cmp A_t
t-component of the elecctromagnetic potential 1-form, divided by .
virtual void del_deriv() const
Deletes all the derived quantities.
double GyroMag() const
Gyromagnetic ratio .
virtual double mass_g() const
Gravitational mass.
bool is_conduct() const
Tells if the star is made of conducting or isolating material.
virtual void sauve(FILE *) const
Save in a file.
Cmp j_t
t-component of the current 4-vector
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
Tenseur E_em
electromagnetic energy density in the Eulerian frame
Tenseur Jp_em
component of the electromagnetic momentum density 3-vector, as measured in the Eulerian frame.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
double Q
In the case of a perfect conductor, the requated baryonic charge.
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
Tenseur Elec() const
Computes the electric field spherical components in Lorene's units.
Cmp B_phi
-component of the magnetic field
Tenseur Magn() const
Computes the magnetic field spherical components in Lorene's units.
Class for isolated rotating stars *** DEPRECATED : use class Star_rot instead ***.
virtual void del_deriv() const
Deletes all the derived quantities.
void operator=(const Etoile_rot &)
Assignment to another Etoile_rot.
virtual void set_der_0x0() const
Sets to 0x0 all the pointers on derived quantities.
virtual const Itbl & l_surf() const
Description of the stellar surface: returns a 2-D Itbl containing the values of the domain index l on...
virtual ostream & operator>>(ostream &) const
Operator >> (virtual function called by the operator <<).
virtual void del_hydro_euler()
Sets to ETATNONDEF (undefined state) the hydrodynamical quantities relative to the Eulerian observer.
virtual void sauve(FILE *) const
Save in a file.
int nzet
Number of domains of *mp occupied by the star.
Map & mp
Mapping associated with the star.
Tenseur press
Fluid pressure.
const Tbl & xi_surf() const
Description of the stellar surface: returns a 2-D Tbl containing the values of the radial coordinat...
Base class for coordinate mappings.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Tensor handling *** DEPRECATED : use class Tensor instead ***.
void sauve(FILE *) const
Save in a file.
Cmp pow(const Cmp &, int)
Power .
int fread_be(int *aa, int size, int nb, FILE *fich)
Reads integer(s) from a binary file according to the big endian convention.
int fwrite_be(const int *aa, int size, int nb, FILE *fich)
Writes integer(s) into a binary file according to the big endian convention.
Standard electro-magnetic units.