Make sure doy uses integer arrays. Fix #123

TUW-GEO · Feb 6, 2018 · eaaa817 · eaaa817
1 parent a83d9c5
commit eaaa817
Showing 1 changed file with 3 additions and 2 deletions.
diff --git a/pytesmo/timedate/julian.py b/pytesmo/timedate/julian.py
@@ -303,8 +303,9 @@ def doy(month, day, year=None):
 
     if year is not None:
         nonleap_years = np.invert(is_leap_year(year))
-        day_of_year = day_of_year - nonleap_years + \
-            np.logical_and(day_of_year < 60, nonleap_years)
+        day_of_year = (day_of_year -
+                       nonleap_years.astype('int') +
+                       np.logical_and(day_of_year < 60, nonleap_years).astype('int'))
 
     return day_of_year