From 6d84eeb1c22beb1ef83d7c11c6bb05f0e07b6c09 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 27 Jul 2021 12:00:56 -0400 Subject: [PATCH] Updated test with successful run --- .../swiftest_vs_swifter.ipynb | 112 ++++++++---------- 1 file changed, 50 insertions(+), 62 deletions(-) diff --git a/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb b/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb index 3016031e3..71a2c4da6 100644 --- a/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb +++ b/examples/symba_swifter_comparison/1pl_1tp_encounter/swiftest_vs_swifter.ipynb @@ -21,9 +21,9 @@ "output_type": "stream", "text": [ "Reading Swifter file param.swifter.in\n", - "Reading in time 1.355e-01\n", + "Reading in time 1.348e-01\n", "Creating Dataset\n", - "Successfully converted 199 output frames.\n", + "Successfully converted 198 output frames.\n", "Swifter simulation data stored as xarray DataSet .ds\n" ] } @@ -75,23 +75,23 @@ }, { "cell_type": "code", - "execution_count": 6, + "execution_count": 11, "metadata": {}, "outputs": [ { "data": { "text/plain": [ - "[,\n", - " ]" + "[,\n", + " ]" ] }, - "execution_count": 6, + "execution_count": 11, "metadata": {}, "output_type": "execute_result" }, { "data": { - "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZAAAAERCAYAAABVU/GxAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAbFElEQVR4nO3dfZQcdZ3v8fenk0hEkgUkgRmGbLImQAJK4M7yfABJskJ0CQ8LlyxiUPbmclQuLuv1ZG+O1+WesxpX3QOrrHtGWG9ULtnlIYaHbFhJyMEFQYMgeRgxrEQzyUiGCEJEEMn3/tGVME56JjPV3VXV3Z/XOXOmq+pX9ftOZ3q++dWv6luKCMzMzEaqlHcAZmbWmJxAzMwsFScQMzNLxQnEzMxScQIxM7NUnEDMzCyVpk0gkv5Z0g5JG2p0vFWSXpJ03yDbvyxpVy36MjNrBE2bQID/C5xXw+N9Abiy0gZJncDBNezLzKzwmjaBRMTDwC/7r5P0rmQk8YSk70o6dgTHWw28MnC9pFGUk8unqo3ZzKyRjM47gIx1AddExGZJpwD/CJxb5TE/DtwTEb2Sqg7QzKxRtEwCkXQQcDpwR78/9Ack2y4G/k+F3bZFxPuGOGY7cClwTk2DNTNrAC2TQCifrnspImYO3BARdwN3pzjmicBU4NkkKR0o6dmImFpNoGZmjaBp50AGioiXgeckXQqgshOqPOb9EXFEREyOiMnAq04eZtYqmjaBSLod+B5wjKQeSVcDVwBXS/oRsBGYN4LjfRe4A5iVHG/QU1tmZq1ALuduZmZpNO0IxMzM6qspJ9EPO+ywmDx5ct5hmJk1jCeeeOKFiJgwkn2aMoFMnjyZdevW5R2GmVnDkPSzke7jU1hmZpaKE4iZmaXiBGJmZqk4gZiZWSpOIGZmlkquCUTSeZKekfSspEUVtkvSPyTbn5Z0Uh5xmpnZvnJLIMlzNG4GzgdmAPMlzRjQ7HxgWvK1EPhqpkGamdmg8rwP5GTg2Yj4KYCkZZRrU23q12Ye8I0o11t5TNLBktoiorceAT3ZdQ2H7XqmHoc2M6ubkuCwqX/MAX/6d5n2m2cCORLY2m+5BzhlGG2OBPZJIJIWUh6lMGnSpFQBrd/2EkfHb1Lta2aWp9+8+GuyLgWeZwKp9Pi+gZUdh9OmvDKii/ITB+ns7ExVIfJDNyxLs5uZWW56XnyVMz//EH834z2ZJ5A8J9F7gKP6LXcA21O0MTNrWaXkCat5VFbPM4H8AJgmaYqktwGXA/cMaHMP8KHkaqxTgV/Va/7DzKwR7XlC9+4cnsyR2ymsiPidpI8DDwCjgH+OiI2Srkm2/xOwEpgLPAu8Cnw4r3jNzIpozwhkdw4jkFyr8UbESspJov+6f+r3OoCPZR2XmVmjyHME4jvRzcwaWKvOgZiZWZX2nsLKYQjiBGJm1sBKPoVlZmZpKMdJdCcQM7MGtmcEkkP+cAIxM2tkeV7G6wRiZtbA3kogOfSdfZdmZlYrb90H4hGImZmNgO8DMTOzVHwZr5mZpeJJdDMzS8W1sMzMLBVJSJ4DMTOzFEqST2GZmdnIleRTWGZmloJaaQQi6VBJ35G0Ofl+SIU2R0l6SFK3pI2SrssjVjOzoiuptWphLQJWR8Q0YHWyPNDvgL+KiOnAqcDHJM3IMEYzs4ZQklrqeSDzgKXJ66XAhQMbRERvRPwwef0K0A0cmVWAZmaNojyJnkO/2XcJwOER0QvlRAFMHKqxpMnAicDjQ7RZKGmdpHV9fX21jNXMrNCkfG4kHF2vA0t6EDiiwqbFIzzOQcBdwCci4uXB2kVEF9AF0NnZmUMuNjPLR0nK5T6QuiWQiJg92DZJz0tqi4heSW3AjkHajaGcPG6LiLvrFKqZWUNrtct47wEWJK8XACsGNlD5OY23At0R8fcZxmZm1lBa7UbCJcAcSZuBOckyktolrUzanAFcCZwr6anka24+4ZqZFZdymkSv2ymsoUTETmBWhfXbgbnJ6/8AlHFoZmYNp+RaWGZmlkarncIyM7MaabVJdDMzq5GWqoVlZma1Uyq1Vi0sMzOrEc+BmJlZKq1WC8vMzGokr1pYTiBmZg0ur1pYTiBmZg2uJNi9O4d+s+/SzMxqyZPoZmaWSl61sJxAzMwanGthmZlZKj6FZWZmqbgWlpmZpeJaWGZmlkp5DiSHfrPv0szMaqml5kAkHSrpO5I2J98PGaLtKElPSrovyxjNzBpFSyUQYBGwOiKmAauT5cFcB3RnEpWZWQNSi02izwOWJq+XAhdWaiSpA3g/cEs2YZmZNZ5Wq4V1eET0AiTfJw7S7kbgU8B+q7xIWihpnaR1fX19NQvUzKzoSqV8RiCj63VgSQ8CR1TYtHiY+38A2BERT0g6Z3/tI6IL6ALo7OzM4a00M8tHXnMgdUsgETF7sG2SnpfUFhG9ktqAHRWanQFcIGkuMBYYL+lbEfHBOoVsZtaQWq0W1j3AguT1AmDFwAYR8dcR0RERk4HLgTVOHmZm+2q1WlhLgDmSNgNzkmUktUtamVNMZmYNqelOYQ0lInYCsyqs3w7MrbB+LbC27oGZmTUgP1DKzMxSkUQeVw45gZiZNbhWmwMxM7MaabVSJmZmViOlFruM18zMaqRcC8sjEDMzG6FyLawc+s2+SzMzq6WSRyBmZpaGJ9HNzCwVSb6R0MzMRs73gZiZWSq+jNfMzFIpP1DKIxAzMxuhVnseiJmZ1YjnQMzMLBVfxmtmZqm01CS6pEMlfUfS5uT7IYO0O1jSnZJ+LKlb0mlZx2pmVnStVgtrEbA6IqYBq5PlSm4CVkXEscAJQHdG8ZmZNYxWq4U1D1iavF4KXDiwgaTxwFnArQAR8duIeCmj+MzMGkar1cI6PCJ6AZLvEyu0+SOgD/i6pCcl3SLpHYMdUNJCSeskrevr66tP1GZmBdR0k+iSHpS0ocLXvGEeYjRwEvDViDgR+DWDn+oiIroiojMiOidMmFCDn8DMrDHkdR/I6HodOCJmD7ZN0vOS2iKiV1IbsKNCsx6gJyIeT5bvZIgEYmbWqlrtPpB7gAXJ6wXAioENIuIXwFZJxySrZgGbsgnPzKxxtNRlvMASYI6kzcCcZBlJ7ZJW9mt3LXCbpKeBmcBnsw7UzKzo8ppEr9sprKFExE7KI4qB67cDc/stPwV0ZheZmVnjUXIZb0QgKbN+fSe6mVmDKyVJI+tBiBOImVmDKyWDjqxPYzmBmJk1uFKSQbKeSHcCMTNrcPIIxMzM0vAciJmZpeI5EDMzS2XPCMQJxMzMRkTyJLqZmaWw5xRW1vWwnEDMzBpcySMQMzNLw5PoZmaWijyJbmZmafg+EDMzS6XQp7AkXT1geZSkz9QnJDMzG4miT6LPkrRSUpuk44HHgHF1jMvMzIZpby2sjDPIsB4oFRF/Lum/AuuBV4H5EfFI2k4lHQr8CzAZ2AJcFhEvVmj3l8BfAJH0/eGIeC1tv2ZmzajQcyCSpgHXAXdR/oN/paQDq+h3EbA6IqYBq5PlgX0eCfwPoDMijgdGAZdX0aeZWVMqJX/JCzkHAtwLfDoi/jtwNvAT4AdV9DsPWJq8XgpcOEi70cDbJY0GDgS2V9GnmVlTKnotrJOBEyTdDdxJ+ZRSNaOBwyOiFyD5PnFgg4jYBnwR+DnQC/wqIv59sANKWihpnaR1fX19VYRmZtZYil4L6xbgOODLwFeA6cCnh9pB0oOSNlT4mjecDiUdQnmkMgVoB94h6YODtY+IrojojIjOCRMmDPPHMjNrfHnVwhrWJDpwTESc0G/5IUk/GmqHiJg92DZJz0tqi4heSW3AjgrNZgPPRURfss/dwOnAt4YZs5lZSyj6ZbxPSjp1z4KkU4DUV2EB9wALktcLgBUV2vwcOFXSgSqPz2YB3VX0aWbWlAp9IyFwCvCopC2StgDfA86WtF7S0yn6XQLMkbQZmJMsI6ld0kqAiHic8nzLDylfwlsCulL0ZWbW1PKqhTXcU1jn1bLTiNhJeUQxcP12YG6/5c8AvuPdzGwIed0HMtwbCX9W70DMzCydop/CMjOzgir6JLqZmRWUPAIxM7M03poDcQIxM7MR8CksMzNLpZRTOXcnEDOzBlf0WlhmZlZQedXCcgIxM2twpZJHIGZmloJvJDQzs1TyqoXlBGJm1uAK/Ux0MzMrLp/CMjOzVHwjoZmZpeJaWGZmlkpL1cKSdKmkjZJ2S+ocot15kp6R9KykRVnGaGbWKFrtFNYG4GLg4cEaSBoF3AycD8wA5kuakU14ZmaNI69J9OE+0ramIqIb3rp2eRAnA89GxE+TtsuAecCmugdoZtZAXAtrX0cCW/st9yTrzMysn7xqYdVtBCLpQeCICpsWR8SK4RyiwrpB3x1JC4GFAJMmTRpWjGZmzaCU053odUsgETG7ykP0AEf1W+4Atg/RXxfQBdDZ2ZnxQM7MLD97E8jujPvNtrsR+QEwTdIUSW8DLgfuyTkmM7PCaan7QCRdJKkHOA24X9IDyfp2SSsBIuJ3wMeBB4Bu4F8jYmMe8ZqZFdmecu5Z18LK6yqs5cDyCuu3A3P7La8EVmYYmplZw3EtLDMzS6XVbiQ0M7Maaak5EDMzq52WqoVlZma141NYZmaWiifRzcwsFdfCMjOzVPKqheUEYmbW4PKqheUEYmbW4DyJbmZmqfg+EDMzS+Wt+0Ay7jfb7szMrNb2Xsab8TksJxAzswbnORAzM0vFcyBmZpaKJCTfB2JmZimUJJ/CMjOzkRMtcgpL0qWSNkraLalzkDZHSXpIUnfS9rqs4zQzaxStNALZAFwMPDxEm98BfxUR04FTgY9JmpFFcGZmjSaPOZC8noneDW9VkBykTS/Qm7x+RVI3cCSwKYsYzcwaSXkE0gKnsEZK0mTgRODxIdoslLRO0rq+vr7MYjMzK4KSsr8PpG4jEEkPAkdU2LQ4IlaM4DgHAXcBn4iIlwdrFxFdQBdAZ2dnxm+jmVm+8hiB1C2BRMTsao8haQzl5HFbRNxdfVRmZs2pPAeSbZ+5zIEMh8oTJLcC3RHx99Ue74033qCnp4fXXnut+uAKZOzYsXR0dDBmzJi8QzGzHJVKTTQCGYqki4AvAxOA+yU9FRHvk9QO3BIRc4EzgCuB9ZKeSnb9XxGxMk2fPT09jBs3jsmTJw85ed9IIoKdO3fS09PDlClT8g7HzHLUVKewhhIRy4HlFdZvB+Ymr/+D8r0xNfHaa681VfKA8lVs73znO/FFA2aWxyR6Q1yFVSvNlDz2aMafycxGTpJrYZmZ2ciVBLt3Z9xntt01r9NPP73i+quuuoo777wz42jMrNX4RsIG9uijj+Ydgpm1sDxqYRX2Mt5Gc9BBB7Fr1y4igmuvvZY1a9YwZcqUzM9Jmllr8vNAmsDy5ct55plnWL9+PV/72tc8MjGzTPgUVhN4+OGHmT9/PqNGjaK9vZ1zzz0375DMrAX4Mt4m4UtrzSxrJYmsT5g7gdTYWWedxbJly3jzzTfp7e3loYceyjskM2sBUvZPJPQkeo1ddNFFrFmzhne/+90cffTRnH322XmHZGYtoJTDjYROIDWya9cuoHz66itf+UrO0ZhZqylJvpHQzMxGLo9TWE4gZmZNII8bCZ1AzMyaQKnkGwnNzCwF30hoZmapqFVOYUm6VNJGSbslde6n7ShJT0q6L6v4zMwaTamFJtE3ABcDDw+j7XVAd33Dqb+tW7fy3ve+l+nTp3Pcccdx00035R2SmTWR8n0gGfeZbXdlEdEdEc/sr52kDuD9wC31j6q+Ro8ezZe+9CW6u7t57LHHuPnmm9m0aVPeYZlZk8hjBFL0GwlvBD4FjNtfQ0kLgYUAkyZNGrLtDfduZNP2l2sQ3ltmtI/nM3963KDb29raaGtrA2DcuHFMnz6dbdu2MWPGjJrGYWatSc00iS7pQUkbKnzNG+b+HwB2RMQTw2kfEV0R0RkRnRMmTKgq9nrbsmULTz75JKecckreoZhZk8ijGm/dRiARMbvKQ5wBXCBpLjAWGC/pWxHxwWpjG2qkUG+7du3ikksu4cYbb2T8+PG5xWFmzaUk8WbGtUwKexlvRPx1RHRExGTgcmBNLZJHnt544w0uueQSrrjiCi6++OK8wzGzJtIyd6JLukhSD3AacL+kB5L17ZJW5hFTvUUEV199NdOnT+f666/POxwzazItUwsrIpYno4sDIuLwiHhfsn57RMyt0H5tRHwg+0hr55FHHuGb3/wma9asYebMmcycOZOVK5syV5pZDvIYgRT9KqymceaZZ2Zep8bMWkdJroVlZmYpuBaWmZmlIj9QyszM0milWlhmZlZDLVMLy8zMaqtU8gjEzMxSaKpaWLavj3zkI0ycOJHjjz9+77pf/vKXzJkzh2nTpjFnzhxefPHFvds+97nPMXXqVI455hgeeOCBPEI2swbhU1hN7qqrrmLVqlW/t27JkiXMmjWLzZs3M2vWLJYsWQLApk2bWLZsGRs3bmTVqlV89KMf5c0338wjbDNrAC7nnpV/WwS/WF/bYx7xbjh/yZBNzjrrLLZs2fJ761asWMHatWsBWLBgAeeccw6f//znWbFiBZdffjkHHHAAU6ZMYerUqXz/+9/ntNNOq23cZtYUWqYWlr3l+eef3/uckLa2Nnbs2AHAtm3bOOqoo/a26+joYNu2bbnEaGbFl0ctrNYcgexnpFAElUoSSMohEjNrBJ4DaUGHH344vb29APT29jJx4kSgPOLYunXr3nY9PT20t7fnEqOZFZ9vJGxBF1xwAUuXLgVg6dKlzJs3b+/6ZcuW8frrr/Pcc8+xefNmTj755DxDNbMCy6MWVmuewsrJ/PnzWbt2LS+88AIdHR3ccMMNLFq0iMsuu4xbb72VSZMmcccddwBw3HHHcdlllzFjxgxGjx7NzTffzKhRo3L+CcysqORy7s3t9ttvr7h+9erVFdcvXryYxYsX1zMkM2sS7+n4g9Yo5y7pUkkbJe2W1DlEu4Ml3Snpx5K6JfkaVjOzCuafPIkll7wn0z7zmgPZAFwMPLyfdjcBqyLiWOAEoLvegZmZ2fDkcgorIrph6MtSJY0HzgKuSvb5LfDbKvttukth/ZRDM8tLka/C+iOgD/i6pCcl3SLpHYM1lrRQ0jpJ6/r6+vbZPnbsWHbu3NlUf3Ajgp07dzJ27Ni8QzGzFlS3EYikB4EjKmxaHBErhnGI0cBJwLUR8bikm4BFwKcrNY6ILqALoLOzc58s0dHRQU9PD5WSSyMbO3YsHR0deYdhZi2obgkkImZXeYgeoCciHk+W76ScQFIZM2YMU6ZMqTIkMzPbo7CnsCLiF8BWScckq2YBm3IMyczM+snrMt6LJPUApwH3S3ogWd8uaWW/ptcCt0l6GpgJfDbzYM3MrKK8rsJaDiyvsH47MLff8lPAoPeJmJlZftRMVyXtIakP+FnK3Q8DXqhhOFlwzNlotJgbLV5wzFmpFPMfRsSEkRykKRNINSSti4iGGvU45mw0WsyNFi845qzUKubCTqKbmVmxOYGYmVkqTiD76so7gBQcczYaLeZGixccc1ZqErPnQMzMLBWPQMzMLBUnEDMzS6VlEoik8yQ9I+lZSfvU1FLZPyTbn5Z00nD3LVrMko6S9FDyEK6Nkq4resz9to9Kqi/f1wgx5/XQsypj/svk92KDpNslZVLOeRgxHyvpe5Jel/TJkexbtJgL/hkc9H1Otg//MxgRTf8FjAL+k3KJ+LcBPwJmDGgzF/g3QMCpwOPD3beAMbcBJyWvxwE/KXrM/bZfD/w/4L6i/24k25YCf5G8fhtwcJFjBo4EngPeniz/K3BVQWKeCPwx8LfAJ0eybwFjLvJnsGLM/bYP+zPYKiOQk4FnI+KnUX4w1TJg3oA284BvRNljwMGS2oa5b6FijojeiPghQES8QvlJjkcWOWYASR3A+4FbMoi16pj11kPPboXyQ88i4qUix5xsGw28XdJo4EBgexFijogdEfED4I2R7lu0mIv8GRzifR7xZ7BVEsiRwNZ+yz3s+485WJvh7FsP1cS8l6TJwInA49RftTHfCHwK2F2n+CqpJuYRPfSshlLHHBHbgC8CPwd6gV9FxL/XMdYh48lg32rUpN8CfgaHciMj+Ay2SgKp9BzbgdcvD9ZmOPvWQzUxlzdKBwF3AZ+IiJdrGNtgUscs6QPAjoh4ovZhDama93nPQ8++GhEnAr+mimfWjEA17/MhlP9HOgVoB94h6YM1jq+Saj5HRf4MDn2AYn4GK++Y4jPYKgmkBziq33IH+w7bB2sznH3roZqYkTSG8i/ubRFxdx3jHFY8w2hzBnCBpC2Uh93nSvpW/ULdbzzDaVPpoWcnUX/VxDwbeC4i+iLiDeBu4PQ6xrq/eOq9bzWq6rfAn8HBjPwzWO9JnSJ8Uf6f4k8p/69rz8TScQPavJ/fn3T8/nD3LWDMAr4B3Ngo7/OANueQ3SR6VTED3wWOSV7/DfCFIscMnAJspDz3IcoXAVxbhJj7tf0bfn9CurCfwSFiLuxncLCYB2wb1mcwsx8s7y/KV6X8hPIVCouTddcA1/T7B7852b4e6Bxq3yLHDJxJedj6NPBU8jW3yDEPOMawfnmLEDPlB52tS97rbwOHNEDMNwA/BjYA3wQOKEjMR1D+H/TLwEvJ6/GD7VvkmAv+GRz0fe53jGF9Bl3KxMzMUmmVORAzM6sxJxAzM0vFCcTMzFJxAjEzs1ScQMzMLBUnELMhJNV2P9pvuV3SnXXq60JJ/3s/bb4o6dx69G82Ur6M12wISR2j+yLi+Az6ehS4ICJeGKLNHwJfi4g/qXc8ZvvjEYjZ0JYA75L0lKQvSJosaQOApKskfVvSvZKek/RxSdcnhRUfk3Ro0u5dklZJekLSdyUdO7ATSUcDr0fEC5LGJccbk2wbL2mLpDER8TPgnZKOyPA9MKvICcRsaIuA/4yImRHxPytsPx74c8pltP8WeDXKhRW/B3woadNFuVzIfwE+CfxjheOcAfQv/72WcjkSgMuBu6Jcu4qk3RlV/lxmVRuddwBmDe6h5A/+K5J+BdybrF8PvCepxno6cIe0t1DqARWO00a5NPwet1Auq/1t4MPAf+u3bQflSrpmuXICMavO6/1e7+63vJvy56sEvBQRM/dznN8Af7BnISIeSU6XnQ2MiogN/dqOTdqb5cqnsMyG9grlR5KmEuVnQDwn6VLY+6zyEyo07QamDlj3DeB24OsD1h9NuRCiWa6cQMyGEBE7gUckbZD0hZSHuQK4WtKPKJdSr/Q41oeBE9XvPBdwG3AI5SQC7H3GxFTKFYDNcuXLeM0KQtJNwL0R8WCy/GfAvIi4sl+bi4CTIuLTOYVptpfnQMyK47OUH/iEpC8D51N+tkN/o4EvZRyXWUUegZiZWSqeAzEzs1ScQMzMLBUnEDMzS8UJxMzMUnECMTOzVP4/uuPo9YaKC4oAAAAASUVORK5CYII=\n", + "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZgAAAEGCAYAAABYV4NmAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjMuNCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8QVMy6AAAACXBIWXMAAAsTAAALEwEAmpwYAAAWF0lEQVR4nO3dfZBddZ3n8fd3kkCGIchjoEMHkzGBSXhYjL0hoIUKZipEN1GZscg4Q/CJQsSHZVg3M9bujrU1mirHXXTNSAWRShxnUg7Kg1aAiYCLhRMkPAiEGJMBJB1aiFGQrMuj3/3j3mRvOreT2337d89N8n5V3ep7zu/3O+fbN33y6d85t8+NzESSpNH2e1UXIEk6MBkwkqQiDBhJUhEGjCSpCANGklTE2KoL6KRjjz02p0yZUnUZkrRfuf/++3+ZmccNd9xBFTBTpkxh3bp1VZchSfuViPj5SMZ5ikySVIQBI0kqwoCRJBVxUF2DkaS9eeWVV+jv7+fFF1+supRKjB8/nt7eXsaNGzcq2zNgJKmuv7+fCRMmMGXKFCKi6nI6KjPZvn07/f39TJ06dVS26SkySap78cUXOeaYYw66cAGICI455phRnb0ZMJLU4GAMl51G+3s3YCRJRRgwklShc845p+n6Sy65hBtuuKHD1YwuA0aSKvSjH/2o6hKK8V1kklShww8/nB07dpCZfPzjH+fOO+9k6tSpHAifNuwMRpK6wI033sjGjRt55JFHuPbaaw+ImY0BI0ld4O6772bRokWMGTOGSZMmcd5551VdUtsMGEnqEgfaW6QNGEnqAueeey6rVq3itddeY2BggLvuuqvqktrmRX5J6gLvec97uPPOOzn99NM5+eSTeetb31p1SW0zYCSpQjt27ABqp8e+8pWvVFzN6PIUmSSpCANGklSEASNJKsKAkSQVYcBIkoowYCRJRRgwktRFtmzZwtvf/nZmzJjBqaeeype+9KU9+mQmn/jEJ5g2bRpnnHEGDzzwQAWV7pt/ByNJXWTs2LF88YtfZNasWbzwwgu86U1vYu7cucycOXNXn1tvvZVNmzaxadMm7r33Xj760Y9y7733Vlh1c5XOYCJiXkRsjIjNEbGkSXtExJfr7Q9HxKxB7WMi4sGI+F7nqpakcnp6epg1q/Zf3YQJE5gxYwZbt27drc/NN9/MxRdfTEQwZ84cnnvuOQYGBqood68qm8FExBhgGTAX6Afui4hbMvOxhm4XANPrj7OAr9a/7vRJYANwREeKlnTQ+Ox31/PY078Z1W3OnHQE/+0/nNpy/yeffJIHH3yQs846a7f1W7duZfLkybuWe3t72bp1Kz09PaNW62iocgYzG9icmY9n5svAKmDhoD4LgZVZsxY4MiJ6ACKiF3gn8LVOFi1JnbBjxw4uvPBCrr76ao44YvffoZt9GFk33om5ymswJwJbGpb72X12MlSfE4EB4Grg08CEve0kIi4FLgU46aST2ipY0sFjODON0fbKK69w4YUX8v73v5/3vve9e7T39vayZcv//6+xv7+fSZMmdbLEllQ5g2kWt4NjuWmfiHgX8Gxm3r+vnWTm8szsy8y+4447biR1SlLHZCYf+tCHmDFjBldeeWXTPgsWLGDlypVkJmvXruV1r3td150eg2pnMP3A5IblXuDpFvv8CbAgIuYD44EjIuIfMvPPC9YrScXdc889fOMb3+D000/nzDPPBOBzn/scTz31FACXXXYZ8+fPZ/Xq1UybNo3DDjuM66+/vsKKh1ZlwNwHTI+IqcBW4CLgzwb1uQW4IiJWUTt99nxmDgB/VX8QEW8DrjJcJB0I3vKWtzS9xtIoIli2bFmHKhq5ygImM1+NiCuA24ExwNczc31EXFZvvwZYDcwHNgO/BT5QVb2SpOGp9A8tM3M1tRBpXHdNw/MEPraPbfwA+EGB8iRJbfBWMZKkIgwYSVIRBowkqQgDRpJUhAEjSV3kgx/8IBMnTuS0007bte5Xv/oVc+fOZfr06cydO5df//rXu9o+//nPM23aNE455RRuv/32ptvc2/iSDBhJ6iKXXHIJt912227rli5dyvnnn8+mTZs4//zzWbp0KQCPPfYYq1atYv369dx2221cfvnlvPbaa3tsc6jxpRkwktRFzj33XI4++ujd1t18880sXrwYgMWLF3PTTTftWn/RRRdx6KGHMnXqVKZNm8aPf/zjPbY51PjS/MAxSWrm1iXwi0dGd5snnA4XDH/28Mwzz+y611hPTw/PPvssULtt/5w5c3b123nb/lbHl+YMRpL2U91+235nMJLUzAhmGqUcf/zxDAwM0NPTw8DAABMnTgRav23/UONLcwYjSV1uwYIFrFixAoAVK1awcOHCXetXrVrFSy+9xBNPPMGmTZuYPXt2y+NLM2AkqYssWrSIs88+m40bN9Lb28t1113HkiVLWLNmDdOnT2fNmjUsWbIEgFNPPZX3ve99zJw5k3nz5rFs2TLGjBkDwIc//GHWrVsHMOT40mJft4U+kPT19eXOF1ySBtuwYQMzZsyouoxKNXsNIuL+zOwb7racwUiSijBgJElFGDCS1OBgumww2Gh/7waMJNWNHz+e7du3H5Qhk5ls376d8ePHj9o2/TsYSarr7e2lv7+fbdu2VV1KJcaPH09vb++obc+AkaS6cePGMXXq1KrLOGB4ikySVIQBI0kqwoCRJBVhwEiSijBgJElFGDCSpCIMGElSEQaMJKkIA0aSVIQBI0kqwoCRJBVRacBExLyI2BgRmyNij8/wjJov19sfjohZ9fWTI+KuiNgQEesj4pOdr16StDeVBUxEjAGWARcAM4FFETFzULcLgOn1x6XAV+vrXwX+MjNnAHOAjzUZK0mqUJUzmNnA5sx8PDNfBlYBCwf1WQiszJq1wJER0ZOZA5n5AEBmvgBsAE7sZPGSpL2rMmBOBLY0LPezZ0jss09ETAHeCNw7+iVKkkaqyoCJJusGf4zcXvtExOHAt4FPZeZvmu4k4tKIWBcR6w7WDxGSpCpUGTD9wOSG5V7g6Vb7RMQ4auHyzcz8zlA7yczlmdmXmX3HHXfcqBQuSdq3KgPmPmB6REyNiEOAi4BbBvW5Bbi4/m6yOcDzmTkQEQFcB2zIzP/R2bIlSa2o7COTM/PViLgCuB0YA3w9M9dHxGX19muA1cB8YDPwW+AD9eFvBv4CeCQiHqqv++vMXN3Bb0GStBeROfiyx4Grr68v161bV3UZkrRfiYj7M7NvuOP8S35JUhEGjCSpCANGklSEASNJKsKAkSQVYcBIkoowYCRJRRgwkqQiDBhJUhEGjCSpCANGklSEASNJKsKAkSQVYcBIkoowYCRJRRgwkqQiDBhJUhEGjCSpCANGklSEASNJKsKAkSQVYcBIkoowYCRJRRgwkqQiDBhJUhEGjCSpCANGklSEASNJKsKAkSQVYcBIkoowYCRJRewzYCJiYpN1p4zGziNiXkRsjIjNEbGkSXtExJfr7Q9HxKxWx0qSqtXKDOaHEfG+nQsR8ZfAje3uOCLGAMuAC4CZwKKImDmo2wXA9PrjUuCrwxgrSarQ2Bb6vA1YHhF/ChwPbABmj8K+ZwObM/NxgIhYBSwEHmvosxBYmZkJrI2IIyOiB5jSwthRs/bvP8KE5zaU2LQkdcQLR85gzuXXdnSf+5zBZOYAcBtwNrX/2Fdm5o5R2PeJwJaG5f76ulb6tDIWgIi4NCLWRcS6bdu2tV20JKk1+5zBRMQaYAA4DegFvh4Rd2fmVW3uO5qsyxb7tDK2tjJzObAcoK+vr2mffel06kvSgaCVazC3An+dmc9l5qPAOcDzo7DvfmByw3Iv8HSLfVoZK0mqUCsBMwG4PSJ+GBEfA47JzP8+Cvu+D5geEVMj4hDgIuCWQX1uAS6uv5tsDvB8/ZRdK2MlSRVq5RrMZzPzVOBjwCTgf0fE99vdcWa+ClwB3E7tjQPfysz1EXFZRFxW77YaeBzYDFwLXL63se3WJEkaPa28i2ynZ4FfANuBPf42ZiQyczW1EGlcd03D86QWbC2NlSR1j1b+0PKjEfED4A7gWOAjmXlG6cIkSfu3VmYwrwc+lZkPFa5FknQA2WfAZKa3YZEkDZs3u5QkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUVUEjARcXRErImITfWvRw3Rb15EbIyIzRGxpGH9FyLipxHxcETcGBFHdqx4SVJLqprBLAHuyMzpwB315d1ExBhgGXABMBNYFBEz681rgNMy8wzgZ8BfdaRqSVLLqgqYhcCK+vMVwLub9JkNbM7MxzPzZWBVfRyZ+S+Z+Wq931qgt2y5kqThqipgjs/MAYD614lN+pwIbGlY7q+vG+yDwK2jXqEkqS1jS204Ir4PnNCk6TOtbqLJuhy0j88ArwLf3EsdlwKXApx00kkt7lqS1K5iAZOZ7xiqLSKeiYiezByIiB7g2Sbd+oHJDcu9wNMN21gMvAs4PzOTIWTmcmA5QF9f35D9JEmjq6pTZLcAi+vPFwM3N+lzHzA9IqZGxCHARfVxRMQ84D8DCzLztx2oV5I0TFUFzFJgbkRsAubWl4mISRGxGqB+Ef8K4HZgA/CtzFxfH/8VYAKwJiIeiohrOv0NSJL2rtgpsr3JzO3A+U3WPw3Mb1heDaxu0m9a0QIlSW3zL/klSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFWHASJKKMGAkSUUYMJKkIgwYSVIRBowkqQgDRpJUhAEjSSrCgJEkFVFJwETE0RGxJiI21b8eNUS/eRGxMSI2R8SSJu1XRURGxLHlq5YkDUdVM5glwB2ZOR24o768m4gYAywDLgBmAosiYmZD+2RgLvBURyqWJA1LVQGzEFhRf74CeHeTPrOBzZn5eGa+DKyqj9vpfwKfBrJgnZKkEaoqYI7PzAGA+teJTfqcCGxpWO6vryMiFgBbM/Mn+9pRRFwaEesiYt22bdvar1yS1JKxpTYcEd8HTmjS9JlWN9FkXUbEYfVt/HErG8nM5cBygL6+Pmc7ktQhxQImM98xVFtEPBMRPZk5EBE9wLNNuvUDkxuWe4GngTcAU4GfRMTO9Q9ExOzM/MWofQOSpLZUdYrsFmBx/fli4OYmfe4DpkfE1Ig4BLgIuCUzH8nMiZk5JTOnUAuiWYaLJHWXqgJmKTA3IjZReyfYUoCImBQRqwEy81XgCuB2YAPwrcxcX1G9kqRhKnaKbG8ycztwfpP1TwPzG5ZXA6v3sa0po12fJKl9/iW/JKkIA0aSVIQBI0kqwoCRJBVhwEiSijBgJElFGDCSpCIMGElSEQaMJKkIA0aSVIQBI0kqwoCRJBVhwEiSijBgJElFGDCSpCIMGElSEQaMJKkIA0aSVIQBI0kqwoCRJBVhwEiSijBgJElFGDCSpCIMGElSEZGZVdfQMRGxDfj5CIcfC/xyFMvpBGvunP2xbmvujAOh5tdn5nHD3chBFTDtiIh1mdlXdR3DYc2dsz/Wbc2dcTDX7CkySVIRBowkqQgDpnXLqy5gBKy5c/bHuq25Mw7amr0GI0kqwhmMJKkIA0aSVIQBA0TEvIjYGBGbI2JJk/aIiC/X2x+OiFmtju22miNickTcFREbImJ9RHyy22tuaB8TEQ9GxPf2h5oj4siIuCEiflp/vc/eD2r+j/Wfi0cj4p8iYnyX1PxHEfGvEfFSRFw1nLHdVnOVx2A7dTe0t34cZuZB/QDGAP8G/CFwCPATYOagPvOBW4EA5gD3tjq2C2vuAWbVn08AftbtNTe0Xwn8I/C9bv/ZqLetAD5cf34IcGQ31wycCDwB/H59+VvAJV1S80Tg3wN/C1w1nLFdWHMlx2C7dTe0t3wcOoOB2cDmzHw8M18GVgELB/VZCKzMmrXAkRHR0+LYrqo5Mwcy8wGAzHwB2EDtP5aurRkgInqBdwJf60CtbdccEUcA5wLXAWTmy5n5XDfXXG8bC/x+RIwFDgOe7oaaM/PZzLwPeGW4Y7ut5gqPQWjvtR72cWjA1P5htzQs97PnP/ZQfVoZW0I7Ne8SEVOANwL3jn6Je2i35quBTwO/K1RfM+3U/IfANuD6+umEr0XEH5Qsdh/17LNPZm4F/g54ChgAns/MfylY617r6cDYdozKfjt8DEL7dV/NMI5DA6Z2mmCwwe/dHqpPK2NLaKfmWmPE4cC3gU9l5m9GsbahjLjmiHgX8Gxm3j/6Ze1VO6/zWGAW8NXMfCPwf4BOXB9o53U+itpvs1OBScAfRMSfj3J9zbRzHHXzMbj3DXT+GIQ26h7JcWjA1BJ8csNyL3ueFhiqTytjS2inZiJiHLUf7G9m5ncK1tlSPS30eTOwICKepDalPy8i/qFcqfusp5U+/UB/Zu78zfQGaoFTWjs1vwN4IjO3ZeYrwHeAcwrWuq96So9tR1v7regYhPbqHv5x2IkLS938oPab5uPUfmvbedHr1EF93snuF0V/3OrYLqw5gJXA1fvL6zyoz9vo3EX+tmoGfgicUn/+N8AXurlm4CxgPbVrL0HtTQof74aaG/r+DbtfMO/aY3AvNVdyDLZb96C2lo7Djn5z3fqg9q6an1F7d8Vn6usuAy5r+IFYVm9/BOjb29hurhl4C7Up8cPAQ/XH/G6uedA2WvrB7oaagTOBdfXX+ibgqP2g5s8CPwUeBb4BHNolNZ9A7bfv3wDP1Z8fMdTYbq65ymOw3de6YRstHYfeKkaSVITXYCRJRRgwkqQiDBhJUhEGjCSpCANGklSEASONUP1uyZc3LE+KiBsK7evdEfFf99Hn7yLivBL7l0bCtylLI1S/j9T3MvO0DuzrR8CCzPzlXvq8Hrg2M/+4dD1SK5zBSCO3FHhDRDwUEV+IiCkR8ShARFwSETdFxHcj4omIuCIirqzf+HJtRBxd7/eGiLgtIu6PiB9GxB8N3klEnAy8lJm/jIgJ9e2Nq7cdERFPRsS4zPw5cExEnNDB10AakgEjjdwS4N8y88zM/E9N2k8D/ozaLdL/Fvht1m58+a/AxfU+y6ndjuVNwFXA3zfZzpuBxtu7/4Da7V4ALgK+nbV7h1Hv9+Y2vy9pVIytugDpAHZXPRBeiIjnge/W1z8CnFG/m+45wD9H7LrJ7aFNttND7db/O32N2i3TbwI+AHykoe1ZandClipnwEjlvNTw/HcNy7+jduz9HvBcZp65j+38X+B1Oxcy85766bi3AmMy89GGvuPr/aXKeYpMGrkXqH3k7Yhk7TNAnoiIPwWImn/XpOsGYNqgdSuBfwKuH7T+ZGo3qpQqZ8BII5SZ24F7IuLRiPjCCDfzfuBDEfETarfKb/Zxv3cDb4yG82jAN4GjqIUMsOszRqZRu4OzVDnfpiztByLiS8B3M/P79eU/ARZm5l809HkPMCsz/0tFZUq78RqMtH/4HLUPBCMi/hdwAbXP9Wg0Fvhih+uShuQMRpJUhNdgJElFGDCSpCIMGElSEQaMJKkIA0aSVMT/AxawhCPk7KksAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] @@ -103,7 +103,7 @@ } ], "source": [ - "swiftdiff['px'].plot.line(x=\"time (y)\")" + "swiftdiff['vx'].plot.line(x=\"time (y)\")" ] }, { @@ -465,41 +465,33 @@ " stroke: currentColor;\n", " fill: currentColor;\n", "}\n", - "
<xarray.DataArray 'vx' (time (y): 199)>\n",
-       "array([ 0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,  0.,\n",
-       "        0.,  0.,  0., nan])\n",
+       "
<xarray.DataArray 'vx' (time (y): 198)>\n",
+       "array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n",
+       "       0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])\n",
        "Coordinates:\n",
-       "    id        int64 100\n",
-       "  * time (y)  (time (y)) float64 0.0 0.0006845 0.001369 ... 0.1342 0.1348 0.1355
  • " ], "text/plain": [ - "\n", - "array([ 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", - " 0., 0., 0., nan])\n", + "\n", + "array([0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,\n", + " 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.])\n", "Coordinates:\n", - " id int64 100\n", - " * time (y) (time (y)) float64 0.0 0.0006845 0.001369 ... 0.1342 0.1348 0.1355" + " id float64 100.0\n", + " * time (y) (time (y)) float64 0.0 0.0006845 0.001369 ... 0.1335 0.1342 0.1348" ] }, "execution_count": 7,