Changeset 9733 for issm/trunk/src/m/model/mesh/rifts/meshprocessoutsiderifts.m

Timestamp:

09/09/11 14:58:09 (14 years ago)

Author:

seroussi

Message:

keep building mesh

File:

• issm/trunk/src/m/model/mesh/rifts/meshprocessoutsiderifts.m (modified) (6 diffs)

issm/trunk/src/m/model/mesh/rifts/meshprocessoutsiderifts.m

@@ -12 +12 @@
         
         %first, flag nodes that belong to the domain outline
-        flags=ContourToMesh(md.elements,md.x,md.y,domainoutline,'node',0);
+        flags=ContourToMesh(md.mesh.elements,md.x,md.y,domainoutline,'node',0);
 
         rift=md.rifts.riftstruct(i);
@@ -43 +43 @@
                 while  flags(B), %as long as B does not belong to the domain outline, keep looking.
                         %detect elements on edge A,B:
-                        edgeelements=ElementsFromEdge(md.elements,A,B);
+                        edgeelements=ElementsFromEdge(md.mesh.elements,A,B);
                         %rule out those we already detected
                         already_detected=ismember(edgeelements,elements);
@@ -50 +50 @@
                         elements=[elements;nextelement];
                         %new B:
-                        B=md.elements(nextelement,find(~ismember(md.elements(nextelement,:),[A B])));
+                        B=md.mesh.elements(nextelement,find(~ismember(md.mesh.elements(nextelement,:),[A B])));
                 end
                 
@@ -61 +61 @@
                 
                 %replace tip in elements
-                newelements=md.elements(elements,:);
+                newelements=md.mesh.elements(elements,:);
                 pos=find(newelements==tip);
                 newelements(pos)=num;
-                md.elements(elements,:)=newelements;
+                md.mesh.elements(elements,:)=newelements;
                 md.rifts.riftstruct(i).tips=[md.rifts.riftstruct(i).tips num];
 
@@ -73 +73 @@
                         pos=find(md.mesh.segments(segment_index,1:2)~=tip);
                         other_node=md.mesh.segments(segment_index,pos);
-                        if ~isconnected(md.elements,other_node,tip),
+                        if ~isconnected(md.mesh.elements,other_node,tip),
                                 pos=find(md.mesh.segments(segment_index,1:2)==tip);
                                 md.mesh.segments(segment_index,pos)=num;
@@ -83 +83 @@
 
 %Fill in rest of fields:
-md.mesh.numberofelements=length(md.elements);
+md.mesh.numberofelements=length(md.mesh.elements);
 md.mesh.numberofvertices=length(md.x);
 md.z=zeros(md.mesh.numberofvertices,1);