Changeset 903 in lmdz_wrf for trunk/tools/nc_var_tools.py

Timestamp:

Jun 19, 2016, 6:56:35 PM (8 years ago)

Author:

lfita

Message:

Adding `soil-layers' dimension and variable dimension in 'ModelChar?'

File:

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

trunk/tools/nc_var_tools.py

@@ -2152 +2152 @@
                 self.sname = var.getncattr('standard_name')
             else:
-#                print '    variable_inf.classpy: variable "' + self.name + '" does not have attribute "standard_name"'
                 self.sname = None
 
@@ -2158 +2157 @@
                 self.lname = var.getncattr('long_name')
             else:
-#                print '    variable_inf.classpy: variable "' + self.name + '" does not have attribute "long_name"'
                 self.lname = None
 
@@ -2164 +2162 @@
                 self.coor = var.getncattr('coordinates')
             else:
-#                print '    variable_inf.classpy: variable "' + self.name + '" does not have attribute "coordinates"'
                 self.coor = None
 
@@ -2170 +2167 @@
                 self.units = var.getncattr('units')
             else:
-#                print '    variable_inf.classpy: variable "' + self.name + '" does not have attribute "units"'
                 self.units = None
 
@@ -2176 +2172 @@
                 self.FillValue = var.getncattr('_FillValue')
             else:
-#                print '    variable_inf.classpy: variable "' + self.name + '" does not have attribute "_FillValue"'
                 self.FillValue = None
              
@@ -15745 +15740 @@
             self.dimzn = 'bottom_top'
             self.dimtn = 'Time'
+            self.dimsn = 'soil_layers_stag'
             self.vardxn = 'XLONG'
             self.vardyn = 'XLAT'
             self.vardzn = 'ZNU'
             self.vardtn = 'Times'
+            self.vardsn = 'ZS'
             if ncfile is not None and ncfile[0:2] == 'geo':
                 self.vardxn = 'XLONG_M'
@@ -15800 +15797 @@
       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 [hPa] (decreassing in hPa)
+        [interplevs]: ':' separated list of pressure values to interpolate [Pa] (decreassing in Pa)
         [linloginterp]: kind of linear interpolation
           1: linear interp in pressure
@@ -15952 +15949 @@
             elif vn == 'WRFrh': 
                 # relative humidity
-                ovar1 = 10.*0.6112*exp(17.67*(temp-273.16)/(temp-29.65))
+                ovar1 = 10.*0.6112*np.exp(17.67*(temp-273.16)/(temp-29.65))
                 ovar2 = 0.622*ovar1/(0.01*pres-(1.-0.622)*ovar1)
                 varin = 100.*qv/ovar2
@@ -16016 +16013 @@
     onewnc.setncattr('version', '1.0')
     onewnc.setncattr('oringal', 'Interpolation using p_interp.F90 subroutines')
-    newattr = set_atributek(onewnc, 'original_subroutines_author', unicode('Cindy ' +\
+    newattr = set_attributek(onewnc, 'original_subroutines_author', unicode('Cindy '+\
       'Bruy' + unichr(232) + 're'), 'U')
     onewnc.setncattr('original_subroutines_date', 'November 2007')