Index: /issm/trunk-jpl/src/c/classes/Elements/Element.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Element.h	(revision 27850)
@@ -257,5 +257,5 @@
 		virtual void       ComputeStressTensor(void)=0;
 		virtual void       ComputeEsaStrainAndVorticity(void)=0;
-		virtual void       ComputeMeanEla(IssmDouble* paltitude,int* pcounter)=0;
+		//virtual void       ComputeMeanEla(IssmDouble* paltitude,int* pcounter)=0;
 		virtual void       Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters,Inputs* inputsin)=0;
 		virtual void       ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp)=0;
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 27850)
@@ -5283,26 +5283,26 @@
 /*}}}*/
 #endif
-void       Penta::ComputeMeanEla(IssmDouble* paltitude, int* pcounter){/*{{{*/
-
-	if(!this->IsOnSurface() & !this->IsIceInElement()) return;
-        int     ipos=0,ineg=0, counter=0;
-        IssmDouble mean_element_surface=0;
+//void       Penta::ComputeMeanEla(IssmDouble* paltitude, int* pcounter){/*{{{*/
+//
+//	if(!this->IsOnSurface() & !this->IsIceInElement()) return;
+//        int     ipos=0,ineg=0, counter=0;
+//        IssmDouble mean_element_surface=0;
         
 	/* input */
-        IssmDouble smb[NUMVERTICES],surface[NUMVERTICES];
-	Element::GetInputListOnVertices(&smb[0],SmbMassBalanceEnum);
-	Element::GetInputListOnVertices(&surface[0],SurfaceEnum);
-
-	for(int iv=0;iv<NUMVERTICES;iv++){
-		if(smb[iv]>0) ipos=ipos+1;
-		if(smb[iv]<0) ineg=ineg+1;
-	}  
+//        IssmDouble smb[NUMVERTICES],surface[NUMVERTICES];
+//	Element::GetInputListOnVertices(&smb[0],SmbMassBalanceEnum);
+//	Element::GetInputListOnVertices(&surface[0],SurfaceEnum);
+
+//	for(int iv=0;iv<NUMVERTICES;iv++){
+//		if(smb[iv]>0) ipos=ipos+1;
+//		if(smb[iv]<0) ineg=ineg+1;
+//	}  
 
 	/* we define ELA if an element has pos and neg smb */
-	if(ipos>0 & ineg>0){
-		for(int iv=0;iv<NUMVERTICES;iv++) mean_element_surface=mean_element_surface+surface[iv]/double(NUMVERTICES);
-		*paltitude+=mean_element_surface;
-		*pcounter+=counter+1;
-	}
-}
-/*}}}*/
+//	if(ipos>0 & ineg>0){
+//		for(int iv=0;iv<NUMVERTICES;iv++) mean_element_surface=mean_element_surface+surface[iv]/double(NUMVERTICES);
+//		*paltitude+=mean_element_surface;
+//		*pcounter+=counter+1;
+//	}
+//}
+/*}}}*/
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.h	(revision 27850)
@@ -64,4 +64,5 @@
 		void           ComputeSigmaNN(){_error_("not implemented yet");};
 		void           ComputeStressTensor();
+		//void           ComputeMeanEla(IssmDouble* paltitude, int* pcounter);
 		void           Configure(Elements* elements,Loads* loads,Nodes* nodes,Vertices* vertices,Materials* materials,Parameters* parameters,Inputs* inputsin);
 		void           ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp);
Index: /issm/trunk-jpl/src/c/classes/Elements/Seg.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Seg.h	(revision 27850)
@@ -50,5 +50,5 @@
 		void        ComputeSigmaNN(){_error_("not implemented yet");};
 		void        ComputeStressTensor(){_error_("not implemented yet");};
-		void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+		//void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
 		void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters,Inputs* inputsin){_error_("not implemented yet");};
 		void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp){_error_("not implemented yet");};
Index: /issm/trunk-jpl/src/c/classes/Elements/Tetra.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tetra.h	(revision 27850)
@@ -46,5 +46,5 @@
 		void        ComputeSigmaNN(){_error_("not implemented yet");};
 		void        ComputeStressTensor(){_error_("not implemented yet");};
-		void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+		//void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
 		void        ComputeDeviatoricStressTensor(){_error_("not implemented yet");};
 		void        ComputeEsaStrainAndVorticity(){_error_("not implemented yet!");};
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.h
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 27849)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.h	(revision 27850)
@@ -72,5 +72,5 @@
 		void        ComputeStressTensor();
 		void        ComputeSurfaceNormalVelocity();
-		void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
+		//void        ComputeMeanEla(IssmDouble* paltitude, int* pcounter){_error_("not implemented yet");};
 		void        Configure(Elements* elements,Loads* loads,Nodes* nodesin,Vertices* verticesin,Materials* materials,Parameters* parameters,Inputs* inputsin);
 		void        ControlInputSetGradient(IssmDouble* gradient,int enum_type,int control_index,int offset,int M,int N,int interp);