Index: ../trunk-jpl/src/m/classes/masstransport.py
===================================================================
--- ../trunk-jpl/src/m/classes/masstransport.py	(revision 24286)
+++ ../trunk-jpl/src/m/classes/masstransport.py	(revision 24287)
@@ -73,6 +73,7 @@
         md = checkfield(md, 'fieldname', 'masstransport.stabilization', 'values', [0, 1, 2, 3, 4])
         md = checkfield(md, 'fieldname', 'masstransport.min_thickness', '>', 0)
         md = checkfield(md, 'fieldname', 'masstransport.requested_outputs', 'stringrow', 1)
+        md = checkfield(md, 'fieldname', 'masstransport.vertex_pairing', '>', 0)
 
         return md
     # }}}
Index: ../trunk-jpl/src/m/classes/stressbalance.py
===================================================================
--- ../trunk-jpl/src/m/classes/stressbalance.py	(revision 24286)
+++ ../trunk-jpl/src/m/classes/stressbalance.py	(revision 24287)
@@ -133,7 +133,7 @@
         md = checkfield(md, 'fieldname', 'stressbalance.referential', 'size', [md.mesh.numberofvertices, 6])
         md = checkfield(md, 'fieldname', 'stressbalance.loadingforce', 'size', [md.mesh.numberofvertices, 3])
         md = checkfield(md, 'fieldname', 'stressbalance.requested_outputs', 'stringrow', 1)
-
+        md = checkfield(md, 'fieldname', 'stressbalance.vertex_pairing', '>', 0)
         #singular solution
         #        if ~any((~isnan(md.stressbalance.spcvx) + ~isnan(md.stressbalance.spcvy)) == 2),
         if not np.any(np.logical_and(np.logical_not(np.isnan(md.stressbalance.spcvx)), np.logical_not(np.isnan(md.stressbalance.spcvy)))):