Changeset 280 in lmdz_wrf

Timestamp:

Feb 25, 2015, 12:12:56 PM (10 years ago)

Author:

lfita

Message:

Adding slices when we are not at a border point

File:

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

trunk/tools/nc_var_tools.py

@@ -10446 +10446 @@
                 yrangeslice2D.append([yinit2D, yend2D])
                 xrangeslice2D.append([xinit2D, xend2D])
+            else:
+                yrangeslice.append([gtrajvals[it,2]-box2, gtrajvals[it,2]+box2 + 1])
+                xrangeslice.append([gtrajvals[it,1]-box2, gtrajvals[it,1]+box2 + 1])
+                yrangeslice2D.append([gtrajvals[it,2]-box2, gtrajvals[it,2]+box2 + 1])
+                xrangeslice2D.append([gtrajvals[it,1]-box2, gtrajvals[it,1]+box2 + 1])
 
 # circle values
@@ -10486 +10491 @@
                 cyrangeslice2D.append([yinit2D, yend2D])
                 cxrangeslice2D.append([xinit2D, xend2D])
+            else:
+                cyrangeslice.append([gtrajvals[it,2]-Nrad, gtrajvals[it,2]+Nrad + 1])
+                cxrangeslice.append([gtrajvals[it,1]-Nrad, gtrajvals[it,1]+Nrad + 1])
+                cyrangeslice2D.append([gtrajvals[it,2]-Nrad, gtrajvals[it,2]+Nrad + 1])
+                cxrangeslice2D.append([gtrajvals[it,1]-Nrad, gtrajvals[it,1]+Nrad + 1])
+
 
         iline = iline + 1
 
     trajobj.close()
+
+    print 'Lluis _________'
+    print yrangeslice
+
+    quit()
 
 # Creation of the netCDF file
@@ -10526 +10542 @@
  
             if not objofile.dimensions.has_key('z'):
-                newdim = objofile.createDimension('z', dimt)
+                newdim = objofile.createDimension('z', dimz)
 
                 varzobj = objfile.variables[zn]
@@ -10560 +10576 @@
             for it in range(dimt):
 # box values
-                slicev.appen(slice(0,dimz))
-                slicev.appen(slice(yrangeslice[it]))
-                slicev.appen(slice(xrangeslice[it]))
-
-                slicevnoT.appen(slice(0,dimz))
-                slicevnoT.appen(slice(yrangeslice[it]))
-                slicevnoT.appen(slice(xrangeslice[it]))
-
-                slice2D.appen(slice(0,dimz))
-                slice2D.appen(slice(yrangeslice2D[it]))
-                slice2D.appen(slice(xrangeslice2D[it]))
+                print it,'/',dimt, yrangeslice
+                quit()
+                slicev.append(slice(0,dimz))
+                slicev.append(slice(yrangeslice[it]))
+                slicev.append(slice(xrangeslice[it]))
+
+                slicevnoT.append(slice(0,dimz))
+                slicevnoT.append(slice(yrangeslice[it]))
+                slicevnoT.append(slice(xrangeslice[it]))
+
+                slice2D.append(slice(0,dimz))
+                slice2D.append(slice(yrangeslice2D[it]))
+                slice2D.append(slice(xrangeslice2D[it]))
 
                 if gtrajvals[it,2]-box2 < 0 or gtrajvals[it,2]+box2 + 1 > dimy + 1   \
@@ -10622 +10640 @@
 
 # Circle values
-                cslicev.appen(slice(0,dimz))
-                cslicev.appen(slice(cyrangeslice[it]))
-                cslicev.appen(slice(cxrangeslice[it]))
-
-                cslicevnoT.appen(slice(0,dimz))
-                cslicevnoT.appen(slice(cyrangeslice[it]))
-                cslicevnoT.appen(slice(cxrangeslice[it]))
-
-                cslice2D.appen(slice(0,dimz))
-                cslice2D.appen(slice(cyrangeslice2D[it]))
-                cslice2D.appen(slice(cxrangeslice2D[it]))
+                cslicev.append(slice(0,dimz))
+                cslicev.append(slice(cyrangeslice[it]))
+                cslicev.append(slice(cxrangeslice[it]))
+
+                cslicevnoT.append(slice(0,dimz))
+                cslicevnoT.append(slice(cyrangeslice[it]))
+                cslicevnoT.append(slice(cxrangeslice[it]))
+
+                cslice2D.append(slice(0,dimz))
+                cslice2D.append(slice(cyrangeslice2D[it]))
+                cslice2D.append(slice(cxrangeslice2D[it]))
 
                 if gtrajvals[it,2]-Nrad < 0 or gtrajvals[it,2]+Nrad + 1 > dimy + 1   \
@@ -10771 +10789 @@
             for it in range(dimt):
 # box values
-                slicev.appen(slice(yrangeslice[it]))
-                slicev.appen(slice(xrangeslice[it]))
-
-                slicevnoT.appen(slice(yrangeslice[it]))
-                slicevnoT.appen(slice(xrangeslice[it]))
-
-                slice2D.appen(slice(yrangeslice2D[it]))
-                slice2D.appen(slice(xrangeslice2D[it]))
+                slicev.append(slice(yrangeslice[it]))
+                slicev.append(slice(xrangeslice[it]))
+
+                slicevnoT.append(slice(yrangeslice[it]))
+                slicevnoT.append(slice(xrangeslice[it]))
+
+                slice2D.append(slice(yrangeslice2D[it]))
+                slice2D.append(slice(xrangeslice2D[it]))
 
                 if gtrajvals[it,2]-box2 < 0 or gtrajvals[it,2]+box2 + 1 > dimy + 1   \
@@ -10839 +10857 @@
 
 # Circle values
-                cslicev.appen(slice(cyrangeslice[it]))
-                cslicev.appen(slice(cxrangeslice[it]))
-
-                cslicevnoT.appen(slice(cyrangeslice[it]))
-                cslicevnoT.appen(slice(cxrangeslice[it]))
-
-                cslice2D.appen(slice(cyrangeslice2D[it]))
-                cslice2D.appen(slice(cxrangeslice2D[it]))
+                cslicev.append(slice(cyrangeslice[it]))
+                cslicev.append(slice(cxrangeslice[it]))
+
+                cslicevnoT.append(slice(cyrangeslice[it]))
+                cslicevnoT.append(slice(cxrangeslice[it]))
+
+                cslice2D.append(slice(cyrangeslice2D[it]))
+                cslice2D.append(slice(cxrangeslice2D[it]))
                 if gtrajvals[it,2]-Nrad < 0 or gtrajvals[it,2]+Nrad + 1 > dimy + 1   \
                   or gtrajvals[it,1]-Nrad < 0 or gtrajvals[it,1]+Nrad + 1 > dimx + 1: