8 #include "../Matrix/Matrix.h"
24 double G(
double alpha);
33 double T(
double alpha);
37 double Y(
double alpha);
46 double f_interp(
double E,
double f1,
double E1,
double f2,
double E2);
52 double mu_calc(
double L,
double pc,
double Alpha);
58 double Jc_calc(
double L,
double pc,
double Alpha);
69 double max(
double v1,
double v2);
88 double Dfe (
double alpha,
double Ke,
double L,
double Kp);
93 double Coulomb(
double L,
double energy);
107 double B (
double Lparam);
118 double Kfunc(
double pc);
122 double pfunc(
double K);
155 double J_L7 (
double K,
double x1,
double y1,
double x2,
double y2);
169 double mu2pc (
double L,
double mu,
double alpha);
174 double pc2mu (
double L,
double pc,
double alpha);
double max(double v1, double v2)
Return maximum.
double Precipitation(double Kp, double L, double energy)
Empirical electron lifetime.
double f_interp(double E, double f1, double E1, double f2, double E2)
Interpolation f on specific energy.
double Jc_calc(double L, double pc, double Alpha)
Caltulating J*c.
double mu2pc(double L, double mu, double alpha)
Convert μ to pc.
double Outer_Boundary_Flux_at_L7(double E7, string J_L7_function)
Call various flux model of the outer boundary.
string dtostr(double n)
Convert double to string.
double J_L7_corrected(double K)
Outer boundary spectrum (accurate)
double pfunc(double K)
Computation of moumentum from Kinetic energy .
Various useful constants, see source code for details.
double Coulomb(double L, double energy)
Coulomb scattering electron lifetime.
double J_L7(double K, double x1, double y1, double x2, double y2)
Interpolate spectrum.
double bounce_time_new(double L, double pc, double alpha)
bounce time
double Dfe(double alpha, double Ke, double L, double Kp)
Electrostatic radial diffusion coefficient.
double pc2mu(double L, double pc, double alpha)
Convert pc to μ.
double mu_calc(double L, double pc, double Alpha)
Caltulating μ.
double B(double Lparam)
Dipole magnetic field.
double T(double alpha)
Function for bounce time.
double alc(double L)
Loss cone calculations.
double J_L7_corrected_0(double K)
Outer boundary spectrum (accurate with 0 end the end)
double Y(double alpha)
Y(α) function.
bool check_time(int iter, int d_iter)
Small routine for checking if it is right iteration for the next output.
double density(double L)
Chorus density model.
double J_L7_old(double K)
Outer boundary spectrum (old)
double Kfunc(double pc)
Computation of Kinetic energy from given momentum .
double G(double alpha)
G-function.