Changeset 824 in lmdz_wrf for trunk

Timestamp:

Jun 12, 2016, 11:03:48 PM (8 years ago)

Author:

lfita

Message:

Adding 2 new errors:

• When there is only one grid in all the sub-basin
• When there is not coincidence (distance zero) with any grid point from the file

File:

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

trunk/tools/nc_var_tools.py

@@ -14927 +14927 @@
 
     Lncfile = len(ncfile)
-    filen = ncfile[Lncfile-10:Lncfile]
-
-    if filen == 'routing.nc':
+
+    if ncfile[Lncfile-10:Lncfile] == 'routing.nc':
         lonname = 'nav_lon'
         latname = 'nav_lat'
@@ -14935 +14934 @@
         Soutflows = ['1: N', '2: NE', '3: E', '4: SE', '5: S', '6: SW', '7: W',      \
           '8: NW', '97: sub-basin', '98: small to sea', '99: large to sea']
+        filen = 'routing.nc'
     else:
         print errormsg
@@ -14947 +14947 @@
     lons = olons[:]
     lats = olats[:]
-    outflows = ooutflows[:]
+    outflows = np.array(ooutflows[:], dtype=int)
 
     difflonlat = np.sqrt((lons - sublon)**2 + (lats - sublat)**2)
     mindifflonlat = np.min(difflonlat)
+    if mindifflonlat != 0.:
+        print errormsg
+        print '  ' + fname + ': wrong initial lon,lat of the basin!!'
+        print '    for lon,lat:', sublon,',',sublat,'there is not grid point ' +     \
+          "coincidence in file '" + ncfile + "'"
+        quit(-1)
+
     isub = gen.index_mat(difflonlat, mindifflonlat)
     print '    ' + fname + "initial point of subbasin '" + subname + "' in the " +   \
-      'matrix:', isub, 'mindiff:', mindifflonlat
+      'matrix:', isub
 
     masksubbasin, subflows, subflowspt = gen.subbasin_point(outflows, isub)
+
+    if len(subflows.keys()) == 1:
+        print errormsg
+        print '  ' + fname + ": Only 1 point for subbasin '" + subname +             \
+          "' at lon,lat:", sublon,',',sublat
+        print ' values around (3x3 centered at: ', isub, ')'
+        lon3x3 = gen.vals_around(lons, isub)
+        lat3x3 = gen.vals_around(lats, isub)
+        out3x3 = gen.vals_around(outflows, isub)
+        print ' lon lat out _______'
+        for i in range(3):
+            print lon3x3[i,:], lat3x3[i,:], out3x3[i,:]
+        print 'out equivalencies:', Soutflows
+        quit(-1)
 
     onewnc = NetCDFFile(ofile, 'w')