Index: /issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateElementsNodesAndMaterialsBalancedthickness.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateElementsNodesAndMaterialsBalancedthickness.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Balancedthickness/CreateElementsNodesAndMaterialsBalancedthickness.cpp	(revision 3696)
@@ -44,4 +44,6 @@
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
+		IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
 		
 		for (i=0;i<iomodel->numberofelements;i++){
@@ -66,4 +68,6 @@
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 
 	}
@@ -81,4 +85,6 @@
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"accumulation_rate");
 	
 		for (i=0;i<iomodel->numberofelements;i++){
@@ -103,4 +109,6 @@
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 
 	} //if (strcmp(meshtype,"2d")==0)
Index: /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Balancedthickness2/CreateElementsNodesAndMaterialsBalancedthickness2.cpp	(revision 3696)
@@ -48,4 +48,6 @@
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"accumulation_rate");
 		
 		for (i=0;i<iomodel->numberofelements;i++){
@@ -69,4 +71,6 @@
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 
 	}
Index: /issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateElementsNodesAndMaterialsBalancedvelocities.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateElementsNodesAndMaterialsBalancedvelocities.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Balancedvelocities/CreateElementsNodesAndMaterialsBalancedvelocities.cpp	(revision 3696)
@@ -44,4 +44,6 @@
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"accumulation_rate");
 		
 		for (i=0;i<iomodel->numberofelements;i++){
@@ -68,4 +70,6 @@
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 
 	}
@@ -83,4 +87,6 @@
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
+		IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
 
 		for (i=0;i<iomodel->numberofelements;i++){
@@ -107,4 +113,6 @@
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 
 	} //if (strcmp(meshtype,"2d")==0)
Index: /issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Prognostic/CreateElementsNodesAndMaterialsPrognostic.cpp	(revision 3696)
@@ -42,4 +42,6 @@
 		IoModelFetchData(&iomodel->elementoniceshelf,NULL,NULL,iomodel_handle,"elementoniceshelf");
 		IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater");
+		IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
@@ -65,4 +67,6 @@
 		xfree((void**)&iomodel->elementoniceshelf);
 		xfree((void**)&iomodel->elementonwater);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
@@ -80,4 +84,8 @@
 		IoModelFetchData(&iomodel->elementonsurface,NULL,NULL,iomodel_handle,"elementonsurface");
 		IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater");
+		IoModelFetchData(&iomodel->pressure,NULL,NULL,iomodel_handle,"pressure");
+		IoModelFetchData(&iomodel->temperature,NULL,NULL,iomodel_handle,"temperature");
+		IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
@@ -102,4 +110,8 @@
 		xfree((void**)&iomodel->elementonsurface);
 		xfree((void**)&iomodel->elementonwater);
+		xfree((void**)&iomodel->pressure);
+		xfree((void**)&iomodel->temperature);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
Index: /issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Prognostic2/CreateElementsNodesAndMaterialsPrognostic2.cpp	(revision 3696)
@@ -46,4 +46,6 @@
 		IoModelFetchData(&iomodel->elementoniceshelf,NULL,NULL,iomodel_handle,"elementoniceshelf");
 		IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater");
+		IoModelFetchData(&iomodel->melting_rate,NULL,NULL,iomodel_handle,"melting_rate");
+		IoModelFetchData(&iomodel->accumulation_rate,NULL,NULL,iomodel_handle,"accumulation_rate");
 		IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 		IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
@@ -67,4 +69,6 @@
 		xfree((void**)&iomodel->elementoniceshelf);
 		xfree((void**)&iomodel->elementonwater);
+		xfree((void**)&iomodel->melting_rate);
+		xfree((void**)&iomodel->accumulation_rate);
 		xfree((void**)&iomodel->vx);
 		xfree((void**)&iomodel->vy);
Index: /issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp
===================================================================
--- /issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp	(revision 3695)
+++ /issm/trunk/src/c/ModelProcessorx/Thermal/CreateElementsNodesAndMaterialsThermal.cpp	(revision 3696)
@@ -48,4 +48,6 @@
 	IoModelFetchData(&iomodel->rheology_n,NULL,NULL,iomodel_handle,"rheology_n");
 	IoModelFetchData(&iomodel->elementonwater,NULL,NULL,iomodel_handle,"elementonwater");
+	IoModelFetchData(&iomodel->pressure,NULL,NULL,iomodel_handle,"pressure");
+	IoModelFetchData(&iomodel->temperature,NULL,NULL,iomodel_handle,"temperature");
 	IoModelFetchData(&iomodel->vx,NULL,NULL,iomodel_handle,"vx");
 	IoModelFetchData(&iomodel->vy,NULL,NULL,iomodel_handle,"vy");
@@ -78,4 +80,6 @@
 	xfree((void**)&iomodel->rheology_B);
 	xfree((void**)&iomodel->elementonwater);
+	xfree((void**)&iomodel->pressure);
+	xfree((void**)&iomodel->temperature);
 	xfree((void**)&iomodel->vx);
 	xfree((void**)&iomodel->vy);