                 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 

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

    while ((c = getopt(argc, argv, "ao:L:l:w:smSd:p")) != -1) {
        switch(c) {
            case 'L':
                printf("L is set\n");
