Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17246)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 17247)
@@ -2138,5 +2138,5 @@
 	ElementMatrix* Ke = element->NewElementMatrix(HOApproximationEnum);
 	IssmDouble*    B  = xNew<IssmDouble>((dim-1)*numdof);
-	IssmDouble*    D  = xNew<IssmDouble>((dim-1)*(dim-1));
+	IssmDouble*    D  = xNewZeroInit<IssmDouble>((dim-1)*(dim-1));
 
 	/*Retrieve all inputs and parameters*/