PLUTO
4.4-patch2
|
Average staggered magnetic field to zone center. More...
#include "pluto.h"
Functions | |
void | CT_AverageStaggeredFields (double ****Vs, double ****UU, RBox *box, Grid *grid) |
Average staggered magnetic field components to form a zone-centered field, e.g., . In cylindrical coordinates the volume-averaged radial component is obtained as
this yields
Simlar expressions hold in spherical coordinates. Check with the following MAPLE script:
The averaging is done inside an arbitrary rectangular box. The box may include boundary cells only during the predictor step of the CT-CTU algorithm, whereas is useless for RK time stepping.
When the CT_EN_CORRECTION flag is enabled, we also redefine the zone total energy using the newly formed cell-centered field, i.e.,