Go to the source code of this file.
void GetAreaFlux |
( |
const State_1D * |
state, |
|
|
double ** |
fA, |
|
|
double ** |
fvA, |
|
|
int |
beg, |
|
|
int |
end, |
|
|
Grid * |
grid |
|
) |
| |
Definition at line 4 of file get_area_flux.c.
11 double A, rp, rp2, rp3,
s, s2;
12 double **flux = state->
flux;
16 double **tc_flux = state->
tc_flux;
22 for (i = beg - 1; i <= end; i++){
25 #if GEOMETRY == CYLINDRICAL
29 IF_DUST(fA[i][MX3_D] = flux[i][MX3_D]*A*A;)
39 NVAR_LOOP(nv) fA[i][nv] = flux[i][nv]*rp2;
43 IF_DUST(fA[i][MX3_D] = flux[i][MX3_D]*rp3;)
53 EXPAND(fvA[i][
MX1] = visc_flux[i][
MX1]*A; ,
58 #elif GEOMETRY == POLAR && COMPONENTS >= 2
60 #elif GEOMETRY == SPHERICAL && COMPONENTS == 3
64 fvA[
i][ENG] = visc_flux[
i][ENG]*A;
71 for (j = beg - 1; j <= end; j++){
73 #if GEOMETRY == SPHERICAL
80 IF_DUST(fA[j][MX3_D] = flux[j][MX3_D]*s2;)
84 EXPAND(fvA[j][
iMR] = visc_flux[j][
iMR]*s; ,
88 fvA[j][ENG] = visc_flux[j][ENG]*s;
97 print1 (
"! GetAreaFlux(): should not be called during x3 sweep\n");
double ** flux
upwind flux computed with the Riemann solver
double ** visc_src
Viscosity source term.
void print1(const char *fmt,...)
double ** res_flux
Resistive flux (current)
double ** visc_flux
Viscosity flux.
int g_dir
Specifies the current sweep or direction of integration.
double ** tc_flux
Thermal conduction flux.
#define QUIT_PLUTO(e_code)
double * A
Right interface area, A[i] = .