     extern char *optarg;
     extern int optind, optopt, opterr;
     int c;
+    char *infilename0 = NULL;
+    char *infilename1 = NULL;
     char *outfilename = NULL;
     char defaultoutfilename[256];
     sprintf(defaultoutfilename,"default_out.dat");

            case 'l':
                log_binning_flag = atoi(optarg);
                printf("Will use log binning.\n");
                break;
            case 's':
                scale_factor = 206264.0; // To convert arcseconds to radians.


    for(int i=0; i<NUM_GALAXIES1; i++)
    {
        //fscanf(infile1, "%f %f %f", &temp0, &temp1, &temp2);
        fscanf(infile1, "%e %e %e %e %e %e", &tempdum, &tempdum, &tempdum, &temp0, &temp1, &temp2);
        //printf("%f\n",tempdum);
        h_x1[i] = temp0; ///scale_factor;
        h_y1[i] = temp1; ///scale_factor;
        h_z1[i] = temp2; ///scale_factor;
        if (i<10)
            printf("%f %f %f\n", h_x1[i], h_y1[i], h_z1[i]);
    }
