Index: /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py	(revision 21353)
+++ /issm/trunk-jpl/src/m/classes/mismipbasalforcings.py	(revision 21354)
@@ -6,5 +6,5 @@
 
 class mismipbasalforcings(object):
-    """
+    """ 
     MISMIP Basal Forcings class definition
 
@@ -13,5 +13,5 @@
     """
 
-    def __init__(self,md): # {{{
+    def __init__(self): # {{{
 
         self.groundedice_melting_rate = float('NaN')
@@ -20,8 +20,4 @@
         self.upperdepth_melt = float('NaN')
         self.geothermalflux = float('NaN')
-
-        if np.all(np.isnan(self.groundedice_melting_rate)):
-            self.groundedice_melting_rate=np.zeros(md.mesh.numberofvertices)
-            print ' no basalforcings.groundedice_melting_rate specified: values set as zero'
 
 	self.setdefaultparameters()
@@ -41,4 +37,10 @@
         self.geothermalflux=project3d(md,'vector',self.geothermalflux,'type','node','layer',1)    #bedrock only gets geothermal flux
 	return self
+    #}}}
+    def initialize(self,md): # {{{
+        if np.all(np.isnan(self.groundedice_melting_rate)):
+            self.groundedice_melting_rate=np.zeros(md.mesh.numberofvertices)
+            print ' no basalforcings.groundedice_melting_rate specified: values set as zero'
+        return self
     #}}}
     def setdefaultparameters(self): # {{{