Index: /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14968)
+++ /issm/trunk-jpl/src/c/classes/objects/Elements/Tria.cpp	(revision 14969)
@@ -1914,8 +1914,9 @@
 		}
 		/*Add input to the element: */
-		this->inputs->AddInput(new TriaP1Input(BasisIntegralEnum,values));
-		
+		this->inputs->AddInput(new TriaP1Input(name,values));
 		/*Free ressources:*/
 		xDelete<int>(doflist);
+		return;
+
 	default:
 	
@@ -6190,5 +6191,4 @@
 		/*Loading term*/
 		water_input->GetInputValue(&water_load,gauss);
-		if(this->id==1)	printf("Qin %e \n ", water_load);
 		scalar = Jdet*gauss->weight*water_load;
 		if(reCast<bool,IssmDouble>(dt)) scalar = scalar*dt;
@@ -6248,5 +6248,4 @@
 		/*Loading term*/
 		residual_input->GetInputValue(&residual_load,gauss);
-		if(this->id==1)	printf("Qres %e \n ", residual_load);
 		scalar = Jdet*gauss->weight*residual_load;
 		if(reCast<bool,IssmDouble>(dt)) scalar = scalar*dt;
@@ -6354,6 +6353,4 @@
 	this->inputs->GetInputValue(&converged,ConvergedEnum);
 	GetInputListOnVertices(&intbasis[0],BasisIntegralEnum);
-
-	if(this->id==1)	printf("val %e \n ", values[1]);
 	
 	if(converged){
@@ -6372,5 +6369,5 @@
 		}
 		if(this->id==1){
-			printf("res  %e val %e h-max %e Stor %e \n ", residual[1], values[1], h_max, sediment_storing);
+			printf("res  %e val %e h-max %e Stor %e intbasis %e\n ", residual[1], values[1], h_max, sediment_storing, intbasis[1]);
 			printf("kappa %e kmax %e pen %e dt %e \n", kappa, kmax,penalty_factor,dt);
 		}
@@ -6463,4 +6460,5 @@
 
 	GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
+	GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
 
 	/* Start looping on the number of gaussian points: */
Index: /issm/trunk-jpl/src/c/solvers/solver_hydro_nonlinear.cpp
===================================================================
--- /issm/trunk-jpl/src/c/solvers/solver_hydro_nonlinear.cpp	(revision 14968)
+++ /issm/trunk-jpl/src/c/solvers/solver_hydro_nonlinear.cpp	(revision 14969)
@@ -33,4 +33,6 @@
 	/*Recover parameters: */
 	femmodel->parameters->FindParam(&isefficientlayer,HydrologydcIsefficientlayerEnum);
+	femmodel->SetCurrentConfiguration(HydrologyDCInefficientAnalysisEnum);//FIXME
+	femmodel->BasisIntegralsx();
 	hydro_maxiter=100;
 	count=1;
@@ -38,4 +40,11 @@
 
 	for(;;){
+		
+		/*Computing the transfer term
+			
+
+		 */
+
+		
 
 		/*First layer*/