Index: /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.py	(revision 26903)
+++ /issm/trunk-jpl/src/m/classes/SMBpddSicopolis.py	(revision 26904)
@@ -129,5 +129,5 @@
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 's0t', 'format', 'DoubleMat', 'mattype', 1)
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'rlaps', 'format', 'Double')
-        WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'monthlytemperatures', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', md.constants.yts)
+        WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'monthlytemperatures', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', yts)
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'precipitation', 'format', 'DoubleMat', 'mattype', 1, 'scale', 1. / yts, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', yts)
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'temperature_anomaly', 'format', 'DoubleMat', 'mattype', 1, 'timeserieslength', md.mesh.numberofvertices + 1, 'yts', yts)
Index: /issm/trunk-jpl/src/m/classes/m1qn3inversion.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/m1qn3inversion.py	(revision 26903)
+++ /issm/trunk-jpl/src/m/classes/m1qn3inversion.py	(revision 26904)
@@ -1,4 +1,3 @@
 import numpy as np
-
 from checkfield import checkfield
 from fielddisplay import fielddisplay
Index: /issm/trunk-jpl/src/m/classes/massfluxatgate.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/massfluxatgate.py	(revision 26903)
+++ /issm/trunk-jpl/src/m/classes/massfluxatgate.py	(revision 26904)
@@ -1,2 +1,3 @@
+import numpy as np
 from fielddisplay import fielddisplay
 from pairoptions import pairoptions
@@ -20,5 +21,5 @@
         self.definitionstring = ''
         self.profilename = ''
-        self.segments = float('NaN')
+        self.segments = np.nan
 
     #set defaults
Index: /issm/trunk-jpl/src/m/classes/regionaloutput.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/regionaloutput.py	(revision 26903)
+++ /issm/trunk-jpl/src/m/classes/regionaloutput.py	(revision 26904)
@@ -25,7 +25,8 @@
 
         self.name = ''
+        self.model= ''
         self.definitionstring = ''
         self.outputnamestring = ''
-        self.mask = float('NaN')
+        self.mask = np.nan
         self.maskexpstring = ''
 
@@ -46,6 +47,6 @@
                 self.setmaskfromexp(modelname)
 
-        if (len(self.mask) <= 1 & np.any(np.isnan(self.mask))):
-            error('regionaloutput error message: ''mask'' field or ''maskexpstring'' and ''model'' fields should be defined!')
+        # if (len(self.mask) <= 1 & np.any(np.isnan(self.mask))):
+        #     raise IOError('regionaloutput error message: ''mask'' field or ''maskexpstring'' and ''model'' fields should be defined!')
 
     #}}}
@@ -76,4 +77,5 @@
         return self
     # }}}
+
     def checkconsistency(self, md, solution, analyses):  # {{{