Changeset 1849 in lmdz_wrf

Timestamp:

Mar 25, 2018, 3:00:18 PM (7 years ago)

Author:

lfita

Message:

Working version with projection information added

File:

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

trunk/tools/nc_var_tools.py

@@ -22689 +22689 @@
                 given axis
           [globattrfilen]: name of ASCII file with the given metadata to be added as global attributes
-            with a line per attribute [attrname] [value] it has to have at least (underparenthesis for description)
-            use '#' for comments (use '[unknown]' for that values set-up according to content from file:
-              institute_id [value, use '!' for spaces] (CV)
+            with a line per attribute [attrname] [value](use '!' for spaces]  (underparenthesis for description)
+            use '#' for comments
+            File has to have at least:
+              institute_id [value] (CV)
               institution  [value]
               model_id [value] (CORDEX starts with <institute_id>-)
@@ -22715 +22716 @@
               basetime [value] (reference time "1949-12-01T00:00:00Z" formatted acording to UNIDATA conventions)
               calendar [value] (CF conventions as in driving model)
-              grid [value] (according to CF conventions
+              grid [value] (according to CF conventions)
           [projectfile]: ASCII file with all the CF-mandatory projection information. A line for each attribute. Name of 
             the variable to be used to refer the projection inside the file will be the same as global attribute `grid'
@@ -22785 +22786 @@
         quit(-1)
 
-    pinf = {}
+    pattr = {}
     opinff = open(projfilen, 'r')
     for line in opinff:
         if line[0:1] != '#' and len(line) > 1:
             linevals = line.replace('\n','').replace('\t','').split(' ')
-            pinf[linevals[0]] = [linevals[1].replace('!',' '), linevals[2]]
+            if len(linevals) != 3:
+                print errormsg
+                print '  ' + fname + ": wrong values in file '" + projfilen + "' !!"
+                print '    each line has to have three words (space separated values)'
+                print "    from line: '" + line + "'"
+                print "    values found:", linevals
+                print "      remember to use '!' for spaces"
+                quit(-1)
+            pattr[linevals[0]] = [linevals[1].replace('!',' '), linevals[2]]
     opinff.close()
 
@@ -22858 +22867 @@
               tunitsval='minutes')
             oaxisv = tvals
-     
+
         # CF values of axis
         CFdimvalues = gen.CFcorValues(axn)
@@ -22894 +22903 @@
                     CFvardimvalues.append(dn)
                     varslice.append(slice(0,len(onc.dimensions[dn])))
-                    if dn == dimn and CFdimvalues['length'] != -1: 
+                    if dn == dimn and CFdimvalues['length'] != -1:
                         CFdimvalues['length'] = len(onc.dimensions[dn])
             axsiv = oaxisv[tuple(varslice)]
@@ -22927 +22936 @@
               'Y-m-dTH:M:SZ')
             secfreq = gen.CFtime_freq(axisv[0], axisv[1], tunits)
-            Sfreq = -1
             # Labelling frequency
-            if np.mod(secfreq, 100 * 365 * 24 * 3600.) == 0.:
-                Sfreq = str(int(secfreq/(100 * 365 * 24 * 3600.))) + 'C'
-            elif Sfreq == -1 and np.mod(secfreq, 365 * 24 * 3600.) == 0.: 
-                Sfreq = str(int(secfreq/(365 * 24 * 3600.))) + 'Y'
-            elif Sfreq == -1 and np.mod(secfreq, 7 * 24 * 3600.) == 0.: 
-                Sfreq = str(int(secfreq/(7 * 24 * 3600.))) + 'W'
-            elif Sfreq == -1 and np.mod(secfreq, 24 * 3600.) == 0.:
-                Sfreq = str(int(secfreq/(24.*3600.))) + 'D'
-            elif Sfreq == -1 and np.mod(secfreq,3600.) == 0.: 
-                Sfreq = str(int(secfreq/3600.)) + 'H'
-            elif Sfreq == -1 and np.mod(secfreq,60.) == 0.:
-                Sfreq = str(int(secfreq/60.)) + 'M'
-            else: 
-                # Getting frequency from the global attributes
+            if gattr['frequency'] == '[unknown]':
+                if np.mod(secfreq, 100 * 365 * 24 * 3600.) == 0.:
+                    Sfreq = str(int(secfreq/(100 * 365 * 24 * 3600.))) + 'C'
+                elif np.mod(secfreq, 365 * 24 * 3600.) == 0.: 
+                    Sfreq = str(int(secfreq/(365 * 24 * 3600.))) + 'Y'
+                elif np.mod(secfreq, 7 * 24 * 3600.) == 0.: 
+                    Sfreq = str(int(secfreq/(7 * 24 * 3600.))) + 'W'
+                elif np.mod(secfreq, 24 * 3600.) == 0.:
+                    Sfreq = str(int(secfreq/(24.*3600.))) + 'D'
+                elif np.mod(secfreq,3600.) == 0.: 
+                    Sfreq = str(int(secfreq/3600.)) + 'H'
+                elif np.mod(secfreq,60.) == 0.:
+                    Sfreq = str(int(secfreq/60.)) + 'M'
+                else:
+                    print errormsg
+                    print '  ' + fname + ": no frequency has been provided '" +      \
+                      gattr['frequency'] + "' but, data from file does not conform "+\
+                      "any know time-units!!"
+                    print '    frequency from file:', secfreq, ' seconds _______'
+                    print '      century:', secfreq/(100 * 365 * 24 * 3600.)
+                    print '      year (365d):', secfreq/(365 * 24 * 3600.)
+                    print '      week:', secfreq/(7 * 24 * 3600.)
+                    print '      day:', secfreq/(24 * 3600.)
+                    print '      hour:', secfreq/(3600.)
+                    print '      minute:', secfreq/(60.)
+                    print '      second:', secfreq
+                    print "  provide a frequency value in the global attributes " + \
+                      "file '" + globalattrfilen + "' as frequency [value]"
+                    quit(-1)
+                gattr['frequency'] = Sfreq
+            else:
                 Sfreq = gattr['frequency']
-            if gattr['frequency'] == '[unknown]': gattr['frequency'] = Sfreq
-            print '    frquency in file:', secfreq, 's', Sfreq
 
     # Renaming dimensions for 2D lon,lat variable-dimensions
@@ -22978 +23001 @@
                 axinf = axisinf[axn1]
                 cfaxinf = CFdimvals[axn1]
-                if axinf[0] == fdn: 
+                if axinf[0] == fdn:
                     cfdimaxes.append(cfaxinf['dimn'])
                     if filecfdimn[fdn] is None: filecfdimn[fdn] = cfaxinf['dimn']
@@ -22990 +23013 @@
 
         CFaxisvardimvals[axn] = cfdimaxes
- 
+
     print 'file --> CF equivalencies of dimensions _______'
     print filecfdimn
@@ -23036 +23059 @@
           '_' + gattr['rcm_version_id'] + '_' + Sfreq + '_' + Sinit + '-' + Sendt +  \
           '.nc'
+
         onewnc = NetCDFFile(filen,'w')
         now = dt.datetime.now()