Changeset 2446 in lmdz_wrf

Timestamp:

Apr 17, 2019, 5:10:04 PM (6 years ago)

Author:

lfita

Message:

Fixing `impose_generic' when yr < yrref

File:

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

trunk/tools/generic_tools.py

@@ -1902 +1902 @@
     return newdate
 
-
 def impose_gregorian(timev, Tunits, cal):
     """ Function to impose gregorian calendar to a series of times with a 
@@ -1912 +1911 @@
         'noleap')
       [ 800. 901.]
+      >>> impose_gregorian(-7277.460938, 'days since 1949-12-01 00:00:00', 'noleap')
+      [-7272.460938]
     """
     import datetime as dt
@@ -1978 +1979 @@
             datetime = refdate + deltaT
             yr = datetime.year
-            Nleapd = leapdays(yrref, yr)
+            if yr > yrref:
+                Nleapd = leapdays(yrref, yr)
+            else:
+                Nleapd = leapdays(yr, yrref)
             if datetime.month <= 2 and days_month(yr, 2) > 28: Nleapd = Nleapd - 1
             newtimes[it] = timevsecs[it] + Nleapd*24*3600.
@@ -2051 +2055 @@
 
     return realdate
+
+import datetime as dt
+print timeref_datetime(dt.datetime(1949,12,1,0,0,0), -7277.460938, 'days', cal='noleap')
+quit()
 
 def timeref_datetime_mat(refd, timeval, tu):