                 break;
             case 'd':
                 cuda_device = atoi(optarg); // Use this CUDA device.
-                conv_factor_angle *= 3600.0; // convert radians to arcseconds.
                 printf("Will attempt to use CUDA device %d\n",cuda_device);
                 break;
             case 'S':


    for(int i=0; i<NUM_GALAXIES0; i++)
    {
        fscanf(infile0, "%f %f", &temp0, &temp1, &temp2);
        h_x0[i] = temp0/scale_factor;
        h_y0[i] = temp1/scale_factor;
        h_z0[i] = temp2/scale_factor;
        //if (i<10)
        //printf("%e %e\n", h_x0[i], h_y0[i], h_y0[i],);
    }



                }

                atomicAdd(&shared_hist[bin_index],1);
                //atomicAdd(&shared_hist[9],1);
                /*
                if (int(dist)<hist_max)
                    atomicAdd(&shared_hist[int(dist)],1);
                    */

            }
        }
