Changeset 9619 for issm/trunk/src/m/model/mesh/meshyams.m

Timestamp:

09/06/11 13:40:47 (14 years ago)

Author:

Mathieu Morlighem

Message:

Added rifts object

File:

• issm/trunk/src/m/model/mesh/meshyams.m (modified) (3 diffs)

issm/trunk/src/m/model/mesh/meshyams.m

@@ -94 +94 @@
         %if we have rifts, we just messed them up, we need to recreate the segments that constitute those 
         %rifts, because the segments are used in YamsCall to freeze the rifts elements during refinement.
-        if md.numrifts, 
+        if md.rifts.numrifts, 
                 md.nodeconnectivity=NodeConnectivity(md.elements,md.numberofnodes);
                 md.elementconnectivity=ElementConnectivity(md.elements,md.nodeconnectivity);
@@ -136 +136 @@
 
 %deal with rifts 
-if md.numrifts,
+if md.rifts.numrifts,
         %first, recreate rift segments
         md=meshyamsrecreateriftsegments(md);
 
         %using the segments, recreate the penaltypairs
-        for j=1:md.numrifts,
-                rift=md.rifts(j);
+        for j=1:md.rifts.numrifts,
+                rift=md.rifts.riftstruct(j);
 
                 %build normals and lengths of segments:
@@ -167 +167 @@
                 rift.penaltypairs(:,6)=rift.penaltypairs(:,6)./norms;
 
-                md.rifts(j)=rift;
+                md.rifts.riftstruct(j)=rift;
         end