Changeset 819 in lmdz_wrf for trunk/tools

Timestamp:

Jun 12, 2016, 3:53:24 PM (9 years ago)

Author:

lfita

Message:

Adding `trip' variable and meaning of its in 'subbasins'

File:

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

trunk/tools/nc_var_tools.py

@@ -14929 +14929 @@
         latname = 'nav_lat'
         tripname = 'trip'
+        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']
     else:
         print errormsg
         print '  ' + fname + ": file '" + ncfile + "' not ready!!"
         quit(-1)
+    Ndiroutflow = len(Soutflows)
 
     olons = onc.variables[lonname]
@@ -14973 +14976 @@
     newdim = onewnc.createDimension('Lstring', Lstr)
     newdim = onewnc.createDimension('Ngridsubflow', Nmaxgridsubflow)
+    newdim = onewnc.createDimension('outflow', Ndiroutflow)
     
     if ncfile == 'routing.nc':
@@ -15006 +15010 @@
     onewnc.sync()
 
+# sub-flow grid point data
     newvar = onewnc.createVariable('Nsubflow', 'i', ('Nsubflow'))
     basicvardef(newvar, 'Nsubflow', 'Number of grid points of the sub-flow', '-')
@@ -15014 +15019 @@
     vardef = basicvardef(newvarlat, 'latsubflow', 'Latitude of the grid point of ' + \
       'the sub-flow','degrees_North')
+    newvartrip = onewnc.createVariable('outsubflow', 'i', ('Nsubflow','Ngridsubflow'))
+    vardef = basicvardef(newvartrip, 'outsubflow', 'out-flow of the grid point of '+ \
+      'the sub-flow','outflow')
+
     isub = 0
     for subf in nsubflows:
@@ -15024 +15033 @@
             lonv = lons[jigrid[0],jigrid[1]]
             latv = lats[jigrid[0],jigrid[1]]
+            tripv = outflows[jigrid[0],jigrid[1]]
 #            print '    jigrid:', jigrid,' lon:', lonv,'lat:',latv
             newvarlon[isub,igrid] = lonv
             newvarlat[isub,igrid] = latv
+            newvartrip[isub,igrid] = tripv
+
         isub = isub + 1
     onewnc.sync()
 
+    newvar = onewnc.createVariable('outflow', 'c', ('outflow', 'Lstring'))
+    vardef = basicvardef(newvar, 'outflow', 'out-flow direction','1')
+    newvals = writing_str_nc(newvar, Soutflows, Lstr)
+    onewnc.sync()
+    
 # Global attributes
     onewnc.setncattr('script',  fname)