Changeset 506 in lmdz_wrf for trunk/tools

Timestamp:

Jun 17, 2015, 10:44:42 AM (10 years ago)

Author:

lfita

Message:

Definning wind direction as [0,360)

File:

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

trunk/tools/validation_sim.py

@@ -996 +996 @@
 
                 theta = np.arctan2(ua, va)
+                theta = np.where(theta < 0., theta + 2.*np.pi, theta)
+                varNOcheckv = 360.*theta/(2.*np.pi)
+
                 dimensions = tuple(['Time','bottom_top','south_north','west_east'])
                 shape = ua.shape
-                varNOcheckv = 360.*theta/(2.*np.pi)
 
             elif varn == 'WRFwds':
 #        print '    ' + main + ': computing surface wind direction from WRF as ATAN2(V,U) ...'
                 theta = np.arctan2(ncobj.variables['V10'][:], ncobj.variables['U10'][:])
+                theta = np.where(theta < 0., theta + 2.*np.pi, theta)
                 
                 varNOcheckv = 360.*theta/(2.*np.pi)
@@ -1011 +1014 @@
 #        print '    ' + main + ': computing surface wind direction from TSs as ATAN2(v,u) ...'
                 theta = np.arctan2(ncobj.variables['v'][:], ncobj.variables['u'][:])
+                theta = np.where(theta < 0., theta + 2.*np.pi, theta)
 
                 varNOcheckv = 360.*theta/(2.*np.pi)