Index: /issm/trunk-jpl/src/m/mesh/roundmesh.m
===================================================================
--- /issm/trunk-jpl/src/m/mesh/roundmesh.m	(revision 27842)
+++ /issm/trunk-jpl/src/m/mesh/roundmesh.m	(revision 27843)
@@ -11,5 +11,5 @@
 
 %First we have to create the domain outline 
-if nargin>=4:
+if nargin>=4
 	expname = varargin{1};
 else
@@ -27,10 +27,10 @@
 expwrite(A,expname);
 
-%Call Bamg
+%Call mesher
 md=triangle(md,expname,resolution);
 %md=bamg(md,'domain','RoundDomainOutline.exp','hmin',resolution);
 
 %move the closest node to the center
-[mini pos]=min(md.mesh.x.^2+md.mesh.y.^2);
+[minimum pos]=min(md.mesh.x.^2+md.mesh.y.^2);
 md.mesh.x(pos)=0.;
 md.mesh.y(pos)=0.;
Index: /issm/trunk-jpl/src/m/mesh/roundmesh.py
===================================================================
--- /issm/trunk-jpl/src/m/mesh/roundmesh.py	(revision 27842)
+++ /issm/trunk-jpl/src/m/mesh/roundmesh.py	(revision 27843)
@@ -37,5 +37,5 @@
     expwrite(A, expname)
 
-    # Call Bamg
+    # Call mesher
     md = triangle(md, expname, resolution)
     # md = bamg(md, 'domain', 'RoundDomainOutline.exp', 'hmin', resolution)
@@ -47,5 +47,6 @@
 
     # Delete domain
-    os.remove(expname)
+    if not len(args):
+        os.remove(expname)
 
     return md