Changeset 1054 in lmdz_wrf for trunk

Timestamp:

Aug 27, 2016, 9:40:33 AM (9 years ago)

Author:

lfita

Message:

Adding `CF' file in 'pinterp' (not tested!)

File:

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

trunk/tools/nc_var_tools.py

@@ -16220 +16220 @@
           'lat', 'time_counter']
         MODvardims = ['lon', 'lat', 'time_counter']
+    else:
+        modname = 'CF'
+        print warnmsg
+        print '  ' + fname + ": gessing that file '" + ncfile + "' comes from '" +   \
+          modname + "' !!"
+
+        # Variables to interpolate
+        MODdims = ['time', 'pres', 'lat', 'lon']
+        newMODdims = ['time', 'pres', 'lat', 'lon']
+        notCHK = []
+        MODvarrequired = ['p', 'ps', 'z', 'orog', 'ta', 'hus', 'lon', 'lat', 'time']
+        MODvardims = ['lon', 'lat', 'time']
 
     varns = gen.str_list(variables, ',')
@@ -16274 +16286 @@
         CFdimvs = {'time': dimt, 'presnivs': dimz, 'lat': dimy, 'lon': dimx}
         unstaggerDIM = 'lon'
+    else:
+        ovar1 = onc.variables['p']
+        pres = ovar1[:].astype('float64')
+        dimx = pres.shape[3]
+        dimy = pres.shape[2]
+        dimz = pres.shape[1]
+        dimt = pres.shape[0]
+        MODdimvs = {'time':dimt, 'pres':dimz, 'lat':dimy, 'lon': dimx}
+        CFdimvs = {'time': dimt, 'pres': dimz, 'lat': dimy, 'lon': dimx}
+        unstaggerDIM = 'lon'
 
     # sfc pressure
@@ -16281 +16303 @@
     elif modname == 'LMDZ':
         ovar1 = onc.variables['psol']
+        psfc = ovar1[:].astype('float64')
+    else:
+        ovar1 = onc.variables['ps']
         psfc = ovar1[:].astype('float64')
 
@@ -16295 +16320 @@
         ovar1 = onc.variables['geop']
         geop = ovar1[:].astype('float64')
+    else:
+        ovar1 = onc.variables['z']
+        geop = ovar1[:].astype('float64')
 
     # terrain height
@@ -16304 +16332 @@
         ovar1 = onc.variables['phis']
         hgt = (ovar1[0,:,:]/grav).astype('float64')
+    else:
+        ovar1 = onc.variables['orog']
+        hgt = ovar1[:].astype('float64')
 
     # water vapour mixing ratio
@@ -16311 +16342 @@
     elif modname == 'LMDZ':
         ovar1 = onc.variables['ovap']
+        qv = ovar1[:].astype('float64')
+    else:
+        ovar1 = onc.variables['hus']
         qv = ovar1[:].astype('float64')
 
@@ -16323 +16357 @@
     elif modname == 'LMDZ':
         ovar1 = onc.variables['temp']
+        temp = ovar1[:].astype('float64')
+    else:
+        ovar1 = onc.variables['ta']
         temp = ovar1[:].astype('float64')
 
@@ -16337 +16374 @@
         newdim = onewnc.createDimension('pres', len(interplevs))
         newdim = onewnc.createDimension('time_counter', None)
+    else:
+        newdim = onewnc.createDimension('lon', dimx)
+        newdim = onewnc.createDimension('lat', dimy)
+        newdim = onewnc.createDimension('pres', len(interplevs))
+        newdim = onewnc.createDimension('time', None)
 
 # Creation of variable dimensions