Changeset 432 in lmdz_wrf

Timestamp:

May 26, 2015, 2:53:12 PM (10 years ago)

Author:

lfita

Message:

Fixing the same in a different way...

File:

• trunk/tools/nc_var_tools.py (modified) (3 diffs)

trunk/tools/nc_var_tools.py

@@ -11156 +11156 @@
       'squared mean value within', 'standard deviation value within',                \
       'accumulated value within']
-    Nsts = len(stsn)+1
+    Nsts = len(stsn)
     statnames = []
     cstatnames = []
@@ -11225 +11225 @@
             rlatvals = np.ones(tuple([Ttraj,dimz,Nrad*2+1,Nrad*2+1]), dtype=np.float)
 
-            statvarvals = np.ones(tuple([Ttraj,dimz,Nsts]), dtype=np.float)
-            rstatvarvals = np.ones(tuple([Ttraj,dimz,Nsts]), dtype=np.float)
+# One dimension plus for the values at the center of the trajectory
+            statvarvals = np.ones(tuple([Ttraj,dimz,Nsts+1]), dtype=np.float)
+            rstatvarvals = np.ones(tuple([Ttraj,dimz,Nsts+1]), dtype=np.float)
 
             for it in range(Ttraj):
@@ -11477 +11478 @@
             rlatvals = np.ones(tuple([Ttraj,Nrad*2+1,Nrad*2+1]), dtype=np.float)
 
-            statvarvals = np.ones(tuple([Ttraj,Nsts]), dtype=np.float)
-            rstatvarvals = np.ones(tuple([Ttraj,Nsts]), dtype=np.float)
+# One dimension plus for the values at the center of the trajectory
+            statvarvals = np.ones(tuple([Ttraj,Nsts+1]), dtype=np.float)
+            rstatvarvals = np.ones(tuple([Ttraj,Nsts+1]), dtype=np.float)
 
             for it in range(Ttraj):