Index: /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 175)
+++ /issm/trunk/src/c/ModelProcessorx/DiagnosticHoriz/CreateElementsNodesAndMaterialsDiagnosticHoriz.cpp	(revision 176)
@@ -228,4 +228,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++){
@@ -269,4 +271,13 @@
 			tria_p=model->p[i];
 			tria_q=model->q[i];
+
+			/*meling and accumulation*/
+			tria_melting[0]=*(model->melting+        ((int)*(model->elements+elements_width*i+0)-1));
+			tria_melting[1]=*(model->melting+        ((int)*(model->elements+elements_width*i+1)-1));
+			tria_melting[2]=*(model->melting+        ((int)*(model->elements+elements_width*i+2)-1));
+
+			tria_accumulation[0]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+0)-1));
+			tria_accumulation[1]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+1)-1));
+			tria_accumulation[2]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+2)-1));
 
 			/*element on iceshelf?:*/
@@ -368,4 +379,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)); 
+				tria_melting[j]=*(model->melting+        ((int)*(model->elements+elements_width*i+j)-1));
+				tria_accumulation[j]=*(model->accumulation+        ((int)*(model->elements+elements_width*i+j)-1));
 			}