Changeset 940 in lmdz_wrf for trunk/tools

Timestamp:

Jun 24, 2016, 9:52:21 AM (8 years ago)

Author:

lfita

Message:

Working of new version of `DataSetSection?' without monotones

File:

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

trunk/tools/nc_var_tools.py

@@ -5708 +5708 @@
 
                 if endv == -1: iendv = len(objdim)
-                print 'Lluis endv, begv, intv:', ibegv, iendv, iintv,':', begv, endv, intv
                 dimsize = (iendv - ibegv + 1)/iintv - 1
             else:
@@ -5716 +5715 @@
                     iendv = len(objdim) - 1
                     iintv = 0
-            print fname + '; Lluis dimn beg end int:', ibegv, iendv, iintv
         else:
             if len(objdim) > 1:
@@ -5729 +5727 @@
                 iintv = 0
 
-        print fname + '; Lluis begv, endv:', ibegv, iendv, iintv, '<>', begv, endv, intv
         if ibegv != iendv:
             dimsecslice[dims] = slice(ibegv,iendv,iintv)
@@ -5745 +5742 @@
 
     ncoobj.sync()
-    print fname + '; Lluis: dimsecslice:', dimsecslice
 # Creating variables
 ## 
@@ -5757 +5753 @@
         if not gen.searchInlist(list(vardims),dimn):
             # removing monotones
-            print fname + '; Lluis looking for monotones!'
             varslice = []
             filevardims = []
             for dimname in varobj.dimensions:
                 lendim = len(nciobj.dimensions[dimname])
-                print fname + '; Lluis monotone?:', dimname,' len;', lendim
                 if lendim > 1:
                     varslice.append(slice(0, lendim))
@@ -5776 +5770 @@
                 if dimname == dimn:
                     varslice.append(dimsecslice[dimn])
-                    filevardims.append(dimname)
+                    if type(dimsecslice[dimn]) == type(slice(1)):
+                        filevardims.append(dimname)
                 else:
                     Ldim = len(nciobj.dimensions[dimname])
@@ -5785 +5780 @@
                         varslice.append(0)
 
-                print fname  + '; Lluis varslice:', varslice,' filevardims:', filevardims
             varvals = varobj[tuple(varslice)]
 
-        print fname + ': Lluis filevardims:', filevardims
 # Adding fill value attribute
 ##