shearingbox.h

Go to the documentation of this file.

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

#if INTERPOLATION == LINEAR || INTERPOLATION == WENO3 || INTERPOLATION == LimO3
 #define SB_ORDER              2    
#else
 #define SB_ORDER              3    
#endif                                 

#ifndef FARGO
 #define SB_SYMMETRIZE_HYDRO  YES  
 #define SB_SYMMETRIZE_EY    (YES  && (DIMENSIONS == 3)) 
 #define SB_SYMMETRIZE_EZ     YES  
 #define SB_FORCE_EMF_PERIODS NO  
#else
 #define SB_SYMMETRIZE_HYDRO  NO  
 #define SB_SYMMETRIZE_EY    (NO  && (DIMENSIONS == 3)) 
 #define SB_SYMMETRIZE_EZ     NO 
 #define SB_FORCE_EMF_PERIODS NO 
#endif

/* ----------------------------
    Global variables
   ---------------------------- */

extern double sb_q;  
extern double sb_Omega; 
extern double sb_vy;  
#define sb_A (-0.5*sb_Omega*sb_q)  
void SB_Boundary (const Data *, int, Grid *) ;
void SB_ShearingInterp (double *, double *, double, int, Grid *);
void SB_CorrectFluxes  (Data_Arr, double, double, Grid *);
#ifdef STAGGERED_MHD
void SB_CorrectEMF (EMF *, Data_Arr, double, Grid *);
#endif
int  SB_JSHIFT (int);
void SB_SaveFluxes (State_1D *, Grid *);
void SB_SetBoundaryVar(double ***, RBox *, int, double, Grid *);
void SB_ShiftBoundaryVar(double ***, RBox *, int, double, Grid *);
void SB_FillBoundaryGhost(double ***, RBox *, int, int, Grid *);

#ifdef PARALLEL
 void ExchangeX (real *, real *, int, Grid *);
#endif