Index: /issm/trunk/src/m/classes/public/modelextract2.m
===================================================================
--- /issm/trunk/src/m/classes/public/modelextract2.m	(revision 1399)
+++ /issm/trunk/src/m/classes/public/modelextract2.m	(revision 1400)
@@ -212,4 +212,14 @@
 	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)