// exterior of the rectangular cylinder
#define BHXYZVSXYZ_EXT_NUM_X_POS    13
#define BHXYZVSXYZ_EXT_NUM_Y_POS    13
#define BHXYZVSXYZ_EXT_NUM_Z_POS    77
#define BHXYZVSXYZ_EXT_X_MIN    -0.0030
#define BHXYZVSXYZ_EXT_Y_MIN    -0.0030
#define BHXYZVSXYZ_EXT_Z_MIN    -0.0760
#define BHXYZVSXYZ_EXT_X_STEP    0.000500
#define BHXYZVSXYZ_EXT_Y_STEP    0.000500
#define BHXYZVSXYZ_EXT_Z_STEP    0.002000


// interior of the rectangular cylinder
#define BHXYZVSXYZ_NUM_X_POS    5
#define BHXYZVSXYZ_NUM_Y_POS    5
#define BHXYZVSXYZ_NUM_Z_POS    11
#define BHXYZVSXYZ_X_MIN    -0.0020
#define BHXYZVSXYZ_Y_MIN    -0.0020
#define BHXYZVSXYZ_Z_MIN    -0.0500
#define BHXYZVSXYZ_X_STEP    0.001000
#define BHXYZVSXYZ_Y_STEP    0.001000
#define BHXYZVSXYZ_Z_STEP    0.010000


double  BHXYZVSXYZ_EXT_X_POS[BHXYZVSXYZ_EXT_NUM_X_POS];
double  BHXYZVSXYZ_EXT_Y_POS[BHXYZVSXYZ_EXT_NUM_Y_POS];
double  BHXYZVSXYZ_EXT_Z_POS[BHXYZVSXYZ_EXT_NUM_Z_POS];
double  BHXYZVSXYZ_X_POS[BHXYZVSXYZ_NUM_X_POS];
double  BHXYZVSXYZ_Y_POS[BHXYZVSXYZ_NUM_Y_POS];
double  BHXYZVSXYZ_Z_POS[BHXYZVSXYZ_NUM_Z_POS];


for (i = 0; i < BHXYZVSXYZ_EXT_NUM_X_POS; i++) BHXYZVSXYZ_EXT_X_POS[i] = (double)BHXYZVSXYZ_EXT_X_MIN + i * (double)BHXYZVSXYZ_EXT_X_STEP;
for (i = 0; i < BHXYZVSXYZ_EXT_NUM_Y_POS; i++) BHXYZVSXYZ_EXT_Y_POS[i] = (double)BHXYZVSXYZ_EXT_Y_MIN + i * (double)BHXYZVSXYZ_EXT_Y_STEP;
for (i = 0; i < BHXYZVSXYZ_EXT_NUM_Z_POS; i++) BHXYZVSXYZ_EXT_Z_POS[i] = (double)BHXYZVSXYZ_EXT_Z_MIN + i * (double)BHXYZVSXYZ_EXT_Z_STEP;
for (i = 0; i < BHXYZVSXYZ_NUM_X_POS; i++) BHXYZVSXYZ_X_POS[i] = (double)BHXYZVSXYZ_X_MIN + i * (double)BHXYZVSXYZ_X_STEP;
for (i = 0; i < BHXYZVSXYZ_NUM_Y_POS; i++) BHXYZVSXYZ_Y_POS[i] = (double)BHXYZVSXYZ_Y_MIN + i * (double)BHXYZVSXYZ_Y_STEP;
for (i = 0; i < BHXYZVSXYZ_NUM_Z_POS; i++) BHXYZVSXYZ_Z_POS[i] = (double)BHXYZVSXYZ_Z_MIN + i * (double)BHXYZVSXYZ_Z_STEP;

// fine scan through the surfaces 
// X=XMIN surface
k=0;
for (l = 0; l < BHXYZVSXYZ_EXT_NUM_Y_POS; l++)
{
            for (z = 0; z < BHXYZVSXYZ_EXT_NUM_Z_POS; z++)
            {
                
            }}   

// Y=YMIN surface
l=0;
for (k = 0; k < BHXYZVSXYZ_EXT_NUM_X_POS; k++)
{
for (z = 0; z < BHXYZVSXYZ_EXT_NUM_Z_POS; z++)
{

            }}
            
                                     
// X=XMAX surface
k=BHXYZVSXYZ_EXT_NUM_X_POS-1;

for (l = 0; l < BHXYZVSXYZ_EXT_NUM_Y_POS; l++)
{
            for (z = 0; z < BHXYZVSXYZ_EXT_NUM_Z_POS; z++)
            {
                
            }}     
            
            

// Y=YMAX surface
l=BHXYZVSXYZ_EXT_NUM_Y_POS-1;

for (k = 0; k < BHXYZVSXYZ_EXT_NUM_X_POS; k++)
{
for (z = 0; z < BHXYZVSXYZ_EXT_NUM_Z_POS; z++)
{

            }}
            
            
            
// coarse scan through the volume 

/* Perform the Hall probe measurements */
/* Measure in the xz plane at each y position.  */
for (l = 0; l < BHXYZVSXYZ_NUM_Y_POS; l++)
{
for (k = 0; k < BHXYZVSXYZ_NUM_X_POS; k++)
{
for (z = 0; z < BHXYZVSXYZ_NUM_Z_POS; z++)
{

            }}}