Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17354)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 17355)
@@ -410,5 +410,6 @@
 	switch(transfermethod){
 	case 0:
-		/*Just keepping the transfer to zero, should be OK with the initial value of transfer*/
+		/*Just keepping the transfer to zero*/
+		transfer=0.0;
 		break;
 	case 1:
@@ -456,5 +457,6 @@
 	switch(transfermethod){
 	case 0:
-		/*Just keepping the transfer to zero, should be OK with the initial value of transfer*/
+		/*Just keepping the transfer to zero*/
+		transfer=0.0;
 		break;
 	case 1:
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17354)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 17355)
@@ -513,5 +513,6 @@
 	switch(transfermethod){
 	case 0:
-		/*Just keepping the transfer to zero, should be OK with the initial value of transfer*/
+		/*Just keepping the transfer to zero*/
+		transfer=0.0;
 		break;
 	case 1:
@@ -560,5 +561,6 @@
 	switch(transfermethod){
 	case 0:
-		/*Just keepping the transfer to zero, should be OK with the initial value of transfer*/
+		/*Just keepping the transfer to zero*/
+		transfer=0.0;
 		break;
 	case 1:
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17354)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 17355)
@@ -4156,17 +4156,4 @@
 #endif
 
-/*FUNCTION Penta::CreateEPLDomainMassMatrix {{{*/
-ElementMatrix* Penta::CreateEPLDomainMassMatrix(void){
-
-	if (!IsOnBed()) return NULL;
-
-	Tria* tria=(Tria*)SpawnTria(0); //lower face is 0, upper face is 1.
-	ElementMatrix* Ke=tria->CreateEPLDomainMassMatrix();
-	delete tria->material; delete tria;
-
-	/*clean up and return*/
-	return Ke;
-}
-/*}}}*/
 /*FUNCTION Penta::GetHydrologyDCInefficientHmax{{{*/
 void  Penta::GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode){
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 17354)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 17355)
@@ -219,7 +219,7 @@
 		Gauss*         NewGauss(void);
 		Gauss*         NewGauss(int order);
-      Gauss*         NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order){_error_("not implemented yet");};
-      Gauss*         NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert);
-      Gauss*         NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order){_error_("not implemented yet");};
+		Gauss*         NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order){_error_("not implemented yet");};
+		Gauss*         NewGauss(IssmDouble* xyz_list, IssmDouble* xyz_list_front,int order_horiz,int order_vert);
+		Gauss*         NewGauss(int point1,IssmDouble fraction1,IssmDouble fraction2,bool mainlyfloating,int order){_error_("not implemented yet");};
 		Gauss*         NewGaussBase(int order);
 		Gauss*         NewGaussLine(int vertex1,int vertex2,int order);
@@ -241,5 +241,4 @@
 		IssmDouble     StabilizationParameter(IssmDouble u, IssmDouble v, IssmDouble w, IssmDouble diameter, IssmDouble kappa);
 
-		ElementMatrix* CreateEPLDomainMassMatrix(void);
 		void           GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode);
 		void           HydrologyEPLGetActive(Vector<IssmDouble>* active_vec);
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17354)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 17355)
@@ -4392,43 +4392,4 @@
 }
 /*}}}*/
-/*FUNCTION Tria::CreatEPLDomainMassMatrix {{{*/
-ElementMatrix* Tria::CreateEPLDomainMassMatrix(void){
-
-	/* Intermediaries */
-	IssmDouble  D,Jdet;
-	IssmDouble  xyz_list[NUMVERTICES][3];
-
-	/*Fetch number of nodes and dof for this finite element*/
-	int numnodes = this->NumberofNodes();
-
-	/*Initialize Element matrix and vectors*/
-	ElementMatrix* Ke    = new ElementMatrix(nodes,numnodes,this->parameters,NoneApproximationEnum);
-	IssmDouble*    basis = xNew<IssmDouble>(numnodes);
-
-	/*Retrieve all inputs and parameters*/
-	::GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-
-	/* Start looping on the number of gaussian points: */
-	GaussTria* gauss=new GaussTria(2);
-	for(int ig=gauss->begin();ig<gauss->end();ig++){
-
-		gauss->GaussPoint(ig);
-
-		GetNodalFunctions(basis,gauss);
-		GetJacobianDeterminant(&Jdet, &xyz_list[0][0],gauss);
-		D=gauss->weight*Jdet;
-
-		TripleMultiply(basis,1,numnodes,1,
-					&D,1,1,0,
-					basis,1,numnodes,0,
-					&Ke->values[0],1);
-	}
-
-	/*Clean up and return*/
-	delete gauss;
-	xDelete<IssmDouble>(basis);
-	return Ke;
-}
-/*}}}*/
 /*FUNCTION Tria::GetSolutionFromInputsOneDof{{{*/
 void  Tria::GetSolutionFromInputsOneDof(Vector<IssmDouble>* solution, int enum_type){
@@ -4588,5 +4549,4 @@
 void  Tria::ComputeEPLThickness(void){
 
-	int         i;
 	const int   numdof         = NDOF1 *NUMVERTICES;
 	bool        isefficientlayer;
@@ -4596,8 +4556,8 @@
 	IssmDouble  gravity,latentheat,EPLgrad2;
 	IssmDouble  EPL_N,epl_conductivity;
-	IssmDouble  activeEpl[numdof],thickness[numdof];
-	IssmDouble  eplhead[numdof],old_eplhead[numdof];
+	IssmDouble  thickness[numdof];
+	IssmDouble  eplhead[numdof];
 	IssmDouble  epl_slopeX[numdof],epl_slopeY[numdof];
-	IssmDouble  preceding_thickness[numdof],old_thickness[numdof];
+	IssmDouble  old_thickness[numdof];
 	IssmDouble  ice_thickness[numdof],bed[numdof];
 
@@ -4637,5 +4597,4 @@
 		}
 		else{
-			GetInputListOnVertices(&preceding_thickness[0],HydrologydcEplThicknessEnum);
 			for(int i=0;i<numdof;i++){
 				
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 17354)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 17355)
@@ -252,5 +252,4 @@
 		void UpdateConstraintsExtrudeFromTop(void);
 
-		ElementMatrix* CreateEPLDomainMassMatrix(void);
 		void           CreateHydrologyWaterVelocityInput(void);
 		void           GetHydrologyDCInefficientHmax(IssmDouble* ph_max, Node* innode);