Index: /issm/trunk-jpl/src/c/classes/FemModel.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17101)
+++ /issm/trunk-jpl/src/c/classes/FemModel.cpp	(revision 17102)
@@ -553,5 +553,5 @@
 						/*Allocate matrices*/
 						int         nbe       = this->elements->NumberOfElements();
-						IssmDouble* values    = xNew<IssmDouble>(nbe*nodesperelement);
+						IssmDouble* values    = xNewZeroInit<IssmDouble>(nbe*nodesperelement);
 						IssmDouble* allvalues = xNew<IssmDouble>(nbe*nodesperelement);
 
@@ -567,4 +567,5 @@
 
 						if(save_results)results->AddResult(new GenericExternalResult<IssmDouble*>(results->Size()+1,output_enum,allvalues,nbe,nodesperelement,step,time));
+						xDelete<IssmDouble>(allvalues);
 
 					}
Index: /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 17101)
+++ /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.cpp	(revision 17102)
@@ -142,4 +142,17 @@
 }
 /*}}}*/
+/*FUNCTION PentaInput::ResultToPatch{{{*/
+void PentaInput::ResultToPatch(IssmDouble* values,int nodesperelement,int sid){
+
+	int numnodes = this->NumberofNodes();
+
+	/*Some checks*/
+	_assert_(values);
+	_assert_(numnodes==nodesperelement);
+
+	/*Fill in arrays*/
+	for(int i=0;i<numnodes;i++) values[sid*numnodes + i] = this->values[i];
+}
+/*}}}*/
 
 /*Object functions*/
Index: /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.h	(revision 17101)
+++ /issm/trunk-jpl/src/c/classes/Inputs/PentaInput.h	(revision 17102)
@@ -41,5 +41,5 @@
 		int  GetResultInterpolation(void);
 		int  GetResultNumberOfNodes(void);
-		void ResultToPatch(IssmDouble* values,int nodesperelement,int sid){_error_("not supported yet");};
+		void ResultToPatch(IssmDouble* values,int nodesperelement,int sid);
 		void AddTimeValues(IssmDouble* values,int step,IssmDouble time){_error_("not supported yet");};
 		void Configure(Parameters* parameters);
Index: /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp	(revision 17101)
+++ /issm/trunk-jpl/src/c/classes/Inputs/TriaInput.cpp	(revision 17102)
@@ -139,10 +139,5 @@
 
 	/*Fill in arrays*/
-	for(int i=0;i<numnodes;i++){
-		values[sid*numnodes + i] = this->values[i];
-	}
-
-
-
+	for(int i=0;i<numnodes;i++) values[sid*numnodes + i] = this->values[i];
 }
 /*}}}*/