Index: /issm/trunk-jpl/src/m/classes/SMBgemb.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 26358)
+++ /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 26359)
@@ -22,13 +22,13 @@
 
         #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
 
Index: /issm/trunk-jpl/src/m/classes/dsl.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/dsl.py	(revision 26358)
+++ /issm/trunk-jpl/src/m/classes/dsl.py	(revision 26359)
@@ -68,6 +68,4 @@
 
     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)
Index: /issm/trunk-jpl/src/m/classes/matdamageice.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/matdamageice.py	(revision 26358)
+++ /issm/trunk-jpl/src/m/classes/matdamageice.py	(revision 26359)
@@ -1,2 +1,4 @@
+import numpy as np
+
 from fielddisplay import fielddisplay
 from project3d import project3d