                 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_GALAXIES1; i++)
    {
        fscanf(infile1, "%f %f", &temp0, &temp1, &temp2);
        h_x1[i] = temp0/scale_factor;
        h_y1[i] = temp1/scale_factor;
        h_z1[i] = temp2/scale_factor;
        //if (i<10)
        //printf("%e %e\n", h_x1[i], h_y1[i], h_z1[i]);
    }

// get device diagnostics 


    float lo = hist_lower_range;
    float hi = 0;
    for(int k=0; k<nbins+2; k++)
    {
        if (k==0)
        {
            fprintf(outfile, "Underflow below %.3e %s %lu \n", lo, ",",  hist_array[k]);
        }
        else
        {
