Index: /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp	(revision 22306)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceSIAAnalysis.cpp	(revision 22307)
@@ -331,5 +331,6 @@
 	IssmDouble  gravity    = element->GetMaterialParameter(ConstantsGEnum);
 	IssmDouble  n          = element->GetMaterialParameter(MaterialsRheologyNEnum);
-	IssmDouble  B          = element->GetMaterialParameter(MaterialsRheologyBbarEnum);
+	IssmDouble B;
+	Input* B_input         = element->GetInput(MaterialsRheologyBbarEnum);_assert_(B_input);
 	Input* slopex_input    = element->GetInput(SurfaceSlopeXEnum);        _assert_(slopex_input);
 	Input* slopey_input    = element->GetInput(SurfaceSlopeYEnum);        _assert_(slopey_input);
@@ -347,4 +348,5 @@
 		connectivity=(IssmDouble)element->VertexConnectivity(iv);
 
+		B_input->GetInputValue(&B,gauss);
 		thickness_input->GetInputValue(&thickness,gauss);
 		surface_input->GetInputValue(&surface,gauss);
@@ -418,5 +420,6 @@
 	IssmDouble  gravity    = element->GetMaterialParameter(ConstantsGEnum);
 	IssmDouble  n          = element->GetMaterialParameter(MaterialsRheologyNEnum);
-	IssmDouble  B          = element->GetMaterialParameter(MaterialsRheologyBEnum);
+	IssmDouble B;
+	Input* B_input         = element->GetInput(MaterialsRheologyBEnum);   _assert_(B_input);
 	Input* surface_input   = element->GetInput(SurfaceEnum);              _assert_(surface_input);
 	Input* slopex_input    = element->GetInput(SurfaceSlopeXEnum);        _assert_(slopex_input);
@@ -448,4 +451,5 @@
 			gauss->GaussPoint(ig);
 
+			B_input->GetInputValue(&B,gauss);
 			slopex_input->GetInputValue(&slope[0],gauss);
 			slopey_input->GetInputValue(&slope[1],gauss);
Index: /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 22306)
+++ /issm/trunk-jpl/src/c/classes/Materials/Matice.cpp	(revision 22307)
@@ -217,4 +217,6 @@
 IssmDouble Matice::GetB(Gauss* gauss){/*{{{*/
 
+	_assert_(gauss); 
+
 	/*Output*/
 	IssmDouble B;
@@ -226,4 +228,6 @@
 /*}}}*/
 IssmDouble Matice::GetBbar(Gauss* gauss){/*{{{*/
+
+	_assert_(gauss); 
 
 	/*Output*/