Changeset 2536 in lmdz_wrf

Timestamp:

May 16, 2019, 5:05:50 PM (6 years ago)

Author:

lfita

Message:

Fixing `sellonlatbox' no mask on numpy arrays !

File:

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

trunk/tools/nc_var_tools.py

@@ -7026 +7026 @@
         for j in range(dimy):
             for i in range(dimx):
-                if malonlat.mask[j,i] == False and i < ilon: 
+#                if malonlat.mask[j,i] == False and i < ilon: 
+                if malonlat[j,i] == False and i < ilon: 
                     ilon = i
-                if malonlat.mask[j,i] == False and i > elon: 
+#                if malonlat.mask[j,i] == False and i > elon: 
+                if malonlat[j,i] == False and i > elon: 
                     elon = i
-                if malonlat.mask[j,i] == False and j < ilat: 
+#                if malonlat.mask[j,i] == False and j < ilat: 
+                if malonlat[j,i] == False and j < ilat: 
                     ilat = j
-                if malonlat.mask[j,i] == False and j > elat: 
+#                if malonlat.mask[j,i] == False and j > elat: 
+                if malonlat[j,i] == False and j > elat: 
                     elat = j
     else:
@@ -29384 +29388 @@
                 getvarybndst = ydimvarbnds2D.transpose()
 
+                print 'ref:', dxref, dyref
+                print 'Lluis looking for:', reflon, 'lat.', reflat
+                for j in range(dyref):
+                    print 'Lluis ref lat.', yslice2D[:,j,0]
+                print 'Lluis get lon.', getlon[0,:]
+                print 'Lluis get lat.', getlat[:,0]
+                for i in range(dxget):
+                    print i,'Lluis get lat.', xdimvarbnds2D[:,0,i]
+                for j in range(dyget):
+                    print j,'Lluis get lat.', ydimvarbnds2D[:,j,0]
+
                 Ngridsint, gridsint, areas2Dt, areast, percenst =                    \
                   fsci.module_scientific.grid_spacepercen(xcavals=reflont,           \