Changeset 1846 in lmdz_wrf

Timestamp:

Mar 23, 2018, 12:57:53 AM (7 years ago)

Author:

lfita

Message:

Working on 'CFmorzization'

Location:

trunk/tools

Files:

• nc_var.py (modified) (1 diff)
• nc_var_tools.py (modified) (9 diffs)

trunk/tools/nc_var.py

@@ -18 +18 @@
 ## e.g. # nc_var.py -o WRF_to_newCF -f wrfout_d01_1995-01-01_00\:00\:00 -S 'XLONG:XLAT:Times:19491201000000:minutes' -v QVAPOR,T2,Q2
 ## e.g. # nc_var.py -o reconstruct_matrix_from_vector -f cruncep_halfdeg_1958.nc -S 'nav_lon:nav_lat:-90.:-7.:-67.:15.:latlon:0.5:0.5:0.05' -v all 
-## e.g. # nc_var.py -o CFmorzization -S 'X|west_east|XLONG,Y|south_north|XLAT,T|Times|WRFtime:GlobalAttr1995.inf:proj1995.inf' -f ~/PY/wrfout_d01_1995-01-01_00\:00\:00 -v QFX
+## e.g. # nc_var.py -o CFmorzization -S 'X|west_east|XLONG,Y|south_north|XLAT,T|Times|WRFtime,Z|bottom_top|ZNU:GlobalAttr1995.inf:proj1995.inf' -f ~/PY/wrfout_d01_1995-01-01_00\:00\:00 -v QFX
 
 ## e.g. ccrc468-17 # ./nc_var.py -v time -f 123/CCRC_NARCliM_Sydney_All_1990-1999_pr10max.nc -o out -S 1:-1

trunk/tools/nc_var_tools.py

@@ -22707 +22707 @@
                 members of the CMIP5 list of forcing agents are not necessarily forcing agents of RCM (besides implicitly
                 in driving_experimnt)
-              rcm_version_ID [value] (<free text string>; valid characters only indicates model modifs during the project
+              rcm_version_id [value] (<free text string>; valid characters only indicates model modifs during the project
                 e.g.: parameterizations, small upgrades)
               project_id [value] (single value)
@@ -22744 +22744 @@
       'experiment', 'contact', 'product', 'Conventions', 'creation_date','frequency',\
       'driving_model_id', 'driving_model_ensemble_member', 'driving_experiment_name',\
-      'driving_experiment', 'rcm_version_ID', 'project_id', 'CORDEX_domain',         \
+      'driving_experiment', 'rcm_version_id', 'project_id', 'CORDEX_domain',         \
       'tracking_id', 'basetime', 'calendar', 'grid']
 
@@ -22855 +22855 @@
             oaxisv = tvals
      
-
         # CF values of axis
         CFdimvalues = gen.CFcorValues(axn)
@@ -22876 +22875 @@
             axisv = oaxisv[:]
             if CFdimvalues['length'] != -1:CFdimvalues['length'] = oaxisv.shape[0]
-            CFvardimvalues.append[dimn]
+            CFvardimvalues.append(dimn)
         elif len(oaxisv.shape) == CFdimvalues['maxrank'] + 1:
             print warnmsg
@@ -22889 +22888 @@
                 if dn == dimTn: varslice.append(0)
                 else: 
-                    CFvardimvalues.append(dn)
+                    CFvardimvalues.append(dimn)
                     varslice.append(slice(0,len(onc.dimensions[dn])))
                     if dn == dimn and CFdimvalues['length'] != -1: 
@@ -22912 +22911 @@
             dimt = len(axisv)
             CFitime = axisv[0]
-            CFetime = axisv[dimt]
-            tunits = oaxisv.units
-            tcal = oaxisv.calendar
-            Sinit=datetimeStr_conversion(str(CFitime),'cftime,'+tunits,'Y-m-dTH:M:SZ')
-            Sendt=datetimeStr_conversion(str(CFetime),'cftime,'+tunits,'Y-m-dTH:M:SZ')
-            secfreq = CFtime_freq(axisv[0], axisv[1], tunits)
+            CFetime = axisv[dimt-1]
+            if vardimn != 'WRFtime': 
+                tunits = oaxisv.units
+                tcal = oaxisv.calendar
+            else: 
+                tunits = urefvals
+                tcal = 'standard'
+            Sinit = gen.datetimeStr_conversion(str(CFitime), 'cfTime,'+tunits,       \
+              'Y-m-dTH:M:SZ')
+            Sendt = gen.datetimeStr_conversion(str(CFetime), 'cfTime,'+tunits,       \
+              'Y-m-dTH:M:SZ')
+            secfreq = gen.CFtime_freq(axisv[0], axisv[1], tunits)
             # Labelling frequency
             if np.mod(secfreq, 100 * 365 * 24 * 3600.) == 0.:
@@ -22933 +22938 @@
             else: Sfreq = 'UNKNOWN'
 
+    print 'Values for axes ________'
+    print CFdimvals
+    print 'Values for variables-axes ________'
+    print CFaxisvardimvals
+
     # Processing variables in file
     ##
-    for vn in varns:
+    for vn in Varns:
         # getting CF information
         varinf = gen.variables_values(vn)
@@ -22958 +22968 @@
             axisv = axisinf[axn]
             dimn = axisv[0]
-            if gen.searchInfile(ovar.dimensions, dimn): 
-                varaxes.append[axn]
+            if gen.searchInlist(ovar.dimensions, dimn): 
+                varaxes.append(axn)
                 CFvardims.append(CFdimvals[axn]['dimn'])
                 
-        filen = gattr['institude_id'] + '_' + cfvarn + '_' + gattr['CORDEX_domain'] +\
+        filen = gattr['institute_id'] + '_' + cfvarn + '_' + gattr['CORDEX_domain'] +\
           '_' + gattr['driving_model_id'] + '_' + gattr['driving_experiment_name'] + \
           '_' + gattr['driving_model_ensemble_member'] + '_' + gattr['model_id'] +   \
@@ -22972 +22982 @@
         for axn in varaxes:
             dvals = CFdimvals[axn]
-            newdim = onewnc.createDimension(dvals['dimn'], dvals['legnth'])
+            newdim = onewnc.createDimension(dvals['dimn'], dvals['length'])
 
         # Definition only of that variables-coordinates from the axes in the variable
         for axn in varaxes:
             dvals = CFdimvals[axn]
+            print 'Lluis: ', CFaxisvardimvals[axn]
             newvar = onewnc.createVariable(dvals['dimn'], 'f8',                      \
-              tuple(CFaxesdimvarvals[axn]))
+              tuple(CFaxisvardimvals[axn]))
             basicvardef(newvar, dvals['stdn'], dvals['longname'], dvals['units'])
             for ivn in dvals.keys():