 
     float scale_factor = 1.0; // For if we need to convert input to arcsec or arcmin
     float conv_factor_angle = 57.2957795; // 180/pi // For if we need to convert arcdistance to arcsec or arcmin
-
+    int radec_input = 1; // are we using ra/dec coords, or x/y/z coords? 
     bool silent_on_GPU_testing = false;
 
     ////////////////////////////////////////////////////////////////////////////////

//#define DEFAULT_NBINS 30 

#define CONV_FACTOR 57.2957795 // 180/pi

////////////////////////////////////////////////////////////////////////
// Kernel to calculate angular distances between galaxies and histogram

                printf("Silent mode - don't run the GPU test (suppresses some output)");
                silent_on_GPU_testing = true;
                break;

            case 'p':
                printf("Using input files in Mpc format");
                radec_input = 0;
                break;
            case '?':
                printf("unknown arg %c\n", optopt);
                break;
