Index: /issm/trunk-jpl/src/c/classes/Elements/Element.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.cpp	(revision 16914)
@@ -104,2 +104,6 @@
 	return shelf;
 }/*}}}*/
+IssmDouble Element::TMeltingPoint(IssmDouble pressure){/*{{{*/
+	_assert_(matpar);
+	return this->matpar->TMeltingPoint(pressure);
+}/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 16914)
@@ -62,4 +62,5 @@
 		IssmDouble GetMaterialParameter(int enum_in);
 		bool       IsFloating(); 
+		IssmDouble TMeltingPoint(IssmDouble pressure);
 
 		/*Virtual functions*/
@@ -179,5 +180,4 @@
 		virtual void   SmbGradients(void)=0;
 	   virtual Element*   SpawnBasalElement(void)=0;
-		virtual IssmDouble TMeltingPoint(IssmDouble pressure)=0;
 		virtual void   ReduceMatrices(ElementMatrix* Ke,ElementVector* pe)=0;
 		virtual void   ResetCoordinateSystem()=0;
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 16914)
@@ -3012,12 +3012,4 @@
 	  /*Update inputs*/
 	  this->inputs->AddInput(new PentaInput(SurfaceforcingsMassBalanceEnum,&smb[0],P1Enum));
-}
-/*}}}*/
-/*FUNCTION Penta::TMeltingPoint{{{*/
-IssmDouble Penta::TMeltingPoint(IssmDouble pressure){
-
-	_assert_(matpar);
-	return this->matpar->TMeltingPoint(pressure);
-
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 16914)
@@ -116,5 +116,4 @@
 		void   SmbGradients();
 	   Element* SpawnBasalElement(void);
-		IssmDouble  TMeltingPoint(IssmDouble pressure);
 		IssmDouble SurfaceArea(void);
 		void   Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finitelement);
Index: /issm/trunk-jpl/src/c/classes/Elements/Seg.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 16914)
@@ -192,5 +192,4 @@
 		void        ReduceMatrices(ElementMatrix* Ke,ElementVector* pe){_error_("not implemented yet");};
 		void	      SmbGradients(){_error_("not implemented yet");};
-		IssmDouble  TMeltingPoint(IssmDouble pressure){_error_("not implemented yet");};
 		IssmDouble  SurfaceArea(void){_error_("not implemented yet");};
 		void        Update(int index, IoModel* iomodel,int analysis_counter,int analysis_type,int finitelement){_error_("not implemented yet");};
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 16914)
@@ -2659,12 +2659,4 @@
 	  /*Update inputs*/
 	  this->inputs->AddInput(new TriaInput(SurfaceforcingsMassBalanceEnum,&smb[0],P1Enum));
-}
-/*}}}*/
-/*FUNCTION Tria::TMeltingPoint{{{*/
-IssmDouble Tria::TMeltingPoint(IssmDouble pressure){
-
-	_assert_(matpar);
-	return this->matpar->TMeltingPoint(pressure);
-
 }
 /*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 16913)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 16914)
@@ -119,5 +119,4 @@
 		void	      SmbGradients();
 	   Element*    SpawnBasalElement(void);
-		IssmDouble  TMeltingPoint(IssmDouble pressure);
 		int         VelocityInterpolation();
 		IssmDouble  PureIceEnthalpy(IssmDouble pressure){_error_("not implemented yet");};
