Index: /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 23507)
+++ /issm/trunk-jpl/src/c/shared/Elements/PrintArrays.cpp	(revision 23508)
@@ -49,4 +49,16 @@
 	_printf_("\n");
 }
+void printarray_matlab(const char* filename,int* array,int lines,int cols){
+	FILE *f = fopen(filename,"w");
+	fprintf(f,"%% Matrix of size %ix%i\n",lines,cols);
+	fprintf(f,"\n");
+	fprintf(f,"A=[...\n");
+	for(int i=0;i<lines;i++){  
+		for(int j=0;j<cols;j++) fprintf(f," %i",array[i*cols+j]);
+		fprintf(f,"\n");
+	}  
+	fprintf(f,"];\n");
+	fclose(f);
+}
 void printarray(bool* array,int lines,int cols){
 	_printf_("\n");
Index: /issm/trunk-jpl/src/c/shared/Elements/elements.h
===================================================================
--- /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 23507)
+++ /issm/trunk-jpl/src/c/shared/Elements/elements.h	(revision 23508)
@@ -49,4 +49,5 @@
 void printarray(IssmDouble* array,int lines,int cols=1);
 #endif
+void printarray_matlab(const char* filename,int* array,int lines,int cols=1);
 void printarray(int* array,int lines,int cols=1);
 void printarray(bool* array,int lines,int cols=1);