Changeset 1400

Timestamp:

07/28/09 16:06:09 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed modelextract2 for Stokes (Yessa)

File:

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

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

@@ -212 +212 @@
         end
 
+        %Recreate border stokes
+        if md1.isstokes,
+                pos=find(md2.gridondirichlet_diag);                 %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)==StokesEnum()); %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.segmentonneumann_diag)