Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 21402)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 21403)
@@ -3898,6 +3898,5 @@
 
 	int         i,dim;
-	IssmDouble  Jdet,water_pressure,bed;
-	IssmDouble	normal[3];
+	IssmDouble  Jdet,water_pressure,base;
 	IssmDouble *xyz_list_base = NULL;
 
@@ -3936,8 +3935,6 @@
 		element->NodalFunctionsVelocity(vbasis,gauss);
 
-		element->NormalBase(&normal[0],xyz_list_base);
-		_assert_(normal[dim-1]<0.);
-		base_input->GetInputValue(&bed, gauss);
-		water_pressure=gravity*rho_water*bed;
+		base_input->GetInputValue(&base, gauss);
+		water_pressure=gravity*rho_water*base;
 
 		for(i=0;i<vnumnodes;i++){
@@ -3951,5 +3948,5 @@
 	if(shelf_dampening) {
 		Input*      mb_input=element->GetInput(BasalforcingsFloatingiceMeltingRateEnum); _assert_(mb_input);
-		IssmDouble dt,mb,normal_b;
+		IssmDouble  dt,mb;
 		element->FindParam(&dt,TimesteppingTimeStepEnum);
 		for(int ig=gauss->begin();ig<gauss->end();ig++){
Index: /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 21402)
+++ /issm/trunk-jpl/src/c/classes/Elements/Penta.cpp	(revision 21403)
@@ -717,5 +717,5 @@
 
 		/*Create gauss point in the middle of the basal edge*/
-		Gauss* gauss=NewGauss();
+		Gauss* gauss=NewGaussBase(1);
 		gauss->GaussPoint(0);