Index: ../trunk-jpl/src/m/classes/SMBgemb.py
===================================================================
--- ../trunk-jpl/src/m/classes/SMBgemb.py	(revision 26358)
+++ ../trunk-jpl/src/m/classes/SMBgemb.py	(revision 26359)
@@ -21,15 +21,15 @@
         """
 
         #solution choices
-        self.isgraingrowth
-        self.isalbedo
-        self.isshortwave
-        self.isthermal
-        self.isaccumulation
-        self.ismelt
-        self.isdensification
-        self.isturbulentflux
-        self.isconstrainsurfaceT
+        self.isgraingrowth = 0
+        self.isalbedo = 0
+        self.isshortwave = 0
+        self.isthermal = 0
+        self.isaccumulation = 0
+        self.ismelt = 0
+        self.isdensification = 0
+        self.isturbulentflux = 0
+        self.isconstrainsurfaceT = 0
         self.isclimatology = np.nan
 
         # Inputs
Index: ../trunk-jpl/src/m/classes/dsl.py
===================================================================
--- ../trunk-jpl/src/m/classes/dsl.py	(revision 26358)
+++ ../trunk-jpl/src/m/classes/dsl.py	(revision 26359)
@@ -67,8 +67,6 @@
     #}}}
 
     def initialize(self, md):  #{{{
-        print('sdasa')
-        exit()
         if np.isnan(self.global_average_thermosteric_sea_level):
             self.global_average_thermosteric_sea_level = np.array([0, 0]).reshape(-1, 1)
             print('      no dsl.global_average_thermosteric_sea_level specified: transient values set at zero')
Index: ../trunk-jpl/src/m/classes/matdamageice.py
===================================================================
--- ../trunk-jpl/src/m/classes/matdamageice.py	(revision 26358)
+++ ../trunk-jpl/src/m/classes/matdamageice.py	(revision 26359)
@@ -1,3 +1,5 @@
+import numpy as np
+
 from fielddisplay import fielddisplay
 from project3d import project3d
 from checkfield import checkfield