Index: /issm/trunk-jpl/src/c/classes/DofIndexing.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/DofIndexing.cpp	(revision 15644)
+++ /issm/trunk-jpl/src/c/classes/DofIndexing.cpp	(revision 15645)
@@ -167,4 +167,6 @@
 	_printf_("DofIndexing:\n");
 	_printf_("   gsize:  " << gsize << "\n");
+	_printf_("   fsize:  " << fsize << "\n");
+	_printf_("   ssize:  " << ssize << "\n");
 	_printf_("   clone:  " << clone << "\n");
 	_printf_("   active: " << active << "\n");
Index: /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp
===================================================================
--- /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15644)
+++ /issm/trunk-jpl/src/c/classes/Elements/Tria.cpp	(revision 15645)
@@ -3166,11 +3166,10 @@
 	IssmDouble  normal[2];
 
-	/*Fetch number of nodes and dof for this finite element*/
+	/*Fetch number of nodes for this finite element*/
 	int numnodes = this->NumberofNodes();
-	int numdof   = numnodes*NDOF2;
 
 	/*Initialize Element vector and other vectors*/
 	GetVerticesCoordinates(&xyz_list[0][0],vertices,NUMVERTICES);
-	ElementVector* pe    = new ElementVector(nodes,numnodes,this->parameters,HOApproximationEnum);
+	ElementVector* pe    = new ElementVector(nodes,numnodes,this->parameters,SSAApproximationEnum);
 	IssmDouble*    basis = xNew<IssmDouble>(numnodes);