Index: /issm/trunk/src/c/Makefile.am
===================================================================
--- /issm/trunk/src/c/Makefile.am	(revision 3323)
+++ /issm/trunk/src/c/Makefile.am	(revision 3324)
@@ -402,4 +402,10 @@
 libpISSM_a_SOURCES = ./objects/objects.h\
 					./objects/Object.h\
+					./objects/BamgGeom.h\
+					./objects/BamgGeom.cpp\
+					./objects/BamgMesh.h\
+					./objects/BamgMesh.cpp\
+					./objects/BamgOpts.h\
+					./objects/BamgOpts.cpp\
 					./objects/Element.h\
 					./objects/Element.cpp\
Index: /issm/trunk/src/c/io/WriteData.cpp
===================================================================
--- /issm/trunk/src/c/io/WriteData.cpp	(revision 3323)
+++ /issm/trunk/src/c/io/WriteData.cpp	(revision 3324)
@@ -189,5 +189,5 @@
 	mxArray*    pfield=NULL;
 	mxArray*    pfield2=NULL;
-	int         numfields=15;
+	int         numfields=16;
 	const char* fnames[numfields];
 	int         fsize[numfields][2];
@@ -203,5 +203,4 @@
 	fnames[++i] = "Segments";                 fsize[i][0]=bm->SegmentsSize[0];                  fsize[i][1]=bm->SegmentsSize[1];                 fpointer[i]=&bm->Segments;
 	fnames[++i] = "ElementEdges";             fsize[i][0]=bm->ElementEdgesSize[0];              fsize[i][1]=bm->ElementEdgesSize[1];             fpointer[i]=&bm->ElementEdges;
-	fnames[++i] = "CrackedEdges";             fsize[i][0]=bm->CrackedEdgesSize[0];              fsize[i][1]=bm->CrackedEdgesSize[1];             fpointer[i]=&bm->CrackedEdges;
 	fnames[++i] = "Quadrilaterals";           fsize[i][0]=bm->QuadrilateralsSize[0];            fsize[i][1]=bm->QuadrilateralsSize[1];           fpointer[i]=&bm->Quadrilaterals;
 	fnames[++i] = "VerticesOnGeometricVertex";fsize[i][0]=bm->VerticesOnGeometricVertexSize[0]; fsize[i][1]=bm->VerticesOnGeometricVertexSize[1];fpointer[i]=&bm->VerticesOnGeometricVertex;
@@ -213,4 +212,6 @@
 	fnames[++i] = "NodalConnectivity";        fsize[i][0]=bm->NodalConnectivitySize[0];         fsize[i][1]=bm->NodalConnectivitySize[1];        fpointer[i]=&bm->NodalConnectivity;
 	fnames[++i] = "NodalElementConnectivity"; fsize[i][0]=bm->NodalElementConnectivitySize[0];  fsize[i][1]=bm->NodalElementConnectivitySize[1]; fpointer[i]=&bm->NodalElementConnectivity;
+	fnames[++i] = "CrackedVertices";          fsize[i][0]=bm->CrackedVerticesSize[0];           fsize[i][1]=bm->CrackedVerticesSize[1];          fpointer[i]=&bm->CrackedVertices;
+	fnames[++i] = "CrackedEdges";             fsize[i][0]=bm->CrackedEdgesSize[0];              fsize[i][1]=bm->CrackedEdgesSize[1];             fpointer[i]=&bm->CrackedEdges;
 	assert(i==numfields-1);
 
Index: /issm/trunk/src/c/objects/BamgMesh.cpp
===================================================================
--- /issm/trunk/src/c/objects/BamgMesh.cpp	(revision 3323)
+++ /issm/trunk/src/c/objects/BamgMesh.cpp	(revision 3324)
@@ -9,5 +9,4 @@
 	bamgmesh->ElementEdgesSize[0]=0,  bamgmesh->ElementEdgesSize[1]=0;   bamgmesh->ElementEdges=NULL;
 	bamgmesh->SegmentsSize[0]=0,      bamgmesh->SegmentsSize[1]=0;       bamgmesh->Segments=NULL;
-	bamgmesh->CrackedEdgesSize[0]=0,  bamgmesh->CrackedEdgesSize[1]=0;   bamgmesh->CrackedEdges=NULL;
 	bamgmesh->QuadrilateralsSize[0]=0,bamgmesh->QuadrilateralsSize[1]=0; bamgmesh->Quadrilaterals=NULL;
 	bamgmesh->VerticesOnGeometricVertexSize[0]=0, bamgmesh->VerticesOnGeometricVertexSize[1]=0;bamgmesh->VerticesOnGeometricVertex=NULL;
@@ -20,4 +19,6 @@
 	bamgmesh->NodalConnectivitySize[0]=0,        bamgmesh->NodalConnectivitySize[1]=0;        bamgmesh->NodalConnectivity=NULL;
 	bamgmesh->NodalElementConnectivitySize[0]=0, bamgmesh->NodalElementConnectivitySize[1]=0; bamgmesh->NodalElementConnectivity=NULL;
+	bamgmesh->CrackedVerticesSize[0]=0, bamgmesh->CrackedVerticesSize[1]=0; bamgmesh->CrackedVertices=NULL;
+	bamgmesh->CrackedEdgesSize[0]=0, bamgmesh->CrackedEdgesSize[1]=0; bamgmesh->CrackedEdges=NULL;
 
 }
Index: /issm/trunk/src/c/objects/BamgMesh.h
===================================================================
--- /issm/trunk/src/c/objects/BamgMesh.h	(revision 3323)
+++ /issm/trunk/src/c/objects/BamgMesh.h	(revision 3324)
@@ -21,7 +21,4 @@
 	int     SegmentsSize[2];
 	double* Segments;
-
-	int     CrackedEdgesSize[2];
-	double* CrackedEdges;
 
 	int     QuadrilateralsSize[2];
@@ -53,4 +50,10 @@
 	int     NodalElementConnectivitySize[2];
 	double* NodalElementConnectivity;
+
+	int     CrackedVerticesSize[2];
+	double* CrackedVertices;
+
+	int     CrackedEdgesSize[2];
+	double* CrackedEdges;
 };