Changeset 2461 in lmdz_wrf

Timestamp:

Apr 23, 2019, 3:26:30 PM (6 years ago)

Author:

lfita

Message:

Adding:

• `juliandate': Function to provide the julian date (in days since -4712 January 1st 00 UTC) for any greogiran date

File:

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

trunk/tools/generic_tools.py

@@ -138 +138 @@
 # IsNumber: Function to determine if a given String is a given type of number
 # join_NOstrlist: Function to join a non-string list of values
-# julday_360d: Function to provide the julian day of a date in a 360 days/yr (or 12 30-days months) calendar
+# julday_360d: Function to provide the julian day of a date in a 360 days/yr (or 12 
+#   30-days months) calendar
+# juliandate: Function to provide the julian date (in days since -4712 January 1st 00 UTC) 
+#   for any greogiran date
 # latex_fig_array: Function to add an array of figures to an existing tex file
 # latex_text: Function to transform a text to LaTeX following style rules
@@ -15842 +15845 @@
 
 def juliandate(yr,mo,dd,hr,mi,ss):
-    """ Function to provide the julian date (in days) for a date between 1801 and 2099
+    """ Function to provide the julian date (in days since -4712 January 1st 00 UTC) 
+          for any greogiran date
       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(1877,8,11,7,30,0)
+    2406843.3125
+    >>> juliandate(877,8,11,7,30,0)
+    2041601.3125
+    >>> juliandate(1970,1,1,0,0,0)
+    2440588.0
     >>> juliandate(1978,1,1,0,0,0)
-
+    2443510.0
     >>> juliandate(1978,7,21,15,0,0)
-
+    2443712.625
     """
     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.)
+    sec1 = dd - 32075.
+    sec2 = 1461.*(yr+4800.+(mo-14)/12.)/4.
+    sec3 = 367.*(mo-2-(mo-14)/12.*12)/12.
+    sec4 = 3.*((yr+4900.+(mo-14.)/12.)/100.)/4.
+ 
+    juliand = int(sec1 + sec2 + sec3 - sec4)
+
+    juliand = juliand*1. + (hr + mi/60. + ss/3600.)/24.
+
+    # Alternative (without full agreement, just kept in case)
+    #K = yr + 0
+    #M = mo + 0
+    #I = dd + 0
+    #UT = hr*1. + mi/60. + ss/3600.
+
+    #sec1 = 367.*K
+    #sec2 = int((7.*(K+int((M+9)/12.)))/4.)
+    #sec3 = int((275.*M)/9.)
+    #sec4 = I*1.
+    #sec5 = 1721013.5
+    #sec6 = UT/24.
+    #sec7 = 0.5*sign(100.*K+M*1.-190002.5)
+    #sec8 = 0.5
+
+    #juliand = sec1 - sec2 + sec3 + sec4 + sec5 + sec6 - sec7 + sec8
+    #print '  juliand:', juliand
 
     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,                                      \