Changeset 2810 in lmdz_wrf

Timestamp:

Apr 8, 2020, 9:45:23 PM (5 years ago)

Author:

lfita

Message:

Adding:

• `truncate': Function to truncate a Float value to a certain decimal

File:

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

trunk/tools/generic_tools.py

@@ -256 +256 @@
 #   of time-steps for a whole period is multiple of a given number
 # timefmt_timelab: Function to transform from C-like time format to general one
+# truncate: Function to truncate a Float value to a certain decimal
 # unitsdsDate: Function to know how many units of time are from a given pair of dates
 # units_lunits: Function to provide LaTeX equivalences from a given units
@@ -14095 +14096 @@
     return values
 
+def truncate(fv, decimals=7):
+    """ Function to truncate a Float value to a certain decimal
+      FROM: https://realpython.com/python-rounding/#truncation
+      fv: float value
+      decimals: decimals to which truncate (7, default)
+    >>> truncate(np.pi)
+    3.1415926
+    >>> truncate(np.pi, decimals=3)
+    3.141
+    >>> truncate(np.pi, decimals=0)
+    3.0
+    >>> truncate(1./3., decimals=7)
+    0.3333333
+    """
+    fname = 'truncate'
+
+    multiplier = 10. ** decimals
+    tfv = int(fv * multiplier) / multiplier
+
+    return int(fv * multiplier) / multiplier
+
 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
+     3, 59, 15.54
      >>> angle_DegMinSec(-23.0025)
-     -23 0 9.0                         
+     -23, 0, 9.0                         
    """
    fname = 'angle_DegMinSec'
@@ -14108 +14130 @@
    minv = int((abs(angle) - degv*1.)*60.)
    secv = (abs(angle) - degv*1. - minv/60.)*3600.
+
+   if secv < 1.e-8: secv = 0.
                                                   
    if angle < 0.:
@@ -14113 +14137 @@
    else:
        return degv, minv, secv
+
+print angle_DegMinSec(49.1)
+print angle_DegMinSec(3.98765)
+quit()
 
 def DegMinSec_angle(degv, minv, secv):
@@ -17568 +17596 @@
     return newlist
 
-#quit()
-
-