Changeset 916 in lmdz_wrf for trunk

Timestamp:

Jun 20, 2016, 8:54:53 PM (8 years ago)

Author:

lfita

Message:

Adding `-9' for the last value in 'SliceVarDict?'

File:

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

trunk/tools/nc_var_tools.py

@@ -8964 +8964 @@
             * [integer]: which value of the dimension
             * -1: all along the dimension
+            * -9: last value along the dimension
             * [beg]:[end] slice from [beg] to [end]
             NOTE: that dimension of the variables which not appear al values will be used
@@ -9126 +9127 @@
             objnewnc.sync()
     
+    # Statistics of the variable dimension along the dimensions operated
+    basicoper = ['min', 'max', 'menan', 'mean2', 'sum']
+
 # operating variable dimension
     print '  ' + fname + ': Opering slice of variable dimension _______'
@@ -9221 +9225 @@
             objnewnc.sync()
 
+            # Statistics of the variable dimension along the dimensions operated
+            for op in basicoper:
+                if op != operkind:
+                    finalvarv, finalds = operdim(varvals,varsliceddims,dimsoper,op)
+                    if varinf.FillValue is not None:
+                        newvar = objnewnc.createVariable(vn+op, nctype(varinf.dtype),\
+                          tuple(finalds), fillValue=varinf.FillValue)
+                    else:
+                        newvar = objnewnc.createVariable(vn+op, nctype(varinf.dtype),\
+                          tuple(finalds))
+                    newNvar = newNvar + 1
+                    newvar[:] = finalvarv[:]
+                    newattr = basicvardef(newvar, stdname + op, lname + ' '+ op +    \
+                      ' along ' + dimsoperS, uname)
+            objnewnc.sync()
+
 # Including operation dimensions on the new slice for the variables dimension
 #   but, which value along the operated dimension should be taken?
@@ -9267 +9287 @@
                 attrv = ov.getncattr(attrn)
                 newattr = set_attribute(newvar, attrn, attrv)
-#            LLUIS
             objnewnc.sync()
-
 
 # global attributes
@@ -13344 +13362 @@
         slicedict[dimn] = [value]
           [value] = -1, all the values
+          [value] = -9, last values
           [value] = int, a single value
           [value] = [beg, end, frq], from a beginning to an end with a given frequency
@@ -13366 +13385 @@
                 if dictv == -1:
                     slicevals.append(slice(0,varshape[idim]))
+                    slicefinaldims.append(dim)
+                elif dictv == -9:
+                    slicevals.append(varshape[idim]-1)
                     slicefinaldims.append(dim)
                 else: