Changeset 17733

Timestamp:

04/14/14 22:44:14 (11 years ago)

Author:

cborstad

Message:

NEW: template for python version of SectionValues

Location:

issm/trunk-jpl/src/m/interp

Files:

• SectionValues.m (modified) (2 diffs)
• SectionValues.py (added)

issm/trunk-jpl/src/m/interp/SectionValues.m

@@ -92 +92 @@
         %vertically extrude mesh
 
-        %Get bed and surface for each 2d point, offset to make sure that it is inside the glacier system
+        %Get base and surface for each 2d point, offset to make sure that it is inside the glacier system
         offset=10^-3;
-        bed=InterpFromMeshToMesh2d(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d,project2d(md,md.geometry.base,1),X,Y)+offset;
+        base=InterpFromMeshToMesh2d(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d,project2d(md,md.geometry.base,1),X,Y)+offset;
         surface=InterpFromMeshToMesh2d(md.mesh.elements2d,md.mesh.x2d,md.mesh.y2d,project2d(md,md.geometry.surface,1),X,Y)-offset;
 
@@ -111 +111 @@
                 X3(i:layers:end)=X;
                 Y3(i:layers:end)=Y;
-                Z3(i:layers:end)=bed+(i-1)*(surface-bed)/(layers-1);
+                Z3(i:layers:end)=base+(i-1)*(surface-base)/(layers-1);
                 S3(i:layers:end)=S;