Changeset 24732

Timestamp:

04/23/20 11:29:16 (5 years ago)

Author:

Mathieu Morlighem

Message:

CHG: time should be in years

File:

• issm/trunk-jpl/scripts/OutbinRead.py (modified) (4 diffs)

issm/trunk-jpl/scripts/OutbinRead.py

@@ -7 +7 @@
 
 def OutbinRead(filin, filout='', verbose=0):  #{{{
+
+    yts  = 365*24*3600.
 
     print("reading binary file.")
@@ -31 +33 @@
         time = struct.unpack('d', f.read(struct.calcsize('d')))[0]
         step = struct.unpack('i', f.read(struct.calcsize('i')))[0]
+        if time!=-9999.: time = time/yts
         print("time = {} step = {}".format(time,step))
 
@@ -40 +43 @@
         M = struct.unpack('i', f.read(struct.calcsize('i')))[0]
         N = 1 #default
-        print("M = {}".format(M))
 
         #Step 2: read the data itself.
@@ -72 +74 @@
             s[1] = struct.unpack('i', f.read(struct.calcsize('i')))[0]
             print("size = [{}x{}]".format(s[0], s[1]))
-            data = np.zeros((s[0], s[1]))
-            for i in range(s[0]):
-                for j in range(s[1]):
-                    data[i][j] = struct.unpack('d', f.read(struct.calcsize('d')))[0]    #get to the "c" convention, hence the transpose
-                    if verbose > 2:
-                        print("data[{}, {}] = {}".format(i, j, data[i][j]))
+            if True:
+               f.read(struct.calcsize('d')*s[0]*s[1])
+            else:
+                                        data = np.zeros((s[0], s[1]))
+                                        for i in range(s[0]):
+                                                 for j in range(s[1]):
+                                                          data[i][j] = struct.unpack('d', f.read(struct.calcsize('d')))[0]    #get to the "c" convention, hence the transpose
+                                                          if verbose > 2:
+                                                                        print("data[{}, {}] = {}".format(i, j, data[i][j]))
 
         else: