Changeset 17687 for issm/trunk-jpl/src/m/boundaryconditions/SetMarineIceSheetBC.m

Timestamp:

04/08/14 21:13:16 (11 years ago)

Author:

Mathieu Morlighem

Message:

CHG: replacing some meshtype by meshdimension

File:

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

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

@@ -52 +52 @@
 
 %First find segments that are not completely on the front
-if ~strcmp(md.mesh.domaintype(),'3D'),
+if strcmp(elementtype(md.mesh),'Penta'),
+        numbernodesfront=4;
+elseif strcmp(elementtype(md.mesh),'Tria'),
         numbernodesfront=2;
-else 
-        numbernodesfront=4;
+else
+        error('mesh type not supported yet');
 end
 segmentsfront=md.mask.ice_levelset(md.mesh.segments(:,1:numbernodesfront))==0;