Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18150)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18151)
@@ -4410,4 +4410,7 @@
 		else{
 			_assert_(tnumnodes==4);
+			d_zz = xNew<IssmDouble>(tnumnodes);
+			d_xz = xNew<IssmDouble>(tnumnodes);
+			d_yz = xNew<IssmDouble>(tnumnodes);
 			Matrix4x4Solve(&d_xx[0],Ke,pe_xx);
 			Matrix4x4Solve(&d_yy[0],Ke,pe_yy);
@@ -4495,4 +4498,9 @@
 		IssmDouble* tau_xz = NULL;
 		IssmDouble* tau_yz = NULL;
+		if(dim==3){
+			tau_zz = xNew<IssmDouble>(tnumnodes);
+			tau_xz = xNew<IssmDouble>(tnumnodes);
+			tau_yz = xNew<IssmDouble>(tnumnodes);
+		}
 		Gauss* gauss = element->NewGauss();
 		for(int ig=0;ig<tnumnodes;ig++){