Index: /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 24390)
+++ /issm/trunk-jpl/src/c/analyses/HydrologyShaktiAnalysis.cpp	(revision 24391)
@@ -412,5 +412,5 @@
 
 	/*Get gap height derivatives at the center of the element*/
-	Gauss* gauss=element->NewGauss(0);
+	Gauss* gauss=element->NewGauss(1);
 	head_input->GetInputDerivativeValue(&dh[0],xyz_list,gauss);
 	delete gauss;
Index: /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp	(revision 24390)
+++ /issm/trunk-jpl/src/c/shared/Numerics/GaussPoints.cpp	(revision 24391)
@@ -20,4 +20,5 @@
 		For degree p, the required number of Gauss-Legendre points is
 		n>=(p+1)/2.*/
+	_assert_(ngaus>0);
 
 	/*Intermediaries*/
@@ -98,4 +99,5 @@
 	  Symmetrical Gaussian Quadrature Rules for the Triangle", IJNME,
 	  Vol. 21, pp. 1129-1148 (1985), as transcribed for Probe rules3.*/
+	_assert_(iord>0);
 
 	/*Intermediaries*/
@@ -1218,4 +1220,5 @@
 		Quadrature Formulas", Computer Methods in Applied Mechanics and
 		Engineering, Vol. 55, pp. 339-348 (1986).*/
+	_assert_(iord>0);
 
 	/*Intermediaries*/
@@ -1484,4 +1487,5 @@
 	  For degree p, the required number of Gauss-Lobatto points is
 	  n>=(p+1)/2+1 (one more than Gauss-Legendre).*/
+	_assert_(ngaus>0);
 
 	int i;