Index: /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp
===================================================================
--- /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18217)
+++ /issm/trunk-jpl/src/c/analyses/StressbalanceAnalysis.cpp	(revision 18218)
@@ -7,5 +7,5 @@
 #include "../cores/cores.h"
 
-//#define FSANALYTICAL 2
+//#define FSANALYTICAL 24
 
 /*Model processing*/
@@ -3188,4 +3188,5 @@
 	IssmDouble  x_coord,y_coord,z_coord;
 	IssmDouble *xyz_list_base = NULL;
+	IssmDouble *xyz_list      = NULL;
 	Gauss*      gauss         = NULL;
 
@@ -3205,4 +3206,5 @@
 	/*Retrieve all inputs and parameters*/
 	element->GetVerticesCoordinatesBase(&xyz_list_base);
+	element->GetVerticesCoordinates(&xyz_list);
 	Input* vx_input         = element->GetInput(VxEnum);      _assert_(vx_input);
 	Input* vy_input         = element->GetInput(VyEnum);      _assert_(vy_input);
@@ -3215,7 +3217,7 @@
 		gauss->GaussPoint(ig);
 
-		x_coord=element->GetXcoord(xyz_list_base,gauss);
-		y_coord=element->GetYcoord(xyz_list_base,gauss);
-		if(dim==3) z_coord=element->GetZcoord(xyz_list_base,gauss);
+		x_coord=element->GetXcoord(xyz_list,gauss);
+		y_coord=element->GetYcoord(xyz_list,gauss);
+		if(dim==3) z_coord=element->GetZcoord(xyz_list,gauss);
 		else z_coord=0.;
 
@@ -3236,4 +3238,5 @@
 	/*Clean up and return*/
 	delete gauss;
+	xDelete<IssmDouble>(xyz_list);
 	xDelete<IssmDouble>(xyz_list_base);
 	xDelete<IssmDouble>(B);