Paste code

coastlines_filename = &#039;./data/coastlines/Global_coastlines_2015_v1_low_res.shp&#039;
reconstructed_coastlines_filename = &#039;./tmp/coastlines.shp&#039;
rotation_filename = &#039;./data/rotations/PALEOMAP_PlateModel.rot&#039;
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 = &#039;./tmp/coastlines&#039;  + str(tid) + &#039;.shp&#039;
    
    pygplates.reconstruct(coastlines_filename, 
                      rotation_filename, 
                      reconstructed_coastlines_filename, 
                      tid)
    
    m = Basemap(projection = &#039;robin&#039;, lon_0 = 0.0, lat_0 = 0.0, resolution = &#039;c&#039;, ax = ax_map)
    coastlines_info = m.readshapefile(&#039;./tmp/coastlines&#039;+ str(tid), &#039;coastlines&#039;, 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][&#039;PLATEID1&#039;]))
        poly = Polygon(coast_geom, facecolor = col, edgecolor = &#039;black&#039;, 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])