Changeset 2688

Timestamp:

11/30/09 15:32:44 (15 years ago)

Author:

Mathieu Morlighem

Message:

Added rift outline

File:

• issm/trunk/src/m/classes/public/mesh/meshyams.m (modified) (2 diffs)

issm/trunk/src/m/classes/public/mesh/meshyams.m

@@ -29 +29 @@
 domainoutline=getfieldvalueerr(options,'domainoutline');
 disp(sprintf('   %-15s: ''%s''','DomainOutline',domainoutline));
+riftoutline=getfieldvalue(options,'riftoutline','N/A');
+disp(sprintf('   %-15s: ''%s''','riftoutline',riftoutline));
 groundeddomain=getfieldvalue(options,'groundeddomain','N/A');
 disp(sprintf('   %-15s: ''%s''','GroundedDomain',groundeddomain));
@@ -48 +50 @@
 %mesh with initial resolution
 disp('Initial mesh generation...');
-md=mesh(md,domainoutline,resolution);
+if strcmpi(riftoutline,'N/A');
+        md=mesh(md,domainoutline,resolution);
+else
+        md=mesh(md,domainoutline,riftoutline,resolution);
+end
 disp(['Initial mesh, number of elements: ' num2str(md.numberofelements)]);