source: issm/oecreview/Archive/23390-24306/ISSM-24110-24111.diff

../trunk-jpl/src/m/parameterization/killberg.m

@@ -16 +16 @@
 mask         = zeros(md.mesh.numberofvertices,1);
 element_flag = zeros(md.mesh.numberofelements,1);
 
-%Step 1, go through all elements and put 1 in local_mask if the element is grounded
 disp('Looking for isolated patches of floating ice (icebergs)');
-isice      = min(md.mask.ice_levelset(md.mesh.elements),[],2)<0;
+
+%do not go through elements that don't have ice, mark flag as 1 (done)
+isice = min(md.mask.ice_levelset(md.mesh.elements),[],2)<0;
+element_flag(find(~isice)) = 1;
+
+%do not go through elements that are grounded, mark flag as 1 (done)
+%and initialize mask as 1 for all vertices of these elements
 isgrounded = max(md.mask.groundedice_levelset(md.mesh.elements),[],2)>0;
 pos = find(isgrounded);
 element_flag(pos) = 1;
@@ -30 +35 @@
         disp(['   -- iteration ' num2str(iter)]);
         more = false;
 
-        pos = find(~element_flag & isice);
-        if numel(pos)==0
-                break;
-        end
-
-        for i=pos'
+        for i=find(~element_flag)'
                 indices = md.mesh.elements(i,:);
                 MAX=max(mask(indices));
                 if(MAX==0)