Changeset 2744

Timestamp:

01/04/10 14:33:53 (15 years ago)

Author:

Eric.Larour

Message:

Allow presence of rifts in YamsCall

File:

• issm/trunk/src/m/utils/Mesh/YamsCall.m (modified) (1 diff)

issm/trunk/src/m/utils/Mesh/YamsCall.m

@@ -55 +55 @@
 fprintf(fid,'%i %i %i %i\n',[md.elements zeros(md.numberofelements,1)]');
 numberofelements1=md.numberofelements;
+        
+%Deal with rifts
+if md.numrifts,
+        
+        %we have the list of triangles that make up the rift. keep those triangles around during refinement.
+        triangles=[];
+        for i=1:md.numrifts,
+                triangles=[triangles md.rifts(i).segments(:,3)'];
+        end
+
+        fprintf(fid,'\n\n%s\n%i\n\n','RequiredTriangles',length(triangles));
+        fprintf(fid,'%i\n',triangles);
+end
 
 %close