Index: /issm/trunk/src/c/modules/GetSolutionFromInputsx/GetSolutionFromInputsx.cpp
===================================================================
--- /issm/trunk/src/c/modules/GetSolutionFromInputsx/GetSolutionFromInputsx.cpp	(revision 4384)
+++ /issm/trunk/src/c/modules/GetSolutionFromInputsx/GetSolutionFromInputsx.cpp	(revision 4385)
@@ -32,4 +32,5 @@
 	/*Get size of vector: */
 	gsize=nodes->NumberOfDofs(analysis_type);
+	if (gsize==0) ISSMERROR("Allocating a Vec of size 0 as gsize=0 for analysis: %s",EnumAsString(analysis_type));
 	
 	/*Initialize solution: */
Index: /issm/trunk/src/c/objects/Elements/Tria.cpp
===================================================================
--- /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4384)
+++ /issm/trunk/src/c/objects/Elements/Tria.cpp	(revision 4385)
@@ -5212,4 +5212,5 @@
 
 	/*Add value to global vector*/
+	printarray(&doflist[0],numdof);
 	VecSetValues(solution,numdof,doflist,(const double*)values,INSERT_VALUES);
 
Index: /issm/trunk/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk/src/c/shared/Elements/elements.h	(revision 4384)
+++ /issm/trunk/src/c/shared/Elements/elements.h	(revision 4385)
@@ -12,5 +12,5 @@
 double Paterson(double temperature);
 int GetVerticesCoordinates(double* xyz,  Node** nodes, int numgrids);
-inline void printarray(double* array,int lines,int cols){
+inline void printarray(double* array,int lines,int cols=1){
 	int i,j;
 	printf("\n");
@@ -24,5 +24,5 @@
 	printf("\n");
 }
-inline void printarray(int* array,int lines,int cols){
+inline void printarray(int* array,int lines,int cols=1){
 	int i,j;
 	printf("\n");