Index: /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 22260)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 22261)
@@ -139,4 +139,7 @@
 ElementMatrix* StressbalanceVerticalAnalysis::CreateKMatrix(Element* element){/*{{{*/
 
+	/* Check if ice in element */
+	if(!element->IsIceInElement()) return NULL;
+
 	bool hack = false;
 
@@ -278,4 +281,7 @@
 }/*}}}*/
 ElementVector* StressbalanceVerticalAnalysis::CreatePVector(Element* element){/*{{{*/
+
+	/* Check if ice in element */
+	if(!element->IsIceInElement()) return NULL;
 
 	bool hack = false;
@@ -617,4 +623,5 @@
 void           StressbalanceVerticalAnalysis::UpdateConstraints(FemModel* femmodel){/*{{{*/
 	/*Default, do nothing*/
+	SetActiveNodesLSMx(femmodel);
 	return;
 }/*}}}*/
Index: /issm/trunk-jpl/src/c/cores/movingfront_core.cpp
===================================================================
--- /issm/trunk-jpl/src/c/cores/movingfront_core.cpp	(revision 22260)
+++ /issm/trunk-jpl/src/c/cores/movingfront_core.cpp	(revision 22261)
@@ -40,5 +40,5 @@
 	/* determine variables for extrapolation */
 	num_extrapol_vars=0;
-	if(isstressbalance) num_extrapol_vars+=2;
+	if(isstressbalance) num_extrapol_vars+=3;
 	if(ismasstransport) num_extrapol_vars+=1;
 	if(isthermal && domaintype==Domain3DEnum) num_extrapol_vars+=1;
@@ -48,4 +48,5 @@
 		extrapol_vars[index]=VxEnum; index++;
 		extrapol_vars[index]=VyEnum; index++;
+		extrapol_vars[index]=VzEnum; index++;
 	}
 	if(ismasstransport){extrapol_vars[index]=ThicknessEnum; index++;}