Exemplo n.º 1
0
    def test_select_dummy(self):
        nutime = self.TestTime(datetime(1950, 1, 1), 366, 24, "hours since 1900-01-01", "standard")

        dates = [datetime(1950, 1, 2, 6), datetime(1950, 1, 3), datetime(1950, 1, 3, 18)]

        t = date2index(dates, nutime, select="before")
        assert_equal(t, [1, 2, 2])

        t = date2index(dates, nutime, select="after")
        assert_equal(t, [2, 2, 3])

        t = date2index(dates, nutime, select="nearest")
        assert_equal(t, [1, 2, 3])
Exemplo n.º 2
0
    def test_nonuniform(self):
        """Check that the fallback mechanism works. """
        nutime = self.TestTime(datetime(1950, 1, 1), 366, 24, "hours since 1900-01-01", "standard")

        # Let's remove the second entry, so that the computed stride is not
        # representative and the bisection method is needed.
        nutime._data = nutime._data[numpy.r_[0, slice(2, 200)]]

        t = date2index(datetime(1950, 2, 1), nutime)
        assert_equal(t, 30)

        t = date2index([datetime(1950, 2, 1), datetime(1950, 2, 3)], nutime)
        assert_equal(t, [30, 32])
Exemplo n.º 3
0
    def test_nonuniform(self):
        """Check that the fallback mechanism works. """
        nutime = self.TestTime(datetime(1950, 1, 1), 366, 24,
                               'hours since 1900-01-01', 'standard')

        # Let's remove the second entry, so that the computed stride is not
        # representative and the bisection method is needed.
        nutime._data = nutime._data[numpy.r_[0, slice(2, 200)]]

        t = date2index(datetime(1950, 2, 1), nutime)
        assert_equal(t, 30)

        t = date2index([datetime(1950, 2, 1), datetime(1950, 2, 3)], nutime)
        assert_equal(t, [30, 32])
Exemplo n.º 4
0
 def test_select_dummy(self):
     nutime = self.TestTime(datetime(1950, 1, 1), 366, 24,
       'hours since 1900-01-01', 'standard')
     
     dates = [datetime(1950,1,2,6), datetime(1950,1,3), datetime(1950,1,3,18)]
     
     t = date2index(dates, nutime, select='before')
     assert_equal(t, [1, 2, 2])
     
     t = date2index(dates, nutime, select='after')
     assert_equal(t, [2, 2, 3])
     
     t = date2index(dates, nutime, select='nearest')
     assert_equal(t, [1,2,3])
Exemplo n.º 5
0
 def test_select_dummy(self):
     nutime = self.TestTime(datetime(1950, 1, 1), 366, 24,
       'hours since 1900-01-01', 'standard')
     
     dates = [datetime(1950,1,2,6), datetime(1950,1,3), datetime(1950,1,3,18)]
     
     t = date2index(dates, nutime, select='before')
     assert_equal(t, [1, 2, 2])
     
     t = date2index(dates, nutime, select='after')
     assert_equal(t, [2, 2, 3])
     
     t = date2index(dates, nutime, select='nearest')
     assert_equal(t, [1,2,3])
Exemplo n.º 6
0
 def test_singletime(self):
     # issue 215 test (date2index with time variable length == 1)
     f = Dataset(self.file)
     time2 = f.variables['time2']
     result_index = date2index(self.first_timestamp,time2, select="exact")
     assert_equal(result_index, 0)
     f.close()
Exemplo n.º 7
0
 def test_singletime(self):
     # issue 215 test (date2index with time variable length == 1)
     f = Dataset(self.file)
     time2 = f.variables['time2']
     result_index = date2index(self.first_timestamp, time2, select="exact")
     assert_equal(result_index, 0)
     f.close()
Exemplo n.º 8
0
 def test_failure(self):
     nutime = self.TestTime(datetime(1950, 1, 1), 366, 24, "hours since 1900-01-01", "standard")
     try:
         t = date2index(datetime(1949, 2, 1), nutime)
     except ValueError:
         pass
     else:
         raise ValueError("This test should have failed.")
Exemplo n.º 9
0
 def test_select_nc(self):
     f = Dataset(self.file, 'r')
     nutime = f.variables['time']
     
     dates = [datetime(1950,1,2,6), datetime(1950,1,3), datetime(1950,1,3,18)]
     
     t = date2index(dates, nutime, select='before')
     assert_equal(t, [1, 2, 2])
     
     t = date2index(dates, nutime, select='after')
     assert_equal(t, [2, 2, 3])
     
     t = date2index(dates, nutime, select='nearest')
     assert_equal(t, [1,2,3])
 
     # Test dates outside the support with select
     t = date2index(datetime(1949,12,1), nutime, select='nearest')
     assert_equal(t, 0)
     
     t = date2index(datetime(1978,1,1), nutime, select='nearest')
     assert_equal(t, 365)
     
     # Test dates outsied the support with before
     self.assertRaises(ValueError, date2index, datetime(1949,12,1), nutime, select='before')
     
     t = date2index(datetime(1978,1,1), nutime, select='before')
     assert_equal(t, 365)
     
     # Test dates outsied the support with after
     t = date2index(datetime(1949,12,1), nutime, select='after')
     assert_equal(t, 0)
     
     self.assertRaises(ValueError, date2index, datetime(1978,1,1), nutime, select='after')
Exemplo n.º 10
0
 def test_failure(self):
     nutime = self.TestTime(datetime(1950, 1, 1), 366, 24,
                            'hours since 1900-01-01', 'standard')
     try:
         t = date2index(datetime(1949, 2, 1), nutime)
     except ValueError:
         pass
     else:
         raise ValueError('This test should have failed.')
Exemplo n.º 11
0
 def test_issue444(self):
     # make sure integer overflow not causing error in
     # calculation of nearest index when sum of adjacent
     # time values won't fit in 32 bits.
     query_time = datetime(2037, 1, 3, 21, 12)
     print(self.time_vars['time3'])
     index = date2index(query_time, self.time_vars['time3'],
                        select='nearest')
     assert(index == 11)
Exemplo n.º 12
0
    def test_select_nc(self):
        f = Dataset(self.file, "r")
        nutime = f.variables["time"]

        dates = [datetime(1950, 1, 2, 6), datetime(1950, 1, 3), datetime(1950, 1, 3, 18)]

        t = date2index(dates, nutime, select="before")
        assert_equal(t, [1, 2, 2])

        t = date2index(dates, nutime, select="after")
        assert_equal(t, [2, 2, 3])

        t = date2index(dates, nutime, select="nearest")
        assert_equal(t, [1, 2, 3])

        # Test dates outside the support with select
        t = date2index(datetime(1949, 12, 1), nutime, select="nearest")
        assert_equal(t, 0)

        t = date2index(datetime(1978, 1, 1), nutime, select="nearest")
        assert_equal(t, 365)

        # Test dates outside the support with before
        self.assertRaises(ValueError, date2index, datetime(1949, 12, 1), nutime, select="before")

        t = date2index(datetime(1978, 1, 1), nutime, select="before")
        assert_equal(t, 365)

        # Test dates outside the support with after
        t = date2index(datetime(1949, 12, 1), nutime, select="after")
        assert_equal(t, 0)

        self.assertRaises(ValueError, date2index, datetime(1978, 1, 1), nutime, select="after")
        # test microsecond and millisecond units
        unix_epoch = "milliseconds since 1970-01-01T00:00:00Z"
        from netCDF4 import date2num

        d = datetime(2038, 1, 19, 3, 14, 7)
        millisecs = int(date2num(d, unix_epoch, calendar="proleptic_gregorian"))
        assert_equal(millisecs, (2 ** 32 / 2 - 1) * 1000)
        unix_epoch = "microseconds since 1970-01-01T00:00:00Z"
        microsecs = int(date2num(d, unix_epoch))
        assert_equal(microsecs, (2 ** 32 / 2 - 1) * 1000000)
        # test microsecond accuracy in date2num/num2date roundtrip
        # note: microsecond accuracy lost for time intervals greater
        # than about 270 years.
        from dateutil.tz import tzutc

        units = "microseconds since 1776-07-04 00:00:00-12:00"
        dates = [
            datetime(1962, 10, 27, 6, 1, 30, 9001),
            datetime(1993, 11, 21, 12, 5, 25, 999),
            datetime(1995, 11, 25, 18, 7, 59, 999999),
        ]
        times2 = date2num(dates, units)
        dates2 = num2date(times2, units)
        for date, date2 in zip(dates, dates2):
            assert_equal(date.replace(tzinfo=tzutc()), date2)
        f.close()
Exemplo n.º 13
0
    def test_select_nc(self):
        f = Dataset(self.file, 'r')
        nutime = f.variables['time']

        dates = [datetime(1950, 1, 2, 6), datetime(
            1950, 1, 3), datetime(1950, 1, 3, 18)]

        t = date2index(dates, nutime, select='before')
        assert_equal(t, [1, 2, 2])

        t = date2index(dates, nutime, select='after')
        assert_equal(t, [2, 2, 3])

        t = date2index(dates, nutime, select='nearest')
        assert_equal(t, [1, 2, 3])

        # Test dates outside the support with select
        t = date2index(datetime(1949, 12, 1), nutime, select='nearest')
        assert_equal(t, 0)

        t = date2index(datetime(1978, 1, 1), nutime, select='nearest')
        assert_equal(t, 365)

        # Test dates outside the support with before
        self.assertRaises(
            ValueError, date2index, datetime(1949, 12, 1), nutime, select='before')

        t = date2index(datetime(1978, 1, 1), nutime, select='before')
        assert_equal(t, 365)

        # Test dates outside the support with after
        t = date2index(datetime(1949, 12, 1), nutime, select='after')
        assert_equal(t, 0)

        self.assertRaises(
            ValueError, date2index, datetime(1978, 1, 1), nutime, select='after')
        # test microsecond and millisecond units
        unix_epoch = "milliseconds since 1970-01-01T00:00:00Z"
        from netCDF4 import date2num
        d = datetime(2038, 1, 19, 3, 14, 7)
        millisecs = int(
            date2num(d, unix_epoch, calendar='proleptic_gregorian'))
        assert_equal(millisecs, (2 ** 32 / 2 - 1) * 1000)
        unix_epoch = "microseconds since 1970-01-01T00:00:00Z"
        microsecs = int(date2num(d, unix_epoch))
        assert_equal(microsecs, (2 ** 32 / 2 - 1) * 1000000)
        # test microsecond accuracy in date2num/num2date roundtrip
        # note: microsecond accuracy lost for time intervals greater
        # than about 270 years.
        from dateutil.tz import tzutc
        units = 'microseconds since 1776-07-04 00:00:00-12:00'
        dates =\
            [datetime(1962, 10, 27, 6, 1, 30, 9001), datetime(
                1993, 11, 21, 12, 5, 25, 999), datetime(1995, 11, 25, 18, 7, 59, 999999)]
        times2 = date2num(dates, units)
        dates2 = num2date(times2, units)
        for date, date2 in zip(dates, dates2):
            assert_equal(date.replace(tzinfo=tzutc()), date2)
        f.close()
Exemplo n.º 14
0
 def test_issue444(self):
     # make sure integer overflow not causing error in
     # calculation of nearest index when sum of adjacent
     # time values won't fit in 32 bits.
     ntimes = 20
     f = Dataset(self.file, 'r')
     query_time = datetime(2037, 1, 3, 21, 12)
     index = date2index(query_time, f.variables['time3'], select='nearest')
     assert (index == 11)
     f.close()
Exemplo n.º 15
0
 def test_simple(self):
     t = date2index(datetime(1950,2,1), self.standardtime)
     assert_equal(t, 31)
Exemplo n.º 16
0
 def test_list(self):
     t = date2index(
         [datetime(1950, 2, 1), datetime(1950, 2, 3)], self.standardtime)
     assert_equal(t, [31, 33])
Exemplo n.º 17
0
 def test_list(self):
     t = date2index([datetime(1950,2,1), datetime(1950,2,3)], self.standardtime)
     assert_equal(t, [31, 33])
Exemplo n.º 18
0
 def test_simple(self):
     t = date2index(datetime(1950, 2, 1), self.standardtime)
     assert_equal(t, 31)
Exemplo n.º 19
0
 def test_singletime(self):
     # issue 215 test (date2index with time variable length == 1)
     time2 = self.time_vars['time2']
     result_index = date2index(self.first_timestamp, time2, select="exact")
     assert_equal(result_index, 0)