Changeset 4348

Timestamp:

06/30/10 12:16:26 (15 years ago)

Author:

Mathieu Morlighem

Message:

Fixed modelextract (no more md.type)

File:

• issm/trunk/src/m/classes/public/modelextract.m (modified) (5 diffs)

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

@@ -68 +68 @@
 elements_2(:,2)=Pgrid(elements_2(:,2));
 elements_2(:,3)=Pgrid(elements_2(:,3));
-if strcmpi(md1.type,'3d'),
+if md1.dim==3,
         elements_2(:,4)=Pgrid(elements_2(:,4));
         elements_2(:,5)=Pgrid(elements_2(:,5));
@@ -107 +107 @@
 
         %uppernodes lowernodes
-        if strcmpi(md1.type,'3d')
+        if md1.dim==3
                 md2.uppergrids=md1.uppergrids(pos_grid);
                 pos=find(~isnan(md2.uppergrids));
@@ -118 +118 @@
 
         %Initial 2d mesh 
-        if strcmpi(md1.type,'3d')
+        if md1.dim==3
                 flag_elem_2d=flag_elem(1:md1.numberofelements2d);
                 pos_elem_2d=find(flag_elem_2d);
@@ -162 +162 @@
 
         %recreate segments
-        if strcmpi(md1.type,'2d')
+        if md1.dim==2
                 md2.nodeconnectivity=NodeConnectivity(md2.elements,md2.numberofgrids);
                 md2.elementconnectivity=ElementConnectivity(md2.elements,md2.nodeconnectivity);
@@ -215 +215 @@
                 md2.pressureload(:,2)=Pgrid(md1.pressureload(:,2)); 
                 md2.pressureload(:,end-1)=Pelem(md1.pressureload(:,end-1));
-                if strcmpi(md1.type,'3d')
+                if md1.dim==3
                         md2.pressureload(:,3)=Pgrid(md1.pressureload(:,3)); 
                         md2.pressureload(:,4)=Pgrid(md1.pressureload(:,4));