18 double ***Ch_dt, ***Cp_dt;
19 static double *cmax, **dcoeff;
24 if (state.
rhs == NULL) {
53 for (nv =
NVAR; nv--; ) state.
v[i][nv] = d->
Vc[nv][k][j][i];
59 Ch_dt[
k][
j][
i] = 0.5*(cmax[
i] + cmax[i-1])*inv_dl[i];
61 inv_dl2 = 0.5*inv_dl[
i]*inv_dl[
i];
62 Cp_dt[
k][
j][
i] = (dcoeff[i-1][ENG] + dcoeff[
i][ENG])*inv_dl2;
76 for (nv =
NVAR; nv--; ) state.
v[j][nv] = d->
Vc[nv][k][j][i];
82 Ch_dt[
k][
j][
i] += 0.5*(cmax[
j] + cmax[j-1])*inv_dl[j];
84 inv_dl2 = 0.5*inv_dl[
j]*inv_dl[
j];
85 Cp_dt[
k][
j][
i] += (dcoeff[j-1][ENG] + dcoeff[
j][ENG])*inv_dl2;
double ** v
Cell-centered primitive varables at the base time level, v[i] = .
void TC_Flux(double ***, const State_1D *, double **, int, int, Grid *)
double ** rhs
Conservative right hand side.
void States(const State_1D *, int, int, 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.
double * GetInverse_dl(const Grid *)
#define ARRAY_3D(nx, ny, nz, type)
void MakeState(State_1D *)
void ComputeUserVar(const Data *d, Grid *grid)
double *** GetUserVar(char *)
#define TOT_LOOP(k, j, i)
int g_dir
Specifies the current sweep or direction of integration.
long int NX3_TOT
Total number of zones in the X3 direction (boundaries included) for the local processor.
long int IEND
Upper grid index of the computational domain in the the X1 direction for the local processor...
void SetIndexes(Index *indx, Grid *grid)
#define ARRAY_1D(nx, type)
long int NMAX_POINT
Maximum number of points among the three directions, boundaries excluded.
Riemann_Solver HLL_Solver
#define ARRAY_2D(nx, ny, type)
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.
static Image image[MAX_IMAGE_NUMBER]
long int IBEG
Lower grid index of the computational domain in the the X1 direction for the local processor...