/** \name Indexing convention for the aggregate grid function
 *  containing MHD arrays in second shell 
 */
#define MHD_POTENTIAL_NORTH   0
#define MHD_POTENTIAL_SOUTH   1
#define MHD_FAC_NORTH         2
#define MHD_FAC_SOUTH         3
#define MHD_SOUND_SPEED_NORTH 4
#define MHD_SOUND_SPEED_SOUTH 5
#define MHD_DENSITY_NORTH     6
#define MHD_DENSITY_SOUTH     7

/** \name Indexing convention for the aggregate grid function
 *  containing Solver arrays
 */
#define SOLVER_POTENTIAL_NORTH   0
#define SOLVER_POTENTIAL_SOUTH   1
#define SOLVER_FAC_NORTH         2
#define SOLVER_FAC_SOUTH         3
#define SOLVER_SOUND_SPEED_NORTH 4
#define SOLVER_SOUND_SPEED_SOUTH 5
#define SOLVER_DENSITY_NORTH     6
#define SOLVER_DENSITY_SOUTH     7
#define SOLVER_SIGMAP_NORTH      8
#define SOLVER_SIGMAP_SOUTH      9
#define SOLVER_SIGMAH_NORTH     10
#define SOLVER_SIGMAH_SOUTH     11
#define SOLVER_AVG_ENG_NORTH    12
#define SOLVER_AVG_ENG_SOUTH    13
#define SOLVER_NUM_FLUX_NORTH   14
#define SOLVER_NUM_FLUX_SOUTH   15
#define SOLVER_WIND_NORTH       16
#define SOLVER_WIND_SOUTH       17

/** \name Indexing convention for the aggregate grid function
 *  containing SM Globe arrays
 */
#define SM_GLOBE_POTENTIAL   0
#define SM_GLOBE_AVG_ENG     1
#define SM_GLOBE_NUM_FLUX    2