Index: /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp
===================================================================
--- /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp	(revision 13241)
+++ /issm/trunk-jpl/src/c/matlab/io/WriteMatlabData.cpp	(revision 13242)
@@ -9,10 +9,8 @@
 #endif
 
+#include <mex.h>
 #include "../../include/include.h"
 #include "../../shared/shared.h"
 #include "./matlabio.h"
-
-#include <mex.h>
-
 
 /*Primitive data types*/
Index: /issm/trunk-jpl/src/c/matlab/io/matlabio.h
===================================================================
--- /issm/trunk-jpl/src/c/matlab/io/matlabio.h	(revision 13241)
+++ /issm/trunk-jpl/src/c/matlab/io/matlabio.h	(revision 13242)
@@ -12,9 +12,8 @@
 #endif 
 
+#include <mex.h>
 #include "../../classes/classes.h"
 #include "../../Container/Container.h"
 #include "../../include/include.h"
-
-#include <mex.h>
 
 void WriteData(mxArray** pdataref,SeqMat<double>* matrix);
Index: /issm/trunk-jpl/src/c/python/io/WritePythonData.cpp
===================================================================
--- /issm/trunk-jpl/src/c/python/io/WritePythonData.cpp	(revision 13241)
+++ /issm/trunk-jpl/src/c/python/io/WritePythonData.cpp	(revision 13242)
@@ -20,5 +20,5 @@
 #include "../../EnumDefinitions/EnumDefinitions.h"
 
-
+/*Primitive data types*/
 /*FUNCTION WriteData(PyObject* py_tuple,int index,int integer){{{*/
 void WriteData(PyObject* py_tuple, int index, int integer){
@@ -33,6 +33,32 @@
 
 }/*}}}*/
-/*FUNCTION WriteData(PyObject* tuple,int index,SeqMat<double>* matrix){{{*/
-void WriteData(PyObject* tuple,int index,SeqMat<double>* matrix){
+/*FUNCTION WriteData(PyObject* py_tuple,int index, double* matrix, int M, int N){{{*/
+void WriteData(PyObject* tuple, int index, double* matrix, int M,int N){
+
+	npy_intp dims[2]={0,0};
+	PyObject* array=NULL;
+
+	dims[0]=(npy_intp)M;
+	dims[1]=(npy_intp)N;
+	array=PyArray_SimpleNewFromData(2,dims,NPY_DOUBLE,matrix);
+
+	PyTuple_SetItem(tuple, index, array);
+}/*}}}*/
+
+/*ISSM objects*/
+/*FUNCTION WriteData(PyObject* py_tuple,int index,BamgGeom* bamggeom){{{*/
+void WriteData(PyObject* py_tuple,int index,BamgGeom* bamggeom){
+
+	_error_("not implemented yes, see ../../matlab/io/WriteMatlabData.cpp");
+}
+/*}}}*/
+/*FUNCTION WriteData(PyObject* py_tuple,int index,BamgMesh* bamgmesh){{{*/
+void WriteData(PyObject* py_tuple,int index,BamgMesh* bamgmesh){
+
+	_error_("not implemented yes, see ../../matlab/io/WriteMatlabData.cpp");
+}
+/*}}}*/
+/*FUNCTION WriteData(PyObject* py_tuple,int index,SeqMat<double>* matrix){{{*/
+void WriteData(PyObject* py_tuple,int index,SeqMat<double>* matrix){
 	
 	int M,N;
@@ -47,5 +73,5 @@
 	array=PyArray_SimpleNewFromData(2,dims,NPY_DOUBLE,buffer);
 	
-	PyTuple_SetItem(tuple, index, array);
+	PyTuple_SetItem(py_tuple, index, array);
 
 }/*}}}*/
@@ -64,19 +90,4 @@
 	
 	PyTuple_SetItem(tuple, index, array);
-
-
 }
 /*}}}*/
-/*FUNCTION WriteData(PyObject* py_tuple,int index, double* matrix, int M, int N){{{*/
-void WriteData(PyObject* tuple, int index, double* matrix, int M,int N){
-	
-	npy_intp dims[2]={0,0};
-	PyObject* array=NULL;
-	
-	dims[0]=(npy_intp)M;
-	dims[1]=(npy_intp)N;
-	array=PyArray_SimpleNewFromData(2,dims,NPY_DOUBLE,matrix);
-	
-	PyTuple_SetItem(tuple, index, array);
-
-}/*}}}*/
Index: /issm/trunk-jpl/src/c/python/io/pythonio.h
===================================================================
--- /issm/trunk-jpl/src/c/python/io/pythonio.h	(revision 13241)
+++ /issm/trunk-jpl/src/c/python/io/pythonio.h	(revision 13242)
@@ -12,37 +12,16 @@
 #endif 
 
-
 #include "../../classes/classes.h"
 #include "../../Container/Container.h"
 #include "../../include/include.h"
 
-class DataSet;
-class Parameters;
+void WriteData(PyObject* py_tuple,int index, double* matrix, int M,int N);
+void WriteData(PyObject* py_tuple,int index, int integer);
+void WriteData(PyObject* py_tuple,int index, char* string);
+void WriteData(PyObject* py_tuple,int index, SeqMat<double>* matrix);
+void WriteData(PyObject* py_tuple,int index, SeqVec<double>* vector);
+void WriteData(PyObject* py_tuple,int index, BamgGeom* bamggeom);
+void WriteData(PyObject* py_tuple,int index, BamgMesh* bamgmesh);
 
-//void WriteData(PyObject* py_tuple,DataSet* dataset);
-//void WriteData(PyObject* py_tuple,int*    matrix, int M,int N);
-//void WriteData(PyObject* py_tuple,double* vector, int M);
-//void WriteData(PyObject* py_tuple,int integer);
-//void WriteData(PyObject* py_tuple,bool boolean);
-//void WriteData(PyObject* py_tuple,double scalar);
-//void WriteData(DataHandle* py_tuple,Parameters* parameters);
-void WriteData(PyObject* py_tuple, int index, double* matrix, int M,int N);
-void WriteData(PyObject* py_tuple, int index, int integer);
-void WriteData(PyObject* py_tuple, int index, char* string);
-void WriteData(PyObject* py_tuple, int index, SeqMat<double>* matrix);
-void WriteData(PyObject* py_tuple, int index, SeqVec<double>* vector);
-
-//void FetchData(DataSet** pdataset,PyObject* py_ref);
-//void FetchData(double** pmatrix,int* pnumel,int* pndims,int** psize,PyObject* py_ref);
-//void FetchData(int** pmatrix,int* pM,int *pN,PyObject* py_ref);
-//void FetchData(bool** pmatrix,int* pM,int *pN,PyObject* py_ref);
-//void FetchData(bool** pmatrix,int* pnumel,int* pndims,int** psize,PyObject* py_ref);
-//void FetchData(Matrix** pmatrix,PyObject* py_ref);
-//void FetchData(int** pvector,int* pM,PyObject* py_ref);
-//void FetchData(float** pvector,int* pM,PyObject* py_ref);
-//void FetchData(bool** pvector,int* pM,PyObject* py_ref);
-//void FetchData(Vector** pvector,PyObject* py_ref);
-//void FetchData(char** pmatrix,int* pnumel,int* pndims,int** psize,PyObject* py_ref);
-//void FetchData(Parameters** pparameters, DataHandle py_ref);
 void FetchData(double** pvector,int* pM,PyObject* py_ref);
 void FetchData(double** pmatrix,int* pM,int *pN,PyObject* py_array);