Changeset 9640 for issm/trunk/src/m/model/mask.m

Timestamp:

09/07/11 08:24:03 (14 years ago)

Author:

Mathieu Morlighem

Message:

temporary commit

File:

• issm/trunk/src/m/model/mask.m (copied) (copied from issm/trunk/src/m/model/geography.m ) (2 diffs)

issm/trunk/src/m/model/mask.m

@@ -1 @@
-function md=geography(md,iceshelfname,icesheetname)
+function md=setmask(md,floatingicename,groundedicename)
 %GEOGRAPHY - establish boundaries between grounded and floating ice.
 %
-%   By default, ice is considered grounded. The contour iceshelfname defines nodes 
-%   for which ice is floating. The contour icesheetname defines nodes inside an iceshelf, 
+%   By default, ice is considered grounded. The contour floatingicename defines nodes 
+%   for which ice is floating. The contour groundedicename defines nodes inside an floatingice, 
 %   that are grounded (ie: ice rises, islands, etc ...)
 %   All input files are in the Argus format (extension .exp).
 %
 %   Usage:
-%      md=geography(md,iceshelfname,icesheetname)
+%      md=setmask(md,floatingicename,groundedicename)
 %
 %   Examples:
-%      md=geography(md,'all','');
-%      md=geography(md,'Iceshelves.exp','Islands.exp');
+%      md=setmask(md,'all','');
+%      md=setmask(md,'Iceshelves.exp','Islands.exp');
 
 %some checks on list of arguments
 if ((nargin~=3) | (nargout~=1)),
-        help geography
-        error('geography error message');
+        help mask
+        error('mask error message');
 end
 
@@ -25 +25 @@
 elements=md.elements;
 
-%Assign elementoniceshelf, elementonicesheet, nodeonicesheet and nodeoniceshelf. Only change at your own peril! This is synchronized heavily with the GroundingLineMigration module. {{{1
-elementoniceshelf=FlagElements(md,iceshelfname);
-elementonicesheet=FlagElements(md,icesheetname);
+%Assign elementonfloatingice, elementongroundedice, vertexongroundedice and vertexonfloatingice. Only change at your own peril! This is synchronized heavily with the GroundingLineMigration module. {{{1
+elementonfloatingice=FlagElements(md,floatingicename);
+elementongroundedice=FlagElements(md,groundedicename);
 
-%Because icesheet nodes and elements can be included into an iceshelf, we need to update. Remember, all the previous 
+%Because groundedice nodes and elements can be included into an floatingice, we need to update. Remember, all the previous 
 %arrays come from domain outlines that can intersect one another: 
-elementoniceshelf=double((elementoniceshelf & ~elementonicesheet));
-elementonicesheet=double(~elementoniceshelf);
+elementonfloatingice=double((elementonfloatingice & ~elementongroundedice));
+elementongroundedice=double(~elementonfloatingice);
 
-%the order here is important. we choose nodeonicesheet as default on the grounding line.
-nodeoniceshelf=zeros(md.numberofnodes,1);
-nodeonicesheet=zeros(md.numberofnodes,1);
-nodeonicesheet(md.elements(find(elementonicesheet),:))=1;
-nodeoniceshelf(find(~nodeonicesheet))=1;
+%the order here is important. we choose vertexongroundedice as default on the grounding line.
+vertexonfloatingice=zeros(md.numberofnodes,1);
+vertexongroundedice=zeros(md.numberofnodes,1);
+vertexongroundedice(md.elements(find(elementongroundedice),:))=1;
+vertexonfloatingice(find(~vertexongroundedice))=1;
 %}}}
 
 %Return: 
-md.elementoniceshelf=elementoniceshelf;
-md.nodeoniceshelf=nodeoniceshelf;
-
-md.elementonicesheet=elementonicesheet;
-md.nodeonicesheet=nodeonicesheet;
-
-md.nodeonwater=zeros(md.numberofnodes,1);
-md.elementonwater=zeros(md.numberofelements,1);
+md.mask.elementonfloatingice=elementonfloatingice;
+md.mask.vertexonfloatingice=vertexonfloatingice;
+md.mask.elementongroundedice=elementongroundedice;
+md.mask.vertexongroundedice=vertexongroundedice;
+md.mask.vertexonwater=zeros(md.numberofnodes,1);
+md.mask.elementonwater=zeros(md.numberofelements,1);