Changeset 564 in lmdz_wrf for trunk/tools

Timestamp:

Jul 6, 2015, 12:41:33 PM (10 years ago)

Author:

lfita

Message:

Adding `TStd', dew-point temperature from TS variables

File:

• trunk/tools/validation_sim.py (modified) (3 diffs)

trunk/tools/validation_sim.py

@@ -829 +829 @@
         'WRFT': CF-time from WRF variables
         'WRFt': temperature from WRF variables
+        'TStd': dew-point temperature from TS variables
         'WRFtd': dew-point temperature from WRF variables
         'WRFwd': wind direction from WRF variables
@@ -961 +962 @@
                 dimensions = ncobj.variables['T'].dimensions
                 shape = ncobj.variables['P'].shape
+
+            elif varn == 'TStd':
+#        print '    ' + main + ': computing dew-point temperature from TS as t and Tetens...'
+# tacking from: http://en.wikipedia.org/wiki/Dew_point
+                p=ncobj.variables['p'][:]
+
+                temp = ncobj.variables['t'][:]
+
+                qv = ncobj.variables['q'][:]
+
+                tk = temp
+                data1 = 10.*0.6112*np.exp(17.67*(tk-273.16)/(tk-29.65))
+                data2 = 0.622*data1/(0.01*p-(1.-0.622)*data1)
+
+                rh = qv/data2
+                
+                pa = rh * data1/100.
+                varNOcheckv = 257.44*np.log(pa/6.1121)/(18.678-np.log(pa/6.1121))
+
+                dimensions = ncobj.variables['t'].dimensions
+                shape = ncobj.variables['t'].shape
 
             elif varn == 'WRFtd':
@@ -1431 +1453 @@
   "'WRFp': pressure from WRF variables; " +                                          \
   "'WRFrh': relative humidty fom WRF variables; " +                                  \
+  "'TSrhs': surface relative humidity from TS variables; " +                         \
   "'WRFrhs': surface relative humidity from WRF variables; " +                       \
   "'WRFT': CF-time from WRF variables; " +                                           \
   "'WRFt': temperature from WRF variables; " +                                       \
+  "'TStd': dew-point temperature from TS variables; " +                              \
   "'WRFtd': dew-point temperature from WRF variables; " +                            \
   "'WRFwd': wind direction from WRF variables; " +                                   \
+  "'TSwds': surface wind direction from TS variables; " +                            \
   "'WRFwds': surface wind direction from WRF variables; " +                          \
   "'WRFws': wind speed from WRF variables; " +                                       \
+  "'TSwss': surface wind speed from TS variables; " +                                \
   "'WRFwss': surface wind speed from WRF variables; " +                              \
   "'WRFz': height from WRF variables"