 	//Copy our texture properties (linear interpolation, texture access) to data array on device
 	checkCudaErrors(cudaBindTextureToArray(dataTex, dataArray, channelDesc));
 
+//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//
+
 	//Declare pointers to arrays with line data (output of integration), one each on device and host
 	float4 *d_lines, *h_lines;
 

	dataTex.addressMode[1] = cudaAddressModeBorder;
	dataTex.addressMode[2] = cudaAddressModeBorder;


	//Copy data to device
	cudaMemcpy3DParms copyParms = {0};
	copyParms.srcPtr = make_cudaPitchedPtr((void *)hostvfield[0][0], extent.width*sizeof(float4), extent.height, extent.depth);
	copyParms.dstArray = dataArray;

		std::cout << "Warning: origin read from file not equal to origin from constants.h" << std::endl;
	}

//%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%//

	//Allocate array on device
	cudaArray* dataArray;
	cudaChannelFormatDesc channelDesc = cudaCreateChannelDesc(32,32,32,32,cudaChannelFormatKindFloat);
