source: issm/oecreview/Archive/24684-25833/ISSM-25515-25516.diff@ 25834

../trunk-jpl/src/m/modules/BamgTriangulate.py

@@ -2 +2 @@
 
 
 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]

../trunk-jpl/test/NightlyRun/test2004.m

@@ -153 +153 @@
                         longAIS=delH{:,1};
                         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;
                         pos=find(long>360);
@@ -262 +263 @@
                 longGIS=delH(:,1);
                 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;
                 pos=find(long>360);
@@ -275 +277 @@
                 longGLA=delH(:,1);
                 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; 
                 pos=find(long>360);