 			//r_radii are still needed, so no memory free just yet.
 
 		//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
-	
+
 			//Allocating arrays to store the alpha and beta data to compute winding numbers.
 			float *d_alpha, *d_beta;
 			checkCudaErrors(cudaMalloc(&d_alpha, dataCount*sizeof(float)));

	//Allocate host array for the winding numbers
	float* h_windingdata = (float*) malloc(BIGnroflines*sizeof(float));

	for (int xindex = 0; xindex < BIGgridSize.x; xindex += gridSizeRK4.x) {
		for (int yindex = 0; yindex < BIGgridSize.y; yindex += gridSizeRK4.y) {

			float4 startloc = BIGstartloc + xindex/BIGgridSize.x * BIGxvec + yindex/BIGgridSize.y * BIGyvec;
			float4 xvec = BIGxvec * (gridSizeRK4.x/BIGgridSize.x);

	dim3 blockSizeRK4(8,8); //gridSizeRK4*blockSizeRK4*steps should not exceed 2^26, to fit on 4GM RAM

	int BIGnroflines = BIGgridSize.x*BIGgridSize.y*blockSizeRK4.x*blockSizeRK4.y;
	std::cout << " BIGnroflines = " << BIGnroflines << std::endl;

	float4 BIGstartloc = make_float4(-2.0,0,-1.0,0); //Location (in Smietcoords) to start the integration, to be varied
	float4 BIGxvec = {2.0,0,0,0};
