Changeset 2169 in lmdz_wrf

Timestamp:

Oct 5, 2018, 10:26:17 PM (7 years ago)

Author:

lfita

Message:

Fixing `retrieve_stations' when data does not have height

File:

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

trunk/tools/nc_var_tools.py

@@ -25834 +25834 @@
     olat = donc.variables[datavlatn]
     lat0 = olat[:]
-    oheight = donc.variables[datavheightn]
-    height0 = oheight[:]
 
     lon, lat = gen.lonlat2D(lon0, lat0)
@@ -25842 +25840 @@
 
     if datavheightn != 'None':
+        oheight = donc.variables[datavheightn]
+        height0 = oheight[:]
         if len(height0.shape) == 3: height = height0[0,:,:]
         elif len(height0.shape) == 2: height = height0[:]
@@ -26100 +26100 @@
 
             stgridvals[0] = ist
-            if ist == 1:
-                print 'Lluis shapes ____'
-                print 'newvar:', newvar.shape
-                print 'vslice:', vslice
-                print 'vstslice:', vstslice
-                print 'stgridvals:', stgridvals
-                print 'varvalues:', varvalues.shape
-                print 'weights:', vwghts.shape
-                print 'stvarv:', stvarv.shape
-                vvv =  varvalues*vwghts
-                print ' Lluis shape prod:', vvv.shape
-
             newvar[tuple(stgridvals)] = stvarv[:]
             onewnc.sync()