Changeset 25311

Timestamp:

07/29/20 11:23:45 (5 years ago)

Author:

jdquinn

Message:

BUG: WriteData matrix addressing (in Python, it's a list of NumPy ndarray)

Location:

issm/trunk-jpl/src/m

Files:

• classes/SMBgemb.py (modified) (1 diff)
• solve/WriteData.py (modified) (2 diffs)

issm/trunk-jpl/src/m/classes/SMBgemb.py

@@ -377 +377 @@
         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)

issm/trunk-jpl/src/m/solve/WriteData.py

@@ -53 +53 @@
     #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 +65 @@
                 #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: