Index: /issm/trunk-jpl/src/m/mesh/ComputeMetric.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/ComputeMetric.py	(revision 14094)
+++ /issm/trunk-jpl/src/m/mesh/ComputeMetric.py	(revision 14095)
@@ -54,5 +54,5 @@
 	#take care of NaNs if any (use Numpy eig in a loop)
 	pos=numpy.nonzero(numpy.isnan(metric))[0]
-	if pos:
+	if numpy.size(pos):
 		print(" %i NaN found in the metric. Use Numpy routine..." % numpy.size(pos))
 		for posi in pos:
Index: /issm/trunk-jpl/src/m/mesh/squaremesh.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/squaremesh.py	(revision 14094)
+++ /issm/trunk-jpl/src/m/mesh/squaremesh.py	(revision 14095)
@@ -22,6 +22,6 @@
 	#initialization
 	index=numpy.zeros((nel,3),int)
-	x=numpy.zeros((nx*ny,1))
-	y=numpy.zeros((nx*ny,1))
+	x=numpy.zeros((nx*ny))
+	y=numpy.zeros((nx*ny))
 
 	#create coordinates
@@ -59,10 +59,10 @@
 	md.mesh.x=x
 	md.mesh.y=y
-	md.mesh.z=numpy.zeros((nods,1))
+	md.mesh.z=numpy.zeros((nods))
 	md.mesh.numberofvertices=nods
-	md.mesh.vertexonboundary=numpy.zeros((nods,1),bool)
+	md.mesh.vertexonboundary=numpy.zeros((nods),bool)
 	md.mesh.vertexonboundary[segments[:,0:2]-1]=True
-	md.mesh.vertexonbed=numpy.ones((nods,1),bool)
-	md.mesh.vertexonsurface=numpy.ones((nods,1),bool)
+	md.mesh.vertexonbed=numpy.ones((nods),bool)
+	md.mesh.vertexonsurface=numpy.ones((nods),bool)
 
 	#plug elements
@@ -70,6 +70,6 @@
 	md.mesh.segments=segments
 	md.mesh.numberofelements=nel
-	md.mesh.elementonbed=numpy.ones((nel,1),bool)
-	md.mesh.elementonsurface=numpy.ones((nel,1),bool)
+	md.mesh.elementonbed=numpy.ones((nel),bool)
+	md.mesh.elementonsurface=numpy.ones((nel),bool)
 
 	#Now, build the connectivity tables for this mesh.
Index: /issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py
===================================================================
--- /issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py	(revision 14094)
+++ /issm/trunk-jpl/src/m/miscellaneous/MatlabFuncs.py	(revision 14095)
@@ -8,4 +8,12 @@
 
 	if 'Windows' in platform.system():
+		return True
+	else:
+		return False
+
+def ismac():
+	import platform
+
+	if 'Darwin' in platform.system():
 		return True
 	else: