Index: /issm/trunk-jpl/src/m/classes/model.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/model.py	(revision 26572)
+++ /issm/trunk-jpl/src/m/classes/model.py	(revision 26573)
@@ -80,5 +80,5 @@
     """
 
-    def __init__(self, *args): #{{{
+    def __init__(self, *args):  #{{{
         self.mesh = None
         self.mask = None
@@ -189,5 +189,5 @@
     #}}}
 
-    def properties(self): #{{{
+    def properties(self):  #{{{
         # ordered list of properties since vars(self) is random
         return [
@@ -239,5 +239,5 @@
     #}}}
 
-    def setdefaultparameters(self, planet): #{{{
+    def setdefaultparameters(self, planet):  #{{{
         self.mesh = mesh2d()
         self.mask = mask()
@@ -286,5 +286,5 @@
     #}}}
 
-    def checkmessage(self, string): #{{{
+    def checkmessage(self, string):  #{{{
         print("model not consistent: {}".format(string))
         self.private.isconsistent = False
@@ -293,5 +293,5 @@
     #@staticmethod
 
-    def extract(self, area): #{{{
+    def extract(self, area):  #{{{
         """EXTRACT - extract a model according to an Argus contour or flag list
 
@@ -371,5 +371,5 @@
             field = getattr(md1, fieldi)
             fieldsize = np.shape(field)
-            if hasattr(field, '__dict__') and fieldi not in ['results']: #recursive call
+            if hasattr(field, '__dict__') and fieldi not in ['results']:  #recursive call
                 object_fields = vars(field)
                 for fieldj in object_fields:
@@ -572,5 +572,5 @@
     #}}}
 
-    def extrude(md, *args): #{{{
+    def extrude(md, *args):  #{{{
         """EXTRUDE - vertically extrude a 2d mesh
 
@@ -759,5 +759,5 @@
     #}}}
 
-    def collapse(md): #{{{
+    def collapse(md):  #{{{
         """COLLAPSE - collapses a 3d mesh into a 2d mesh
 
@@ -976,5 +976,5 @@
             mesh.scale_factor = project2d(md, md.mesh.scale_factor, 1)
         if hasattr(md.mesh, 'vertexonboundary') and not np.isnan(md.mesh.vertexonboundary).all():
-            mesh.vertexonboundary= project2d(md, md.mesh.vertexonboundary, 1)
+            mesh.vertexonboundary = project2d(md, md.mesh.vertexonboundary, 1)
         if hasattr(md.mesh, 'elementonboundary') and not np.isnan(md.mesh.elementonboundary).all():
             mesh.elementonboundary = project2d(md, md.mesh.elementonboundary, 1)