 	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};

//	float4write("../datadir/linedata.bin", BIGdataCount, h_lines);
	floatwrite("../datadir/windings.bin", BIGnroflines, h_windingdata);

	free(h_windingdata);

	return 0;
}

		//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
		   
			//Free host pointers
			free(h_origins);
			free(h_lines);
			free(h_lengths);
