Changeset 2173 in lmdz_wrf for trunk/tools

Timestamp:

Oct 9, 2018, 8:19:21 PM (7 years ago)

Author:

lfita

Message:

Using right indices!!!

File:

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

trunk/tools/nc_var_tools.py

@@ -25954 +25954 @@
           str(stlats[ist]))
         stdataij[stns[ist]] = [ijst, dst]
-        dlon = lon[ijst[0],ijst[1]]
-        dlat = lat[ijst[0],ijst[1]]
+        dlon = lon[ijst[1],ijst[0]]
+        dlat = lat[ijst[1],ijst[0]]
 
         # Writing in file station values from grid
-        newvari[ist] = ijst[1]
-        newvarj[ist] = ijst[0]
+        newvari[ist] = ijst[0]
+        newvarj[ist] = ijst[1]
         newvarlon[ist] = dlon
         newvarlat[ist] = dlat
@@ -25971 +25971 @@
             weights = np.ones((1), dtype=np.float)
         elif method == 'mean' or method == 'dist':
-            ddx = lon[ijst[0],ijst[1]+1] - dlon
-            ddy = lat[ijst[0]+1,ijst[1]] - dlat
+            ddx = lon[ijst[1],ijst[0]+1] - dlon
+            ddy = lat[ijst[1]+1,ijst[0]] - dlat
             signx = np.abs(ddx)/ddx
             signy = np.abs(ddy)/ddy
@@ -25982 +25982 @@
  
             # Avoiding outside shape
-            xij = np.min([ijst[1]+incx, ddimx-1])
-            yij = np.min([ijst[0]+incy, ddimy-1])
+            xij = np.min([ijst[0]+incx, ddimx-1])
+            yij = np.min([ijst[1]+incy, ddimy-1])
             xij = int(np.max([0, xij]))
             yij = int(np.max([0, yij]))
 
             ptvslice.append(ijst)
-            ptvslice.append([ijst[0],xij])
-            ptvslice.append([yij,ijst[1]])
+            ptvslice.append([ijst[1],xij])
+            ptvslice.append([yij,ijst[0]])
             ptvslice.append([yij,xij])