Index: /issm/trunk-jpl/src/c/bamg/Mesh.cpp
===================================================================
--- /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 21298)
+++ /issm/trunk-jpl/src/c/bamg/Mesh.cpp	(revision 21299)
@@ -500,7 +500,7 @@
 			subdomains = new SubDomain [ nbsubdomains ];
 			for (i=0;i<nbsubdomains;i++) {
-				i3  =(int)bamgmesh->SubDomains[i*3+0];
-				head=(int)bamgmesh->SubDomains[i*3+1]-1;//C indexing
-				direction=(int)bamgmesh->SubDomains[i*3+2];
+				i3  =(int)bamgmesh->SubDomains[i*4+0];
+				head=(int)bamgmesh->SubDomains[i*4+1]-1;//C indexing
+				direction=(int)bamgmesh->SubDomains[i*4+2];
 				if (i3!=3) _error_("Bad Subdomain definition: first number should be 3");
 				if (head<0 || head>=nbt) _error_("Bad Subdomain definition: head should in [1 " << nbt << "] (triangle number)");