Changeset 2808 in lmdz_wrf

Timestamp:

Apr 6, 2020, 3:44:48 PM (5 years ago)

Author:

lfita

Message:

Fixing `interp_curve' indices referred to the center of the 3x3 box, thus adding 'ij3x3 = [-1, 0, 1]'

File:

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

trunk/tools/nc_var_tools.py

@@ -22489 +22489 @@
     curvevar = np.zeros(tuple(dimref+[Ncurv]), dtype=np.float)
 
+    ij3x3 = [-1, 0, 1]
+
     for icv in range(Ncurv):
         # Getting all `useful' weights
@@ -22497 +22499 @@
 
         for isp in range(Nwgt[icv]):
-            dictslice[ydn] = iy + int(ijcurv[icv,0,isp])
-            dictslice[xdn] = ix + int(ijcurv[icv,1,isp])
+            dictslice[ydn] = iy + ij3x3[int(ijcurv[icv,0,isp])]
+            dictslice[xdn] = ix + ij3x3[int(ijcurv[icv,1,isp])]
             varslice, dimvarvals = SliceVarDict(ov, dictslice)
             iwgt = sortedwgt[isp]