Changeset 7314

Timestamp:

02/03/11 16:21:53 (14 years ago)

Author:

Eric.Larour

Message:

improved determination of ice shelf and ice sheet related fields. do not touch under any circunstance...

File:

• issm/trunk/src/m/model/geography.m (modified) (2 diffs)

issm/trunk/src/m/model/geography.m

@@ -32 +32 @@
 elements=md.elements;
 
+%Assign elementoniceshelf, elementonicesheet, gridonicesheet and gridoniceshelf. Only change at your own peril! This is synchronized heavily with the GroundingLineMigration module. {{{1
+
 elementoniceshelf=FlagElements(md,iceshelfname);
 elementonicesheet=FlagElements(md,icesheetname);
@@ -41 +43 @@
 elementoniceshelf=double((elementoniceshelf & ~elementonicesheet));
 elementonicesheet=double(~elementoniceshelf);
+
 %the order here is important. we choose gridoniceshelf as default on the grounding line.
 gridonicesheet(md.elements(find(elementonicesheet),:))=1;
 gridoniceshelf(md.elements(find(elementoniceshelf),:))=1;
+
+%gridoniceshelf exclusively not on ice sheet: 
+gridonicesheet(md.elements(find(elementoniceshelf),:))=0;
+
+%finally, an ice sheet that has one node on an ice shelf, becomes an iceshelf: 
+elementoniceshelf=double(sum(gridoniceshelf(elements),2)>0);
+pos=find(elementoniceshelf); elementonicesheet(pos)=0;
+
+%}}}
 
 %Return: