splits-and-merges.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'0.36.0'"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pygplates\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "#import pandas as pd\n",
    "import glob\n",
    "\n",
    "import splits_and_merges as snm\n",
    "# reload(snm)\n",
    "\n",
    "%matplotlib inline\n",
    "\n",
    "#######################################################\n",
    "# Define Input files for Muller 2016 AREPS model\n",
    "#######################################################\n",
    "\n",
    "\n",
    "MODELDIR = '/Applications/GPlates-2.0.0/SampleData/FeatureCollections/'\n",
    "RotFile_List = ['%s/Rotations/Matthews_etal_GPC_2016_410-0Ma_GK07.rot' % MODELDIR]\n",
    "GPML_List = ['%s/DynamicPolygons/Matthews_etal_GPC_2016_MesozoicCenozoic_PlateTopologies.gpmlz' % MODELDIR,\\\n",
    "             '%s/DynamicPolygons/Matthews_etal_GPC_2016_Paleozoic_PlateTopologies.gpmlz' % MODELDIR,\\\n",
    "             '%s/DynamicPolygons/Matthews_etal_GPC_2016_TopologyBuildingBlocks.gpmlz' % MODELDIR]\n",
    "\n",
    "\n",
    "#MODELDIR = '/Users/Simon/Desktop/Global_Model_WD_Internal_Release_2016_v3_RM17/'\n",
    "#RotFile_List = glob.glob('%s/*.rot' % MODELDIR)\n",
    "#GPML_List = glob.glob('%s/*.gpml' % MODELDIR)\n",
    "\n",
    "\n",
    "#####################################\n",
    "rotation_model = pygplates.RotationModel(RotFile_List)\n",
    "topology_features = pygplates.FeatureCollection()\n",
    "for file in GPML_List:\n",
    "    topology_feature = pygplates.FeatureCollection(file)\n",
    "    topology_features.add(topology_feature)\n",
    "#'''\n",
    "\n",
    "pygplates.__version__\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 51.00 Ma and 50.00 Ma are: \n",
      " {901, 902, 903, 904, 406, 918, 409, 922, 801, 802, 804, 301, 304, 307, 308, 701, 317, 201, 604, 608, 609, 738, 224, 101, 102, 501}\n",
      "\n",
      "plates that appeared between 51.00 Ma and 50.00 Ma are: \n",
      " {2035}\n",
      "\n",
      "plates that disappeared between 51.00 Ma and 50.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 2035 mapped to plate 101 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 50.00 Ma and 49.00 Ma are: \n",
      " {901, 902, 903, 904, 406, 918, 409, 922, 801, 802, 804, 301, 304, 307, 308, 317, 701, 201, 604, 608, 609, 738, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 50.00 Ma and 49.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 50.00 Ma and 49.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 49.00 Ma and 48.00 Ma are: \n",
      " {901, 902, 903, 904, 406, 918, 409, 922, 801, 802, 804, 301, 304, 307, 308, 317, 701, 201, 604, 608, 609, 738, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 49.00 Ma and 48.00 Ma are: \n",
      " {861}\n",
      "\n",
      "plates that disappeared between 49.00 Ma and 48.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 861 mapped to plate 801 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 48.00 Ma and 47.00 Ma are: \n",
      " {901, 902, 903, 904, 406, 918, 409, 922, 801, 802, 804, 301, 304, 307, 308, 317, 701, 201, 604, 861, 608, 609, 738, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 48.00 Ma and 47.00 Ma are: \n",
      " {645, 665, 926, 688}\n",
      "\n",
      "plates that disappeared between 48.00 Ma and 47.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 645 mapped to plate 608 at delta time\n",
      "Centroid for plate 665 mapped to plate 608 at delta time\n",
      "Centroid for plate 926 mapped to plate 608 at delta time\n",
      "Centroid for plate 688 mapped to plate 609 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 47.00 Ma and 46.00 Ma are: \n",
      " {901, 645, 903, 904, 902, 406, 918, 665, 922, 926, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 604, 861, 608, 609, 738, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 47.00 Ma and 46.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 47.00 Ma and 46.00 Ma are: \n",
      " {409}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 409 mapped to plate 101 at delta time\n",
      ">>>>>>>> [(409, 101)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 46.00 Ma and 45.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 918, 665, 922, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 604, 861, 608, 609, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 46.00 Ma and 45.00 Ma are: \n",
      " {603, 735}\n",
      "\n",
      "plates that disappeared between 46.00 Ma and 45.00 Ma are: \n",
      " {738, 926}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 603 mapped to plate 604 at delta time\n",
      "Centroid for plate 735 mapped to plate 738 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 738 mapped to plate 735 at delta time\n",
      "Centroid for plate 926 mapped to plate 735 at delta time\n",
      ">>>>>>>> [(738, 735), (926, 735)]\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 45.00 Ma and 44.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 918, 665, 922, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 603, 604, 861, 735, 608, 609, 224, 101, 102, 2035, 501}\n",
      "\n",
      "plates that appeared between 45.00 Ma and 44.00 Ma are: \n",
      " {855}\n",
      "\n",
      "plates that disappeared between 45.00 Ma and 44.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 855 mapped to plate 901 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 44.00 Ma and 43.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 918, 406, 665, 922, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 855, 603, 604, 861, 735, 608, 609, 224, 101, 102, 2035}\n",
      "\n",
      "plates that appeared between 44.00 Ma and 43.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 44.00 Ma and 43.00 Ma are: \n",
      " {501}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 501 mapped to plate 801 at delta time\n",
      ">>>>>>>> [(501, 801)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 43.00 Ma and 42.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 918, 406, 665, 922, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 855, 603, 604, 861, 735, 608, 609, 224, 101, 102, 2035}\n",
      "\n",
      "plates that appeared between 43.00 Ma and 42.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 43.00 Ma and 42.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 42.00 Ma and 41.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 918, 665, 922, 801, 802, 804, 301, 688, 304, 307, 308, 317, 701, 201, 855, 603, 604, 861, 735, 224, 609, 608, 101, 102, 2035}\n",
      "\n",
      "plates that appeared between 42.00 Ma and 41.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 42.00 Ma and 41.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 41.00 Ma and 40.00 Ma are: \n",
      " {901, 902, 645, 904, 903, 406, 665, 922, 801, 802, 804, 301, 688, 304, 308, 317, 701, 201, 855, 603, 604, 861, 224, 609, 608, 101, 102, 2035}\n",
      "\n",
      "plates that appeared between 41.00 Ma and 40.00 Ma are: \n",
      " {663, 673, 818, 819, 9923, 328, 617}\n",
      "\n",
      "plates that disappeared between 41.00 Ma and 40.00 Ma are: \n",
      " {307, 918, 735}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 663 mapped to plate 735 at delta time\n",
      "Centroid for plate 673 mapped to plate 604 at delta time\n",
      "Centroid for plate 818 mapped to plate 802 at delta time\n",
      "Centroid for plate 819 mapped to plate 802 at delta time\n",
      "Centroid for plate 9923 mapped to plate 801 at delta time\n",
      "Centroid for plate 328 mapped to plate 301 at delta time\n",
      "Centroid for plate 617 mapped to plate 301 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 307 mapped to plate 701 at delta time\n",
      "Centroid for plate 918 mapped to plate 901 at delta time\n",
      "Centroid for plate 735 mapped to plate 663 at delta time\n",
      ">>>>>>>> [(307, 701), (918, 901), (735, 663)]\n",
      "Only works for new splits with two plates\n",
      "(802, [818, 819])\n",
      ">>>>>>>>>> Great Circle Normal is (-0.100146, 0.874088, -0.475332)\n",
      "(301, [328, 617])\n",
      "Only works for new splits with one geometry so far\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 40.00 Ma and 39.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 663, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 701, 317, 9923, 328, 201, 855, 603, 604, 861, 608, 609, 224, 101, 102, 617, 2035}\n",
      "\n",
      "plates that appeared between 40.00 Ma and 39.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 40.00 Ma and 39.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 39.00 Ma and 38.00 Ma are: \n",
      " {901, 645, 903, 904, 902, 406, 663, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 701, 317, 9923, 328, 201, 603, 604, 861, 608, 609, 224, 101, 102, 617, 2035}\n",
      "\n",
      "plates that appeared between 39.00 Ma and 38.00 Ma are: \n",
      " {857}\n",
      "\n",
      "plates that disappeared between 39.00 Ma and 38.00 Ma are: \n",
      " {855}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 857 mapped to plate 855 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 855 mapped to plate 857 at delta time\n",
      ">>>>>>>> [(855, 857)]\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 38.00 Ma and 37.00 Ma are: \n",
      " {901, 645, 903, 904, 902, 406, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 701, 317, 9923, 328, 201, 857, 603, 861, 608, 609, 224, 101, 102, 2035}\n",
      "\n",
      "plates that appeared between 38.00 Ma and 37.00 Ma are: \n",
      " {613}\n",
      "\n",
      "plates that disappeared between 38.00 Ma and 37.00 Ma are: \n",
      " {663, 617, 604}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 613 mapped to plate 617 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 663 mapped to plate 613 at delta time\n",
      "Centroid for plate 617 mapped to plate 613 at delta time\n",
      "Centroid for plate 604 mapped to plate 613 at delta time\n",
      ">>>>>>>> [(663, 613), (617, 613), (604, 613)]\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 37.00 Ma and 36.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 701, 317, 9923, 328, 201, 857, 603, 861, 224, 609, 608, 613, 102, 101, 2035}\n",
      "\n",
      "plates that appeared between 37.00 Ma and 36.00 Ma are: \n",
      " {838, 839, 855}\n",
      "\n",
      "plates that disappeared between 37.00 Ma and 36.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 838 mapped to plate 801 at delta time\n",
      "Centroid for plate 839 mapped to plate 857 at delta time\n",
      "Centroid for plate 855 mapped to plate 857 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "(857, [839, 855])\n",
      "Only works for new splits with one geometry so far\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 36.00 Ma and 35.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 701, 317, 838, 839, 328, 201, 855, 857, 603, 861, 224, 608, 609, 613, 102, 101, 2035}\n",
      "\n",
      "plates that appeared between 36.00 Ma and 35.00 Ma are: \n",
      " {853, 698}\n",
      "\n",
      "plates that disappeared between 36.00 Ma and 35.00 Ma are: \n",
      " {9923}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 853 mapped to plate 9923 at delta time\n",
      "Centroid for plate 698 mapped to plate 9923 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 9923 mapped to plate 853 at delta time\n",
      ">>>>>>>> [(9923, 853)]\n",
      "(9923, [853, 698])\n",
      ">>>>>>>>>> Great Circle Normal is (0.515384, 0.0307048, -0.856409)\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 35.00 Ma and 34.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 665, 922, 801, 802, 673, 804, 301, 688, 304, 818, 819, 308, 698, 317, 701, 838, 839, 328, 201, 853, 855, 857, 603, 861, 224, 608, 613, 102, 101, 2035}\n",
      "\n",
      "plates that appeared between 35.00 Ma and 34.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 35.00 Ma and 34.00 Ma are: \n",
      " {609}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 609 mapped to plate 608 at delta time\n",
      ">>>>>>>> [(609, 608)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 34.00 Ma and 33.00 Ma are: \n",
      " {901, 645, 902, 904, 903, 406, 665, 922, 801, 802, 673, 804, 301, 688, 818, 819, 698, 701, 838, 839, 201, 853, 857, 603, 861, 224, 608, 613, 102, 101, 2035}\n",
      "\n",
      "plates that appeared between 34.00 Ma and 33.00 Ma are: \n",
      " {3201, 522, 689, 306, 708}\n",
      "\n",
      "plates that disappeared between 34.00 Ma and 33.00 Ma are: \n",
      " {304, 308, 317, 328, 855}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 3201 mapped to plate 304 at delta time\n",
      "Centroid for plate 522 mapped to plate 328 at delta time\n",
      "Centroid for plate 689 mapped to plate 688 at delta time\n",
      "Centroid for plate 306 mapped to plate 304 at delta time\n",
      "Centroid for plate 708 mapped to plate 304 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 304 mapped to plate 301 at delta time\n",
      "Centroid for plate 308 mapped to plate 701 at delta time\n",
      "Centroid for plate 317 mapped to plate 301 at delta time\n",
      "Centroid for plate 328 mapped to plate 301 at delta time\n",
      "Centroid for plate 855 mapped to plate 839 at delta time\n",
      ">>>>>>>> [(304, 301), (308, 701), (317, 301), (328, 301), (855, 839)]\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 33.00 Ma and 32.00 Ma are: \n",
      " {3201, 901, 645, 902, 903, 904, 522, 406, 665, 922, 801, 802, 673, 804, 301, 688, 689, 306, 819, 818, 698, 701, 708, 838, 839, 201, 853, 857, 603, 861, 224, 608, 613, 101, 2035}\n",
      "\n",
      "plates that appeared between 33.00 Ma and 32.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 33.00 Ma and 32.00 Ma are: \n",
      " {102}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 102 mapped to plate 101 at delta time\n",
      ">>>>>>>> [(102, 101)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 32.00 Ma and 31.00 Ma are: \n",
      " {3201, 901, 645, 902, 904, 903, 522, 406, 922, 801, 802, 673, 804, 301, 688, 689, 306, 819, 818, 698, 701, 708, 838, 839, 201, 853, 857, 603, 861, 224, 608, 613, 101, 2035}\n",
      "\n",
      "plates that appeared between 32.00 Ma and 31.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 32.00 Ma and 31.00 Ma are: \n",
      " {665}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 665 mapped to plate 645 at delta time\n",
      ">>>>>>>> [(665, 645)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 31.00 Ma and 30.00 Ma are: \n",
      " {3201, 645, 902, 903, 904, 901, 522, 406, 922, 801, 673, 802, 804, 301, 688, 689, 306, 819, 818, 698, 701, 708, 839, 201, 853, 603, 861, 224, 608, 101}\n",
      "\n",
      "plates that appeared between 31.00 Ma and 30.00 Ma are: \n",
      " {653, 659, 686, 835, 709, 604, 607}\n",
      "\n",
      "plates that disappeared between 31.00 Ma and 30.00 Ma are: \n",
      " {838, 857, 613, 2035}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 653 mapped to plate 608 at delta time\n",
      "Centroid for plate 659 mapped to plate 608 at delta time\n",
      "Centroid for plate 686 mapped to plate 603 at delta time\n",
      "Centroid for plate 835 mapped to plate 838 at delta time\n",
      "Centroid for plate 709 mapped to plate 701 at delta time\n",
      "Centroid for plate 604 mapped to plate 613 at delta time\n",
      "Centroid for plate 607 mapped to plate 603 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 838 mapped to plate 835 at delta time\n",
      "Centroid for plate 857 mapped to plate 861 at delta time\n",
      "Centroid for plate 613 mapped to plate 604 at delta time\n",
      "Centroid for plate 2035 mapped to plate 224 at delta time\n",
      ">>>>>>>> [(838, 835), (857, 861), (613, 604), (2035, 224)]\n",
      "(608, [653, 659])\n",
      "Only works for new splits with one geometry so far\n",
      "Only works for new splits with two plates\n",
      "Only works for new splits with two plates\n",
      "(603, [686, 607])\n",
      ">>>>>>>>>> Great Circle Normal is (0.774884, 0.0421399, 0.630697)\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 30.00 Ma and 29.00 Ma are: \n",
      " {3201, 645, 902, 901, 904, 903, 522, 653, 659, 406, 922, 801, 673, 802, 804, 301, 686, 688, 689, 306, 819, 818, 698, 701, 835, 708, 709, 839, 201, 853, 603, 604, 861, 607, 224, 608, 101}\n",
      "\n",
      "plates that appeared between 30.00 Ma and 29.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 30.00 Ma and 29.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 29.00 Ma and 28.00 Ma are: \n",
      " {645, 901, 902, 904, 903, 522, 653, 659, 406, 922, 673, 801, 802, 804, 301, 686, 688, 689, 306, 819, 818, 698, 701, 835, 708, 709, 839, 201, 853, 603, 604, 861, 607, 224, 608, 101}\n",
      "\n",
      "plates that appeared between 29.00 Ma and 28.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 29.00 Ma and 28.00 Ma are: \n",
      " {3201}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 3201 mapped to plate 301 at delta time\n",
      ">>>>>>>> [(3201, 301)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 28.00 Ma and 27.00 Ma are: \n",
      " {645, 902, 901, 904, 903, 522, 653, 659, 406, 922, 673, 801, 802, 804, 301, 686, 688, 689, 306, 818, 819, 698, 701, 835, 708, 709, 839, 201, 853, 603, 604, 861, 607, 224, 608, 101}\n",
      "\n",
      "plates that appeared between 28.00 Ma and 27.00 Ma are: \n",
      " {503}\n",
      "\n",
      "plates that disappeared between 28.00 Ma and 27.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 503 mapped to plate 701 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 27.00 Ma and 26.00 Ma are: \n",
      " {645, 902, 901, 904, 903, 522, 653, 659, 406, 922, 801, 802, 673, 301, 686, 688, 689, 306, 818, 819, 698, 701, 835, 708, 709, 839, 201, 853, 603, 604, 861, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 27.00 Ma and 26.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 27.00 Ma and 26.00 Ma are: \n",
      " {804}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 804 mapped to plate 802 at delta time\n",
      ">>>>>>>> [(804, 802)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 26.00 Ma and 25.00 Ma are: \n",
      " {645, 902, 901, 904, 903, 522, 653, 659, 406, 922, 801, 802, 673, 301, 686, 688, 689, 306, 819, 818, 701, 708, 709, 839, 201, 853, 603, 604, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 26.00 Ma and 25.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 26.00 Ma and 25.00 Ma are: \n",
      " {835, 698, 861}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 835 mapped to plate 853 at delta time\n",
      "Centroid for plate 698 mapped to plate 801 at delta time\n",
      "Centroid for plate 861 mapped to plate 839 at delta time\n",
      ">>>>>>>> [(835, 853), (698, 801), (861, 839)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 25.00 Ma and 24.00 Ma are: \n",
      " {645, 901, 902, 904, 903, 522, 653, 659, 406, 922, 801, 802, 673, 301, 686, 688, 306, 818, 819, 701, 708, 709, 839, 201, 853, 603, 604, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 25.00 Ma and 24.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 25.00 Ma and 24.00 Ma are: \n",
      " {689}\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 689 mapped to plate 688 at delta time\n",
      ">>>>>>>> [(689, 688)]\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 24.00 Ma and 23.00 Ma are: \n",
      " {645, 901, 903, 904, 902, 522, 653, 659, 406, 922, 801, 802, 673, 301, 686, 306, 818, 819, 701, 708, 709, 839, 201, 853, 603, 604, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 24.00 Ma and 23.00 Ma are: \n",
      " {911}\n",
      "\n",
      "plates that disappeared between 24.00 Ma and 23.00 Ma are: \n",
      " {688}\n",
      "\n",
      "Working on appearing plates\n",
      "Centroid for plate 911 mapped to plate 902 at delta time\n",
      "\n",
      "Working on disappearing plates\n",
      "Centroid for plate 688 mapped to plate 653 at delta time\n",
      ">>>>>>>> [(688, 653)]\n",
      "Only works for new splits with two plates\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 23.00 Ma and 22.00 Ma are: \n",
      " {645, 901, 903, 904, 902, 522, 653, 911, 659, 406, 922, 801, 802, 673, 301, 686, 306, 818, 819, 701, 708, 709, 839, 201, 853, 603, 604, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 23.00 Ma and 22.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 23.00 Ma and 22.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n",
      "\n",
      "##############################\n",
      "\n",
      "\n",
      "plates that persist between 22.00 Ma and 21.00 Ma are: \n",
      " {645, 901, 903, 904, 902, 522, 653, 911, 659, 406, 922, 801, 802, 673, 301, 686, 306, 819, 818, 701, 708, 709, 839, 201, 853, 603, 604, 607, 224, 608, 101, 503}\n",
      "\n",
      "plates that appeared between 22.00 Ma and 21.00 Ma are: \n",
      " set()\n",
      "\n",
      "plates that disappeared between 22.00 Ma and 21.00 Ma are: \n",
      " set()\n",
      "\n",
      "Working on appearing plates\n",
      "\n",
      "Working on disappearing plates\n",
      ">>>>>>>> []\n"
     ]
    }
   ],
   "source": [
    "time_list = np.arange(50,20,-1)\n",
    "\n",
    "for time in time_list:\n",
    "    \n",
    "    print('\\n##############################\\n\\n')\n",
    "        \n",
    "    time_from = time+1\n",
    "    time_to = time\n",
    "    \n",
    "    (appearing_plate_map,\n",
    "     disappearing_plate_map) = snm.find_plate_change_for_time_step(topology_features,rotation_model,time_from,time_to)\n",
    "                \n",
    "    print('>>>>>>>> %s' % disappearing_plate_map)\n",
    "    \n",
    "    new_plates_grouped_by_origin_plate = snm.match_splits_to_origin(appearing_plate_map)\n",
    "        \n",
    "    resolved_topologies = []\n",
    "    shared_boundary_sections = []\n",
    "    pygplates.resolve_topologies(topology_features, rotation_model, resolved_topologies, time_to, shared_boundary_sections)\n",
    "    \n",
    "    for plate_split in new_plates_grouped_by_origin_plate:\n",
    "        \n",
    "        if len(plate_split[1])!=2:\n",
    "            print('Only works for new splits with two plates')\n",
    "        else:\n",
    "            gca = snm.get_great_circle_along_plate_split(plate_split,\n",
    "                                                     shared_boundary_sections)\n",
    "\n",
    "            if gca is not None:\n",
    "                print('>>>>>>>>>> Great Circle Normal is %s' % gca.get_great_circle_normal())\n",
    "\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Splitting plate still exists.....\n",
      "Time is 50, Plates are (101, [2035]), Area Ratio 1.00165\n",
      "Splitting plate still exists.....\n",
      "Time is 48, Plates are (801, [861]), Area Ratio 0.99711\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 47, Plates are (608, [645, 665, 926]), Area Ratio 2.48827\n",
      "Splitting plate still exists.....\n",
      "Time is 47, Plates are (609, [688]), Area Ratio 1.06554\n",
      "Time is 45, Plates are (738, [735]), Area Ratio 0.97248\n",
      "Splitting plate still exists.....\n",
      "Time is 45, Plates are (604, [603]), Area Ratio 1.00763\n",
      "Splitting plate still exists.....\n",
      "Time is 44, Plates are (901, [855]), Area Ratio 0.99751\n",
      "Splitting plate still exists.....\n",
      "Time is 40, Plates are (801, [9923]), Area Ratio 1.00144\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 40, Plates are (802, [818, 819]), Area Ratio 1.98210\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 40, Plates are (301, [328, 617]), Area Ratio 1.98032\n",
      "Splitting plate still exists.....\n",
      "Time is 40, Plates are (604, [673]), Area Ratio 1.00723\n",
      "Time is 40, Plates are (735, [663]), Area Ratio 0.96455\n",
      "Time is 38, Plates are (855, [857]), Area Ratio 0.99910\n",
      "Time is 37, Plates are (617, [613]), Area Ratio 2.07222\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 36, Plates are (857, [839, 855]), Area Ratio 2.29989\n",
      "Splitting plate still exists.....\n",
      "Time is 36, Plates are (801, [838]), Area Ratio 0.99297\n",
      "Time is 35, Plates are (9923, [853, 698]), Area Ratio 1.06357\n",
      "Time is 33, Plates are (304, [3201, 306, 708]), Area Ratio 0.18619\n",
      "Time is 33, Plates are (328, [522]), Area Ratio 0.52025\n",
      "Splitting plate still exists.....\n",
      "Time is 33, Plates are (688, [689]), Area Ratio 1.10128\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 30, Plates are (608, [653, 659]), Area Ratio 1.86821\n",
      "Time is 30, Plates are (613, [604]), Area Ratio 0.95792\n",
      "Time is 30, Plates are (838, [835]), Area Ratio 0.72483\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 30, Plates are (603, [686, 607]), Area Ratio 1.31990\n",
      "Splitting plate still exists.....\n",
      "Time is 30, Plates are (701, [709]), Area Ratio 1.00266\n",
      "Splitting plate still exists.....\n",
      "Time is 27, Plates are (701, [503]), Area Ratio 1.00221\n",
      "Splitting plate still exists.....\n",
      "Time is 23, Plates are (902, [911]), Area Ratio 0.98758\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 20, Plates are (801, [501, 511]), Area Ratio 1.59808\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 20, Plates are (701, [328, 334, 337, 340]), Area Ratio 3.98562\n",
      "Splitting plate still exists.....\n",
      "Time is 17, Plates are (818, [815]), Area Ratio 1.07341\n",
      "Splitting plate still exists.....\n",
      "Time is 15, Plates are (911, [971]), Area Ratio 1.04914\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 12, Plates are (801, [821, 615]), Area Ratio 1.90667\n",
      "Time is 12, Plates are (922, [849]), Area Ratio 0.86577\n",
      "Splitting plate still exists.....\n",
      "Time is 10, Plates are (801, [824]), Area Ratio 1.00398\n",
      "Splitting plate still exists.....\n",
      "Time is 10, Plates are (645, [646]), Area Ratio 0.79106\n",
      "Time is 10, Plates are (902, [969, 970, 909]), Area Ratio 1.00615\n",
      "Time is 10, Plates are (971, [972]), Area Ratio 1.18501\n",
      "Splitting plate still exists.....\n",
      "Time is 10, Plates are (301, [677]), Area Ratio 1.09029\n",
      "Splitting plate still exists.....\n",
      "Time is 9, Plates are (801, [566]), Area Ratio 0.99785\n",
      "Splitting plate still exists.....\n",
      "Time is 7, Plates are (849, [829]), Area Ratio 0.81766\n",
      "Time is 7, Plates are (615, [678, 669]), Area Ratio 0.97178\n",
      "Splitting plate still exists.....\n",
      "Time is 5, Plates are (659, [699]), Area Ratio 0.49595\n",
      "Splitting plate still exists.....\n",
      "Time is 5, Plates are (678, [564]), Area Ratio 1.02686\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 5, Plates are (566, [676, 664]), Area Ratio 1.71574\n",
      "Time is 4, Plates are (564, [674]), Area Ratio 0.88573\n",
      "Splitting plate still exists.....\n",
      "Splitting plate still exists.....\n",
      "Time is 4, Plates are (653, [697, 847]), Area Ratio 1.86779\n"
     ]
    }
   ],
   "source": [
    "time_list = np.arange(50,0,-1)\n",
    "verbose=False\n",
    "\n",
    "for time in time_list:\n",
    "    \n",
    "    if verbose:\n",
    "        print('\\n##############################\\n')\n",
    "        \n",
    "    time_from = time+1\n",
    "    time_to = time\n",
    "    \n",
    "    (appearing_plate_map,\n",
    "     disappearing_plate_map) = snm.find_plate_change_for_time_step(topology_features,rotation_model,time_from,time_to,verbose)\n",
    "                \n",
    "    new_plates_grouped_by_origin_plate = snm.match_splits_to_origin(appearing_plate_map)\n",
    "        \n",
    "    resolved_topologies_time_to = []\n",
    "    shared_boundary_sections_time_to = []\n",
    "    pygplates.resolve_topologies(topology_features, \n",
    "                                 rotation_model, \n",
    "                                 resolved_topologies_time_to, \n",
    "                                 time_to, \n",
    "                                 shared_boundary_sections_time_to)\n",
    "    \n",
    "    resolved_topologies_time_from = []\n",
    "    pygplates.resolve_topologies(topology_features, \n",
    "                                 rotation_model, \n",
    "                                 resolved_topologies_time_from, \n",
    "                                 time_from)\n",
    "    \n",
    "    \n",
    "    for plate_split in new_plates_grouped_by_origin_plate:\n",
    "        \n",
    "        #print plate_split\n",
    "        \n",
    "        area_after_split = 0\n",
    "        for topology in resolved_topologies_time_to:\n",
    "            for plate_id in plate_split[1]:\n",
    "                if topology.get_feature().get_reconstruction_plate_id()==plate_id:\n",
    "                    area_after_split += topology.get_resolved_geometry().get_area()\n",
    "                if topology.get_feature().get_reconstruction_plate_id()==plate_split[0]:\n",
    "                    print('Splitting plate still exists.....')\n",
    "                    area_after_split += topology.get_resolved_geometry().get_area()\n",
    "\n",
    "        for topology in resolved_topologies_time_from:\n",
    "            if topology.get_feature().get_reconstruction_plate_id()==plate_split[0]:\n",
    "                area_before_split = topology.get_resolved_geometry().get_area()\n",
    "\n",
    "        area_change_ratio = area_after_split/area_before_split\n",
    "        print('Time is %d, Plates are %s, Area Ratio %0.5f' % (time_to, plate_split, area_change_ratio))\n",
    "\n",
    "        #gca = get_great_circle_along_plate_split(plate_split,\n",
    "        #                                         shared_boundary_sections_time_to)\n",
    "        #if gca is not None:\n",
    "        #    print '>>>>>>>>>> Great Circle Normal is %s' % gca.get_great_circle_normal()\n",
    "\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[(328, 304, 308, 855, 317), (301, 301, 701, 839, 301)]\n"
     ]
    }
   ],
   "source": [
    "tmp = [(328, 301), (304, 301), (308, 701), (855, 839), (317, 301)]\n",
    "print(list(zip(*tmp)))\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "collapsed": false,
    "jupyter": {
     "outputs_hidden": false
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "(409,)\n",
      "(738, 926)\n",
      "(501,)\n",
      "(307, 918, 735)\n",
      "(855,)\n",
      "(663, 617, 604)\n",
      "(9923,)\n",
      "(609,)\n",
      "(304, 308, 317, 328, 855)\n",
      "(102,)\n",
      "(665,)\n",
      "(838, 857, 613, 2035)\n",
      "(3201,)\n",
      "(804,)\n",
      "(835, 698, 861)\n",
      "(689,)\n",
      "(688,)\n",
      "(522, 306, 708)\n",
      "(853, 839)\n",
      "(673,)\n",
      "(922,)\n",
      "(902, 971, 603, 604)\n",
      "(615,)\n",
      "(819, 972)\n",
      "(564,)\n",
      "(904, 969)\n",
      "[(409, 46), (738, 45), (926, 45), (501, 43), (307, 40), (918, 40), (735, 40), (855, 38), (663, 37), (617, 37), (604, 37), (9923, 35), (609, 34), (304, 33), (308, 33), (317, 33), (328, 33), (855, 33), (102, 32), (665, 31), (838, 30), (857, 30), (613, 30), (2035, 30), (3201, 28), (804, 26), (835, 25), (698, 25), (861, 25), (689, 24), (688, 23), (522, 20), (306, 20), (708, 20), (853, 18), (839, 18), (673, 14), (922, 12), (902, 10), (971, 10), (603, 10), (604, 10), (615, 7), (819, 5), (972, 5), (564, 4), (904, 3), (969, 3)]\n",
      "20\n"
     ]
    }
   ],
   "source": [
    "time_list = np.arange(50,0,-1)\n",
    "verbose=False\n",
    "\n",
    "plate_disappearance_time_lut = []\n",
    "\n",
    "for time in time_list:\n",
    "    \n",
    "    if verbose:\n",
    "        print('\\n##############################\\n')\n",
    "        \n",
    "    time_from = time+1\n",
    "    time_to = time\n",
    "    \n",
    "    (appearing_plate_map,\n",
    "     disappearing_plate_map) = snm.find_plate_change_for_time_step(topology_features,rotation_model,time_from,time_to,verbose)\n",
    "        \n",
    "    if len(list(zip(*disappearing_plate_map))) > 0:\n",
    "        print(list(zip(*disappearing_plate_map))[0])\n",
    "        for plate_id in list(zip(*disappearing_plate_map))[0]:\n",
    "            plate_id\n",
    "            plate_disappearance_time_lut.append((plate_id,time))\n",
    "    \n",
    "print(plate_disappearance_time_lut)\n",
    "\n",
    "for plate_disappearance_time in plate_disappearance_time_lut:\n",
    "    if plate_disappearance_time[0]==522:\n",
    "        print(plate_disappearance_time[1])\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.9.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}