Changeset 1181

Timestamp:

06/30/09 11:00:19 (16 years ago)

Author:

Mathieu Morlighem

Message:

added default value (NaN)

File:

• issm/trunk/src/m/classes/public/SectionValues.m (modified) (3 diffs)

issm/trunk/src/m/classes/public/SectionValues.m

@@ -72 +72 @@
 
         %Interpolation of data on specified points
-        data_interp=InterpFromMesh2d(md.elements,md.x,md.y,data,X,Y);
+        data_interp=InterpFromMesh2d(md.elements,md.x,md.y,data,X,Y,NaN);
 
         %Compute index
@@ -83 +83 @@
         %Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system
         offset=10^-10;
-        bed=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y)+offset;
-        surface=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y)-offset;
+        bed=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.bed,1),X,Y,NaN)+offset;
+        surface=InterpFromMesh2d(md.elements2d,md.x2d,md.y2d,project2d(md,md.surface,1),X,Y,NaN)-offset;
 
         %Some useful parameters
@@ -109 +109 @@
 
         %Interpolation of data on specified points
-        data_interp=InterpFromMesh3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3);
+        data_interp=InterpFromMesh3d(md.elements,md.x,md.y,md.z,data,X3,Y3,Z3,NaN);
 
         %build outputs