Changeset 2196 in lmdz_wrf for trunk/tools

Timestamp:

Oct 18, 2018, 7:17:53 PM (6 years ago)

Author:

lfita

Message:

Adding number of grid points used in `compute_slice2Dstats'

File:

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

trunk/tools/nc_var_tools.py

@@ -26338 +26338 @@
 
         newvard = ['slice_'+varn2,'slice_'+varn1]
+        newvarN = onewnc.createVariable(varn+'Nsliced', 'f', tuple(newvard),         \
+          fill_value=gen.fillValueF)
+        basicvardef(newvarN, varn+'Nsliced','number of values in slice of ' + varn + \
+          'sliced by '+varn1+' & '+varn2, vu)
+        add_varattrs(onc, onewnc, [varn], [varn+'Nsliced'])
+
         newvarn = onewnc.createVariable(varn+'minsliced', 'f', tuple(newvard),        \
           fill_value=gen.fillValueF)
@@ -26362 +26368 @@
         add_varattrs(onc, onewnc, [varn], [varn+'stdsliced'])
 
+        newvarN[:] = np.sum(~manewvar.mask, axis=tuple(np.arange(vrank) + 2))
         newvarn[:] = np.min(manewvar, axis=tuple(np.arange(vrank) + 2))
         newvarx[:] = np.max(manewvar, axis=tuple(np.arange(vrank) + 2))
@@ -26394 +26401 @@
 #compute_slice2Dstats(values, '/home/lluis/PY/wrfout_d01_1995-01-01_00:00:00', 'T2,Q2')
 
+def same_delat_sign(values, ncfile, variable):
+    """ Function to determine if a given series of 1D values share the same sign of 
+         increase/decrease between consecutive values
+      values= [slicedims], values to the dimensions to perform a slice of variables
+        [slicedims]: '|' separated list of [dimname]:[valdim]
+        [valdim]: 
+          * [integer]: which value of the dimension
+          * -1: all along the dimension
+          * -9: last value of the dimension
+          * [beg]@[end]@[freq] slice from [beg] to [end] every [freq]
+          * NOTE, no dim name all the dimension size
+      ncfile= netCDF file to use
+      variable: ',' list of variables ('all' for all variables)
+    """
+    fname = 'same_delat_sign'
+
+    if values == 'h':
+        print fname + '_____________________________________________________________'
+        print same_delat_sign.__doc__
+        quit()
+
+    expectargs = 'varn1,minvar1,maxvar1,slcevar1,varn2,minvar2,maxvar2,slcevar2,' +  \
+      'dimvars'
+    gen.check_arguments(fname, values, expectargs, ',')
+
+    varn1 = values.split(',')[0]
+    minvar1 = np.float(values.split(',')[1])
+
+    return
+
 #quit()