 //	dimension of problem:		/*	  0    		 1				*/
 unsigned int		multiGPU=0;	/*	{ singleGPU, multiGPU } 	*/
 
-//	STAT:						/*	  0    1    2    3    4		*/
-unsigned int		STAT	= 2;/*	{ sum, min, max, std, mean }*/
+//	STAT:						/*	  0    1    2    3     4		*/
+unsigned int		STAT	= 0;/*	{ sum, min, max, mean, std }*/
 //	input:
-unsigned int		Nmaps	= 5;
+unsigned int		Nmaps	= 21;
 // 	---SMALL---
 const char			*iDIR 	= "/home/giuliano/work/Projects/LIFE_Project/run_clime_daily/run#2/maps/";
-const char			*clPAR	= "temp_min_h24";
+const char			*clPAR	= "rain_cum_h24";
 // 	---LARGE---
 //const char			*iDIR 	= "/home/giuliano/git/cuda/weatherprog-cudac/data/";
 //const char			*clPAR	= "L5_temp_min_h24";
 //const char			*clPAR	= "L1_temp_min_h24";
 const char			*YEAR	= "2011";
 const char			*MONTH	= "01";
-const char			*SPATMOD= "gprism";
+const char			*SPATMOD= "idw2";
 const char			*RES	= "80";
 const char			*EXT	= ".tif";
 //	output:

		CUDA_CHECK_RETURN( cudaMemcpy(oGRID, dev_oGRID, iMap_bytes, cudaMemcpyDeviceToHost) );
		end_t = clock();
		break;
	}
	// elapsed time [ms]:
	elapsed_time += (int)( (double)(end_t  - start_t ) / (double)CLOCKS_PER_SEC * 1000 );

	// call 2d reduction to computing mean & std for scatterplot:
	dim3 block_2d( BLOCKSIZE*BLOCKSIZE,1,1);
	dim3 grid_2d ( Nmaps,1,1);

	// memset:
	CUDA_CHECK_RETURN( cudaMemset(dev_oPLOT, 0,  (size_t)oPLOT_bytes) 		);

	/*
	 * 		****3-D aggregation****
	 */
	// kernel size:
	unsigned int 	BLOCKSIZE, GRIDSIZE;
	BLOCKSIZE	= 32;//floor(sqrt( devProp.maxThreadsPerBlock ));
