44 static int first_call = 1;
45 double ***pr, ***dn, dp;
63 }
else if (dir ==
JDIR){
66 }
else if (dir ==
KDIR){
87 MPI_Allreduce (&n, &n_glob, 1, MPI_INT, MPI_MAX, MPI_COMM_WORLD);
109 }
else if (dir ==
JDIR){
117 }
else if (dir ==
KDIR){
212 dp = q[
k][
j][i+1] - q[
k][
j][i-1];
213 dp /= q[
k][
j][i+1] + q[
k][
j][i-1];
215 if (fabs(dp) > 1.e-5)
return(i);
226 dp = q[
k][j+1][
i] - q[
k][j-1][
i];
227 dp /= q[
k][j+1][
i] + q[
k][j-1][
i];
229 if (fabs(dp) > 1.e-5)
return(j);
240 dp = q[k+1][
j][
i] - q[k-1][
j][
i];
241 dp /= q[k+1][
j][
i] + q[k-1][
j][
i];
243 if (fabs(dp) > 1.e-5)
return(k);
long int NX1
Number of interior zones in the X1 directions (boundaries excluded) for the local processor...
int rbound
Same as lbound, but for the right edge of the grid.
long int NX2_TOT
Total number of zones in the X2 direction (boundaries included) for the local processor.
double **** Vc
The main four-index data array used for cell-centered primitive variables.
int lend
Local end index for the local array.
int lbeg
Local start index for the local array.
long int g_stepNumber
Gives the current integration step number.
long int NX2
Number of interior zones in the X2 directions (boundaries excluded) for the local processor...
void UnsetJetDomain(const Data *d, int dir, Grid *grid)
static int GetRightmostIndex(int, double ***)
long int NX3_TOT
Total number of zones in the X3 direction (boundaries included) for the local processor.
int nghost
Number of ghost zones.
long int IEND
Upper grid index of the computational domain in the the X1 direction for the local processor...
void SetJetDomain(const Data *d, int dir, int log_freq, Grid *grid)
long int KBEG
Lower grid index of the computational domain in the the X3 direction for the local processor...
long int KEND
Upper grid index of the computational domain in the the X3 direction for the local processor...
long int JBEG
Lower grid index of the computational domain in the the X2 direction for the local processor...
long int JEND
Upper grid index of the computational domain in the the X2 direction for the local processor...
long int NX1_TOT
Total number of zones in the X1 direction (boundaries included) for the local processor.
long int NX3
Number of interior zones in the X3 directions (boundaries excluded) for the local processor...
long int IBEG
Lower grid index of the computational domain in the the X1 direction for the local processor...