Index: ../trunk-jpl/src/c/classes/Inputs2/IntInput2.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Inputs2/IntInput2.cpp	(revision 24938)
+++ ../trunk-jpl/src/c/classes/Inputs2/IntInput2.cpp	(revision 24939)
@@ -75,10 +75,6 @@
 /*IntInput2 management*/
 void IntInput2::GetInput(int* pvalue,int index){/*{{{*/
 
-	if(index<0){
-		printf("-------------- file: IntInput2.cpp line: %i\n",__LINE__); 
-		int* temp = xNew<int>(3);
-	}
 	_assert_(index>=0); 
 	_assert_(index<this->size); 
 
Index: ../trunk-jpl/src/c/classes/Inputs2/DoubleInput2.cpp
===================================================================
--- ../trunk-jpl/src/c/classes/Inputs2/DoubleInput2.cpp	(revision 24938)
+++ ../trunk-jpl/src/c/classes/Inputs2/DoubleInput2.cpp	(revision 24939)
@@ -75,9 +75,6 @@
 /*DoubleInput2 management*/
 void DoubleInput2::GetInput(IssmDouble* pvalue,int index){/*{{{*/
 
-	if(index<0){
-		printf("-------------- file: DoubleInput2.cpp line: %g\n",__LINE__); 
-	}
 	_assert_(index>=0); 
 	_assert_(index<this->size); 
 
Index: ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp
===================================================================
--- ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 24938)
+++ ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.cpp	(revision 24939)
@@ -1286,64 +1286,6 @@
 	xDelete<IssmDouble>(enthalpies);
 	return kappa;
 }/*}}}*/
-void           EnthalpyAnalysis::GetBAdvec(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
-	/*Compute B  matrix. B=[B1 B2 B3 B4 B5 B6] where Bi is of size 5*1.
-	 * For node i, Bi' can be expressed in the actual coordinate system
-	 * by:
-	 *       Bi_advec =[ h ]
-	 *                 [ h ]
-	 *                 [ h ]
-	 * where h is the interpolation function for node i.
-	 *
-	 * We assume B has been allocated already, of size: 3x(1*NUMNODESP1)
-	 */
-
-	/*Fetch number of nodes for this finite element*/
-	int numnodes = element->GetNumberOfNodes();
-
-	/*Get nodal functions*/
-	IssmDouble* basis=xNew<IssmDouble>(numnodes);
-	element->NodalFunctions(basis,gauss);
-
-	/*Build B: */
-	for(int i=0;i<numnodes;i++){
-		B[numnodes*0+i] = basis[i];
-		B[numnodes*1+i] = basis[i];
-		B[numnodes*2+i] = basis[i];
-	}
-
-	/*Clean-up*/
-	xDelete<IssmDouble>(basis);
-}/*}}}*/
-void           EnthalpyAnalysis::GetBAdvecprime(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss){/*{{{*/
-	/*Compute B  matrix. B=[B1 B2 B3 B4 B5 B6] where Bi is of size 5*1.
-	 * For node i, Bi' can be expressed in the actual coordinate system
-	 * by:
-	 *       Biprime_advec=[ dh/dx ]
-	 *                     [ dh/dy ]
-	 *                     [ dh/dz ]
-	 * where h is the interpolation function for node i.
-	 *
-	 * We assume B has been allocated already, of size: 3x(1*numnodes)
-	 */
-
-	/*Fetch number of nodes for this finite element*/
-	int numnodes = element->GetNumberOfNodes();
-
-	/*Get nodal functions derivatives*/
-	IssmDouble* dbasis=xNew<IssmDouble>(3*numnodes);
-	element->NodalFunctionsDerivatives(dbasis,xyz_list,gauss);
-
-	/*Build B: */
-	for(int i=0;i<numnodes;i++){
-		B[numnodes*0+i] = dbasis[0*numnodes+i];
-		B[numnodes*1+i] = dbasis[1*numnodes+i];
-		B[numnodes*2+i] = dbasis[2*numnodes+i];
-	}
-
-	/*Clean-up*/
-	xDelete<IssmDouble>(dbasis);
-}/*}}}*/
 void           EnthalpyAnalysis::GetBasalConstraints(Vector<IssmDouble>* vec_spc,Element* element){/*{{{*/
 
 	/*Intermediary*/
Index: ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.h
===================================================================
--- ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.h	(revision 24938)
+++ ../trunk-jpl/src/c/analyses/EnthalpyAnalysis.h	(revision 24939)
@@ -40,8 +40,6 @@
  		static void       DrainageUpdateEnthalpy(FemModel* femmodel);
 		static IssmDouble EnthalpyDiffusionParameter(Element* element,IssmDouble enthalpy,IssmDouble pressure);
 		static IssmDouble EnthalpyDiffusionParameterVolume(Element* element,int enthalpy_enum);
-		void              GetBAdvec(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
-		void              GetBAdvecprime(IssmDouble* B,Element* element,IssmDouble* xyz_list,Gauss* gauss);
 		static void       GetBasalConstraints(Vector<IssmDouble>* vec_spc,Element* element);
 		static void       GetBasalConstraintsSteadystate(Vector<IssmDouble>* vec_spc,Element* element);
 		static void       GetBasalConstraintsTransient(Vector<IssmDouble>* vec_spc,Element* element);