Changeset 1848 in lmdz_wrf

Timestamp:

Mar 23, 2018, 2:48:42 PM (7 years ago)

Author:

lfita

Message:

Getting almost done!! with `CFmorzization'

File:

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

trunk/tools/nc_var_tools.py

@@ -22690 +22690 @@
           [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 '#' for comments (use '[unknown]' for that values set-up according to content from file:
               institute_id [value, use '!' for spaces] (CV)
               institution  [value]
@@ -22699 +22699 @@
               product [value] (single value: "output")
               Conventions [value] (later than CF1.4 if present / single value: "CF-1.4")
-              creation_date [YYYY-MM-DDTHH:MM:SSZ] (format according to UNIDATA conventions, same in CMIP5 & CORDEX)
+              creation_date [value] ([YYYY-MM-DDTHH:MM:SSZ] format according to UNIDATA conventions, same in CMIP5 & CORDEX)
               frequency [value] (same CV as CMIP5 with "sem" in addition; 3hr highest freq. in CORDEX; no "yr")
               driving_model_id [value] (CMIP5 institute_id, <CMIP5 institute_id>-<member in CMIP5 CV list of model_id>)
@@ -22726 +22726 @@
           NOTE: CF values will be taken from 'variables_values.dat'
     """
+    import datetime as dt
     fname = 'CFmorzization'
 
@@ -22926 +22927 @@
               '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 np.mod(secfreq, 365 * 24 * 3600.) == 0.: 
+            elif Sfreq == -1 and np.mod(secfreq, 365 * 24 * 3600.) == 0.: 
                 Sfreq = str(int(secfreq/(365 * 24 * 3600.))) + 'Y'
-            elif np.mod(secfreq, 7 * 24 * 3600.) == 0.: 
+            elif Sfreq == -1 and np.mod(secfreq, 7 * 24 * 3600.) == 0.: 
                 Sfreq = str(int(secfreq/(7 * 24 * 3600.))) + 'W'
-            elif np.mod(secfreq, 24 * 3600.) == 0.:
+            elif Sfreq == -1 and np.mod(secfreq, 24 * 3600.) == 0.:
                 Sfreq = str(int(secfreq/(24.*3600.))) + 'D'
-            elif np.mod(secfreq,3600.) == 0.: 
+            elif Sfreq == -1 and np.mod(secfreq,3600.) == 0.: 
                 Sfreq = str(int(secfreq/3600.)) + 'H'
-            elif np.mod(secfreq,60.) == 0.:
+            elif Sfreq == -1 and np.mod(secfreq,60.) == 0.:
                 Sfreq = str(int(secfreq/60.)) + 'M'
-            else: Sfreq = 'UNKNOWN'
+            else: 
+                # Getting frequency from the global attributes
+                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
+    if len(filedimvals['X'].shape) == 2 or len(filedimvals['Y'].shape) == 2:
+        print infmsg
+        print '    because Xvals have rank 2:', filedimvals['X'].shape,              \
+          'renaming lon,lat dimensions to x,y'
+        axv = CFdimvals['X']
+        axv['dimn'] = 'x'
+        CFdimvals['X'] = axv
+        axv = CFdimvals['Y']
+        axv['dimn'] = 'y'
+        CFdimvals['Y'] = axv
+        #axiv = axisinf['X']
+        #filecfdimn[axiv[0]] = 'x'
+        #axiv = axisinf['Y']
+        #filecfdimn[axiv[0]] = 'y'
 
     print 'Values for variables-axes from file ________'
@@ -22968 +22990 @@
 
         CFaxisvardimvals[axn] = cfdimaxes
+ 
     print 'file --> CF equivalencies of dimensions _______'
     print filecfdimn
@@ -23013 +23036 @@
           '_' + gattr['rcm_version_id'] + '_' + Sfreq + '_' + Sinit + '-' + Sendt +  \
           '.nc'
-
         onewnc = NetCDFFile(filen,'w')
+        now = dt.datetime.now()
+        if gattr['creation_date'] == '[unknown]': 
+            gattr['creation_date'] = now.strftime("%Y-%m-%dT%H:%M:%SZ")
         # Definition only of that coordinates from the variable
         for axn in varaxes: