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

    // Read in the second file
    ////////////////////////////////////////////////////////////////////////////

    fscanf(infile1, "%d", &NUM_GALAXIES1);

    int size_of_galaxy_array1 = NUM_GALAXIES1 * sizeof(float);    
    printf("SIZE 1 # GALAXIES: %d\n",NUM_GALAXIES1);


    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");
