import numpy as np
d_array = np.arrange(np.datetime64('date_begin'),
np.datetime64('date_end'),
3,
dtype = 'datetime64[h]') # indicates the time step unit
# here every 3hr