Changeset 440

Timestamp:

05/15/09 08:42:32 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed diagnostic if setelements type done before extrusion

Location:

issm/trunk/src/m/classes/public

Files:

• extrude.m (modified) (1 diff)
• ismodelselfconsistent.m (modified) (1 diff)

issm/trunk/src/m/classes/public/extrude.m

@@ -134 +134 @@
         md.gridonpattyn=project3d(md,md.gridonpattyn,'node');
         md.gridonstokes=project3d(md,md.gridonstokes,'node');
+
+        %dead grids
+        md.deadgrids=ones(md.numberofgrids,1);
+        md.deadgrids(md.elements(md.elements_type(:,1)~=macayealenum,:))=0;%non macayeal grids are not dead
+        md.deadgrids(find(md.gridonbed))=0;%grids from elements on bed are not dead
 end
 

issm/trunk/src/m/classes/public/ismodelselfconsistent.m

@@ -115 +115 @@
 
 %SIZE NUMBEROFGRIDS
-fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface'};
+fields={'x','y','z','B','drag','gridondirichlet_diag','melting','accumulation','surface','thickness','bed','gridonbed','gridonsurface','deadgrids'};
 for i=1:length(fields),
         if find(length(eval(['md.' char(fields(i))]))~=md.numberofgrids),