Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 21430)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCEfficientAnalysis.cpp	(revision 21431)
@@ -405,5 +405,5 @@
 	bool     active_element;
 	int      domaintype;
-	Element* basalelement;
+	Element* basalelement=NULL;
 
 	/*Get basal element*/
@@ -414,5 +414,5 @@
 			break;
 		case Domain3DEnum:
-			if(!element->IsOnBase()) return NULL;
+			if(!element->IsOnBase()) return;
 			basalelement = element->SpawnBasalElement();
 			break;
@@ -437,5 +437,5 @@
 
 	/*Use the dof list to index into the solution vector: */
-	for(i=0;i<numnodes;i++){
+	for(int i=0;i<numnodes;i++){
 		eplHeads[i]=solution[doflist[i]];
 		if(xIsNan<IssmDouble>(eplHeads[i])) _error_("NaN found in solution vector");
@@ -653,9 +653,8 @@
 
 	bool        active_element;
-	int         i,j;
 	int         domaintype;
 	IssmDouble  h_max;
 	IssmDouble  sedheadmin;
-	Element*    basalelement;
+	Element*    basalelement=NULL;
 
 	/*Get basal element*/
@@ -666,5 +665,5 @@
 			break;
 		case Domain3DEnum:
-			if(!element->IsOnBase()) return NULL;
+			if(!element->IsOnBase()) return;
 			basalelement = element->SpawnBasalElement();
 			break;
@@ -693,6 +692,6 @@
 	/*Get minimum sediment head of the element*/
 	sedheadmin=sedhead[0];
-	for(i=1;i<numnodes;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i];
-	for(i=0;i<numnodes;i++){
+	for(int i=1;i<numnodes;i++) if(sedhead[i]<=sedheadmin)sedheadmin=sedhead[i];
+	for(int i=0;i<numnodes;i++){
 		/*If node is now closed bring its thickness back to initial*/
 		if (old_active[i]==0.){
@@ -718,5 +717,5 @@
 		GetHydrologyDCInefficientHmax(&h_max,basalelement,basalelement->nodes[i]);
 		if(eplhead[i]>=h_max && active_element){
-			for(j=0;j<numnodes;j++){
+			for(int j=0;j<numnodes;j++){
 				/*Increase of the domain is on the downstream node in term of sediment head*/
 				if(sedhead[j] == sedheadmin){
@@ -743,5 +742,5 @@
 	/*Constants*/
 	int      domaintype;
-	Element*   basalelement;
+	Element*   basalelement=NULL;
 
 	/*Get basal element*/
@@ -752,5 +751,5 @@
 			break;
 		case Domain3DEnum:
-			if(!element->IsOnBase()) return NULL;
+			if(!element->IsOnBase()) return;
 			basalelement = element->SpawnBasalElement();
 			break;
Index: /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 21430)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyDCInefficientAnalysis.cpp	(revision 21431)
@@ -425,5 +425,5 @@
 	/*Intermediaries*/
 	int						domaintype;
-	Element*			basalelement;
+	Element*			basalelement=NULL;
 	bool					converged;
 	int*					doflist	=	NULL;
@@ -436,5 +436,5 @@
 			break;
 		case Domain3DEnum:
-			if(!element->IsOnBase()) return NULL;
+			if(!element->IsOnBase()) return;
 			basalelement = element->SpawnBasalElement();
 			break;