Changeset 1702 in lmdz_wrf

Timestamp:

Dec 8, 2017, 6:24:45 PM (7 years ago)

Author:

lfita

Message:

Fixing and working copy of the script

Location:

trunk/tools

Files:

• diag_tools.py (modified) (1 diff)
• get_stations.py (modified) (2 diffs)

trunk/tools/diag_tools.py

@@ -1334 +1334 @@
 
             if len(ta.shape) != len(p.shape):
-                p = fill_Narray(p, ta*0., filldim=[0,2,3])
+                p = gen.fill_Narray(p, ta*0., filldim=[0,2,3])
 
             self.values = var_td(ta, p, hur)

trunk/tools/get_stations.py

@@ -393 +393 @@
     if len(Xvarvals.shape) == 1:
         newXvarvals, newYvarvals = np.meshgrid(Xvarvals, Yvarvals) 
-        dimvarvalues[dimvariables[dimvariables[xn]] = newXvarvals
+        dimvarvalues[dimvariables[xn]] = newXvarvals
+        if len(Yvarvals.shape) == 1:
+            dimvarvalues[dimvariables['lat']] = newYvarvals
+
 for yn in axesvars['Y']:
     Xvarvals = dimvarvalues[CFdims['lon']]
@@ -399 +402 @@
     if len(Yvarvals.shape) == 1:
         newXvarvals, newYvarvals = np.meshgrid(Xvarvals, Yvarvals) 
-        dimvarvalues[dimvariables[dimvariables[yn]] = newYvarvals
+        dimvarvalues[dimvariables[yn]] = newYvarvals
 
 # Retrieving surface data