Changeset 410 in lmdz_wrf

Timestamp:

May 5, 2015, 11:33:03 AM (10 years ago)

Author:

lfita

Message:

Fixing indices to reach all the domain

File:

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

trunk/tools/nc_var_tools.py

@@ -7430 +7430 @@
             j10 = firstmappairs[0,jnew,inew]
 
-            oldlon2d2 = oldlon2d[slice(j10-dy2,j10+dy2), slice(i10-dx2,i10+dx2)]
-            oldlat2d2 = oldlat2d[slice(j10-dy2,j10+dy2), slice(i10-dx2,i10+dx2)]
+            imin = np.max([i10-dx2, 0])
+            imax = np.min([i10+dx2, olddimx+1])
+            jmin = np.max([j10-dy2, 0])
+            jmax = np.min([j10+dy2, olddimy+1])
+
+            oldlon2d2 = oldlon2d[slice(jmin,jmax), slice(imin,imax)]
+            oldlat2d2 = oldlat2d[slice(jmin,jmax), slice(imin,imax)]
 
             dist=np.sqrt((newlon2d[jnew, inew] - oldlon2d2)**2. +                    \
@@ -7526 +7531 @@
 #                       '% done'
 
+#                    print 'slices new:',newslice0,'old:',oldslice0
                     newvarval[tuple(newslice0)] = ovar[tuple(oldslice0)]