Index: /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 16778)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceVerticalAnalysis.cpp	(revision 16779)
@@ -110,4 +110,13 @@
 	IssmDouble   rho_ice,g;
 
+	/*Get the approximation and do nothing if the element in FS or None*/
+	element->GetInputValue(&approximation,ApproximationEnum);
+	if(approximation==FSApproximationEnum || approximation==NoneApproximationEnum){
+		return;
+	}
+
+	/*Get dof list and vertices coordinates: */
+	element->GetVerticesCoordinates(&xyz_list);
+	element->GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
 	IssmDouble*  values    = xNew<IssmDouble>(numdof);
 	IssmDouble*  vx        = xNew<IssmDouble>(numnodes);
@@ -120,14 +129,4 @@
 	IssmDouble*  pressure  = xNew<IssmDouble>(numnodes);
 	IssmDouble*  surface   = xNew<IssmDouble>(numnodes);
-
-	/*Get the approximation and do nothing if the element in FS or None*/
-	element->GetInputValue(&approximation,ApproximationEnum);
-	if(approximation==FSApproximationEnum || approximation==NoneApproximationEnum){
-		return;
-	}
-
-	/*Get dof list and vertices coordinates: */
-	element->GetVerticesCoordinates(&xyz_list);
-	element->GetDofList(&doflist,NoneApproximationEnum,GsetEnum);
 
 	/*Use the dof list to index into the solution vector vz: */