Index: /issm/trunk/src/m/model/ismodelselfconsistent.m
===================================================================
--- /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 5918)
+++ /issm/trunk/src/m/model/ismodelselfconsistent.m	(revision 5919)
@@ -71,5 +71,5 @@
 checksize(md,fields,[md.numberofgrids 1]);
 %Check the values of vertices_type
-checkvalues(md,{'vertices_type'},[MacAyealApproximationEnum() HutterApproximationEnum() PattynApproximationEnum() MacAyealPattynApproximationEnum() StokesApproximationEnum() NoneApproximationEnum()]);
+checkvalues(md,{'vertices_type'},[MacAyealApproximationEnum() HutterApproximationEnum() PattynApproximationEnum() MacAyealPattynApproximationEnum() StokesApproximationEnum() PattynStokesApproximationEnum() NoneApproximationEnum()]);
 if (md.dim==2),
 	checkvalues(md,{'vertices_type'},[MacAyealApproximationEnum() HutterApproximationEnum()]);
Index: /issm/trunk/src/m/model/setelementstype.m
===================================================================
--- /issm/trunk/src/m/model/setelementstype.m	(revision 5918)
+++ /issm/trunk/src/m/model/setelementstype.m	(revision 5919)
@@ -195,6 +195,4 @@
 pos=find(gridonmacayealpattyn);
 md.vertices_type(pos)=MacAyealPattynApproximationEnum();
-pos=find(gridonpattynstokes);
-md.vertices_type(pos)=PattynStokesApproximationEnum();
 pos=find(gridonstokes);
 md.vertices_type(pos)=StokesApproximationEnum();
@@ -203,4 +201,6 @@
 	md.vertices_type(pos)=NoneApproximationEnum();
 end
+pos=find(gridonpattynstokes);
+md.vertices_type(pos)=PattynStokesApproximationEnum();
 
 %figure out solution types