ct.h

Go to the documentation of this file.

/* ///////////////////////////////////////////////////////////////////// */
/* ///////////////////////////////////////////////////////////////////// */

#define STAGGERED_MHD

/* ---- set labels for CT_EMF_AVERAGE ---- */

#define ARITHMETIC     1
#define UCT0           2
#define UCT_CONTACT    3
#define UCT_HLL        4
#define RIEMANN_2D     5

/* ---- staggered component labels ---- */

#define BX1s  0
#define BX2s  1
#define BX3s  2

/* ---- backward compatibility ---- */

#define BXs BX1s
#define BYs BX2s
#define BZs BX3s


#define FACE_EMF   11
#define EDGE_EMF   12

/*       Now define more convenient and user-friendly 
         pointer labels for geometry setting      */

#if GEOMETRY == CYLINDRICAL 
 #define iBRs    BX1s
 #define iBZs    BX2s
 #define iBPHIs  BX3s
#endif

#if GEOMETRY == POLAR 
 #define iBRs    BX1s
 #define iBPHIs  BX2s
 #define iBZs    BX3s
#endif

#if GEOMETRY == SPHERICAL 
 #define iBRs    BX1s
 #define iBTHs   BX2s
 #define iBPHIs  BX3s
#endif

  
/* ********************************************************************* */
typedef struct ElectroMotiveForce{

  double ***exj; 
  double ***exk; 
  double ***eyi; 
  double ***eyk; 
  double ***ezi; 
  double ***ezj; 
  signed char ***svx, ***svy, ***svz;

  int  ibeg, jbeg, kbeg;
  int  iend, jend, kend;
  double ***SxL;
  double ***SxR;
  double ***SyL;
  double ***SyR;
  double ***SzL;
  double ***SzR;
  double ***ex;
  double ***ey;
  double ***ez;
  double ***dbx_dy, ***dbx_dz;  
  double ***dby_dx, ***dby_dz;
  double ***dbz_dx, ***dbz_dy;

  double ***dvx_dx, ***dvy_dx, ***dvz_dx;
  double ***dvx_dy, ***dvy_dy, ***dvz_dy;
  double ***dvx_dz, ***dvy_dz, ***dvz_dz;
} EMF;

void CT_EMF_ArithmeticAverage (const EMF *, const double);
void CT_EMF_IntegrateToCorner (const Data *, const EMF *, Grid *);
void CT_AverageMagneticField (double ****bf, double ***UU[], Grid *);
void CT_AverageNormalMagField (const Data *, int, Grid *);
void CT_AverageTransverseMagFiled (const Data *, int, Grid *);
void CT_Update(const Data *, Grid *);
void CT_CheckDivB (double ***b[], Grid *);
void EMF_BOUNDARY (EMF *, Grid *);
void CT_StoreEMF (const State_1D *, int, int, Grid *);
void CT_EMF_Init (const State_1D *, int, int, Grid *);

void CT_StoreResistiveEMF (double **, int, int, Grid *);
EMF *CT_GetEMF (const Data *, Grid *);
void FillMagneticField (const Data *, int, Grid *); 

void MULTI_D_Reconstruct (double ***vv[], double ***bf[], double ***, 
                          double ***, double **, double **, Grid *);
/*
void EMF_BOUNDARY_STATES (Data_Arr, Data_Arr, 
                          Data_Arr, Data_Arr, 
                          Data_Arr, Data_Arr, Grid *);
*/
void EMF_USERDEF_BOUNDARY (EMF *, int, int, Grid *);

void CT_GetStagSlopes (const Data_Arr, EMF *, Grid *);


void CT_EMF_HLL_Solver (const Data *, const EMF *, Grid *);
void CT_StoreVelSlopes (EMF *, const State_1D *, int, int);
void CT_EMF_CMUSCL_Average (const Data *, const EMF *, Grid *);