Changeset 2460 in lmdz_wrf for trunk/tools/generic_tools.py

Timestamp:

Apr 23, 2019, 2:59:43 PM (6 years ago)

Author:

lfita

Message:

Ading:

• `sign': Function to provide the sign of a number

File:

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

trunk/tools/generic_tools.py

@@ -192 +192 @@
 # shrinkarray_dim: Function to shrink a given array with a vector of True/False along a given dimension
 # significant_decomposition: Function to decompose a given number by its signifcant potencies
+# sign: Function to provide the sign of a number
 # singleline_printing_class: Function to print all the values of a given class in a single line to be parseavel
 # stagger_unstagger: Function to de-stagger a variable
@@ -15826 +15827 @@
     return isleap
 
+def sign(number):
+    """ Function to provide the sign of a number
+      num: number to provide the sign
+    >>> sign(7109)
+    1
+    >>> sign(-3.43512)
+    -1.0
+    """
+    fname = 'sign'
+
+    sign = np.abs(number)/number
+
+    return sign
+
+def juliandate(yr,mo,dd,hr,mi,ss):
+    """ Function to provide the julian date (in days) for a date between 1801 and 2099
+      refdate: 1 January 4713 BC (= -4712 January 1), Greenwich mean noon (= 12h UT).
+    FROM: https://aa.usno.navy.mil/faq/docs/JD_Formula.php
+      after 1990 edition of the U.S. Naval Observatory's Almanac for Computers
+    >>> juliandate(1978,1,1,0,0,0)
+
+    >>> juliandate(1978,7,21,15,0,0)
+
+    """
+    fname = 'juliandate'
+
+    K = dd
+    I = yr
+    J = mo
+    UT = hr + mi/60. + ss/3600.
+
+    sec1 = 367.*K
+    sec2 = int((7.*(K+int((M+9)/12.)))/4)
+    sec3 = int((275.*M)/9.)
+    sec4 = I
+    sec5 = 1721013.5
+    sec6 = UT/24.
+    sec7 = 0.5*sign(100.*K+M-190002.5)
+
+    juliand = int(dd - 32075. + 1461.*(yr+4800.+(mo-14)/12.)/4. +               \
+      367.*(mo-2-(mo-14)/12.*12)/12. - 3*((yr+4900+(mo-14)/12.)/100.)/4.)
+
+    JD= K-32075+1461*(I+4800+(J-14)/12)/4+367*(J-2-(J-14)/12*12)                \
+      /12-3*((I+4900+(J-14)/12)/100)/4
+
+    print yr,mo,dd,hr,mi,ss,' JD:', JD
+
+#    juliand = juliand + hr/24. + mi/(24.*60.) + ss/(24.*3600.)
+
+    return juliand
+print juliandate(1877,8,11,7,30,0)
+
+print juliandate(1970,1,1,0,0,0)
+
+print juliandate(1978,1,1,0,0,0)
+
+print juliandate(1978,7,21,15,0,0)
+quit()
+
 def fix_CFdates(timevals, origcftimeu, origcal,                                      \
   newcftimeu= 'hours since 1949-12-01 00:00:00', newcal='gregorian'):