Changeset 306 in lmdz_wrf

Timestamp:

Feb 25, 2015, 8:01:31 PM (10 years ago)

Author:

lfita

Message:

Wrong time slice for circdist

File:

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

trunk/tools/nc_var_tools.py

@@ -10976 +10976 @@
                 cslice2D = []
                 cslicevnoT = []
+                cslicevC = []
 
                 slicev.append(gtrajvals[it,0])
@@ -11087 +11088 @@
                     cslice2D.append(slice(gtrajvals[it,1]-Nrad,gtrajvals[it,1]+Nrad+1))
 
+                    cslicevC.append(it)
+                    cslicevC.append(slice(gtrajvals[it,2]-Nrad,gtrajvals[it,2]+Nrad+1))
+                    cslicevC.append(slice(gtrajvals[it,1]-Nrad,gtrajvals[it,1]+Nrad+1))
+
                     rvarvalst = varobj[tuple(cslicev)]
-                    cdist = circdist[tuple(cslicev)]
+                    cdist = circdist[tuple(cslicevC)]
 # circle values
                     rvarvalst = np.where(cdist > np.float(Nrad),fillValue,rvarvalst)