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