Index: /issm/trunk-jpl/test/NightlyRun/python_skipped_tests.txt
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/python_skipped_tests.txt	(revision 13683)
+++ /issm/trunk-jpl/test/NightlyRun/python_skipped_tests.txt	(revision 13684)
@@ -1,2 +1,5 @@
+%Get full list with:
+%  for i in test*.m; do NAME=$(echo $i | sed -e "s/\.m//"); %if [ ! -f $NAME.py ]; then echo $NAME; fi; done
+
 test216    needs TriMeshProcessRifts module in meshprocessrifts.py (skip, per Eric, 10/03/12)
 test218    needs Dakota
@@ -4,2 +7,13 @@
 test235    needs Dakota
 test274    needs TriMeshProcessRifts module in meshprocessrifts.py
+test328    John is working on it
+test329    John is working on it
+test412    needs Dakota
+test413    needs Dakota
+test414    needs Dakota
+test417    needs Dakota
+test418    needs Dakota
+test420    needs Dakota
+test423    needs roundmesh (and not straightforward)
+test511    needs extract
+test613    needs extract
Index: /issm/trunk-jpl/test/NightlyRun/test516.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test516.py	(revision 13684)
+++ /issm/trunk-jpl/test/NightlyRun/test516.py	(revision 13684)
@@ -0,0 +1,28 @@
+import numpy
+from model import *
+from EnumDefinitions import *
+from MatlabFuncs import *
+from triangle import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from solve import *
+
+md=triangle(model(),'../Exp/Pig.exp',30000.)
+md=setmask(md,'../Exp/PigShelves.exp','../Exp/PigIslands.exp')
+md=parameterize(md,'../Par/Pig.py')
+md.extrude(3,1.)
+md=setflowequation(md,'pattyn','all')
+md.thermal.stabilization=2
+md.cluster=generic('name',oshostname(),'np',3)
+md.timestepping.time_step=0
+md.thermal.penalty_threshold=40
+md=solve(md,ThermalSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Temperature','BasalforcingsMeltingRate']
+field_tolerances=[1e-11,1e-11]
+field_values=[\
+	md.results['ThermalSolution'][1]['Temperature'],\
+	md.results['ThermalSolution'][1]['BasalforcingsMeltingRate'],\
+	]
