Changeset 7335

Timestamp:

02/04/11 09:48:26 (14 years ago)

Author:

seroussi

Message:

fixed SetMarineIceSheetBC for the new geography (do not use gridoniceshelf)

File:

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

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

@@ -27 +27 @@
 else
         %Guess where the ice front is
-        gridonicefront=double(md.gridonboundary & md.gridoniceshelf);
+        gridoniceshelf=zeros(md.numberofgrids,1);
+        gridoniceshelf(md.elements(find(md.elementoniceshelf),:))=1;
+        gridonicefront=double(md.gridonboundary & gridoniceshelf);
 end
 pos=find(md.gridonboundary & ~gridonicefront);