Index: /issm/trunk-jpl/src/m/mesh/meshconvert.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/meshconvert.py	(revision 13278)
+++ /issm/trunk-jpl/src/m/mesh/meshconvert.py	(revision 13279)
@@ -33,10 +33,10 @@
 	md.private.bamg['mesh']=bamgmesh(bamgmesh_out)
 	md.private.bamg['geometry']=bamggeom(bamggeom_out)
-	md.mesh.x=bamgmesh_out['Vertices'][:,0]
-	md.mesh.y=bamgmesh_out['Vertices'][:,1]
-	md.mesh.elements=bamgmesh_out['Triangles'][:,0:3]
-	md.mesh.edges=bamgmesh_out['IssmEdges']
-	md.mesh.segments=bamgmesh_out['IssmSegments'][:,0:3]
-	md.mesh.segmentmarkers=bamgmesh_out['IssmSegments'][:,3]
+	md.mesh.x=bamgmesh_out['Vertices'][:,0].copy()
+	md.mesh.y=bamgmesh_out['Vertices'][:,1].copy()
+	md.mesh.elements=bamgmesh_out['Triangles'][:,0:3].copy()
+	md.mesh.edges=bamgmesh_out['IssmEdges'].copy()
+	md.mesh.segments=bamgmesh_out['IssmSegments'][:,0:3].copy()
+	md.mesh.segmentmarkers=bamgmesh_out['IssmSegments'][:,3].copy()
 
 	#Fill in rest of fields:
Index: /issm/trunk-jpl/test/NightlyRun/test106.py
===================================================================
--- /issm/trunk-jpl/test/NightlyRun/test106.py	(revision 13279)
+++ /issm/trunk-jpl/test/NightlyRun/test106.py	(revision 13279)
@@ -0,0 +1,26 @@
+from model import *
+from triangle import *
+from meshconvert import *
+from setmask import *
+from parameterize import *
+from setflowequation import *
+from EnumDefinitions import *
+from solve import *
+from MatlabFuncs import *
+
+md=triangle(model(),'../Exp/Square.exp',150000)
+md=meshconvert(md)
+md=setmask(md,'all','')
+md=parameterize(md,'../Par/SquareShelfConstrained.py')
+md=setflowequation(md,'macayeal','all')
+md.cluster=generic('name',oshostname(),'np',3)
+md.prognostic.stabilization=3
+md.prognostic.spcthickness=md.geometry.thickness
+md=solve(md,PrognosticSolutionEnum())
+
+#Fields and tolerances to track changes
+field_names     =['Thickness']
+field_tolerances=[1e-13]
+field_values=[\
+	md.results['PrognosticSolution'][1]['Thickness'],\
+	]
