Changeset 435

Timestamp:

05/15/09 07:07:09 (16 years ago)

Author:

Mathieu Morlighem

Message:

oops...

File:

• issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m (modified) (1 diff)

issm/trunk/src/m/utils/BC/SetMarineIceSheetBC.m

@@ -25 +25 @@
         gridinsideicefront=ContourToMesh(md.elements,md.x,md.y,expread(icefrontfile,1),'node',2);
         gridonicefront=double(md.gridonboundary & gridinsideicefront);
-
 else
         %Guess where the ice front is
-        pos=find(md.gridonboundary & md.gridoniceshelf);
-        md.gridondirichlet_diag=zeros(md.numberofgrids,1);
-        md.gridondirichlet_diag(pos)=1;
+        gridonicefront=double(md.gridonboundary & md.gridoniceshelf);
 end
 pos=find(md.gridonboundary & ~gridonicefront);