Index: /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 21505)
+++ /issm/trunk-jpl/src/c/analyses/HydrologySommersAnalysis.cpp	(revision 21506)
@@ -392,8 +392,4 @@
 	element->AddInput(HydrologyReynoldsEnum,&reynolds,P0Enum);
 
-	/*Calculate basal flux for output*/
-	IssmDouble q = conductivity*sqrt(dh[0]*dh[0]+dh[1]*dh[1]);
-	element->AddInput(HydrologyBasalFluxEnum,&q,P1Enum);
-
 	/*Free resources:*/
 	xDelete<IssmDouble>(values);
@@ -454,4 +450,5 @@
 	IssmDouble* xyz_list = NULL;
    IssmDouble  dpressure_water[2],dbed[2],PMPheat;
+	IssmDouble q = 0.;
 
 	/*Retrieve all inputs and parameters*/
@@ -540,4 +537,7 @@
 					));
 		totalweights +=gauss->weight*Jdet;
+
+		/* Compute basal water flux */
+      q += gauss->weight*Jdet*(conductivity*sqrt(dh[0]*dh[0]+dh[1]*dh[1]));
 	}
 
@@ -550,8 +550,11 @@
 	if(newgap>thickness)
 	 newgap = thickness;
-	
 	 
 	/*Add new gap as an input*/
 	element->AddInput(HydrologyGapHeightEnum,&newgap,P0Enum);
+ 
+	/*Divide by connectivity, add basal flux as an input*/
+	q = q/totalweights;
+	element->AddInput(HydrologyBasalFluxEnum,&q,P0Enum);
 
 	/*Clean up and return*/