Changeset 5492

Timestamp:

08/23/10 08:29:43 (15 years ago)

Author:

Mathieu Morlighem

Message:

Fixed modelextract for Stokes

File:

• issm/trunk/src/m/classes/public/modelextract.m (modified) (1 diff)

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

@@ -219 +219 @@
         end
 
-        %Recreate border stokes
-        if md1.isstokes,
-                pos=find(sum(md2.spcvelocity(:,1:3),2));                 %find all the grids on the boundary of the domain without icefront
-                md2.gridonstokes(pos)=0;                               %we have to constrain all the boundary of the domain without icefront
-                stokes_elements=find(md2.elements_type(:,2)==StokesApproximationEnum()); %find the elements on the stokes domain
-                borderflags=zeros(numberofgrids2,1); 
-                borderflags(md2.elements(stokes_elements,:))=1; %find all the grids of the elements on stokes domain, ie stokes grids and borderstokes
-                md2.borderstokes=borderflags-md2.gridonstokes;   %remove stokes grids from this list
-        end
-
         %Diagnostic
         if ~isnan(md2.pressureload)