Index: /issm/trunk-jpl/src/m/classes/SMBgemb.py
===================================================================
--- /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 25310)
+++ /issm/trunk-jpl/src/m/classes/SMBgemb.py	(revision 25311)
@@ -377,5 +377,8 @@
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'isconstrainsurfaceT', 'format', 'Boolean')
 
+        print(md.smb.Ta[-1])
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'Ta', 'format', 'DoubleMat', 'mattype', 2, 'timeserieslength', md.mesh.numberofelements + 1, 'yts', yts)
+        print(md.smb.Ta[-1])
+        exit()
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'V', 'format', 'DoubleMat', 'mattype', 2, 'timeserieslength', md.mesh.numberofelements + 1, 'yts', yts)
         WriteData(fid, prefix, 'object', self, 'class', 'smb', 'fieldname', 'dswrf', 'format', 'DoubleMat', 'mattype', 2, 'timeserieslength', md.mesh.numberofelements + 1, 'yts', yts)
Index: /issm/trunk-jpl/src/m/solve/WriteData.py
===================================================================
--- /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 25310)
+++ /issm/trunk-jpl/src/m/solve/WriteData.py	(revision 25311)
@@ -53,6 +53,6 @@
     #Scale data if necessary
     if datatype == 'MatArray':
-        #if it is a matrix array we loop over the matrixes
-        for i in range(np.size(data)):
+        #if it is a matrix array we loop over the matrices
+        for i in range(len(data)):
             if options.exist('scale'):
                 scale = options.getfieldvalue('scale')
@@ -65,5 +65,5 @@
                 #no scaling given, use yts instead
                 yts = options.getfieldvalue('yts')
-                #We scale everything but the last line that holds time
+                #We scale only the last line that holds time
                 data[i][-1, :] = yts * data[i][-1, :]
     else: