Changeset 1370

Timestamp:

07/24/09 13:56:13 (16 years ago)

Author:

Mathieu Morlighem

Message:

fixed extraction 2 for Stokes

File:

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

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

@@ -223 +223 @@
                 md2.segmentonneumann_diag=md2.segmentonneumann_diag(find(md2.segmentonneumann_diag(:,1) & md2.segmentonneumann_diag(:,2)),:);
         end
+        if ~isnan(md2.segmentonneumann_diag_stokes)
+                md2.segmentonneumann_diag_stokes(:,1)=Pgrid(md1.segmentonneumann_diag_stokes(:,1)); 
+                md2.segmentonneumann_diag_stokes(:,2)=Pgrid(md1.segmentonneumann_diag_stokes(:,2)); 
+                md2.segmentonneumann_diag_stokes(:,end)=Pelem(md1.segmentonneumann_diag_stokes(:,end)); 
+                if strcmpi(md1.type,'3d')
+                        md2.segmentonneumann_diag_stokes(:,3)=Pgrid(md1.segmentonneumann_diag_stokes(:,3)); 
+                        md2.segmentonneumann_diag_stokes(:,4)=Pgrid(md1.segmentonneumann_diag_stokes(:,4)); 
+                end
+                md2.segmentonneumann_diag_stokes=md2.segmentonneumann_diag_stokes(find(md2.segmentonneumann_diag_stokes(:,1) & md2.segmentonneumann_diag_stokes(:,2)),:);
+        end
         md2.neumannvalues_diag=NaN;