Changeset 1810 in lmdz_wrf

Timestamp:

Mar 16, 2018, 9:11:22 PM (7 years ago)

Author:

lfita

Message:

Adding WRFw' and WRFhus' on 'pinterp'

File:

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

trunk/tools/nc_var_tools.py

@@ -18051 +18051 @@
 def pinterp(values, ncfile, variables):
     """ Function to vertically interpolate using subroutines from the p_interp.F90 NCAR program
-      Using fortran codes: module_generic.F90
-      foudre: f2py -m module_ForInterpolate --f90exec=/usr/bin/gfortran-4.7 -c module_generic.F90 module_ForInterpolate.F90  >& run_f2py.log
-      ciclad: f2py --f90flags="-fPIC" -L/opt/canopy-1.3.0/Canopy_64bit/System/lib/ -L/usr/lib64/ -L/opt/canopy-1.3.0/Canopy_64bit/System/lib/ -m module_ForInterpolate -c module_generic.F90 module_ForInterpolate.F90 
       values= [interplevs],[linloginterp],[extrap]
         [interplevs]: ':' separated list of pressure values to interpolate [Pa] (decreassing in Pa)
@@ -18063 +18060 @@
       variables = ',' separated list of names of variables to interpolate ('all', fo all 4D-atmospheric variables)
         'WRFght': for WRF geopotential height
+        'WRFhus': for WRF specific humidity
         'WRFrh': for WRF relative humidity
         'WRFt': for WRF temperature
         'WRFu': for WRF x-wind de-staggered
         'WRFv': for WRF y-wind de-staggered
+        'WRFw': for WRF z-wind de-staggered
     """
     import module_ForInt as fin
@@ -18352 +18351 @@
                 isgeop = True
                 varattrs = gen.variables_values('WRFght')
+                CFvn = varattrs[0]
+                newvarattr['standard_name'] = varattrs[1]
+                newvarattr['long_name'] = varattrs[4].replace('|',' ')
+                newvarattr['units'] = varattrs[5]
+            elif vn == 'WRFhus': 
+                # relative humidity
+                varin = qv/(1.+qv)
+                isgeop = False
+                varattrs = gen.variables_values('hus')
                 CFvn = varattrs[0]
                 newvarattr['standard_name'] = varattrs[1]
@@ -18401 +18409 @@
                 newvarattr['long_name'] = varattrs[4].replace('|',' ')
                 newvarattr['units'] = varattrs[5]
+            elif vn == 'WRFw':
+                ovarin = onc.variables['W']
+                print infmsg
+                print '  ' + fname + ': De-staggering z-wind variable !!'
+                print '    from:', ovarin.shape, 'to', (dimt, dimz, dimy, dimx)
+                varin = np.zeros((dimt, dimz, dimy, dimx), dtype=np.float64)
+                varin[:] = 0.5*(ovarin[:,0:dimz,:,:] + ovarin[:,1:dimz+1,:,:])
+                isgeop = False
+                varattrs = gen.variables_values('wa')
+                CFvn = varattrs[0]
+                newvarattr['standard_name'] = varattrs[1]
+                newvarattr['long_name'] = varattrs[4].replace('|',' ')
+                newvarattr['units'] = varattrs[5]
             elif vn == 'pres':
                 varin = pres