Changeset 1855 in lmdz_wrf

Timestamp:

Mar 25, 2018, 6:56:17 PM (7 years ago)

Author:

lfita

Message:

Working version with `time_bnds' of 'CFmorziation'

Location:

trunk/tools

Files:

• nc_var.py (modified) (1 diff)
• nc_var_tools.py (modified) (6 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|Time|WRFtime,Z|bottom_top|ZNU:GlobalAttr1995.inf:proj1995.inf' -f ~/PY/wrfout_d01_1995-01-01_00\:00\:00 -v QFX
+## e.g. # nc_var.py -o nc_var.py -o CFmorzization -S 'X|west_east|XLONG,Y|south_north|XLAT,T|Time|WRFtime,Z|bottom_top|ZNU:GlobalAttr1995.inf:proj1995.inf' -f ~/PY/wrfout_d01_1995-01-01_00\:00\:00 -v QFX@instantaneous@None
 
 ## 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

@@ -22618 +22618 @@
         quit(-1)
 
-    Nbnds = int(dFCt / dtbnds)
+    Nbnds = int((dFCt + dtbnds)/ dtbnds)
     if Nbnds < 2:
         print errormsg
@@ -22660 +22660 @@
 #otime = onc.variables['Times']
 #wtimev = otime[:]
-#tbnds, tc, tu = compute_WRFtime_bnds(wtimev, 'i,1', kindWRFt='begperiod',            \
+#tbnds, tc, tu = compute_WRFtime_bnds(wtimev, 'd,1', kindWRFt='begperiod',            \
 #  refdate='19491201000000', tunitsval='minutes')
 
@@ -22727 +22727 @@
                proj_units [value] (units of the resoltion of the projection)
         ncfile= netCDF file to use
-        variable= ',' list of variables to CF transform ('all' for all variables). A separated file will be created
-          for each variable
+        variable= ',' list of [varfn]@[method]@[Tbnds] to CF transform. A separated file will be 
+          created for each variable
+          use [varfn]@instantaneous@None for instantaneous values
           NOTE: CF values will be taken from 'variables_values.dat'
+            [varfn]: name of the variable inside the file
+            [method]: cell_method of the variable ('!' for spaces) uing, standard CF description: 
+                http://cfconventions.org/cf-conventions/cf-conventions.html#_data_representative_of_cells
+            [var_time_bnds]: variable with the time_bnds values ('WRFtime_bnds'|[tunit]|[quantity], to be computed
+              directly from WRF)
+                'WRFtime_bnds'|[tunit]|[quantity] period to compute the bounds as:
+                  [tunit]: unit of time: 'c' century, 'y' year, 'm' month, 'w' week, 'd' day, 
+                    'h' hour, 'i' minute, 's' second
+                  [quantity] = amount (integer) of [tunit] to cover a period of time_bnds
     """
     import datetime as dt
@@ -23049 +23059 @@
         print '  ', axn
         print '  ', CFdimvals[axn]
-    print 'Values for variables-axes ________'
+    print 'Dimensions for variables-axes ________'
     for axn in axes:
         print '  ', axn, ':', CFaxisvardimvals[axn]
@@ -23055 +23065 @@
     # Processing variables in file
     ##
-    for vn in Varns:
+    for vtmn in Varns:
+        vn = vtmn.split('@')[0]
+        vm = vtmn.split('@')[1].replace('!', ' ')
+        vtb = vtmn.split('@')[2]
+
         # getting CF information
         varinf = gen.variables_values(vn)
@@ -23163 +23177 @@
             onewnc.sync()
 
+            # time_bnds
+            if vm != 'instantaneous':
+                print infmsg
+                print '  ' + fname + ": Adding 'time_bnds' variable ..."
+                set_attribute(newvar, 'cell_methods', vm)
+                ovtime = onewnc.variables['time']
+                set_attribute(ovtime, 'bounds', 'time_bnds')
+                if vtb[0:12] == 'WRFtime_bnds':
+                    vtbvals = vtb.split('|')
+                    if len(vtbvals) != 3:
+                        print errormsg
+                        print '  ' + fname + ": to compute 'time_bnds' from WRF " +  \
+                          "'Times', one needs to provide the period of the " +       \
+                          "bounds as: 'WRFtime_bnds'|[tunit]|[quantity]"
+                        print "    values passed: '" + vtb
+                        onc.close()
+                        quit(-1)
+                    period = vtbvals[1] + ',' + vtbvals[2]
+                    print infmsg
+                    print '      ' + fname + ": creation of variable 'time_bnds' " + \
+                      "from WRF 'Times' using period '" + period + "' !!"
+                    odimvar = onc.variables['Times']
+                    timewrfv = odimvar[:]
+                    refT = gen.datetimeStr_conversion(gattr['basetime'],             \
+                      'Y-m-dTH:M:SZ', 'YmdHMS')
+                    tbvals, tvals, urefvals = compute_WRFtime_bnds(timewrfv, period, \
+                      kindWRFt='begperiod', refdate=refT, tunitsval='minutes')
+                else:
+                    if not onc.variables.has_key(vardimn):
+                        print errormsg
+                        print '  '+fname+"': netcdf file '"+ncfile + "' does not " + \
+                          "contain variable '" + vardimn+ "' with time_bnds values !!"
+                        print '    available ones:', oncvars
+                        onc.close()
+                        quit(-1)
+                    tbvals = onc.variables[vardimn]
+                # Re-setting time values to the center of the bounds
+                ovtime[:] = tvals[:]
+                onewnc.sync()
+
+                # creation of time_bnds variable
+                if not gen.searchInlist(onewnc.dimensions, 'nv'):
+                    onewnc.createDimension('nv', 2)
+                newvar = onewnc.createVariable('time_bnds', 'f8', ('time', 'nv'))
+                newvar[:] = tbvals[:]
+
         # Adding projection variable
         print infmsg