Changeset 1737 in lmdz_wrf

Timestamp:

Dec 16, 2017, 7:15:39 PM (7 years ago)

Author:

lfita

Message:

Adding user providing resolution on `reconstruct_matrix_from_vector'

File:

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

trunk/tools/nc_var_tools.py

@@ -20697 +20697 @@
     """ Function to reconstruct a 2D matrix from a pair of syncronized vectors with 
       the positions on x and y coordinates
-        values=[varXpos]:[varYpos]:[Xmin]:[Xmax]:[Ymin]:[Ymax]:[matProj]:[maxdiff]
+        values=[varXpos]:[varYpos]:[Xmin]:[Xmax]:[Ymin]:[Ymax]:[matProj]:[resX]:[resY]:[maxdiff]
           [varXpos]: variable with the vector of X positions
           [varYpos]: variable with the vector of Y positions
@@ -20706 +20706 @@
           [matProj]: assumed projection of the final matrix
             'latlon': regular lon/lat projection
+          [resX]: resolution along the x-axis (to be interpreted according to the projection)
+          [resY]: resolution along the y-axis (to be interpreted according to the projection)
           [maxdiff]: maximum allowed distance between matrix grid point and vector location
         ncfile= netCDF file to use
@@ -20721 +20723 @@
 
     expectargs = '[varXpos]:[varYpos]:[Xmin]:[Xmax]:[Ymin]:[Ymax]:[matProj]:' +      \
-      '[maxdiff]'
+      '[resX]:[resY]:[maxdiff]'
     gen.check_arguments(fname, values, expectargs, ':')
 
@@ -20734 +20736 @@
         Xmin = values.split(':')[2]
     matProj = values.split(':')[6]
-    maxdiff = np.float(values.split(':')[7])
+    resX = np.float(values.split(':')[7])
+    resY = np.float(values.split(':')[8])
+    maxdiff = np.float(values.split(':')[9])
 
     if not gen.searchInlist(availProj, matProj):
@@ -20799 +20803 @@
     # Matrix values
     if matProj == 'latlon':
-        ddx = xposv[1]- xposv[0]
-        ddy = yposv[1]- yposv[0]
-        dimx = (Xmax - Xmin+ddx)/ddx
-        dimy = (Ymax - Ymin+ddy)/ddy
+        dimx = (Xmax - Xmin+resX)/resX
+        dimy = (Ymax - Ymin+resY)/resY
 
     print dir(Sci)