Changeset 744 in lmdz_wrf

Timestamp:

May 3, 2016, 5:36:37 PM (9 years ago)

Author:

lfita

Message:

Fixing right initial point!

File:

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

trunk/tools/nc_var_tools.py

@@ -20419 +20419 @@
 
 # Reducing the searching points
-    newvarinvals = newvarinpt[:]
+    Ninpts=np.mini([499999, Ninpts-ipoint])
+    newvarinvals = newvarinpt[ipoint:ipoint+Ninpts]
     maskpt = np.where(newvarinvals.mask == True, False, True)
-    points = np.arange(Ninpts)
+    points = np.arange(ipoint,ipoint+Ninpts,1)
     mapoints = ma.array(points, mask=maskpt)
     ptsf = mapoints.compressed()
@@ -20427 +20428 @@
     Nptsf = len(ptsf)
     print Ninpts,'Npoints to find:', len(ptsf), ptsf[0:10], newvarindiff[ptsf[0:10]]
-# Error at 150024, 150025, 151709, 153421
-    print '  ' + fname + ': from:', Ninpts,'re-locating:',Nptsf,'points...'
+    print '  ' + fname + ': from:', Ninpts,'re-locating:',Nptsf,'points starting at',\
+      ipoint,'...'
     if kind == 'Goode':
 #        newvar,newvarin,newvarinpt,newvarindiff =                                    \
 #          fin.module_forinterpolate.coarseinterpolate(projlon, projlat, lonvs,       \
 #          latvs, percen, mindiff, ivar, dimx, dimy, ninpts)
-        Ninpts=499999
-        for ir in range(ptsf[0],Ninpts,fracs):
+        for ir in range(ipoint,ipoint+Ninpts,fracs):
             iri = ir
             ire = ir + fracs