Changeset 2358 in lmdz_wrf for trunk

Timestamp:

Feb 20, 2019, 11:06:22 PM (6 years ago)

Author:

lfita

Message:

Fixing in all places when dimy slice becomes zero sized ! Ndy == 0 in `pinterp'

File:

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

trunk/tools/nc_var_tools.py

@@ -19071 +19071 @@
 
                 Ndy = len(range(yini,yend))
-                if Ndy == 0: continue
+                if Ndy == 0: 
+                    print '      empty dimy slice: yini == yend. Skipping it !!'
+                    continue
                 if gen.searchInlist(notCHK, vn):
                     if vn == 'WRFght': 
@@ -19261 +19263 @@
                     tempt = temp[tini:tend,:,yini:yend,:].transpose()
                     qvt = qv[tini:tend,:,yini:yend,:].transpose()
-                    print 'Lluis shapes :', yini, yend, 'hgt', hgt[yini:yend,:].shape, 'Ndy:', Ndy, dimy
                     varinterp = fin.module_forinterpolate.interp( data_in=varint,    \
                       pres_field=prest, interp_levels=interplevs, psfc=psfct,        \
@@ -19282 +19283 @@
                 yend = iyy + dimyfrac
                 Ndy = len(range(yini,yend))
-
+                if Ndy == 0: 
+                    print '        empty dimy slice: yini == yend. Skipping it !!'
+                    continue
                 if gen.searchInlist(notCHK, vn):
                     if vn == 'WRFght': 
@@ -19489 +19492 @@
 
                 Ndy = len(range(yini,yend))
+                if Ndy == 0: 
+                    print '        empty dimy slice: yini == yend. Skipping it !!'
+                    continue
                 if gen.searchInlist(notCHK, vn):
                     if vn == 'WRFght':