coastlines_filename = './data/coastlines/Global_coastlines_2015_v1_low_res.shp'
reconstructed_coastlines_filename = './tmp/coastlines.shp'
rotation_filename = './data/rotations/PALEOMAP_PlateModel.rot'
reconstruction_time = 550

fig = plt.figure(figsize = (40, 40))

for posisjon, tid in enumerate(range(0,600,50)):
    ax_map=plt.subplot(6,2,posisjon+1)
    plt.title(str(tid), fontsize=40)

    reconstructed_coastlines_filename = './tmp/coastlines'  + str(tid) + '.shp'
    
    pygplates.reconstruct(coastlines_filename, 
                      rotation_filename, 
                      reconstructed_coastlines_filename, 
                      tid)
    
    m = Basemap(projection = 'robin', lon_0 = 0.0, lat_0 = 0.0, resolution = 'c', ax = ax_map)
    coastlines_info = m.readshapefile('./tmp/coastlines'+ str(tid), 'coastlines', drawbounds = False, color = None)

    for index, coast_geom in enumerate(m.coastlines):
        col = plate_tectonic_utils.get_colour_by_plateid(int(m.coastlines_info[index]['PLATEID1']))
        poly = Polygon(coast_geom, facecolor = col, edgecolor = 'black', alpha = 0.9, zorder = 3)
        ax_map.add_patch(poly)
        

   
    m.drawparallels(np.arange(-90., 91., 30.), labels = [True, True, False, False])
    m.drawmeridians(np.arange(-180., 181., 45.), labels = [False, False, False, True])