Index: /issm/trunk-jpl/src/m/parameterization/setflowequation.m
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setflowequation.m	(revision 19534)
+++ /issm/trunk-jpl/src/m/parameterization/setflowequation.m	(revision 19535)
@@ -258,6 +258,4 @@
 %Create vertices_type
 md.flowequation.vertex_equation=zeros(md.mesh.numberofvertices,1);
-pos=find(nodeonSIA);
-md.flowequation.vertex_equation(pos)=1;
 pos=find(nodeonSSA);
 md.flowequation.vertex_equation(pos)=2;
@@ -268,4 +266,7 @@
 pos=find(nodeonFS);
 md.flowequation.vertex_equation(pos)=5;
+%DO SIA LAST! Otherwise spcs might not be set up correctly (SIA should have priority)
+pos=find(nodeonSIA);
+md.flowequation.vertex_equation(pos)=1;
 if any(FSflag),
 	pos=find(~nodeonFS);
Index: /issm/trunk-jpl/src/m/parameterization/setflowequation.py
===================================================================
--- /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 19534)
+++ /issm/trunk-jpl/src/m/parameterization/setflowequation.py	(revision 19535)
@@ -246,6 +246,4 @@
 	#Create vertices_type
 	md.flowequation.vertex_equation=numpy.zeros(md.mesh.numberofvertices,int)
-	pos=numpy.nonzero(nodeonSIA)
-	md.flowequation.vertex_equation[pos]=1
 	pos=numpy.nonzero(nodeonSSA)
 	md.flowequation.vertex_equation[pos]=2
@@ -256,4 +254,7 @@
 	pos=numpy.nonzero(nodeonFS)
 	md.flowequation.vertex_equation[pos]=5
+	#DO SIA LAST! Otherwise spcs might not be set up correctly (SIA should have priority)
+	pos=numpy.nonzero(nodeonSIA)
+	md.flowequation.vertex_equation[pos]=1
 	if any(FSflag):
 		pos=numpy.nonzero(numpy.logical_not(nodeonFS))