Changeset 741 in lmdz_wrf for trunk

Timestamp:

May 3, 2016, 9:56:09 AM (9 years ago)

Author:

lfita

Message:

Fixing `sellonlatbox'
Replacing global institutions' attribute for the unicode version with accents

File:

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

trunk/tools/nc_var_tools.py

@@ -9330 +9330 @@
         quit(-1)
 
-
     lonobj = objfile.variables[lonn]
     latobj = objfile.variables[latn]
@@ -9383 +9382 @@
 
     if len(lonv.shape) == 1:
+        dlonn = loninf.dimns[0]
+        dlatn = latinf.dimns[0]
         dimx = lonv.shape[0]
         dimy = latv.shape[0]
@@ -9389 +9390 @@
             malonlat[:,i] = malon[i] + malat
     elif len(lonv.shape) == 2:
+        dlonn = loninf.dimns[1]
+        dlatn = loninf.dimns[0]
         dimx = lonv.shape[1]
         dimy = lonv.shape[0]
@@ -9398 +9401 @@
         print '  ' + fname + ': assuming that lon/lat is constant to the first dimension!!'
         print '    shape lon:', lonv.shape, 'lat:', latv.shape
+        dlonn = loninf.dimns[2]
+        dlatn = loninf.dimns[1]
         dimx = lonv.shape[2]
         dimy = lonv.shape[1]
@@ -9442 +9447 @@
     newdimx = elon - ilon + 1
     newdimy = elat - ilat + 1
+
+    if len(malon.shape) == 1:
+        datalonSW = lonv[ilon]
+        datalatSW = latv[ilat]
+        datalonNE = lonv[elon]
+        datalatNE = latv[elat]
+    elif len(malon.shape) == 2:
+        datalonSW = lonv[ilat,ilon]
+        datalatSW = latv[ilat,ilon]
+        datalonNE = lonv[elat,elon]
+        datalatNE = latv[elat,elon]
+    elif len(malon.shape) == 3:
+        datalonSW = lonv[0,ilat,ilon]
+        datalatSW = latv[0,ilat,ilon]
+        datalonNE = lonv[0,elat,elon]
+        datalatNE = latv[0,elat,elon]
+    
+    print '  ' + fname + ': final SW:', datalonSW, ',', datalatSW, 'NE:', datalonNE, \
+      ',', datalatNE
     print '  ' + fname + ': new lon size:', ilon, ',', elon, 'new dimx:', newdimx
     print '  ' + fname + ': new lat size:', ilat, ',', elat, 'new dimy:', newdimy
@@ -9448 +9472 @@
 
 # Creation of lon,lat related dims
-    newdim = newncobj.createDimension(lonn, newdimx)
-    newdim = newncobj.createDimension(latn, newdimy)
+    newdim = newncobj.createDimension(dlonn, newdimx)
+    newdim = newncobj.createDimension(dlatn, newdimy)
 
 # Adding lonn/latn to the output file
@@ -9476 +9500 @@
             if not searchInlist(newncobj.dimensions, dimn):
                 newncobj.createDimension(dimn, len(objfile.dimensions[dimn]))
-            if dimn == lonn:
+            if dimn == dlonn:
                 varslice.append(slice(ilon,elon+1))
-            elif dimn == latn:
+            elif dimn == dlatn:
                 varslice.append(slice(ilat,elat+1))
             else:
@@ -9503 +9527 @@
 # global attributes
     newncobj.setncattr('author', 'L. Fita')
-    newncobj.setncattr('institution', 'Laboratire de Meteorologie Dynamique')
+    newattr = set_attributek(newncobj, 'institution', unicode('Laboratoire de M' +       \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
     newncobj.setncattr('university', 'Pierre Marie Curie - Jussieu')
     newncobj.setncattr('center', 'Centre National de Recherches Scientifiques')
@@ -9511 +9536 @@
     newncobj.setncattr('function', 'sellonlatbox')
     newncobj.setncattr('version', '1.0')
+    newncobj.setncattr('original_file', ncfile)
+    newnattr = set_attributek(newncobj, 'lonSW', lonSW, 'R')
+    newnattr = set_attributek(newncobj, 'latSW', latSW, 'R')
+    newnattr = set_attributek(newncobj, 'lonNE', lonNE, 'R')
+    newnattr = set_attributek(newncobj, 'latNE', latNE, 'R')
+    newnattr = set_attributek(newncobj, 'data_lonSW', datalonSW, 'R')
+    newnattr = set_attributek(newncobj, 'data_latSW', datalatSW, 'R')
+    newnattr = set_attributek(newncobj, 'data_lonNE', datalonNE, 'R')
+    newnattr = set_attributek(newncobj, 'data_latNE', datalatNE, 'R')
+    newnattr = set_attributek(newncobj, 'data_pt_lonSW', ilon, 'I')
+    newnattr = set_attributek(newncobj, 'data_pt_latSW', ilat, 'I')
+    newnattr = set_attributek(newncobj, 'data_pt_lonNE', elon, 'I')
+    newnattr = set_attributek(newncobj, 'data_pt_latNE', elat, 'I')
 
     objfile.close()
@@ -17442 +17480 @@
     onewfile.setncattr('version',  '1.0')
     onewfile.setncattr('author',  'L. Fita')
-    onewfile.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+    newattr = set_attributek(onewnfile, 'institution', unicode('Laboratoire de M' +  \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
     onewfile.setncattr('university',  'Pierre et Marie Curie')
     onewfile.setncattr('country',  'France')
@@ -18446 +18485 @@
     onewnc.setncattr('version',  '1.0')
     onewnc.setncattr('author',  'L. Fita')
-    onewnc.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+    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')
@@ -18792 +18832 @@
         onewnc.setncattr('version',  '1.0')
         onewnc.setncattr('author',  'L. Fita')
-        onewnc.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+        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')
@@ -19686 +19727 @@
     newnc.setncattr('version',  '1.0')
     newnc.setncattr('author',  'L. Fita')
-    newnc.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+    newattr = set_attributek(newnc, 'institution', unicode('Laboratoire de M' +      \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
     newnc.setncattr('university',  'Pierre et Marie Curie')
     newnc.setncattr('country',  'France')
@@ -20083 +20125 @@
     newnc.setncattr('version',  '1.0')
     newnc.setncattr('author',  'L. Fita')
-    newnc.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+    newattr = set_attributek(newnc, 'institution', unicode('Laboratoire de M' +      \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
     newnc.setncattr('university',  'Pierre et Marie Curie')
     newnc.setncattr('country',  'France')
@@ -20484 +20527 @@
     newnc.setncattr('version',  '1.0')
     newnc.setncattr('author',  'L. Fita')
-    newnc.setncattr('institution',  'Laboratoire de Meteorology Dynamique')
+    newattr = set_attributek(newnc, 'institution', unicode('Laboratoire de M' +      \
+      unichr(233) + 't' + unichr(233) + 'orologie Dynamique'), 'U')
     newnc.setncattr('university',  'Pierre et Marie Curie')
     newnc.setncattr('country',  'France')