 }
 
 __global__ 
-void InterpolateCudaKernel(const double* __restrict__ source, 
-							double* __restrict__ target,
+void InterpolateCudaKernel(const double* __restrict__ d_source, 
+							double* __restrict__ d_target,
 							const point* __restrict__ d_grid,
 							himan::info_simple sourceInfo, 
 							himan::info_simple targetInfo)

{
	cudaStream_t stream;
	CUDA_CHECK(cudaStreamCreate(&stream));

	/* Determine all grid point coordinates that need to be interpolated.
	 * This is done with newbase by explicitly looping through the grid.

};

__host__ __device__ 
int Index(int x, int y, int sx)
{
	return y * sx + x;
}

__host__ __device__ 
int Index(point p, int sx)
{
	return Index(static_cast<int> (p.x),static_cast<int> (p.y), sx);
}

__global__
