Index: /issm/trunk/src/c/objects/Elements/Penta.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 8489)
+++ /issm/trunk/src/c/objects/Elements/Penta.cpp	(revision 8490)
@@ -3685,5 +3685,5 @@
 	double     mixed_layer_capacity,thermal_exchange_velocity;
 	double     rho_ice,rho_water,pressure,dt,scalar_ocean;
-	double     meltingpoint,beta,heatcapacity,t_pmp;
+	double     heatcapacity,t_pmp;
 	double     xyz_list[NUMVERTICES][3];
 	double     xyz_list_tria[NUMVERTICES2D][3];
@@ -3705,6 +3705,4 @@
 	rho_ice=matpar->GetRhoIce();
 	heatcapacity=matpar->GetHeatCapacity();
-	beta=matpar->GetBeta();
-	meltingpoint=matpar->GetMeltingPoint();
 	this->parameters->FindParam(&dt,DtEnum);
 	Input* pressure_input=inputs->GetInput(PressureEnum); _assert_(pressure_input);
@@ -3720,5 +3718,5 @@
 
 		pressure_input->GetParameterValue(&pressure,gauss);
-		t_pmp=meltingpoint-beta*pressure;
+		t_pmp=matpar->TMeltingPoint(pressure);
 
 		scalar_ocean=gauss->weight*Jdet2d*rho_water*mixed_layer_capacity*thermal_exchange_velocity*(t_pmp)/(heatcapacity*rho_ice);
Index: /issm/trunk/src/c/objects/Loads/Pengrid.cpp
===================================================================
--- /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 8489)
+++ /issm/trunk/src/c/objects/Loads/Pengrid.cpp	(revision 8490)
@@ -465,6 +465,5 @@
 	double pressure;
 	double temperature;
-	double beta,t_pmp;
-	double meltingpoint;
+	double t_pmp;
 	int    new_active;
 	int    unstable=0;
@@ -487,10 +486,8 @@
 
 	//Recover our data:
-	meltingpoint=matpar->GetMeltingPoint();
-	beta=matpar->GetBeta();
 	parameters->FindParam(&stabilize_constraints,StabilizeConstraintsEnum);
 	
 	//Compute pressure melting point
-	t_pmp=meltingpoint-beta*pressure;
+	t_pmp=matpar->TMeltingPoint(pressure);
 
 	//Figure out if temperature is over melting_point, in which case, this penalty needs to be activated.
Index: /issm/trunk/src/c/objects/Materials/Matpar.cpp
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matpar.cpp	(revision 8489)
+++ /issm/trunk/src/c/objects/Materials/Matpar.cpp	(revision 8490)
@@ -385,2 +385,7 @@
 }
 /*}}}1*/
+/*FUNCTION Matpar::TMeltingPoint {{{1*/
+double Matpar::TMeltingPoint(double pressure){
+	return meltingpoint-beta*pressure;
+}
+/*}}}1*/
Index: /issm/trunk/src/c/objects/Materials/Matpar.h
===================================================================
--- /issm/trunk/src/c/objects/Materials/Matpar.h	(revision 8489)
+++ /issm/trunk/src/c/objects/Materials/Matpar.h	(revision 8490)
@@ -79,4 +79,5 @@
 		double GetGamma();
 		double GetKn();
+		double TMeltingPoint(double pressure);
 		/*}}}*/