     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;
+    int cuda_device = 0;
 
     ////////////////////////////////////////////////////////////////////////////////
     ////////////////////////////////////////////////////////////////////////////////
 
-    while ((c = getopt(argc, argv, "ao:L:l:w:smS")) != -1) {
+    while ((c = getopt(argc, argv, "ao:L:l:w:smSd:")) != -1) {
         switch(c) {
             case 'L':
                 printf("L is set\n");

////////////////////////////////////////////////////////////////////////

int doCalcRaDec(FILE *infile0, FILE *infile1, FILE *outfile, bool silent_on_GPU_testing, float scale_factor, int nbins, float hist_lower_range, float hist_upper_range, float hist_bin_width, int log_binning_flag, bool two_different_files, float conv_factor_angle){
   



    
    float *d_alpha0, *d_delta0;
    float *h_alpha0, *h_delta0;

    cudaFree(d_delta0);  
    cudaFree(d_alpha1);
    cudaFree(d_delta1);  
    cudaFree(dev_hist);

    return 0;
}  
//////////////////////////////////////////////////////////////////////




