 __constant__ uint nodeCountPerECM;
 __constant__ uint nodeCountPerCell;
 
-__constant__ double* nodeLocXAddress;
-__constant__ double* nodeLocYAddress;
-__constant__ double* nodeLocZAddress;
+__constant__ double* nodeLocXAddressConst;
+__constant__ double* nodeLocYAddressConst;
+__constant__ double* nodeLocZAddressConst;
 
 // This template method expands an input sequence by
 // replicating each element a variable number of times. For example,

	infoVecs.nodeCellType.resize(maxTotalNodeCount);
	infoVecs.nodeCellRank.resize(maxTotalNodeCount);
	infoVecs.nodeIsActive.resize(maxTotalNodeCount);
}

void SceNodes::initNodeAllocPara(uint totalBdryNodeCount,

			5 * sizeof(double));
	cudaMemcpyToSymbol(sceECMPara, mechPara.sceECMParaCPU, 5 * sizeof(double));
	//std::cout << "finished SceNodes:" << std::endl;

}

void SceNodes::initDimension(double domainMinX, double domainMaxX,
