Changeset 15993

Timestamp:

08/28/13 10:32:57 (12 years ago)

Author:

Mathieu Morlighem

Message:

BUG: ice front needs to be extended so that it covers the elements that are partially floating

File:

• issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m (modified) (1 diff)

issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m

@@ -23 +23 @@
                 error(['SetMarineIceSheetBC error message: ice front file ' icefrontfile ' not found']);
         end
-        nodeinsideicefront=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,icefrontfile,'node',2);
-        vertexonicefront=double(md.mesh.vertexonboundary & nodeinsideicefront);
+        incontour=ContourToMesh(md.mesh.elements,md.mesh.x,md.mesh.y,icefrontfile,'node',2);
+        vertexonicefront=double(md.mesh.vertexonboundary & incontour);
 else
         %Guess where the ice front is
-        vertexonfloatingice=(md.mask.groundedice_levelset<0.);
+        pos=find(sum(md.mask.groundedice_levelset(md.mesh.elements)<0.,2) >0.);
+        vertexonfloatingice=zeros(md.mesh.numberofvertices,1);
+        vertexonfloatingice(md.mesh.elements(pos,:))=1.;
         vertexonicefront=double(md.mesh.vertexonboundary & vertexonfloatingice);
 end