num_rowcol = math.ceil(math.sqrt(num_t))
    fig = plt.figure()
    fig.clf()
    plot_count = 1
    for i in range(data_t.shape[0]):
       data = data_t[i]
       ax = fig.add_subplot(num_rowcol,num_rowcol,plot_count)
       plt.contourf(lat, lon, data, 64)
       plot_count += 1
    fig.tight_layout()
    plt.show()

# In[4]:

g = GDF()
g.debug = False


# In[5]:

start_date = dt2secs(date(year=2010,month=1,day=1))
end_date = dt2secs(date(year=2010, month=1, day=18))
data_request_descriptors = []
data_request_descriptors.append( {'storage_type': 'LS5TM',
                               'variables': ('B30',),
                               'dimensions': {'X': {'range': (140.0, 141.0)},
                                              'Y': {'range': (-36.0, -35.0)},
                                              'T': {'range': (start_date, end_date),
                                              #      'array_range': (0, 4)
                                                    #'crs': 'SSE', # Seconds since epoch
                                                    #'grouping_function': g.null_grouping