In [1]:

    

url = 'http://opendap.nationaleregenradar.nl/thredds/dodsC/radar/TF2400_A/2000/01/01/RAD_TF2400_A_20000101000000.h5'


    

In [2]:

    

import netCDF4


    

In [3]:

    

ds = netCDF4.Dataset(url)
ds.variables


    

    
Out[3]:





OrderedDict([(u'available', <netCDF4.Variable object at 0x3aac3b0>), (u'east', <netCDF4.Variable object at 0x3aac7a0>), (u'north', <netCDF4.Variable object at 0x3aac830>), (u'time', <netCDF4.Variable object at 0x3aac8c0>), (u'precipitation', <netCDF4.Variable object at 0x3aac950>)])

In [4]:

    

time = netCDF4.num2date(ds.variables['time'][:], ds.variables['time'].units)
idx = np.where(time < datetime.datetime.now())[0]


    

In [26]:

    

P = ds.variables['precipitation'][:,:,(idx.max()-90):(idx.max()-2)]
#p = ds.variables['precipitation'][:,:,time.shape[0]-1]


    

In [32]:

    

time_selected = time[idx.max()-90:idx.max()-2]
maxidx = np.argmax(P.sum(0).sum(0))

plt.imshow(P[:,:,maxidx], cmap='Blues')
plt.title(time_selected[maxidx])


    

    
Out[32]:





<matplotlib.text.Text at 0x61f6a90>






    

In [ ]: