Changeset 1083 in lmdz_wrf for trunk

Timestamp:

Sep 5, 2016, 7:37:48 PM (9 years ago)

Author:

lfita

Message:

Starting to introduce turb' as turbulence of a field X' = X - <X> in operdim'

File:

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

trunk/tools/nc_var_tools.py

@@ -8913 +8913 @@
     """ 
     fname = 'operdim'
-    readyop = ['max', 'mean', 'mean2', 'min', 'std', 'sum']
+    readyop = ['max', 'mean', 'mean2', 'min', 'std', 'sum', 'turb', 'var']
 
     operslice = []
@@ -8925 +8925 @@
     noper = 0
     varoper = varv.copy()
+
     for idim in range(len(varv.shape)):
         if operslice[idim]:
@@ -8940 +8941 @@
             elif oper == 'sum':
                 varoper = np.sum(varoper, axis=idim-noper)
+            elif oper == 'turb':
+                varmean = np.mean(varoper, axis=idim-noper)
+                opereddims.append(vardims[idim])
             elif oper == 'var':
                 varoper = np.var(varoper, axis=idim-noper)