4 void Init (
double *us,
double x1,
double x2,
double x3)
11 static int first_call = 1;
12 double Bx, By, Bz, B0;
13 double vx, vy, vz, vA;
24 vA = B0*B0/(w + B0*B0*(1.0 + eta*
eta));
25 vA /= 0.5*(1.0 + sqrt(1.0 - 4.0*eta*eta*vA*vA));
35 kx = 2.0*
CONST_PI*sqrt(2.0)*(c*x1 + s*x2);
54 us[
VX1] = vx*c - vy*
s;
55 us[
VX2] = vx*s + vy*
c;
58 us[
BX1] = Bx*c - By*
s;
59 us[
BX2] = Bx*s + By*
c;
64 us[
AX3] = B0*((-x1*s + x2*
c) - eta*sin(kx)/(2.0*
CONST_PI*sqrt(2.0)));
67 print1 (
"vA = %18.12e\n",vA);
68 print1 (
"T = %18.12e\n",c/vA);
void UserDefBoundary(const Data *d, RBox *box, int side, Grid *grid)
void print1(const char *fmt,...)
void BACKGROUND_FIELD(real x1, real x2, real x3, real *B0)
void Analysis(const Data *d, Grid *grid)
void Init(double *v, double x1, double x2, double x3)