 		photonVals[y][x] = 0;
 		lambdaVals[y][x] = 0;
 	      }
-	      if ((y - 220) * (y-220) + (x-200) * (x-200) > 20000)
+	      if ((y - 230) * (y-230) + (x-200) * (x-200) < 10000)
 	      {
 		photonVals[y][x] = 0;
 		lambdaVals[y][x] = 0;

    float q = sqrt(x * x + y * y);
float val;
    float den = (2 * ((float) CUDART_PI_F) * q * r);
    if (r > 0.5e-4)
{
    if ( r < 1e-4)
    val = 3 * den;
else
    val = 3 * (sinpif(2 * q * r)  - 2 * ((float) CUDART_PI_F) * q * r * cospif(2 * q * r));

    den = den * den * den;
    
    val /= den;
    val = val * val;
}
else
val = 1.0f;
    
    return val;
  }
};

  char* taskid = getenv("SLURM_ARRAY_TASK_ID");
  int tid;
  sscanf(taskid, "%d", &tid);
  spherer.lfactor = 1.0 /** pow(1.01, (tid) - 64)*/;
  spherer.tid = tid;

  for (int img = 0; img < fullsize[0]; img++)
