Changeset 817 in lmdz_wrf for trunk

Timestamp:

Jun 12, 2016, 2:26:48 PM (8 years ago)

Author:

lfita

Message:

Working version 1.0 of `bassins'

File:

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

trunk/tools/nc_var_tools.py

@@ -14854 +14854 @@
 
 def lonlatvarsFile(lonvals,latvals,dnx,dny,oLlnc):
-    """ Functino to CF-define the longitudes and latitudes variables within a file
+    """ Function to CF-define the longitudes and latitudes variables within a file
       lonvals= 2D matrix with longitudes
       latvals= 2D matrix with latitudes
@@ -14942 +14942 @@
     outflows = ooutflows[:]
 
-    onc.close()
-
     difflonlat = np.sqrt((lons - sublon)**2 + (lats - sublat)**2)
     mindifflonlat = np.min(difflonlat)
@@ -15021 +15019 @@
         Ngrid = len(grids)
         newvar[isub] = Ngrid
-        print '  ' + fname + 'Lluis; isub:', isub,'subf:',subf,'grids:',grids,'Ngrid',Ngrid
+#        print '  ' + fname + ' Lluis; isub:', isub,'subf:',subf,'grids:',grids,'Ngrid',Ngrid
         for igrid in range(Ngrid):
             jigrid = subflowspt[grids[igrid]]
             lonv = lons[jigrid[0],jigrid[1]]
             latv = lats[jigrid[0],jigrid[1]]
-            newvarlon[igrid,isub] = lonv
-            newvarlat[igrid,isub] = latv
+#            print '    jigrid:', jigrid,' lon:', lonv,'lat:',latv
+            newvarlon[isub,igrid] = lonv
+            newvarlat[isub,igrid] = latv
         isub = isub + 1
+    onewnc.sync()
+
+# Global attributes
+    onewnc.setncattr('script',  fname)
+    onewnc.setncattr('version',  '1.0')
+    onewnc.setncattr('author',  'L. Fita')
+    newattr = set_attributek(onewnc, 'institution', unicode('Laboratoire de M' +     \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
+    onewnc.setncattr('university',  'Pierre et Marie Curie')
+    onewnc.setncattr('country',  'France')
+    for attrs in onc.ncattrs():
+        attrv = onc.getncattr(attrs)
+        attr = set_attribute(onewnc, attrs, attrv)
+
+    onc.close()
     onewnc.sync()
     onewnc.close()