Changeset 25515 for issm/trunk-jpl/test/NightlyRun/test2004.m

Timestamp:

09/02/20 14:10:08 (5 years ago)

Author:

jdquinn

Message:

CHG: Using different implementation of delaunay

File:

• issm/trunk-jpl/test/NightlyRun/test2004.m (modified) (5 diffs)

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

@@ -154 +154 @@
                         latAIS=delH{:,2};
                         delHAIS=delH{:,3};
-                        index=delaunay(longAIS,latAIS);
+                        points=[longAIS,latAIS];
+                        index=delaunayn(points);
                         lat=md.mesh.lat; 
                         long=md.mesh.long+360;
@@ -190 +191 @@
 
         %recover lat,long: 
-        [md.mesh.long,md.mesh.lat]=gdaltransform(md.mesh.x,md.mesh.y,md.mesh.proj,'EPSG:4326'); 
+        [md.mesh.long,md.mesh.lat]=gdaltransform(md.mesh.x,md.mesh.y,md.mesh.proj,'EPSG:4326');
 
         %mask:  %{{{
@@ -262 +263 @@
                 latGIS=delH(:,2);
                 delHGIS=delH(:,3);
-                index=delaunay(longGIS,latGIS);
+                points=[longGIS,latGIS];
+                index=delaunayn(points);
                 lat=md.mesh.lat;
                 long=md.mesh.long+360;
@@ -274 +276 @@
                 latGLA=delH(:,2);
                 delHGLA=sum(delH(:,3:end),2);
-                index=delaunay(longGLA,latGLA);
+                points=[longGLA,latGLA];
+                index=delaunayn(points);
                 lat=md.mesh.lat; 
                 long=md.mesh.long+360; 
@@ -409 +412 @@
 md=solve(md,'Sealevelrise');
 Srotation=md.results.SealevelriseSolution.Sealevel;
-
 %}}}