Changeset 1396 in lmdz_wrf for trunk/tools

Timestamp:

Dec 20, 2016, 4:25:00 PM (8 years ago)

Author:

lfita

Message:

Fixing `Nwgts' in 'curvelocalize_2D'

File:

• trunk/tools/generic_tools.py (modified) (3 diffs)

trunk/tools/generic_tools.py

@@ -10972 +10972 @@
 
     curvloc = np.ones((curve.shape), dtype=int)*(-1)
+    Nwgts = np.zeros((Npts), dtype=int)
     curveweights = np.zeros((Npts,3,3), dtype=float)
     ijaround3x3 = np.zeros((Npts,2,9), dtype=int)
@@ -10992 +10993 @@
         yaroundvls = vals_around(yvals,ijloc)
 
-        curveweights[icv,:,:], Nwgt, ijaround3x3[icv,:,:,] = linearint_3x3weights(    \
-          xaroundvls, yaroundvls, curve[icv,:])
+        curveweights[icv,:,:], Nwgts[icv], ijaround3x3[icv,:,:,] =                   \
+          linearint_3x3weights(xaroundvls, yaroundvls, curve[icv,:])
         
-    return curvloc, curveweights, Nwgt, ijaround3x3
+    return curvloc, curveweights, Nwgts, ijaround3x3
 
 #lon1d = np.arange(10)
@@ -11035 +11036 @@
 #for ipt in range(Ncurve):
 #    print 'pos:', ipt, 'curve:', curve[ipt,:],'loccurve:', loccurve[ipt,:], 'wgts:', wgts[ipt,:,:], 'Nwgts:', Nwgts, 'loc3x3:', loc3x3[ipt,:]
+
 #quit()