Index: /issm/trunk-jpl/src/m/modules/BamgTriangulate.py
===================================================================
--- /issm/trunk-jpl/src/m/modules/BamgTriangulate.py	(revision 25515)
+++ /issm/trunk-jpl/src/m/modules/BamgTriangulate.py	(revision 25516)
@@ -3,16 +3,15 @@
 
 def BamgTriangulate(x, y):
-    """
-    BAMGTRIANGULATE
+    """BAMGTRIANGULATE
 
     Usage:
         index = BamgTriangulate(x, y)
 
-    index: index of the triangulation
-    x, y: coordinates of the nodes
+        index   : index of the triangulation
+        x, y    : coordinates of the nodes
     """
 
-    # Call mex module
+    # Call Python module
     index = BamgTriangulate_python(x, y)
-    # return
-    return index
+
+    return index[0]
Index: /issm/trunk-jpl/test/NightlyRun/test2004.m
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test2004.m	(revision 25515)
+++ /issm/trunk-jpl/test/NightlyRun/test2004.m	(revision 25516)
@@ -154,6 +154,7 @@
 			latAIS=delH{:,2};
 			delHAIS=delH{:,3};
-			points=[longAIS,latAIS];
-			index=delaunayn(points);
+			% points=[longAIS,latAIS];
+			% index=delaunayn(points);
+			index=BamgTriangulate(longAIS, latAIS);
 			lat=md.mesh.lat; 
 			long=md.mesh.long+360;
@@ -263,6 +264,7 @@
 		latGIS=delH(:,2);
 		delHGIS=delH(:,3);
-		points=[longGIS,latGIS];
-		index=delaunayn(points);
+		% points=[longGIS,latGIS];
+		% index=delaunayn(points);
+		index=BamgTriangulate(longGIS, latGIS);
 		lat=md.mesh.lat;
 		long=md.mesh.long+360;
@@ -276,6 +278,7 @@
 		latGLA=delH(:,2);
 		delHGLA=sum(delH(:,3:end),2);
-		points=[longGLA,latGLA];
-		index=delaunayn(points);
+		% points=[longGLA,latGLA];
+		% index=delaunayn(points);
+		index=BamgTriangulate(longGLA, latGLA);
 		lat=md.mesh.lat; 
 		long=md.mesh.long+360;