Changeset 17696

Timestamp:

04/09/14 16:33:21 (11 years ago)

Author:

cborstad

Message:

BUG: create copy so that original data is not overwritten by processdata

File:

• issm/trunk-jpl/src/m/plot/processdata.py (modified) (4 diffs)

issm/trunk-jpl/src/m/plot/processdata.py ¶

@@ -29 +29 @@
                 numberofelements2d=npy.nan
 
+        procdata=npy.copy(data)
+
         #process patch
 
@@ -35 +37 @@
 
         #get datasize
-        if npy.ndim(data)==1:
-                datasize=npy.array([len(data),1])
+        if npy.ndim(procdata)==1:
+                datasize=npy.array([len(procdata),1])
         else:
-                datasize=npy.shape(data)
+                datasize=npy.shape(procdata)
         
         #process NaN's if any
-        nanfill=options.getfieldvalue('nan',-999)
-        if npy.any(npy.isnan(data)):
-                data[npy.isnan(data)]=nanfill
-                print "WARNING: nan's treated as -999 by default.  Change using pairoption 'nan',nan_fill_value in plotmodel call"
+        nanfill=options.getfieldvalue('nan',-9999)
+        if npy.any(npy.isnan(procdata)):
+                procdata[npy.isnan(procdata)]=nanfill
+                print "WARNING: nan's treated as -9999 by default.  Change using pairoption 'nan',nan_fill_value in plotmodel call"
 
         #non-patch processing 
@@ -72 +74 @@
 
         #convert rank-2 array to rank-1
-        if npy.ndim(data)==2 and npy.shape(data)[1]==1:
-                data=data.reshape(-1,)
+        if npy.ndim(procdata)==2 and npy.shape(procdata)[1]==1:
+                procdata=procdata.reshape(-1,)
 
         #if datatype is still zero, error out
@@ -79 +81 @@
                 raise ValueError("processdata error: data provided not recognized or not supported")
         else:
-                return data, datatype
+                return procdata, datatype