Compute slopes of staggered magnetic fields components: dBx/dy, dBx/dz, dBy/dx, dBy/dz, dBz/dx, dBz/dy,
Exclude normal derivatives (i.e. dbx_dx).
81 double ***bx, ***by, ***bz;
98 for (k = KOFFSET; k <
NX3_TOT - KOFFSET; k++){
99 for (j = JOFFSET; j <
NX2_TOT - JOFFSET; j++){
100 for (i = IOFFSET; i <
NX1_TOT - IOFFSET; i++){
103 bx[k][j][i] - bx[k][j-1][i]);
105 by[k][j][i] - by[k][j][i-1]);
108 bx[k][j][i] - bx[k-1][j][i]);
110 by[k][j][i] - by[k-1][j][i]);
112 bz[k][j][i] - bz[k][j][i-1]);
114 bz[k][j][i] - bz[k][j-1][i]);
long int NX2_TOT
Total number of zones in the X2 direction (boundaries included) for the local processor.
#define ARRAY_3D(nx, ny, nz, type)
long int NX3_TOT
Total number of zones in the X3 direction (boundaries included) for the local processor.
D_EXPAND(tot/[n]=(double) grid[IDIR].np_int_glob;, tot/[n]=(double) grid[JDIR].np_int_glob;, tot/[n]=(double) grid[KDIR].np_int_glob;)
static double MC_LIM2(double dp, double dm)
long int NX1_TOT
Total number of zones in the X1 direction (boundaries included) for the local processor.