Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 238)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 239)
@@ -359,4 +359,6 @@
 		ModelFetchData((void**)&model->B,NULL,NULL,model_handle,"B","Matrix","Mat");
 		ModelFetchData((void**)&model->n,NULL,NULL,model_handle,"n","Matrix","Mat");
+		ModelFetchData((void**)&model->accumulation,NULL,NULL,model_handle,"accumulation","Matrix","Mat");
+		ModelFetchData((void**)&model->melting,NULL,NULL,model_handle,"melting","Matrix","Mat");
 		
 		for (i=0;i<model->numberofelements;i++){
@@ -379,6 +381,6 @@
 				penta_b[j]=*(model->bed+    ((int)*(model->elements+elements_width*i+j)-1)); 
 				penta_k[j]=*(model->drag+        ((int)*(model->elements+elements_width*i+j)-1)); 
-				//penta_melting[j]=*(model->melting+        ((int)*(model->elements+elements_width*i+j)-1));
-				//penta_accumulation[j]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+j)-1));
+				penta_melting[j]=*(model->melting+        ((int)*(model->elements+elements_width*i+j)-1));
+				penta_accumulation[j]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+j)-1));
 			}