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

            {
	        // this is a way simpler calculation. We already have the x,y,z coodis in co-moving distance, so we can simply do the distance

		xdiff = x0[idx] - x1[idx];
		ydiff = y0[idx] - y1[idx];
		zdiff = z0[idx] - z1[idx];

		dist = sqrt( (xdiff*xdiff) + (ydiff*ydiff) + (zdiff*zdiff));    

                if(dist < hist_min)
                    bin_index = 0; 


    for(int i=0; i<NUM_GALAXIES0; i++)
    {
        fscanf(infile0, "%f %f %f", &temp0, &temp1, &temp2);
        h_x0[i] = temp0/scale_factor;
        h_y0[i] = temp1/scale_factor;
        h_z0[i] = temp2/scale_factor;
        if (i<10)
            printf("%f %f %f\n", h_x0[i], h_y0[i], h_z0[i]);
    }


