continents_filename ='./data/continents/Matthews_etal_GPC_2016_ContinentalPolygons.shp'
rotation_filename = './data/rotations/PALEOMAP_PlateModel.rot'
reconstructed_continents_filename = './tmp/continents.shp'
reconstruction_time = 0


fig = plt.figure(figsize = (40, 40))
for bildeserie, tid in enumerate(range(0,600,50)):
   
    plt.title(str(tid-50) + ' MA', fontsize = 50)
    ax_map = plt.subplot(6, 2, bildeserie+1)

    reconstruction_time = tid
    reconstructed_continents_filename = './tmp/continents' + str(tid) + '.shp'

    
    pygplates.reconstruct(continents_filename, 
                      rotation_filename, 
                      reconstructed_continents_filename, 
                      reconstruction_time)
    ## opprett et basemap-objekt
    m = Basemap(projection = 'robin', lon_0 = 0.0, lat_0 = 0.0, resolution = 'c', ax = ax_map)
    #---------------- Legg til kontinentpolygoner -----------------------

    # last shape-fila for rekonstruerte kontinenter inn i kartet
    continents_info = m.readshapefile('./tmp/continents' + str(tid), 'continents', drawbounds = False, color = 'w')
    # gå igjennom for-løkka for kontinentgeometrier og plott som polygoner
    for index, cont_geom in enumerate(m.continents):
        poly = Polygon(cont_geom, facecolor = None, edgecolor = 'black', alpha = 0.2, zorder = 1)
        ax_map.add_patch(poly)

    # forskjønn basemap med linjer og merkelapper   
    m.drawparallels(np.arange(-90., 91., 30.), labels = [True, True, False, False])
    m.drawmeridians(np.arange(-180., 181., 45.), labels = [False, False, False, True])
    
plt.show()