Changeset 1266

Timestamp:

07/08/09 20:05:56 (16 years ago)

Author:

Eric.Larour

Message:

Added threshold

File:

• issm/trunk/src/m/classes/public/mesh/meshadaptation.m (modified) (4 diffs)

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

@@ -1 @@
-function md2=meshadaptation(md,field,scale,epsilon)
+function md2=meshadaptation(md,field,scale,epsilon,threshold)
 %MESHADAPTATION - remesh a model for a given field
 %
@@ -5 +5 @@
 %   scale: 1.2 large scale leads to higher refining.
 %   epsilon: .5  large epsilon will refine fewer elements
+%   threshold: value under which no refinement is done.
 %
 %   Usage:
-%      md2=meshadaptation(md,field,scale,epsilon)
+%      md2=meshadaptation(md,field,scale,epsilon,threshold)
 %
 %   Example:
@@ -29 +30 @@
 grady=zeros(numberofgrids,1);
 metric=zeros(numberofelements,1);
+
+%take the threshold, a characteristic length, and get the area equivalent threshold:
+threshold=threshold^2/2;
 
 %build some usefull variables
@@ -83 +87 @@
 
 %ok, metric can't go under 1km resolution.
-metric(find(metric<10^6))=10^6;
+metric(find(metric<threshold))=threshold;
 
 %Remesh with this new metric