Changeset 1770 in lmdz_wrf

Timestamp:

Feb 9, 2018, 6:58:18 PM (7 years ago)

Author:

lfita

Message:

Adding functions to transform from angle <==> Deg, Minute, Second

File:

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

trunk/tools/generic_tools.py

@@ -49 +49 @@
 
 ####### Content
+# angle_DegMinSec: Function to transform an angle to Degrees Minutes Seconds
 # ASCII_LaTeX: Function to transform from an ASCII character to LaTeX codification
 # ASCII_to: Function to provide the equivalence text to a given lenguage omitting certain characters
@@ -73 +74 @@
 # DateTimeStr_date: Function to transform a string date-time ([YYYY][MM][DD][HH][MI][SS] format) to a date object
 # datetimeStr_datetime: Function to transform a string date-time ([YYYY]-[MM]-[DD]_[HH]:[MI]:[SS] format) to a date object
+# DegMinSec_angle: Function to transform Degrees Minutes Seconds to an angle
 # dictionary_key: Function to provide the first key in a dictionay with a given value
 # dictionary_key_list: Function to provide the first key in a dictionary of lists with a given value
@@ -12610 +12612 @@
     return values
 
+def angle_DegMinSec(angle):
+   """ Function to transform an angle to Degrees Minutes Seconds
+     angle= decimal angle to transform
+     >>> angle_DegMinSec(3.98765)
+     3 59 15.54
+     >>> angle_DegMinSec(-23.0025)
+     -23 0 9.0                         
+   """
+   fname = 'angle_DegMinSec'
+                                      
+   degv = int(abs(angle))
+   minv = int((abs(angle) - degv*1.)*60.)
+   secv = (abs(angle) - degv*1. - minv/60.)*3600.
+                                                  
+   if angle < 0.:
+       return -degv, minv, secv
+   else:
+       return degv, minv, secv
+
+def DegMinSec_angle(degv, minv, secv):
+   """ Function to transform Degrees Minutes Seconds to an angle
+     degv: value in degrees
+     minv: minutes of arc
+     secv: seconds of arc
+     >>> DegMinSec_angle(-23, 0, 9.0)
+     -23.0025
+   """
+   fname = 'DegMinSec_angle'
+                                      
+   angle = np.abs(degv) + np.abs(minv)/60. + np.abs(secv)/3600.
+                                                  
+   if degv < 0.:
+       return -angle
+   else:
+       return angle
+
+d, m, s = angle_DegMinSec(-23.0025)
+print d, m, s
+print DegMinSec_angle(d,m,s)
+
 #quit()
 
-