 
     for(int i=0; i<NUM_GALAXIES1; i++)
     {
-        fscanf(infile1, "%f %f %f", &temp0, &temp1, &temp2);
-        h_x1[i] = temp0/scale_factor;
-        h_y1[i] = temp1/scale_factor;
-        h_z1[i] = temp2/scale_factor;
+        //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]);
     }

            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.

    printf("hist_upper_range: %f\n",hist_upper_range);

    FILE *infile0, *infile1, *outfile ;

    infile0 = fopen(argv[optind],"r");
    infile1 = fopen(argv[optind+1],"r");

    printf("Opening input file 0: %s\n",argv[optind]);
    printf("Opening input file 1: %s\n",argv[optind+1]);

    outfile = fopen(outfilename, "w");

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