diff --git a/analysis/tritium/multi_block_model.ipynb b/analysis/tritium/multi_block_model.ipynb new file mode 100644 index 0000000..4185850 --- /dev/null +++ b/analysis/tritium/multi_block_model.ipynb @@ -0,0 +1,23354 @@ +{ + "cells": [ + { + "cell_type": "code", + "execution_count": 1, + "metadata": {}, + "outputs": [ + { + "name": "stderr", + "output_type": "stream", + "text": [ + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-1-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-1-2 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-1-4 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-3-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-5-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-5-2 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-6-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-6-2 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-7-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-8-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-8-2 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-10-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-11-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-13-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-14-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-15-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-IV_3-16-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-1-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-1-3 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-1-4 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-2-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-3-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-3-4 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-4-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-5-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-7-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-8-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-9-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-10-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-11-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "/home/remidm/anaconda3/envs/baby-1l-run-3/lib/python3.12/site-packages/libra_toolbox/tritium/lsc_measurements.py:140: UserWarning: Activity of 1L-OV_3-12-1 is negative after substracting background. Setting to zero.\n", + " warnings.warn(\n", + "2025-04-30 21:33:23,857 - INFO - LOGGING (log: True)\n", + "2025-04-30 21:33:23,858 - INFO - SOLVER -> RKBS32 (adaptive: True, implicit: False)\n", + "2025-04-30 21:33:23,859 - INFO - ALGEBRAIC PATH DFS (path_length: 1)\n", + "2025-04-30 21:33:23,861 - INFO - SET 'iterations_min' -> 1\n", + "2025-04-30 21:33:23,861 - INFO - STARTING -> TRANSIENT (Duration: 4320000.00s)\n", + "2025-04-30 21:33:23,862 - INFO - TRANSIENT: 0% | elapsed: 00:00:00 (eta: --:--:--) | 0 steps (N/A steps/s)\n", + "2025-04-30 21:33:23,882 - INFO - TRANSIENT: 21% | elapsed: 00:00:00 (eta: --:--:--) | 55 steps (2830.6 steps/s)\n", + "2025-04-30 21:33:23,895 - INFO - TRANSIENT: 40% | elapsed: 00:00:00 (eta: --:--:--) | 96 steps (3138.7 steps/s)\n", + "2025-04-30 21:33:23,943 - INFO - TRANSIENT: 60% | elapsed: 00:00:00 (eta: --:--:--) | 252 steps (3205.6 steps/s)\n", + "2025-04-30 21:33:23,968 - INFO - TRANSIENT: 82% | elapsed: 00:00:00 (eta: 00:00:00) | 336 steps (3362.6 steps/s)\n", + "2025-04-30 21:33:23,972 - INFO - TRANSIENT: 100% | elapsed: 00:00:00 (eta: 00:00:00) | 345 steps (2329.2 steps/s)\n", + "2025-04-30 21:33:23,973 - INFO - TRANSIENT: 100% | elapsed: 00:00:00 (eta: 00:00:00) | 345 steps (3095.3 avg steps/s)\n", + "2025-04-30 21:33:23,974 - INFO - FINISHED -> TRANSIENT (total steps: 345, successful: 323, runtime: 111.45 ms)\n" + ] + }, + { + "name": "stdout", + "output_type": "stream", + "text": [ + "Processed data stored in ../../data/processed_data.json\n" + ] + } + ], + "source": [ + "from multi_system_tritium_model import *\n", + "\n", + "Sim.run(total_duration)\n", + "\n", + "# plot total plant inventory\n", + "time, data = Sco.read()\n", + "\n", + "s_to_days = 1 / (24 * 3600)" + ] + }, + { + "cell_type": "code", + "execution_count": 2, + "metadata": {}, + "outputs": [ + { + "data": { + "text/html": [ + " \n", + " \n", + " " + ] + }, + "metadata": {}, + "output_type": "display_data" + }, + { + "data": { + "application/vnd.plotly.v1+json": { + "config": { + "plotlyServerURL": "https://plot.ly" + }, + "data": [ + { + "mode": "lines", + "name": "salt", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAAAzYMVAX9BOPgW2WNA2w2Y+N0w8mk2rlT5aDIUjCj3KPonCXuiUWAA/hea79AdtND+SjjOsE4hpPwGIic4R6p8/VYaKtEry0z+yuRtU3e4IQHN1cKeUKj9AsAeL6Jx6c0CE46giRFmoQCLaPyuVb95A0imSOr0FE0HSOneDLMdHQVGme1z3uH1Brbu+2JiTskE4wMchaTjnQd48QF1xBR1Coj9p2LccUkJo4VuMz+SBQg+sSuwDHZhCdW5HJE7mn0LSU0ReYrydQlImYngei5lC1xOk4z/ulULOluXhKtSSQiJQL45mKpBCrtnUVSDCi0IAgn5BCNWHQvXA5GQZdoRCvSNzRjiRgUJX7Z34Qip+QohdRGQF5nlCor8qf0Q8dkKQZyuAJhdzQru/N87oY3BC25jQMOAkbEKiag6E0CloQnxa8MrjvmRCxaHgA7fPYUIfwCS2kpVeQhBT64MnQlpCeKhWll6LVkIDnYJoEFtTQkceQrc6AlJCdxfv1xTpTkI5a7Fv04lKQolgh1HnyEZCjK6bEOWPQ0LAeDRdk8tAQsJBxA3i1jxCVX5n9KTCOEK34I1zGkI1QpDLH5heQDJCvj1gswNXL0JdrCM7POgqQnIWqbf1GSdCIkp1hHzVI0L3eHf9UgchQgNgvXZ6PR1C4P1xg7oaGULYJtqWuo0VQtgArmxMgRJCd6J9XIHGD0IbiuLx9EcLQsrU53gkbAdCvM9OigscBEK95IIr50MBQpo9p9p/pf1BOSbxbQl0+UFka1PDZ9r1QXpb9Dwhw/JB+PLm0sUb8EHymn8vAqnrQdYy9pZ3v+dByI3hkpVj5EG2UX7bUoHhQeQr6kv3Dt5Be4ynDpbO2UFNLBC2JSjWQZRgqj7gBdNBDK4mCRRV0EHGJWWvaAvMQcmd5iHzE8hBkDckGx6sxEGGVRMMmb/BQeuH+e7leb5BgX/Vz2QqukFAt3M5+Ha2Qf7Q8rKMSbNB6/FJHC6PsEFLBNM9LW+sQYiJPzibaahBGjJ2rKj1pEHT55LGvP6hQbDsk/pQ5p5BVCVwK3qHmkESfF0l48aWQepJhOYpjpNB5FON4RbKkEHqhhq4VNSMQUrbLgd0wIhBY7jL3DhAhUF8RxsfwT6CQcNEBLg9VH9BQllZq9rlekFdAllf6hd3QXRy0jG703NB3mpBONEFcUFgK9sM5DptQWIKvsqBGGlBY5zaTtKLZUEEZ740qX9iQe/DY4Oxw19Bsk+Y6YpFW0FMi87aEWpXQUi7N/lDGlRByaXuCWBCUUGiJUA84KJNQdTzBs7IcUlBvaFHsnjYRUF/3KgxeMFCQSBh7uVYGkBBzfqSkI+mO0GYVTSZXb03QeGpH63HYTRBzpp5SsZ/MUGRtD5YTgwuQblaaWtNzClBj2jUwy8mJkG8SUlLMQQjQYMv9wmiUyBBnRVIW+0IHEFYbECq0REYQXCUMcpJqhRB469H+Aa+EUEUB9WEM3cOQUPZwAwUKApBWq+NTft0BkF7BXjC10cDQb0/2Pi2jQBBn2iJFals/EC1kRoscmf4QAdhe9nN8/RAnbJjHCX98UAk2Uj4lOPuQBBsmishhepAQxT5Jt/E5kAILJ7jbozjQCE5i4eayOBA8I80nMfR3EAQjYxLQ77YQG+OmnBXPtVA86KgyiM90kAdUmL7d1HPQOYaZlF548pAktY+Nd8Vx0BsmNgG+tHDQMGeFJVPBMFAlElj3U04vUDE4ylESRa5QKqO3DHqibVA4LmzIQZ+skAPNJXp4cCvQFHVlhchQ6tAzLBBa/9np0C03QiQfBikQFLgjgrZQKFA4HR32ECgnUDdI1NgiG+ZQL2XOsyJ1pVA+84pS8+/kkCXUHAY7BiQQBnBeCcdpItA0ghtyUO7h0DvfZzu+V+EQH+V69o5foFAvj+TnaUJfkCePqX4BMp5QMBqwfo5JHZAR5/IeoICc0ALQEMoMFJwQBcI3zhyBmxApYBaXLAPaEDqaSucdahkQEfriwF1vGFAMJa/SoF0XkAzBBdxwyVaQKogp4H+clZA+P+d6yJGU0AWRo/pP4xQQG4vLQwlakxAkfXSNkllSECO1HEW8/FEQCpQNHyN+0FAeurI/9jgPkCc2swsyII6QGkXFyLbwjZAsBzF07OKM0DeqgobHscwQPU6zFE6zyxAPys7WRK8KED5mbXGdTwlQLKdazKGOyJAOdDyrLFOH0B5JcNcF+EaQEZ01GjTExdARDICMzjQE0AnC2xDzQIRQDiETEe2NQ1ADBGtTg8UCUBAKNieAIgFQHawhpJhfAJARFQETQ++/z/aT9E5tED7PwwyNejpZfc/38O2DLIW9D8nrljrTj/xP8NXOCubne0/xsDBoEFt6T8gIBmNlNTlPx+9ZAUgvuI/mPZa5ngX4D+dnA7snaHbP8o1fx8dudc/rio4Tx9e1D8in7WEoHzRP70YrQvjBs4/QDwjpqLHyT8gJ6JKKiLGP+SO9ry5AMM/s6PcU6RQwD9AldUnwwO8P7OW6RdbDbg/2+HQcG2mtD9gKKsHr7qxP+8DawBnca4/BKCpvAojqj9MfbaVmXCmP7Lihe4FRKM/ILPBrmCKoD9v1Nai0GacP8UFB9pPYpg/C4aT5UfvlD8q1vNoJfmRP2vclxh83I4/1OR1OM5+ij/vRWEiNr+GP8ELPMVXh4M/tiZVrQDEgD+Aa4ACa8l8P5lJDmKftng/icZ+GlI3dT9CC/+cpjZyPxbXpnNnRW8/IOaQkjHYaj8WFWHiQwtnP2XFhcfyx2M/xkwphMf6YD9SJ147GCZdP/bPU1zPBFk/5itenhF5VT+BCwg/uG1SP4IvZLU0oU8/bCgKM0EkSz/+ymxd0ElHPwQAs/7l+kM/yi7hKMYjQT/3pTX9/mY9P2GxqEQMNzk/RO4t4rmeNT9ifHsGlogyP1OqI9tuxC8/LEr+Fp83Kz+++TqEmE8nP8OTp9oN9SM/WXXfxf0TIT+bRO8ddzYdPweoYFAR+Bg/iGjObnlTFT9vMOck9DISP0faHKuzBw8//k9xlU5sCj/uPYdKPHgGP6kU76PlEwM/aVuyMgwrAD8WfvWTolj7PoxsE0GcEfc+jQ2/NHVn8z6LASFGQUTwPhG1OjBpKus+aC3AiWSU5j7KvjWc1ariPtMw+xwgrN4+UkOFB5ID2T58C/oyMznUPnvkhVQkL9A+Vfsmkm+XyT4sOj+cmvHDPhLLOl4Ljb4+1JS2tt7ptj4vnDYSfr+wPvI4fIYdtac+V5sD24caoD4oFB8TDbmUPhOCfHONvYg+Z9+PF7Zlej6gbDyuBxdnPkBBf0AURkU+AKev/C8FzD1Ildte86Wivey89RJsarE9c2RHq7xuwL0ZqXGFQgPPPe7QV9DHQ9298YsawIWd6z3i1MNmEA/6vd3xMl0Qlwg+6DA7l0E0F75mKjfuqGYSPsDjv5VXBLc9", + "dtype": "f8" + } + }, + { + "mode": "lines", + "name": "piping", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEDAABwNX30QQYBAHA1ffRBMGYAcDV99EFx6idwNX30Qecxl381ffRBXTCQhTt99EFawf8ut370Qehmk9h7iPRB4wym80CR9EFDhhFOna30Qc8Rk86Q5PRBqluBQesT9UFZIb9ykzz1QVvx3IB7X/VBcdRvo3N99UG6rXOfLpf1QQ8Q4PBFrfVB39cyXj3A9UFTf30KhtD1QQV7RBiB3vVBbknq7IHq9UFZfSYi0PT1QQkwIDGp/fVB85Af4UEF9kF072GCxwv2QUk5ZvxgEfZBN0n+tS8W9kHgU41cUBr2Qcr0F4/bHfZBKUMiceYg9kF7AsgogyP2QcXs3aWdJPZBZ8kW/jn49UGsP4FcR8z1Qeg3VRCooPVB6lFk51t19UFD+bzZYkr1QZ1Ki8m8H/VBXa/ZhWn19EFSP3PNaMv0QR6KXlG6ofRBSz3+tl149EG3ZOSZUk/0QY44Y42YJvRBe9/kHS/+80FXNxLSFdbzQWCSzytMrvNBB1wVqdGG80EJwanEpV/zQYm9v/bHOPNB+FR/tTcS80F/L3p19OvyQVRjD6r9xfJB/cvAxVKg8kFG+3s683ryQX2F2HneVfJBCStN9RMx8kHIKl0ekwzyQcHZvmZb6PFBc3R8QGzE8UF4+w8exaDxQafPenJlffFBJahZsUxa8UFCZ/VOejfxQf8/UcDtFPFBG442e6by8EGdtD72o9DwQfBK26jlrvBBi9ZcC2uN8EEtRviWM2zwQVlcy8U+S/BBRDDgEowq8EHB5i/6GgrwQR2ASfHV0+9BpCY3GPeT70HqmclmmFTvQW8Eity4Fe9BIFPqelfX7kGwCkVFc5nuQYiY3EALXO5BPjPadB4f7kFrW0zqq+LtQY8LJayypu1BfaM3xzFr7UELmjZKKDDtQT7+sEWV9exB3c8PzHe77EEkNpPxzoHsQWOaT8yZSOxBpqsqdNcP7EFxT9gCh9frQXSD15Onn+tBITNvRDho60H5A6szODHrQc8aWIKm+upB7dsBU4LE6kHTp+7Jyo7qQQOWHA1/WepBAC8+RJ4k6kGcJreYJ/DpQXYXmTUavOlBVkCgR3WI6UEpRDD9N1XpQRbtUIZhIulBG/OqFPHv6EGoxoTb5b3oQXxfvw8/jOhBERDT5/ta6EHKXcybGyroQSLeSGWd+edB7Rh0f4DJ50HkbwQnxJnnQX0LOJpnaudBQM3RGGo750GKRxbkygznQeO6yD6J3uZB4BgobaSw5kGlDOy0G4PmQfsIQl3uVeZBBlzKrhsp5kGaSJXzovzlQR4lIHeD0OVBGYBShryk5UFBSntvTXnlQR4GToI1TuVBN/3fD3Qj5UG8eqVqCPnkQaILb+bxzuRBQMRm2C+l5EFHiw2XwXvkQR5qOHqmUuRBm+IN290p5EEFSgMUZwHkQV0p2oBB2eNB4aKdfmyx40HJ3J9r54njQSRxd6exYuNB3OL8kso740HJF0iQMRXjQdXXrQLm7uJBJVG9TufI4kE0oT3aNKPiQeNiKwzOfeJBb0G2TLJY4kE7kD4F4TPiQXnnUqBZD+JBkcWtiRvr4UFTNTMuJsfhQcx47vt4o+FB37gPYhOA4UFvuenQ9FzhQSmS77kcOuFB42uyj4oX4UF8Qt/FPfXgQTKrPNE10+BBg5+oJ3Kx4EFiTBZA8o/gQdvli5K1buBBCX8gmLtN4EFe5vnKAy3gQTSGSqaNDOBBJJOeTLHY30FZCpuQyJjfQRTJIRVgWd9Bn7bS2HYa30EqNVLcC9zeQawURSIent5BuI1Mr6xg3kFBRAKKtiPeQUxS9Lo6591BTlqhTDir3UFLoXRLrm/dQYUwwsWbNN1Btv7Cy//53EGjIJFv2b/cQQUBJMUnhtxBlZ9M4ulM3EE317HeHhTcQQmrzNPF29tBSprk3N2j20HT+gsXZmzbQSFaHKFdNdtBmuSym8P+2kH90iwpl8jaQaXdo23XktpBb7XqjoNd2kEHgom0mijaQTJlugcc9NlB1QJmswbA2UFJDSDkWYzZQYDVI8gUWdlBjt5QjzYm2UHucydrvvPYQexBxY6rwdhBXO/hLv2P2EHdt8uBsl7YQY8FZL/KLdhBMwkcIUX910FPT/HhIM3XQetQaj5dnddBIv6SdPlt10GUQPnD9D7XQWJzqW1OENdBJs4qtAXi1kHSv3vbGbTWQfc1DimKhtZBdszD41VZ1kHm4OlTfCzWQWKENcP8/9VBeUW/fNbT1UEqy/7MCKjVQak4xgGTfNVBUFA9anRR1UHdStxWrCbVQSBWZhk6/NRBhKzjBB3S1EFwNJttVKjUQeWTC6nfftRBxKHjDb5V1EF8Gvrz7izUQf55RLRxBNRBA9jMqEXc00FKnqYsarTTQTzq4ZvejNNBHGR9U6Jl00F4TVaxtD7TQRh/FhQVGNNBOgMg28Lx0kE37HVmvcvSQZj4oRYEptJBjYSWTJaA0kGsNIxpc1vSQequ2s6aNtJBLZzL3QsS0kEsDmf3xe3RQblBOHzIydFBtYwJzBKm0UGrF5ZFpILRQTvML0Z8X9FBlaJXKZo80UFXL0ZI/RnRQQX/Yfmk99BBUuugj5DV0EFYJtFZv7PQQcE2xqEwktBBs4plq+Nw0EGNm42z10/QQVTW0e4LL9BBX5gEiH8O0EHkEQ89Y9zPQU/TsohCnM9Bvjau95pcz0FtCUhqah3PQfp1WJau3s5BCXxhAGWgzkFHcZfzimLOQTABsHgdJc5BFedHSxnozUH6fajNeqvNQcYZrvo9b81BjMiGVV4zzUHWbPbW1vfMQQ3mwNehvMxBpPDM97iBzEFxOoIBFUfMQQ9q0sitDMxBsPVGBXrSy0EzzVMmb5jLQVffDyGBXstB86tONqIky0GN3OavwurKQZWAtZPQsMpBX3eeSrd2ykGFg2A5XzzKQYM+gkitAcpBqJrSV4LGyUElMtmYuorJQZDp+MksTslBxMe3SakQyUHx+0D2+NHIQbNPQsjbkchBb6fnEAdQyEGP0JQ3IwzIQUcw1sLIxcdBqjnoXHx8x0FX7RJYqS/HQcGheOyZ3sZBwOIN5WuIxkE7RmN7/yvGQUKp8zfcx8VBy6TjiQNaxUGtserZnt/EQa90HVJbVMRB0qitqeeww0ECTZzISObCQf45JivzwsFB64mKPAVdwEGwBTZZ2su9QVT57ye1l7lBvsoaYA35tUGEVL80cd2yQc8Org1dMrBB8hLRrczPq0Ets1iCxeCnQdKeOZUtgKRBUthShd+ZoUFIr79mZiKfQT9MOLIrt5xB", + "dtype": "f8" + } + }, + { + "mode": "lines", + "name": "piping_wall", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BAAAAAGXN3UEAAAAAZc3dQQAAAABlzd1BZ/bbNGKd3UFtA+LyzG3dQdsbpW+EPt1BtfZU5YcP3UGddjTZ1uDcQYzI9tFwstxBTooTV1WE3EEJBsTwg1bcQX73ASj8KNxB1FaGhr3720FPJciWx87bQQI8++MZottBYRwP+rN120HBwq1llUnbQbx6OrS9HdtBcbXQcyzy2kGm4UIz4cbaQc1FGYLbm9pBzduQ8Bpx2kG1LpoPn0baQTo62HBnHNpBAkyfpnPy2UG/5vNDw8jZQRSnidxVn9lBOyrCBCt22UF39qtRQk3ZQUFlAVmbJNlBO48nsTX82EHROS3xENTYQafGybAsrNhBqCRciIiE2EHUwukQJF3YQb2EHeT+NdhBsrhGnBgP2EGTD1jUcOjXQUuW5icHwtdB+bAoM9ub10G3F/WS7HXXQfnUweQ6UNdBmkWjxsUq10F4GkvXjAXXQaxbB7aP4NZBXW3BAs671kEeFv1dR5fWQeaG12j7ctZBjWQGxelO1kHg0tYUEivWQTKBLPtzB9ZBf7iAGw/k1UERa+EZ48DVQaJF8JrvndVBA8LhQzR71UE9O3y6sFjVQTEDF6VkNtVBrXmZqk8U1UH7JHpycfLUQdrLvaTJ0NRB+pD26Vev1EHTD0PrG47UQf16TVIVbdRB4btKyUNM1EHdk/n6pivUQcy+oZI+C9RB7RYTPArr00EzuqSjCcvTQfEwNHY8q9NB4ZUkYaKL00GKv10SO2zTQf1qSzgGTdNB5WfcgQMu00HyxYGeMg/TQZQDLj6T8NJBAz5UESXS0kGbYufI57PSQX5hWRbbldJBg2Gaq/530kFv9Rc7UlrSQXBSvHfVPNJB2YftFIgf0kEsuIzGaQLSQVRT9UB65dFBKlL8OLnI0UEsc+9jJqzRQXV4lHfBj9FB5mYoKopz0UGMxl4ygFfRQTPkYEejO9FBMRTNIPMf0UFd9rV2bwTRQTO7oQEY6dBBK2qJeuzN0EExKdia7LLQQUyFahwYmNBBa7yNuW590EFVCP8s8GLQQbzq6jGcSNBBcHrsg3Iu0EG1sQzfchTQQWR7g/859c9B/53cReHBz0GZ1sEL247PQTyApswmXM9BRe3SBMQpz0E+EWMxsvfOQcssRdDwxc5B0ns4YH+UzkGm5ctgXWPOQWWvXFKKMs5BYDAVtgUCzkGbiOsNz9HNQWFZoNzloc1B4n+9pUlyzUHW0ZTt+ULNQTTcPjn2E81B36OZDj7lzEFhaEf00LbMQY1orXGuiMxBNqnyDtZazEHBvf5URy3MQbySeM0BAMxBXjrFAgXTy0HzugaAUKbLQSXfGtHjectBRAiagr5Ny0FRAtYh4CHLQfrZ2DxI9spBaLRjYvbKykHkqO0h6p/KQVCcogsjdcpBZh5isKBKykHNSL6hYiDKQeaf+nFo9slBafUKtLHMyUGpTJL7PaPJQaTA4dwMeslBuGv37B1RyUECUX3BcCjJQXBHyPAEAMlBX+bWEdrXyEHbc1C876/IQWvUg4hFiMhBZHxmD9tgyEG2YpPqrznIQSr1SbTDEshBBA5tBxbsx0H+6oF/psXHQYYlr7h0n8dBOKy7T4B5x0FuvQ3iyFPHQebiqQ1OLsdBUe4xcQ8Jx0G+9uOrDOTGQatWmV1Fv8ZBqarFJrmaxkFf0HWoZ3bGQbTlToRQUsZB8UeNXHMuxkGPkgPUzwrGQWOdGY5l58VB2HnLLjTExUHIbqhaO6HFQYHy0bZ6fsVBbKL66PFbxUG5N2WXoDnFQVx442iGF8VBgiPVBKP1xEGU2CYT9tPEQab3UDx/ssRBE3pWKT6RxEHMwcODMnDEQaRdrfVbT8RBrsCuKbouxEFV6ujKTA7EQZ78AIUT7sNBhL0eBA7Ow0Hm/+r0O67DQf3vjQSdjsNBnD6t4DBvw0HQJWo390/DQY9AX7fvMMNBIi6eDxoSw0Hw+KzvdfPCQcY2gwcD1cJBZteGB8G2wkFClIigr5jCQTHyv4POesJBsMPGYh1dwkFQF5Tvmz/CQeR6dtxJIsJBQ3gN3CYFwkEpLEKhMujBQdXSPt9sy8FBUiBlSdWuwUGSM0OTa5LBQRHshnAvdsFBomDvlCBawUHfKzy0Pj7BQcY2GoKJIsFBVZwNsgAHwUEAMVj3o+vAQXom3ARz0MBB7y76jG21wEE9ampBk5rAQQNLD9Pjf8BBhYDB8V5lwEErzBNMBEvAQQ99DY/TMMBBLhnaZcwWwEENAdny3Pm/QWckJ+Byxr9B/xj82VmTv0EvC1YfkWC/QfDhh+YXLr9BPFe5W+37vkFA+CqfEMq+QQGiNMOAmL5BKpfzyTxnvkH6mZuiQza+QYGAXCaUBb5Be3jLFC3VvUFilrwPDaW9QS9NdpYydb1Bz/kkAJxFvUG8tXB1Rxa9QeYUE+ky57xBpDRFD1y4vEHZe9hTwIm8QVyBxc5cW7xBeJv0Ni4tvEFxefrSMP+7QbOleGdg0btBNrDFIriju0HhuHKFMna7QWe8M0fJSLtBCtycN3Ubu0HlHRAaLu66QTjLG3zqwLpBThpohZ+TukGpASnAQGa6QVKn09i/OLpBoAWRUwwLukFD9Yk1E925QVJJqZ++rrlBRYK/V/V/uUFj2vk7mlC5QQ0QOpuLILlBBpnaaqLvuEFcLH9Psb24QUPiRWqDirhBO8ZH1dpVuEHmY6uwbh+4Qf2Wn5Lo5rdBc3z6FOGrt0ESbrkU2223QcXudvY8LLdBxiuiz0bmtkE8xxp9A5u2Qd2FSAMxSbZB9vLi9RjvtUFSnwP3SIq1QcMFIFMDF7VBrSnZ6e2OtEFsrsiPAOWzQYtgcdIm7rJBPxlu8Xy6sUHTP4E7XXKwQRUaGXsrIK1B3jKfWujHqUF5iBbg9dGmQcc3ba8OM6RBsam7wzPhoUGJHKrfF6efQZ/uaFOEBJxB2L5IVc/MmEGVSGmaz3qWQaxmY4tEEpVB", + "dtype": "f8" + } + }, + { + "mode": "lines", + "name": "neutron_rate", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAAABAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAQAAQAjWmEEBAABACNaYQQEAAEAI1phBAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAEAAEAI1phB", + "dtype": "f8" + } + }, + { + "mode": "lines", + "name": "h2", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/AAAAAAAA8D8AAAAAAADwPwAAAAAAAPA/", + "dtype": "f8" + } + } + ], + "layout": { + "font": { + "size": 10 + }, + "legend": { + "title": { + "text": "Components" + } + }, + "template": { + "data": { + "bar": [ + { + "error_x": { + "color": "#2a3f5f" + }, + "error_y": { + "color": "#2a3f5f" + }, + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "bar" + } + ], + "barpolar": [ + { + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "barpolar" + } + ], + "carpet": [ + { + "aaxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "baxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "type": "carpet" + } + ], + "choropleth": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "choropleth" + } + ], + "contour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "contour" + } + ], + "contourcarpet": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "contourcarpet" + } + ], + "heatmap": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "heatmap" + } + ], + "histogram": [ + { + "marker": { + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "histogram" + } + ], + "histogram2d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2d" + } + ], + "histogram2dcontour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2dcontour" + } + ], + "mesh3d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "mesh3d" + } + ], + "parcoords": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "parcoords" + } + ], + "pie": [ + { + "automargin": true, + "type": "pie" + } + ], + "scatter": [ + { + "fillpattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + }, + "type": "scatter" + } + ], + "scatter3d": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatter3d" + } + ], + "scattercarpet": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattercarpet" + } + ], + "scattergeo": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergeo" + } + ], + "scattergl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergl" + } + ], + "scattermap": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermap" + } + ], + "scattermapbox": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermapbox" + } + ], + "scatterpolar": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolar" + } + ], + "scatterpolargl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolargl" + } + ], + "scatterternary": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterternary" + } + ], + "surface": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "surface" + } + ], + "table": [ + { + "cells": { + "fill": { + "color": "#EBF0F8" + }, + "line": { + "color": "white" + } + }, + "header": { + "fill": { + "color": "#C8D4E3" + }, + "line": { + "color": "white" + } + }, + "type": "table" + } + ] + }, + "layout": { + "annotationdefaults": { + "arrowcolor": "#2a3f5f", + "arrowhead": 0, + "arrowwidth": 1 + }, + "autotypenumbers": "strict", + "coloraxis": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "colorscale": { + "diverging": [ + [ + 0, + "#8e0152" + ], + [ + 0.1, + "#c51b7d" + ], + [ + 0.2, + "#de77ae" + ], + [ + 0.3, + "#f1b6da" + ], + [ + 0.4, + "#fde0ef" + ], + [ + 0.5, + "#f7f7f7" + ], + [ + 0.6, + "#e6f5d0" + ], + [ + 0.7, + "#b8e186" + ], + [ + 0.8, + "#7fbc41" + ], + [ + 0.9, + "#4d9221" + ], + [ + 1, + "#276419" + ] + ], + "sequential": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "sequentialminus": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ] + }, + "colorway": [ + "#636efa", + "#EF553B", + "#00cc96", + "#ab63fa", + "#FFA15A", + "#19d3f3", + "#FF6692", + "#B6E880", + "#FF97FF", + "#FECB52" + ], + "font": { + "color": "#2a3f5f" + }, + "geo": { + "bgcolor": "white", + "lakecolor": "white", + "landcolor": "#E5ECF6", + "showlakes": true, + "showland": true, + "subunitcolor": "white" + }, + "hoverlabel": { + "align": "left" + }, + "hovermode": "closest", + "mapbox": { + "style": "light" + }, + "paper_bgcolor": "white", + "plot_bgcolor": "#E5ECF6", + "polar": { + "angularaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "radialaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "scene": { + "xaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "yaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "zaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + } + }, + "shapedefaults": { + "line": { + "color": "#2a3f5f" + } + }, + "ternary": { + "aaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "baxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "caxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "title": { + "x": 0.05 + }, + "xaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + }, + "yaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + } + } + }, + "xaxis": { + "exponentformat": "power", + "title": { + "text": "Time [day]" + } + }, + "yaxis": { + "exponentformat": "power", + "title": { + "text": "Block value [#]" + } + } + } + }, + "text/html": [ + "
\n", + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "import plotly.graph_objects as go\n", + "\n", + "fig = go.Figure()\n", + "for p, d in enumerate(data):\n", + " lb = Sco.labels[p] if p < len(Sco.labels) else f\"port {p}\"\n", + " fig.add_trace(go.Scatter(x=time * s_to_days, y=d, mode=\"lines\", name=lb))\n", + "\n", + "fig.update_layout(\n", + " legend_title=\"Components\",\n", + " font=dict(size=10),\n", + " # yaxis_type=\"log\",\n", + " # xaxis_type=\"log\",\n", + " xaxis_title=\"Time [day]\",\n", + " yaxis_title=\"Block value [#]\",\n", + ")\n", + "fig.update_xaxes(exponentformat=\"power\")\n", + "fig.update_yaxes(exponentformat=\"power\")\n", + "fig.write_html(\"baby.html\")\n", + "fig.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 3, + "metadata": {}, + "outputs": [], + "source": [ + "from libra_toolbox.tritium.plotting import quantity_to_activity\n", + "from libra_toolbox.tritium.model import ureg\n", + "\n", + "from scipy.integrate import cumulative_trapezoid\n", + "\n", + "# compute cumulative release\n", + "c_salt = data[0] # salt concentration\n", + "c_piping = data[1] # piping concentration\n", + "c_piping_OV = data[2] # piping OV concentration\n", + "h2_in = data[4] # H2 concentration\n", + "\n", + "direct_release_from_salt = (1 - f) * c_salt * k_top * A_top\n", + "release_from_isotopic_exchange = c_piping * h2_in * k_exchange\n", + "top_release = direct_release_from_salt + release_from_isotopic_exchange\n", + "cumulative_top = cumulative_trapezoid(top_release, time, initial=0)\n", + "\n", + "wall_release = c_salt * k_wall * A_wall * h2_in + c_piping_OV * h2_in * k_exchange_OV\n", + "cumulative_wall = cumulative_trapezoid(wall_release, time, initial=0)\n", + "\n", + "cumulative_top = quantity_to_activity(cumulative_top * ureg.particle).to(\"Bq\")\n", + "cumulative_wall = quantity_to_activity(cumulative_wall * ureg.particle).to(\"Bq\")\n", + "\n", + "cumulative_direct_top = cumulative_trapezoid(direct_release_from_salt, time, initial=0)\n", + "cumulative_direct_top = quantity_to_activity(cumulative_direct_top * ureg.particle).to(\"Bq\")\n", + "\n", + "cumulative_release_from_isotopic_exchange = cumulative_trapezoid(\n", + " release_from_isotopic_exchange, time, initial=0\n", + ")\n", + "cumulative_release_from_isotopic_exchange = quantity_to_activity(\n", + " cumulative_release_from_isotopic_exchange * ureg.particle\n", + ").to(\"Bq\")" + ] + }, + { + "cell_type": "code", + "execution_count": 4, + "metadata": {}, + "outputs": [ + { + "data": { + "application/vnd.plotly.v1+json": { + "config": { + "plotlyServerURL": "https://plot.ly" + }, + "data": [ + { + "line": { + "color": "#023047" + }, + "mode": "lines", + "name": "Top release", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAABuPaKwlj5RNxSwOxz+qnc3GBJga9DW0zeH2dvPDAw9OFPbxzzWi6Y4InpivFKaETl4Nn9YrIB7ORiQpYJ9fOU5obJfV0HJUDqKxfn7dTq6Oks0oCqsfSQ7Z3cggS4CkDviWbqpaAP5O+gJjsSpimM8NOtxQqmIzjzZxqk2xNo3PTXEtiHpoqI94AEvxIkeDT5TVnUVyL92PjF+lqA7xeE+XcUL6d26Sz/RrD4LXWerPzO2WoGsntk/YNJxlwzN5j/eT9LDD5L7P2mP3e7UeA1AtyhTW7Z+FUBNM4ofU0wbQCeW3+XnIyBAG1oBVHJHIkAXRIr7ih0kQNLOZAcmsSVA2fwGA6sLJ0AQ5KvzLDUoQKUoS4KaNClAwJCbVecPKkArNsuSL8wqQKBhtFjWbStA6edh65/4K0CnRP4ryG8sQGpDNOMV1ixANRcTUestLUBaN0FlVHktQC4gUPQSui1Ax9QDM6nxLUCFIxS2YiEuQMaE983/tC5ASgFY2c3iLkBcHtOpuQ8vQK9cpmXyOy9A/e3nRIZnL0DWSB9XgZIvQN3+LQ3uvC9A7MhRdNXmL0CSByO0HwgwQAOyJl+ZHDBAb4IStdowMEAMX99O5kQwQA8yP3K+WDBAbKdOHWVsMEBOpp8Q3H8wQONj2NckkzBAgHUa0UCmMEDyAl4zMbkwQBj85hP3yzBAAdn0apPeMEDs0cgXB/EwQGeJGuRSAzFAlr4Ph3cVMUA+tMindScxQDd5kN9NOTFAhhm+uwBLMUCa51G/jlwxQB54WGT4bTFAJI8bHT5/MUDnEihVYJAxQPIXL3JfoTFAlz7H1DuyMUCK3BLZ9cIxQEzLTteN0zFA4ChNJATkMUCB4N4RWfQxQKFqLu+MBDJAKdsNCaAUMkBaCTqqkiQyQLJblBtlNDJAe4pUpBdEMkCiezOKqlMyQJUxkBEeYzJAU6OPfXJyMkC4NTgQqIEyQAp0iQq/kDJAU4+QrLefMkAFGXo1kq4yQOhcoeNOvTJAFLCd9O3LMkDL/k2lb9oyQFfY4jHU6DJAXC/n1Rv3MkA1/UfMRgUzQHDwWk9VEzNAx0jkmEchM0Af/hviHS8zQPBLsmPYPDNAx7bTVXdKM0Cunizw+lczQGNu7GljZTNAT3XI+bByM0D0d/7V438zQAYBVzT8jDNA6XsnSvqZM0AWIFRM3qYzQM2zUW+oszNAoisn51jAM0CsK2/n78wzQE1uWaNt2TNAQBSsTdLlM0DE4MQYHvIzQJhkmjZR/jNA6hm92GsKNEA+c1gwbhY0QMreM25YIjRA3L+zwiouNEBbX9pd5Tk0QIPUSG+IRTRAt+Y/JhRRNEA96aCxiFw0QHiR7j/mZzRASshN/yxzNEAAd4YdXX40QF1QBMh2iTRA8JXXK3qUNEAy2rV1Z580QJa/+tE+qjRA1bSobAC1NEClrmlxrL80QP/ejwtDyjRAOGoWZsTUNEDsGaKrMN80QO0NggaI6TRAS2uwoMrzNECRCdOj+P00QDoePDkSCDVAh+bqiRcSNUCqT4y+CBw1QHide//lJTVAfw/DdK8vNUCvhBxGZTk1QJod8poHQzVARN1empZMNUCQSC9rElY1QGIE4jN7XzVAY3KoGtFoNUB/TGdFFHI1QBY/t9lEezVA8YHl/GKENUDrb/TTbo01QHodnINoljVA3u1KMFCfNUBBJyb+Jag1QI+FChHqsDVAJ8yMjJy5NUBwVvqTPcI1QDinWUrNyjVA+PZq0kvTNUD3wKhOuds1QFJPSOEV5DVA30U6rGHsNUAILCvRnPQ1QIP1g3HH/DVA/IlqruEENkCvS8Ko6ww2QPCcLIHlFDZArWQJWM8cNkDjkXdNqSQ2QAWeVYFzLDZAbQ5CEy40NkDC9Jsi2Ts2QFhug850QzZAniLaNQFLNkCIwEN3flI2QAV7JrHsWTZAfoSrAUxhNkBXib+GnGg2QJMpE17ebzZAcHEbpRF3NkAyURJ5Nn42QOgT9/ZMhTZAYNWOO1WMNkA192RjT5M2QOuUy4o7mjZASvfbzRmhNkDHBndI6qc2QCe9RRatrjZATpa5UmK1NkA6AA0ZCrw2QEfKQ4SkwjZAnZMrrzHJNkD2OFy0sc82QKNBOK4k1jZA50vttorcNkCueHTo4+I2QK7Wklww6TZA7szZLHDvNkDLhKdyo/U2QJFTJ0fK+zZAoyNSw+QBN0BI3e7/8gc3QFDPkhX1DTdAbxeiHOsTN0CjClAt1Rk3QJ6dn1+zHzdAXM1jy4UlN0ArCECITCs3QBmXqK0HMTdAPAjjUrc2N0DpmQaPWzw3QCen/Hj0QTdAvBWBJ4JHN0AdxiKxBE03QMkFRCx8UjdAYQQbr+hXN0A7TLJPSl03QO8+6SOhYjdAuZZ0Qe1nN0Bm7d69Lm03QARKia5lcjdAN7arKJJ3N0C43FVBtHw3QHGybw3MgTdA7yu6odmGN0BPAdAS3Ys3QOOCJnXWkDdAWYEO3cWVN0B7TLVeq5o3QDvMJQ6HnzdAHrhJ/1ikN0AC8upFIak3QKsJtfXfrTdAwO82IpWyN0B73+TeQLc3QNaHGj/juzdAE34dVnzAN0AzBSA3DMU3QJg3RPWSyTdAS6KfoxDON0CTYz9VhdI3QG7hLB3x1jdApC9zDlTbN0DQQCU8rt83QBMCZbn/4zdA+oVrmUjoN0DeaJLviOw3QHifXs/A8DdAKuiMTPD0N0D5HyB7F/k3QHjGcW82/TdAjgdFPk0BOEA6r9z8WwU4QFV7FMFiCThAyVB+oWENOEC87oO1WBE4QJbRjRVIFThAfRMw2y8ZOEA9N14hEB04QPXupgTpIDhAnRl5o7okOEC4YnMehSg4QG0mwJhILDhADXx/OAUwOEB0kUEnuzM4QGDVk5JqNzhAXc6jrBM7OEBY6fqstj44QCgIV9FTQjhAvyekXutFOECkGRyifUk4QEoIkvIKTThAiVXwsZNQOED8WvBOGFQ4QMyzFEeZVzhAnwfwKBdbOEBv28OWkl44QLvIhEkMYjhA9cBTFIVlOEAt2H7o/Wg4QB+xH9p3bDhAjmNyJfRvOECYEgY1dHM4QJQg8aj5djhAVxQ/X4Z6OEBkdN99HH44QHdxdX6+gThAS7qLPG+FOEA9iuMGMok4QDgW4rQKjThATrCQwf2QOEA4BU9uEJU4QHJgXvBImThADigNrq6dOEAFz++SSqI4QJXgDocnpzhA9ILcHFOsOEDzZoiY3rE4QH/isZDgtzhAhkMGqHe+OEDImwZ4z8U4QMidd2wqzjhAr5yx8vjXOEDkJRJAIeQ4QEHLwv6c9ThASWlwERkLOUAZ/KLSnSE5QMy6vCr0QTlAjWWJnMxdOUD2fhz4tHU5QOSyZqk7ijlAFgvAINubOUC2M3d9/Ko5QIw1ggD6tzlAaO3cJiHDOUCs4Y6i8so5QG1V6geYzzlA", + "dtype": "f8" + } + }, + { + "line": { + "color": "green" + }, + "mode": "lines", + "name": "Wall release", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAh0c5QJFdM/4BR/Y8Qm2T+6i0qzroreP6xtVKVjruE/1HGBPLPY4z+N1YFc/szlP+uCZb/Bkuc/7mjNQ2sw6T+9VP47gKvqP2z9xVG+COw/WIp2wzdM7T/oLGqia3nuPyofraNak+8/cfJ7f0xO8D9A7knjrsvwP+kh2AnIQvE/Ju3dKmm08T8dBx0ERiHyPz3pfwT5ifI/13dJ4Abv8j9g/a6j4VDzPxqSMFbrr/M/lWJsPngM9D/4fe/T0Gb0P8SjnWozv/Q/eYelodUV9T8UaI6e5Wr1P9lAthuLvvU/+QqOUOgQ9j/oUf24GmL2P65tksA7svY/vX58VWEB9z+L2btmnk/3P7WQeVADnfc/6lYMOZ7p9z+UitVgezX4P0TE0malgPg/buB8gyXL+D8+pVO7AxX5P+96QwpHXvk/rQnoiPWm+T9o9ImMFO/5P4V5lsKoNvo/q8Y0SLZ9+j+LoYW+QMT6P2wjBlxLCvs/sjF9+9hP+z9ts8wo7JT7Pxnu8iuH2fs/J6R9Eqwd/D8MSae3XGH8P7mmTsuapPw/U3ny12fn/D/7ptRHxSn9P0auY2m0a/0/PJIEczat/T8c01OGTO79P03S8LL3Lv4/lELk+Dhv/j+L669KEa/+P1kCFI+B7v4/4J2Uooot/z+BT8hYLWz/P0qid31qqv8/vimU1ULo/z8VbAaQ2xIAQIdCQgtkMQBAiG33NrtPAEDfpYtr4W0AQDmSiv/WiwBAw7fXR5ypAEDvddmXMccAQDkJnkGX5ABA42z7lc0BAUBd1ark1B4BQHViYHytOwFAOZLfqldYAUDa6Qy903QBQH05/f4hkQFAadECvEKtAUD88rg+NskBQC+9DdH85AFAYMtKvJYAAkCItRxJBBwCQCKamb9FNwJAjNRGZ1tSAkCB/R2HRW0CQGhPkWUEiAJAQYSPSJiiAkAlQId1Ab0CQHcYajFA1wJArUWvwFTxAkCjC1ZnPwsDQLLi52gAJQNAU2p6CJg+A0DcLLGIBlgDQOc6vytMcQNAyaNoM2mKA0AG0APhXaMDQODBenUqvANAXT9MMc/UA0AB6YxUTO0DQKk/6B6iBQRA4Juhz9AdBECPGJWl2DUEQLByON+5TQRAWt+bunRlBECU2Wp1CX0EQLTo7Ex4lARAZGAGfsGrBED+GjlF5cIEQPAupd7j2QRArJ8Jhr3wBEC2CsV2cgcFQDdR1usCHgVAWz7dH280BUDpKhtNt0oFQFqec63bYAVAfO1setx2BUAe1zDtuYwFQK0ejT50ogVAESX0pgu4BUDXf31egM0FQOSO5pzS4gVAjhCTmQL4BUB6tI2LEA0GQA6tiKn8IQZAwj/eKcc2BkArVJFCcEsGQPkBTin4XwZA0h1qE190BkAjxeU1pYgGQPboa8XKnAZAwddS9s+wBkBNxpz8tMQGQKFX+At62AZAGyTBVx/sBkCZPwATpf8GQMu+bHALEwdArjtsolImB0BAWRPbejkHQGRGJkyETAdA+j8ZJ29fB0BMEhGdO3IHQK+Z497phAdAgkIYHXqXB0CAiOiH7KkHQG11QE9BvAdAJx+/onjOB0A6JbexkuAHQOMtL6uP8gdAn2LivW8ECEBU7EAYMxYIQCNvcOjZJwhA/IVMXGQ5CED1PWeh0koIQJuRCeUkXAhAPeQzVFttCEBkfZ4bdn4IQIsEumd1jwhAZv2vZFmgCECvRGM+IrEIQPWNcCDQwQhAcOIuNmPSCEBaIbCq2+IIQP2Bwag58whA8xfsWn0DCUALWnXrphMJQEGsX4S2IwlAku1qT6wzCUAuChV2iEMJQAKTmiFLUwlAdVv3evRiCUByHueqhHIJQAgs5tn7gQlAEyEyMFqRCUCVqsrVn6AJQMVWcvLMrwlABHavreG+CUBnDs0u3s0JQMLk25zC3AlA3Z6zHo/rCUC0AvTaQ/oJQIpXBvjgCApAUu4enGYXCkBq1z7t1CUKQFbNNREsNApAi1ykLWxCCkAtUv5nlVAKQAt9jeWnXgpAzM10y6NsCkBu5bM+iXoKQJ4jK2RYiApAAUmgYBGWCkDdxMNYtKMKQFzJNnFBsQpAbkWSzri+CkDh6G6VGswKQHFcbupm2QpAfN5F8p3mCkB+fMrRv/MKQHcq/63MAAtAdwMlrMQNC0BqCs7xpxoLQAHQ8qR2JwtAX3IL7DA0C0ByfSzu1kALQFZIKNNoTQtARYS1w+ZZC0BTzpvpUGYLQLk152+ncgtAoM4jg+p+C0DtlKJRGosLQHcVyQs3lwtAg45t5ECjC0AVfEEROK8LQKTSTMscuwtANIV8T+/GC0C7W0ffr9ILQOOWbMFe3gtA6mjQQvzpC0DK6nq3iPULQNvsvnsEAQxAANaO9W8MDEDbtwaWyxcMQPjWM9sXIwxAticjUlUuDEDiskGZhDkMQPV9G2OmRAxAy3qHebtPDEAxNFHBxFoMQMtxcz7DZQxAzvb5GLhwDEDh3KOipHsMQNb3Y12KhgxAgEjhAmuRDEAQ+R6NSJwMQInteUAlpwxA5voxtwOyDEB14L7u5rwMQK7wPVfSxwxAEhpW5cnSDEAObAYn0t0MQGbN8Vvw6AxAWB7lkSr0DEAe8ZDGh/8MQImKxQ8QCw1A/yb7y8wWDUD39KPdyCINQMQizfQQLw1AAlcd7LM7DUCQQJFAw0gNQJnP665TVg1AZ+t5B35kDUDH1gZSYHMNQPANgmsfgw1AAX2DYumTDUBddoAL+aUNQNW6YqybuQ1A6WF3fTrPDUBD10DTa+cNQLJON4oVAw5AcYu4870jDkA4T0PuhUwOQOnDTqDFhw5A+ft5uaDRDkDOPmdKZyAPQCRj8wJFlA9A3UEOmyD7D0BunEZrFisQQJS+hSdiUxBA2895NQ13EEAAUECin5YQQOOhHtuRshBAEYPaeE7LEEB0NXdNI90QQGfmAt/05xBA", + "dtype": "f8" + } + }, + { + "marker": { + "color": "#023047" + }, + "mode": "markers", + "name": "Samples (Top)", + "type": "scatter", + "x": [ + 0.5166666666666666, + 1.2180555555555552, + 2.161805555555555, + 4.0361111111111105, + 6.997916666666665, + 10.122916666666665, + 14.02083333333333, + 18.071527777777774, + 19.26180555555555, + 21.059027777777775, + 24.13333333333333, + 28.05694444444444, + 32.130555555555546, + 35.103472222222216, + 39.0611111111111, + 42.03263888888888 + ], + "y": { + "bdata": "rUfhehSu1z8730+Nl24JQFTjpZvEoBVAbOf7qfESIUBmZmZmZqYmQGQ730+NFypA+FPjpZuELUCKQWDl0OIvQEa28/3U2DJAGQRWDi0SNED2KFyPwnU1QKrx0k1iEDZA8KfGSzeJOEDo+6nx0u04QB1aZDvfbzlAAAAAAACAOUA=", + "dtype": "f8" + } + }, + { + "marker": { + "color": "green" + }, + "mode": "markers", + "name": "Samples (Walls)", + "type": "scatter", + "x": [ + 0.5166666666666666, + 6.990972222222221, + 14.02083333333333, + 18.071527777777774, + 19.26180555555555, + 21.059027777777775, + 24.13333333333333, + 28.05694444444444, + 32.130555555555546, + 35.103472222222216, + 39.0611111111111, + 42.03263888888888 + ], + "y": { + "bdata": "cBSuR+F6lD/81HjpJjHIP4PrUbgehdM/////////1z8oXI/C9SjgPw4tsp3vp+Y/MAisHFpk9T+CwMqhRbYBQEbhehSuRwdA1aNwPQrXCkC+yqFFtvMNQNJNYhBYORBA", + "dtype": "f8" + } + } + ], + "layout": { + "font": { + "size": 10 + }, + "height": 600, + "legend": { + "title": { + "text": "Components" + } + }, + "template": { + "data": { + "bar": [ + { + "error_x": { + "color": "#2a3f5f" + }, + "error_y": { + "color": "#2a3f5f" + }, + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "bar" + } + ], + "barpolar": [ + { + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "barpolar" + } + ], + "carpet": [ + { + "aaxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "baxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "type": "carpet" + } + ], + "choropleth": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "choropleth" + } + ], + "contour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "contour" + } + ], + "contourcarpet": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "contourcarpet" + } + ], + "heatmap": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "heatmap" + } + ], + "histogram": [ + { + "marker": { + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "histogram" + } + ], + "histogram2d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2d" + } + ], + "histogram2dcontour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2dcontour" + } + ], + "mesh3d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "mesh3d" + } + ], + "parcoords": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "parcoords" + } + ], + "pie": [ + { + "automargin": true, + "type": "pie" + } + ], + "scatter": [ + { + "fillpattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + }, + "type": "scatter" + } + ], + "scatter3d": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatter3d" + } + ], + "scattercarpet": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattercarpet" + } + ], + "scattergeo": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergeo" + } + ], + "scattergl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergl" + } + ], + "scattermap": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermap" + } + ], + "scattermapbox": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermapbox" + } + ], + "scatterpolar": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolar" + } + ], + "scatterpolargl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolargl" + } + ], + "scatterternary": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterternary" + } + ], + "surface": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "surface" + } + ], + "table": [ + { + "cells": { + "fill": { + "color": "#EBF0F8" + }, + "line": { + "color": "white" + } + }, + "header": { + "fill": { + "color": "#C8D4E3" + }, + "line": { + "color": "white" + } + }, + "type": "table" + } + ] + }, + "layout": { + "annotationdefaults": { + "arrowcolor": "#2a3f5f", + "arrowhead": 0, + "arrowwidth": 1 + }, + "autotypenumbers": "strict", + "coloraxis": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "colorscale": { + "diverging": [ + [ + 0, + "#8e0152" + ], + [ + 0.1, + "#c51b7d" + ], + [ + 0.2, + "#de77ae" + ], + [ + 0.3, + "#f1b6da" + ], + [ + 0.4, + "#fde0ef" + ], + [ + 0.5, + "#f7f7f7" + ], + [ + 0.6, + "#e6f5d0" + ], + [ + 0.7, + "#b8e186" + ], + [ + 0.8, + "#7fbc41" + ], + [ + 0.9, + "#4d9221" + ], + [ + 1, + "#276419" + ] + ], + "sequential": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "sequentialminus": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ] + }, + "colorway": [ + "#636efa", + "#EF553B", + "#00cc96", + "#ab63fa", + "#FFA15A", + "#19d3f3", + "#FF6692", + "#B6E880", + "#FF97FF", + "#FECB52" + ], + "font": { + "color": "#2a3f5f" + }, + "geo": { + "bgcolor": "white", + "lakecolor": "white", + "landcolor": "#E5ECF6", + "showlakes": true, + "showland": true, + "subunitcolor": "white" + }, + "hoverlabel": { + "align": "left" + }, + "hovermode": "closest", + "mapbox": { + "style": "light" + }, + "paper_bgcolor": "white", + "plot_bgcolor": "#E5ECF6", + "polar": { + "angularaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "radialaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "scene": { + "xaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "yaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "zaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + } + }, + "shapedefaults": { + "line": { + "color": "#2a3f5f" + } + }, + "ternary": { + "aaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "baxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "caxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "title": { + "x": 0.05 + }, + "xaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + }, + "yaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + } + } + }, + "width": 800, + "xaxis": { + "exponentformat": "power", + "title": { + "text": "Time [day]" + } + }, + "yaxis": { + "exponentformat": "power", + "title": { + "text": "Cumulative release [Bq]" + } + } + } + }, + "text/html": [ + "
\n", + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "from tritium_model import (\n", + " IV_stream,\n", + " OV_stream,\n", + " replacement_times_top,\n", + " replacement_times_walls,\n", + ")\n", + "\n", + "# plot cumulative release\n", + "fig = go.Figure()\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=time * s_to_days,\n", + " y=cumulative_top,\n", + " mode=\"lines\",\n", + " name=\"Top release\",\n", + " line=dict(color=\"#023047\"),\n", + " )\n", + ")\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=time * s_to_days,\n", + " y=cumulative_wall,\n", + " mode=\"lines\",\n", + " name=\"Wall release\",\n", + " line=dict(color=\"green\"),\n", + " )\n", + ")\n", + "\n", + "\n", + "# Scatter points for replacement times (Top)\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=replacement_times_top,\n", + " y=IV_stream.get_cumulative_activity(\"total\"),\n", + " mode=\"markers\",\n", + " marker=dict(color=\"#023047\"),\n", + " name=\"Samples (Top)\",\n", + " )\n", + ")\n", + "\n", + "# Scatter points for replacement times (Walls)\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=replacement_times_walls,\n", + " y=OV_stream.get_cumulative_activity(\"total\"),\n", + " mode=\"markers\",\n", + " marker=dict(color=\"green\"),\n", + " name=\"Samples (Walls)\",\n", + " )\n", + ")\n", + "\n", + "fig.update_layout(\n", + " legend_title=\"Components\",\n", + " font=dict(size=10),\n", + " # yaxis_type=\"log\",\n", + " # xaxis_type=\"log\",\n", + " xaxis_title=\"Time [day]\",\n", + " yaxis_title=\"Cumulative release [Bq]\",\n", + " height=600,\n", + " width=800,\n", + ")\n", + "fig.update_xaxes(exponentformat=\"power\")\n", + "fig.update_yaxes(exponentformat=\"power\")\n", + "fig.write_html(\"cumulative_release.html\")\n", + "fig.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 5, + "metadata": {}, + "outputs": [ + { + "data": { + "application/vnd.plotly.v1+json": { + "config": { + "plotlyServerURL": "https://plot.ly" + }, + "data": [ + { + "fill": "tozeroy", + "mode": "lines", + "name": "Direct release from salt", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAABuPaKwlj5RNxSwOxz+qnc3GBJga9DW0zeH2dvPDAw9OFPbxzzWi6Y4InpivFKaETl4Nn9YrIB7ORiQpYJ9fOU5obJfV0HJUDqKxfn7dTq6Oks0oCqsfSQ7Z3cggS4CkDviWbqpaAP5O+gJjsSpimM8NOtxQqmIzjzZxqk2xNo3PTXEtiHpoqI94AEvxIkeDT5TVnUVyL92PjF+lqA7xeE+XcUL6d26Sz/RrD4LXWerPzO2WoGsntk/YNJxlwzN5j/eT9LDD5L7P2mP3e7UeA1AtyhTW7Z+FUBNM4ofU0wbQCeW3+XnIyBAG1oBVHJHIkAXRIr7ih0kQNLOZAcmsSVA2fwGA6sLJ0AQ5KvzLDUoQKUoS4KaNClAwJCbVecPKkArNsuSL8wqQKBhtFjWbStA6edh65/4K0CnRP4ryG8sQGpDNOMV1ixANRcTUestLUBaN0FlVHktQC4gUPQSui1Ax9QDM6nxLUCFIxS2YiEuQMA0mdOENS5AkJ4YMY44LkBxIK6xJzsuQCzRaQ1jPS5AKGjvmE0/LkDiJ23C8kAuQPZeKFpcQi5AeBIQzZJDLkAZXv9WnUQuQO/x4y2CRS5AJCTKpkZGLkBLlKtV70YuQL4bviiARy5AvvDlf/xHLkCS19ZAZ0guQON0XOjCSC5A59gwmRFJLkD4xbkoVUkuQEmw+CmPSS5ARLXu9sBJLkATjqy460kuQFOVPm4QSi5A6Sme8i9KLkBp4ssBS0ouQHsAMT1iSi5ABUZiL3ZKLkA9m1pPh0ouQNLIQAOWSi5A0M/IoqJKLkDVET95rUouQH56Sce2Si5AKCFpxL5KLkCEX0WgxUouQM4TyIPLSi5AZ64RktBKLkCNukrp1EouQBjEV6PYSi5AZstz1ttKLkCu4LSV3kouQND9ffHgSi5Av8Xg9+JKLkCgcPG05EouQCTZDjPmSi5A+Fgge+dKLkBU5cqU6EouQGmonobpSi5AsCc+VupKLkC84H8I60ouQOYki6HrSi5AGODwJOxKLkB73sCV7EouQN4PnPbsSi5A3jXESe1KLkA+WimR7UouQNJddc7tSi5A5uQVA+5KLkBY3EQw7kouQDDKD1fuSi5AURVeeO5KLkDeafaU7kouQEdag63uSi5AyViXwu5KLkD8ILDU7kouQN2lOeTuSi5ArZWQ8e5KLkDUggT97kouQIy+2QbvSi5AWPBKD+9KLkD0c4oW70ouQK6FwxzvSi5AVkQbIu9KLkCijrEm70ouQFXBoSrvSi5AbFoDLu9KLkBFheow70ouQPGSaDPvSi5AfmGMNe9KLkChtGI370ouQOSB9jjvSi5A+zFROu9KLkDv2Ho770ouQERmejzvSi5AX85VPe9KLkAILhI+70ouQOzosz7vSi5AwsM+P+9KLkDW+rU/70ouQEdVHEDvSi5AojV0QO9KLkAYqL9A70ouQLVuAEHvSi5A3ws4Qe9KLkBSy2dB70ouQOHJkEHvSi5AE/yzQe9KLkDgM9JB70ouQIMl7EHvSi5AwmsCQu9KLkB5ixVC70ouQL32JULvSi5Afg80Qu9KLkDRKUBC70ouQOqNSkLvSi5AxHlTQu9KLkCXIltC70ouQBO2YULvSi5Ac1tnQu9KLkBiNGxC70ouQM1dcELvSi5Ah/BzQu9KLkDgAXdC70ouQCSkeULvSi5ACOd7Qu9KLkAN2H1C70ouQMaCf0LvSi5AI/GAQu9KLkCuK4JC70ouQL45g0LvSi5AmSGEQu9KLkCq6IRC70ouQJKThULvSi5ATyaGQu9KLkBJpIZC70ouQHMQh0LvSi5AUW2HQu9KLkALvYdC70ouQH4BiELvSi5ARDyIQu9KLkC5bohC70ouQAuaiELvSi5APr+IQu9KLkAt34hC70ouQJf6iELvSi5AIhKJQu9KLkBXJolC70ouQLE3iULvSi5Al0aJQu9KLkBgU4lC70ouQFteiULvSi5AyWeJQu9KLkDib4lC70ouQNR2iULvSi5Ay3yJQu9KLkDqgYlC70ouQFCGiULvSi5AGIqJQu9KLkBVjYlC70ouQB6QiULvSi5AgZKJQu9KLkCNlIlC70ouQFKWiULvSi5A1JeJQu9KLkAfmYlC70ouQD2aiULvSi5AMpuJQu9KLkADnIlC70ouQLmciULvSi5AU52JQu9KLkDanYlC70ouQEueiULvSi5Arp6JQu9KLkABn4lC70ouQEqfiULvSi5AiJ+JQu9KLkC/n4lC70ouQOyfiULvSi5AE6CJQu9KLkA0oIlC70ouQFGgiULvSi5Aa6CJQu9KLkCAoIlC70ouQJKgiULvSi5AoaCJQu9KLkCwoIlC70ouQLygiULvSi5AxaCJQu9KLkDOoIlC70ouQNWgiULvSi5A26CJQu9KLkDhoIlC70ouQOagiULvSi5A6qCJQu9KLkDtoIlC70ouQPCgiULvSi5A86CJQu9KLkD1oIlC70ouQPigiULvSi5A+aCJQu9KLkD5oIlC70ouQPugiULvSi5A/KCJQu9KLkD8oIlC70ouQP6giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A/6CJQu9KLkD/oIlC70ouQP+giULvSi5A", + "dtype": "f8" + } + }, + { + "fill": "tonexty", + "mode": "lines", + "name": "Release from isotopic exchange", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAAAAAAABuPaKwlj5RNxSwOxz+qnc3GBJga9DW0zeH2dvPDAw9OFPbxzzWi6Y4InpivFKaETl4Nn9YrIB7ORiQpYJ9fOU5obJfV0HJUDqKxfn7dTq6Oks0oCqsfSQ7Z3cggS4CkDviWbqpaAP5O+gJjsSpimM8NOtxQqmIzjzZxqk2xNo3PTXEtiHpoqI94AEvxIkeDT5TVnUVyL92PjF+lqA7xeE+XcUL6d26Sz/RrD4LXWerPzO2WoGsntk/YNJxlwzN5j/eT9LDD5L7P2mP3e7UeA1AtyhTW7Z+FUBNM4ofU0wbQCeW3+XnIyBAG1oBVHJHIkAXRIr7ih0kQNLOZAcmsSVA2fwGA6sLJ0AQ5KvzLDUoQKUoS4KaNClAwJCbVecPKkArNsuSL8wqQKBhtFjWbStA6edh65/4K0CnRP4ryG8sQGpDNOMV1ixANRcTUestLUBaN0FlVHktQC4gUPQSui1Ax9QDM6nxLUCFIxS2YiEuQMeE983/tC5ASgFY2c3iLkBdHtOpuQ8vQK5cpmXyOy9AAO7nRIZnL0DYSB9XgZIvQN3+LQ3uvC9A7chRdNXmL0CSByO0HwgwQASyJl+ZHDBAb4IStdowMEAMX99O5kQwQA8yP3K+WDBAbKdOHWVsMEBQpp8Q3H8wQOVj2NckkzBAgXUa0UCmMED0Al4zMbkwQBn85hP3yzBABNn0apPeMEDu0cgXB/EwQGiJGuRSAzFAl74Ph3cVMUBAtMindScxQDh5kN9NOTFAiBm+uwBLMUCd51G/jlwxQCB4WGT4bTFAJo8bHT5/MUDpEihVYJAxQPYXL3JfoTFAmz7H1DuyMUCO3BLZ9cIxQFDLTteN0zFA5ChNJATkMUCE4N4RWfQxQKRqLu+MBDJALNsNCaAUMkBcCTqqkiQyQLVblBtlNDJAfYpUpBdEMkClezOKqlMyQJoxkBEeYzJAV6OPfXJyMkC8NTgQqIEyQA50iQq/kDJAVo+QrLefMkAIGXo1kq4yQOtcoeNOvTJAF7Cd9O3LMkDO/k2lb9oyQFrY4jHU6DJAXi/n1Rv3MkA4/UfMRgUzQHTwWk9VEzNAykjkmEchM0Ah/hviHS8zQPJLsmPYPDNAyrbTVXdKM0Cwnizw+lczQGVu7GljZTNAUnXI+bByM0D3d/7V438zQAkBVzT8jDNA7HsnSvqZM0AZIFRM3qYzQNCzUW+oszNApisn51jAM0CvK2/n78wzQFBuWaNt2TNARBSsTdLlM0DI4MQYHvIzQJxkmjZR/jNA7hm92GsKNEBBc1gwbhY0QM3eM25YIjRA4L+zwiouNEBfX9pd5Tk0QIbUSG+IRTRAuuY/JhRRNEBA6aCxiFw0QHyR7j/mZzRATchN/yxzNEAEd4YdXX40QGFQBMh2iTRA85XXK3qUNEA12rV1Z580QJm/+tE+qjRA2rSobAC1NECprmlxrL80QAPfjwtDyjRAPWoWZsTUNEDxGaKrMN80QPENggaI6TRATmuwoMrzNECUCdOj+P00QDwePDkSCDVAiebqiRcSNUCtT4y+CBw1QHyde//lJTVAgg/DdK8vNUCyhBxGZTk1QJ0d8poHQzVARt1empZMNUCSSC9rElY1QGQE4jN7XzVAZnKoGtFoNUCCTGdFFHI1QBk/t9lEezVA8oHl/GKENUDub/TTbo01QHwdnINoljVA4O1KMFCfNUBDJyb+Jag1QJCFChHqsDVAKsyMjJy5NUByVvqTPcI1QDynWUrNyjVA+/Zq0kvTNUD6wKhOuds1QFRPSOEV5DVA4kU6rGHsNUAKLCvRnPQ1QIb1g3HH/DVA/4lqruEENkCyS8Ko6ww2QPOcLIHlFDZAsWQJWM8cNkDmkXdNqSQ2QAieVYFzLDZAcg5CEy40NkDG9Jsi2Ts2QFxug850QzZAoiLaNQFLNkCMwEN3flI2QAp7JrHsWTZAgoSrAUxhNkBbib+GnGg2QJYpE17ebzZAdXEbpRF3NkA2URJ5Nn42QOsT9/ZMhTZAZdWOO1WMNkA692RjT5M2QPCUy4o7mjZAUPfbzRmhNkDNBndI6qc2QC69RRatrjZAVJa5UmK1NkBAAA0ZCrw2QE3KQ4SkwjZAo5MrrzHJNkD8OFy0sc82QKhBOK4k1jZA7EvttorcNkCzeHTo4+I2QLTWklww6TZA8szZLHDvNkDQhKdyo/U2QJZTJ0fK+zZApyNSw+QBN0BN3e7/8gc3QFXPkhX1DTdAdBeiHOsTN0CpClAt1Rk3QKKdn1+zHzdAYs1jy4UlN0AvCECITCs3QByXqK0HMTdAQAjjUrc2N0DumQaPWzw3QCyn/Hj0QTdAwBWBJ4JHN0AkxiKxBE03QM4FRCx8UjdAZgQbr+hXN0BATLJPSl03QPQ+6SOhYjdAvZZ0Qe1nN0Br7d69Lm03QApKia5lcjdAPLarKJJ3N0C+3FVBtHw3QHaybw3MgTdA9Cu6odmGN0BUAdAS3Ys3QOmCJnXWkDdAXoEO3cWVN0CATLVeq5o3QEDMJQ6HnzdAI7hJ/1ikN0AF8upFIak3QLAJtfXfrTdAxO82IpWyN0CA3+TeQLc3QNqHGj/juzdAGH4dVnzAN0A4BSA3DMU3QJs3RPWSyTdATqKfoxDON0CWYz9VhdI3QHHhLB3x1jdAqC9zDlTbN0DSQCU8rt83QBUCZbn/4zdA/IVrmUjoN0DgaJLviOw3QHqfXs/A8DdALOiMTPD0N0D7HyB7F/k3QHvGcW82/TdAjwdFPk0BOEA8r9z8WwU4QFZ7FMFiCThAzVB+oWENOEDA7oO1WBE4QJjRjRVIFThAfxMw2y8ZOEBBN14hEB04QPjupgTpIDhAoBl5o7okOEC6YnMehSg4QG4mwJhILDhAD3x/OAUwOEB3kUEnuzM4QGTVk5JqNzhAYM6jrBM7OEBa6fqstj44QCoIV9FTQjhAwCekXutFOECmGRyifUk4QEsIkvIKTThAi1XwsZNQOED+WvBOGFQ4QM2zFEeZVzhAogfwKBdbOEBx28OWkl44QL3IhEkMYjhA+MBTFIVlOEAw2H7o/Wg4QCSxH9p3bDhAkWNyJfRvOECeEgY1dHM4QJgg8aj5djhAWxQ/X4Z6OEBodN99HH44QHtxdX6+gThAT7qLPG+FOEBCiuMGMok4QD4W4rQKjThAUrCQwf2QOEA8BU9uEJU4QHdgXvBImThAEygNrq6dOEAKz++SSqI4QJrgDocnpzhA+oLcHFOsOED4ZoiY3rE4QIPisZDgtzhAikMGqHe+OEDOmwZ4z8U4QM6dd2wqzjhAtJyx8vjXOEDqJRJAIeQ4QEbLwv6c9ThATmlwERkLOUAe/KLSnSE5QNG6vCr0QTlAk2WJnMxdOUD6fhz4tHU5QOmyZqk7ijlAGwvAINubOUC6M3d9/Ko5QJA1ggD6tzlAbe3cJiHDOUCw4Y6i8so5QHJV6geYzzlA", + "dtype": "f8" + } + }, + { + "marker": { + "color": "black" + }, + "mode": "markers", + "name": "Samples (Top)", + "type": "scatter", + "x": [ + 0.5166666666666666, + 1.2180555555555552, + 2.161805555555555, + 4.0361111111111105, + 6.997916666666665, + 10.122916666666665, + 14.02083333333333, + 18.071527777777774, + 19.26180555555555, + 21.059027777777775, + 24.13333333333333, + 28.05694444444444, + 32.130555555555546, + 35.103472222222216, + 39.0611111111111, + 42.03263888888888 + ], + "y": { + "bdata": "rUfhehSu1z8730+Nl24JQFTjpZvEoBVAbOf7qfESIUBmZmZmZqYmQGQ730+NFypA+FPjpZuELUCKQWDl0OIvQEa28/3U2DJAGQRWDi0SNED2KFyPwnU1QKrx0k1iEDZA8KfGSzeJOEDo+6nx0u04QB1aZDvfbzlAAAAAAACAOUA=", + "dtype": "f8" + } + } + ], + "layout": { + "font": { + "size": 10 + }, + "height": 600, + "legend": { + "title": { + "text": "Components" + } + }, + "template": { + "data": { + "bar": [ + { + "error_x": { + "color": "#2a3f5f" + }, + "error_y": { + "color": "#2a3f5f" + }, + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "bar" + } + ], + "barpolar": [ + { + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "barpolar" + } + ], + "carpet": [ + { + "aaxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "baxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "type": "carpet" + } + ], + "choropleth": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "choropleth" + } + ], + "contour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "contour" + } + ], + "contourcarpet": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "contourcarpet" + } + ], + "heatmap": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "heatmap" + } + ], + "histogram": [ + { + "marker": { + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "histogram" + } + ], + "histogram2d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2d" + } + ], + "histogram2dcontour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2dcontour" + } + ], + "mesh3d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "mesh3d" + } + ], + "parcoords": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "parcoords" + } + ], + "pie": [ + { + "automargin": true, + "type": "pie" + } + ], + "scatter": [ + { + "fillpattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + }, + "type": "scatter" + } + ], + "scatter3d": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatter3d" + } + ], + "scattercarpet": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattercarpet" + } + ], + "scattergeo": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergeo" + } + ], + "scattergl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergl" + } + ], + "scattermap": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermap" + } + ], + "scattermapbox": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermapbox" + } + ], + "scatterpolar": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolar" + } + ], + "scatterpolargl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolargl" + } + ], + "scatterternary": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterternary" + } + ], + "surface": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "surface" + } + ], + "table": [ + { + "cells": { + "fill": { + "color": "#EBF0F8" + }, + "line": { + "color": "white" + } + }, + "header": { + "fill": { + "color": "#C8D4E3" + }, + "line": { + "color": "white" + } + }, + "type": "table" + } + ] + }, + "layout": { + "annotationdefaults": { + "arrowcolor": "#2a3f5f", + "arrowhead": 0, + "arrowwidth": 1 + }, + "autotypenumbers": "strict", + "coloraxis": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "colorscale": { + "diverging": [ + [ + 0, + "#8e0152" + ], + [ + 0.1, + "#c51b7d" + ], + [ + 0.2, + "#de77ae" + ], + [ + 0.3, + "#f1b6da" + ], + [ + 0.4, + "#fde0ef" + ], + [ + 0.5, + "#f7f7f7" + ], + [ + 0.6, + "#e6f5d0" + ], + [ + 0.7, + "#b8e186" + ], + [ + 0.8, + "#7fbc41" + ], + [ + 0.9, + "#4d9221" + ], + [ + 1, + "#276419" + ] + ], + "sequential": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "sequentialminus": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ] + }, + "colorway": [ + "#636efa", + "#EF553B", + "#00cc96", + "#ab63fa", + "#FFA15A", + "#19d3f3", + "#FF6692", + "#B6E880", + "#FF97FF", + "#FECB52" + ], + "font": { + "color": "#2a3f5f" + }, + "geo": { + "bgcolor": "white", + "lakecolor": "white", + "landcolor": "#E5ECF6", + "showlakes": true, + "showland": true, + "subunitcolor": "white" + }, + "hoverlabel": { + "align": "left" + }, + "hovermode": "closest", + "mapbox": { + "style": "light" + }, + "paper_bgcolor": "white", + "plot_bgcolor": "#E5ECF6", + "polar": { + "angularaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "radialaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "scene": { + "xaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "yaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "zaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + } + }, + "shapedefaults": { + "line": { + "color": "#2a3f5f" + } + }, + "ternary": { + "aaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "baxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "caxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "title": { + "x": 0.05 + }, + "xaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + }, + "yaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + } + } + }, + "width": 800, + "xaxis": { + "exponentformat": "power", + "title": { + "text": "Time [day]" + } + }, + "yaxis": { + "exponentformat": "power", + "title": { + "text": "Cumulative release [Bq]" + } + } + } + }, + "text/html": [ + "
\n", + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "fig = go.Figure()\n", + "\n", + "# First trace (bottom of the stack)\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=time * s_to_days,\n", + " y=cumulative_direct_top,\n", + " mode=\"lines\",\n", + " name=\"Direct release from salt\",\n", + " fill=\"tozeroy\", # Fill from the line to the x-axis\n", + " )\n", + ")\n", + "\n", + "# Second trace (top of the stack)\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=time * s_to_days,\n", + " y=cumulative_direct_top\n", + " + cumulative_release_from_isotopic_exchange, # Add values for stacking\n", + " mode=\"lines\",\n", + " name=\"Release from isotopic exchange\",\n", + " fill=\"tonexty\", # Fill between this trace and the previous one\n", + " )\n", + ")\n", + "\n", + "# Scatter points for replacement times (Top)\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=replacement_times_top,\n", + " y=IV_stream.get_cumulative_activity(\"total\"),\n", + " mode=\"markers\",\n", + " name=\"Samples (Top)\",\n", + " marker=dict(color=\"black\"),\n", + " )\n", + ")\n", + "\n", + "# Update layout\n", + "fig.update_layout(\n", + " legend_title=\"Components\",\n", + " font=dict(size=10),\n", + " # yaxis_type=\"log\", # Uncomment if needed\n", + " # xaxis_type=\"log\", # Uncomment if needed\n", + " xaxis_title=\"Time [day]\",\n", + " yaxis_title=\"Cumulative release [Bq]\",\n", + " height=600,\n", + " width=800,\n", + ")\n", + "\n", + "# Format axes\n", + "fig.update_xaxes(exponentformat=\"power\")\n", + "fig.update_yaxes(exponentformat=\"power\")\n", + "\n", + "# Show figure\n", + "fig.show()" + ] + }, + { + "cell_type": "code", + "execution_count": 6, + "metadata": {}, + "outputs": [ + { + "data": { + "application/vnd.plotly.v1+json": { + "config": { + "plotlyServerURL": "https://plot.ly" + }, + "data": [ + { + "mode": "lines", + "name": "Piping concentration", + "type": "scatter", + "x": { + "bdata": "AAAAAAAAAACAvxagAGWZO6ZGqDTVHKs7dSlZxIRo1Tv2TcA9NTMJPLd8OnYUTj88py916K2NczzIj2SGtXCoPED2d2zWjN48EJ4SPAUYEz0ZOiByBt5HPTFniQuI1X09aPIDB3Wlsj0+sb5I0k7nPdaV7dqGIh0+KnHUSJQ1Uj7kiwlb+cKGPqzuy7G3c7w+KHUfz1LI8T5xUueCZzomPw0noWMByVs/aLhE3qBdkT8gfrZqWFnBP4+utMPy59c/AAAAAAAA4D9STQc6mufrP0XYynRB1Po/G77pZv/fA0CW8NA061UKQKcKq5DrZRBASk8Rh+GgE0Au53d919sWQOx/3nPNFhpAnBhFasNRHUCl2FWwXEYgQAMliavX4yFAW3G8plKBI0C4ve+hzR4lQBMKI51IvCZAblZWmMNZKEDEoomTPvcpQB/vvI65lCtAfzvwiTQyLUDmhyOFr88uQCNqK0CVNjBATRDFvVIFMUB2tl47ENQxQN3d3d3dNTJAwahdTVZGMkDYwd+Ww1YyQIlKWMYwZzJAOByU9Z13MkBTYM8kC4gyQCWjClR4mDJA9OVFg+WoMkDCKIGyUrkyQJBrvOG/yTJAX673EC3aMkAt8TJAmuoyQPozbm8H+zJAyXapnnQLM0CXueTN4RszQGX8H/1OLDNAND9bLLw8M0ADgpZbKU0zQNDE0YqWXTNAngcNugNuM0BsSkjpcH4zQDuNgxjejjNACdC+R0ufM0DXEvp2uK8zQKVVNaYlwDNAdJhw1ZLQM0BC26sEAOEzQBAe5zNt8TNA3mAiY9oBNECso12SRxI0QHvmmMG0IjRASSnU8CEzNEAXbA8gj0M0QOauSk/8UzRAtPGFfmlkNECBNMGt1nQ0QFB3/NxDhTRAHro3DLGVNEDs/HI7HqY0QLs/rmqLtjRAiYLpmfjGNEBXxSTJZdc0QCYIYPjS5zRA9EqbJ0D4NEDCjdZWrQg1QJHQEYYaGTVAXxNNtYcpNUAsVojk9Dk1QPuYwxNiSjVAydv+Qs9aNUCXHjpyPGs1QGZhdaGpezVANaSw0BaMNUAD5+v/g5w1QNIpJy/xrDVAoWxiXl69NUBwr52Ny801QD/y2Lw43jVADjUU7KXuNUDed08bE/81QK26ikqADzZAff3Fee0fNkBOQAGpWjA2QB6DPNjHQDZA78V3BzVRNkDACLM2omE2QJJL7mUPcjZAZY4plXyCNkA60WTE6ZI2QA4UoPNWozZA41bbIsSzNkC6mRZSMcQ2QJLcUYGe1DZAbB+NsAvlNkBIYsjfePU2QCalAw/mBTdAB+g+PlMWN0DqKnptwCY3QNFttZwtNzdAvLDwy5pHN0Cr8yv7B1g3QKE2Zyp1aDdAnXmiWeJ4N0CgvN2IT4k3QKz/GLi8mTdAwkJU5ymqN0DkhY8Wl7o3QBPJykUEyzdAUwwGdXHbN0ClT0Gk3us3QA6TfNNL/DdAj9a3ArkMOEAuGvMxJh04QO9dLmGTLThA2aFpkAA+OEDx5aS/bU44QEAq4O7aXjhAz24bHkhvOECms1ZNtX84QNP4kXwikDhAYz7Nq4+gOEBohAjb/LA4QPXKQwpqwThAIBJ/OdfROEADWrpoROI4QL6i9Zex8jhAcewwxx4DOUBIN2z2ixM5QHGDpyX5IzlAJNHiVGY0OUCjIB6E00Q5QDpyWbNAVTlAQMaU4q1lOUAbHdARG3Y5QEV3C0GIhjlASdVGcPWWOUDGN4KfYqc5QH2fvc7PtzlARw35/TzIOUAmgjQtqtg5QEX/b1wX6TlAAIari4T5OUDqF+e68Qk6QN62IupeGjpAAGVeGcwqOkDRJJpIOTs6QDr51XemSzpAn+URpxNcOkD27U3WgGw6QNYWigXufDpAnmXGNFuNOkCK4AJkyJ06QOGOP5M1rjpAG3l8wqK+OkAYqbnxD886QFQq9yB93zpAMQo1UOrvOkBIWHN/VwA7QMEmsq7EEDtAwYrx3TEhO0DrnDENnzE7QPB5cjwMQjtAPEO0a3lSO0C7H/ea5mI7QMM8O8pTcztAIs+A+cCDO0BZFMgoLpQ7QA1UEVibpDtAreFchwi1O0BoHqu2dcU7QG97/OXi1TtAnXxRFVDmO0CCu6pEvfY7QPjqCHQqBzxAVNtso5cXPEA8f9fSBCg8QE/xSQJyODxAu3rFMd9IPEDqmkthTFk8QG4Q3pC5aTxAbuN+wCZ6PEC+cTDwk4o8QAF99R8BmzxAHjvRT26rPEBbacd/27s8QLhi3K9IzDxA4jkV4LXcPEB613cQI+08QFAdC0GQ/TxAdg/Xcf0NPUAJBOWiah49QOraP9TXLj1Alz70BUU/PUC67xA4sk89QCwdp2ofYD1Ab8rKnYxwPUAWR5PR+YA9QM25GwZnkT1AUcKDO9ShPUAUNvBxQbI9QPP7i6muwj1AIQyJ4hvTPUA1miEdieM9QF5wmVn28z1AuoM/mGMEPkBXy2/Z0BQ+QMhklR0+JT5AEBMtZas1PkDcJsiwGEY+QFPhDwGGVj5Aq2XJVvNmPkAXUdqyYHc+QGgUThbOhz5AVi9cgjuYPkCGcm/4qKg+QLpzLnoWuT5AcWaFCYTJPkDGk7Go8dk+QPS1Tlpf6j5AH4dmIc36PkAn4IIBOws/QKPTwv6oGz9Az0LzHRcsP0D8fqtkhTw/QEuibtnzTD9AXWbSg2JdP0CiYKxs0W0/QMixRp5Afj9AXWKdJLCOP0Ay2qUNIJ8/QFscommQrz9AgreBSwHAP0Dyq1LJctA/QGLnxPzk4D9ALGXDA1jxP0D6QJMA5gBAQKzTwY4gCUBAHpEPxFsRQEBH4AK8lxlAQHpfqpbUIUBAI9BbeRIqQEC1Y5KPUTJAQO+O8QuSOkBAzl9yKdRCQECQTMIsGEtAQKqU22VeU0BAh57gMadbQECVRUX98mNAQETVUkZCbEBAgYEVoJV0QECsmcS17XxAQK2HuU5LhUBAxf0MU6+NQECLnPXQGpZAQKr3BwOPnkBAsl99Vw2nQEAPdK14l69AQAaN7lYvuEBAPv0aNNfAQEAmowWxkclAQF4zON1h0kBAm2pqSUvbQEDw7TwdUuRAQFn55zB77UBARGjBK8z2QEA4iMWoSwBBQD4MtGIBCkFAGLjkafYTQUA+08tmNR5BQEAXX+3KKEFADIVe58UzQUA8WkkeOD9BQNGDEvI2S0FAHB98UNxXQUBOcTQMSGVBQLMR28Whc0FA1SPtuhuDQUCujPUN95NBQGXq15mKpkFAXBCDc067QUBSMy7N8NJBQDdgwQCA7kFA2uE/ec0PQkCJ/wEnnjpCQJLZ6hdue0JAkVcVpP7QQkDipS8w5DJDQEWA/qeU0UNADQoQ37NwRECZEfMX1A9FQFxFL1P0rkVA8fJwjhRORkBKrbLJNO1GQMNn9ARVjEdAOiI2QHUrSEBth0o3matIQAAAAAAAAElA", + "dtype": "f8" + }, + "y": { + "bdata": "AAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEAAABwNX30QQAAAHA1ffRBAAAAcDV99EEDAABwNX30QQYBAHA1ffRBMGYAcDV99EFx6idwNX30Qecxl381ffRBXTCQhTt99EFawf8ut370Qehmk9h7iPRB4wym80CR9EFDhhFOna30Qc8Rk86Q5PRBqluBQesT9UFZIb9ykzz1QVvx3IB7X/VBcdRvo3N99UG6rXOfLpf1QQ8Q4PBFrfVB39cyXj3A9UFTf30KhtD1QQV7RBiB3vVBbknq7IHq9UFZfSYi0PT1QQkwIDGp/fVB85Af4UEF9kF072GCxwv2QUk5ZvxgEfZBN0n+tS8W9kHgU41cUBr2Qcr0F4/bHfZBKUMiceYg9kF7AsgogyP2QcXs3aWdJPZBZ8kW/jn49UGsP4FcR8z1Qeg3VRCooPVB6lFk51t19UFD+bzZYkr1QZ1Ki8m8H/VBXa/ZhWn19EFSP3PNaMv0QR6KXlG6ofRBSz3+tl149EG3ZOSZUk/0QY44Y42YJvRBe9/kHS/+80FXNxLSFdbzQWCSzytMrvNBB1wVqdGG80EJwanEpV/zQYm9v/bHOPNB+FR/tTcS80F/L3p19OvyQVRjD6r9xfJB/cvAxVKg8kFG+3s683ryQX2F2HneVfJBCStN9RMx8kHIKl0ekwzyQcHZvmZb6PFBc3R8QGzE8UF4+w8exaDxQafPenJlffFBJahZsUxa8UFCZ/VOejfxQf8/UcDtFPFBG442e6by8EGdtD72o9DwQfBK26jlrvBBi9ZcC2uN8EEtRviWM2zwQVlcy8U+S/BBRDDgEowq8EHB5i/6GgrwQR2ASfHV0+9BpCY3GPeT70HqmclmmFTvQW8Eity4Fe9BIFPqelfX7kGwCkVFc5nuQYiY3EALXO5BPjPadB4f7kFrW0zqq+LtQY8LJayypu1BfaM3xzFr7UELmjZKKDDtQT7+sEWV9exB3c8PzHe77EEkNpPxzoHsQWOaT8yZSOxBpqsqdNcP7EFxT9gCh9frQXSD15Onn+tBITNvRDho60H5A6szODHrQc8aWIKm+upB7dsBU4LE6kHTp+7Jyo7qQQOWHA1/WepBAC8+RJ4k6kGcJreYJ/DpQXYXmTUavOlBVkCgR3WI6UEpRDD9N1XpQRbtUIZhIulBG/OqFPHv6EGoxoTb5b3oQXxfvw8/jOhBERDT5/ta6EHKXcybGyroQSLeSGWd+edB7Rh0f4DJ50HkbwQnxJnnQX0LOJpnaudBQM3RGGo750GKRxbkygznQeO6yD6J3uZB4BgobaSw5kGlDOy0G4PmQfsIQl3uVeZBBlzKrhsp5kGaSJXzovzlQR4lIHeD0OVBGYBShryk5UFBSntvTXnlQR4GToI1TuVBN/3fD3Qj5UG8eqVqCPnkQaILb+bxzuRBQMRm2C+l5EFHiw2XwXvkQR5qOHqmUuRBm+IN290p5EEFSgMUZwHkQV0p2oBB2eNB4aKdfmyx40HJ3J9r54njQSRxd6exYuNB3OL8kso740HJF0iQMRXjQdXXrQLm7uJBJVG9TufI4kE0oT3aNKPiQeNiKwzOfeJBb0G2TLJY4kE7kD4F4TPiQXnnUqBZD+JBkcWtiRvr4UFTNTMuJsfhQcx47vt4o+FB37gPYhOA4UFvuenQ9FzhQSmS77kcOuFB42uyj4oX4UF8Qt/FPfXgQTKrPNE10+BBg5+oJ3Kx4EFiTBZA8o/gQdvli5K1buBBCX8gmLtN4EFe5vnKAy3gQTSGSqaNDOBBJJOeTLHY30FZCpuQyJjfQRTJIRVgWd9Bn7bS2HYa30EqNVLcC9zeQawURSIent5BuI1Mr6xg3kFBRAKKtiPeQUxS9Lo6591BTlqhTDir3UFLoXRLrm/dQYUwwsWbNN1Btv7Cy//53EGjIJFv2b/cQQUBJMUnhtxBlZ9M4ulM3EE317HeHhTcQQmrzNPF29tBSprk3N2j20HT+gsXZmzbQSFaHKFdNdtBmuSym8P+2kH90iwpl8jaQaXdo23XktpBb7XqjoNd2kEHgom0mijaQTJlugcc9NlB1QJmswbA2UFJDSDkWYzZQYDVI8gUWdlBjt5QjzYm2UHucydrvvPYQexBxY6rwdhBXO/hLv2P2EHdt8uBsl7YQY8FZL/KLdhBMwkcIUX910FPT/HhIM3XQetQaj5dnddBIv6SdPlt10GUQPnD9D7XQWJzqW1OENdBJs4qtAXi1kHSv3vbGbTWQfc1DimKhtZBdszD41VZ1kHm4OlTfCzWQWKENcP8/9VBeUW/fNbT1UEqy/7MCKjVQak4xgGTfNVBUFA9anRR1UHdStxWrCbVQSBWZhk6/NRBhKzjBB3S1EFwNJttVKjUQeWTC6nfftRBxKHjDb5V1EF8Gvrz7izUQf55RLRxBNRBA9jMqEXc00FKnqYsarTTQTzq4ZvejNNBHGR9U6Jl00F4TVaxtD7TQRh/FhQVGNNBOgMg28Lx0kE37HVmvcvSQZj4oRYEptJBjYSWTJaA0kGsNIxpc1vSQequ2s6aNtJBLZzL3QsS0kEsDmf3xe3RQblBOHzIydFBtYwJzBKm0UGrF5ZFpILRQTvML0Z8X9FBlaJXKZo80UFXL0ZI/RnRQQX/Yfmk99BBUuugj5DV0EFYJtFZv7PQQcE2xqEwktBBs4plq+Nw0EGNm42z10/QQVTW0e4LL9BBX5gEiH8O0EHkEQ89Y9zPQU/TsohCnM9Bvjau95pcz0FtCUhqah3PQfp1WJau3s5BCXxhAGWgzkFHcZfzimLOQTABsHgdJc5BFedHSxnozUH6fajNeqvNQcYZrvo9b81BjMiGVV4zzUHWbPbW1vfMQQ3mwNehvMxBpPDM97iBzEFxOoIBFUfMQQ9q0sitDMxBsPVGBXrSy0EzzVMmb5jLQVffDyGBXstB86tONqIky0GN3OavwurKQZWAtZPQsMpBX3eeSrd2ykGFg2A5XzzKQYM+gkitAcpBqJrSV4LGyUElMtmYuorJQZDp+MksTslBxMe3SakQyUHx+0D2+NHIQbNPQsjbkchBb6fnEAdQyEGP0JQ3IwzIQUcw1sLIxcdBqjnoXHx8x0FX7RJYqS/HQcGheOyZ3sZBwOIN5WuIxkE7RmN7/yvGQUKp8zfcx8VBy6TjiQNaxUGtserZnt/EQa90HVJbVMRB0qitqeeww0ECTZzISObCQf45JivzwsFB64mKPAVdwEGwBTZZ2su9QVT57ye1l7lBvsoaYA35tUGEVL80cd2yQc8Org1dMrBB8hLRrczPq0Ets1iCxeCnQdKeOZUtgKRBUthShd+ZoUFIr79mZiKfQT9MOLIrt5xB", + "dtype": "f8" + } + } + ], + "layout": { + "font": { + "size": 10 + }, + "legend": { + "title": { + "text": "Components" + } + }, + "template": { + "data": { + "bar": [ + { + "error_x": { + "color": "#2a3f5f" + }, + "error_y": { + "color": "#2a3f5f" + }, + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "bar" + } + ], + "barpolar": [ + { + "marker": { + "line": { + "color": "#E5ECF6", + "width": 0.5 + }, + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "barpolar" + } + ], + "carpet": [ + { + "aaxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "baxis": { + "endlinecolor": "#2a3f5f", + "gridcolor": "white", + "linecolor": "white", + "minorgridcolor": "white", + "startlinecolor": "#2a3f5f" + }, + "type": "carpet" + } + ], + "choropleth": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "choropleth" + } + ], + "contour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "contour" + } + ], + "contourcarpet": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "contourcarpet" + } + ], + "heatmap": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "heatmap" + } + ], + "histogram": [ + { + "marker": { + "pattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + } + }, + "type": "histogram" + } + ], + "histogram2d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2d" + } + ], + "histogram2dcontour": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "histogram2dcontour" + } + ], + "mesh3d": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "type": "mesh3d" + } + ], + "parcoords": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "parcoords" + } + ], + "pie": [ + { + "automargin": true, + "type": "pie" + } + ], + "scatter": [ + { + "fillpattern": { + "fillmode": "overlay", + "size": 10, + "solidity": 0.2 + }, + "type": "scatter" + } + ], + "scatter3d": [ + { + "line": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatter3d" + } + ], + "scattercarpet": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattercarpet" + } + ], + "scattergeo": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergeo" + } + ], + "scattergl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattergl" + } + ], + "scattermap": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermap" + } + ], + "scattermapbox": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scattermapbox" + } + ], + "scatterpolar": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolar" + } + ], + "scatterpolargl": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterpolargl" + } + ], + "scatterternary": [ + { + "marker": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "type": "scatterternary" + } + ], + "surface": [ + { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + }, + "colorscale": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "type": "surface" + } + ], + "table": [ + { + "cells": { + "fill": { + "color": "#EBF0F8" + }, + "line": { + "color": "white" + } + }, + "header": { + "fill": { + "color": "#C8D4E3" + }, + "line": { + "color": "white" + } + }, + "type": "table" + } + ] + }, + "layout": { + "annotationdefaults": { + "arrowcolor": "#2a3f5f", + "arrowhead": 0, + "arrowwidth": 1 + }, + "autotypenumbers": "strict", + "coloraxis": { + "colorbar": { + "outlinewidth": 0, + "ticks": "" + } + }, + "colorscale": { + "diverging": [ + [ + 0, + "#8e0152" + ], + [ + 0.1, + "#c51b7d" + ], + [ + 0.2, + "#de77ae" + ], + [ + 0.3, + "#f1b6da" + ], + [ + 0.4, + "#fde0ef" + ], + [ + 0.5, + "#f7f7f7" + ], + [ + 0.6, + "#e6f5d0" + ], + [ + 0.7, + "#b8e186" + ], + [ + 0.8, + "#7fbc41" + ], + [ + 0.9, + "#4d9221" + ], + [ + 1, + "#276419" + ] + ], + "sequential": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ], + "sequentialminus": [ + [ + 0, + "#0d0887" + ], + [ + 0.1111111111111111, + "#46039f" + ], + [ + 0.2222222222222222, + "#7201a8" + ], + [ + 0.3333333333333333, + "#9c179e" + ], + [ + 0.4444444444444444, + "#bd3786" + ], + [ + 0.5555555555555556, + "#d8576b" + ], + [ + 0.6666666666666666, + "#ed7953" + ], + [ + 0.7777777777777778, + "#fb9f3a" + ], + [ + 0.8888888888888888, + "#fdca26" + ], + [ + 1, + "#f0f921" + ] + ] + }, + "colorway": [ + "#636efa", + "#EF553B", + "#00cc96", + "#ab63fa", + "#FFA15A", + "#19d3f3", + "#FF6692", + "#B6E880", + "#FF97FF", + "#FECB52" + ], + "font": { + "color": "#2a3f5f" + }, + "geo": { + "bgcolor": "white", + "lakecolor": "white", + "landcolor": "#E5ECF6", + "showlakes": true, + "showland": true, + "subunitcolor": "white" + }, + "hoverlabel": { + "align": "left" + }, + "hovermode": "closest", + "mapbox": { + "style": "light" + }, + "paper_bgcolor": "white", + "plot_bgcolor": "#E5ECF6", + "polar": { + "angularaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "radialaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "scene": { + "xaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "yaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + }, + "zaxis": { + "backgroundcolor": "#E5ECF6", + "gridcolor": "white", + "gridwidth": 2, + "linecolor": "white", + "showbackground": true, + "ticks": "", + "zerolinecolor": "white" + } + }, + "shapedefaults": { + "line": { + "color": "#2a3f5f" + } + }, + "ternary": { + "aaxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "baxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + }, + "bgcolor": "#E5ECF6", + "caxis": { + "gridcolor": "white", + "linecolor": "white", + "ticks": "" + } + }, + "title": { + "x": 0.05 + }, + "xaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + }, + "yaxis": { + "automargin": true, + "gridcolor": "white", + "linecolor": "white", + "ticks": "", + "title": { + "standoff": 15 + }, + "zerolinecolor": "white", + "zerolinewidth": 2 + } + } + }, + "xaxis": { + "exponentformat": "power", + "title": { + "text": "Time [day]" + } + }, + "yaxis": { + "exponentformat": "power", + "title": { + "text": "Block value [Bq]" + } + } + } + }, + "text/html": [ + "
\n", + "
" + ] + }, + "metadata": {}, + "output_type": "display_data" + } + ], + "source": [ + "fig = go.Figure()\n", + "\n", + "fig.add_trace(\n", + " go.Scatter(\n", + " x=time * s_to_days,\n", + " y=data[1],\n", + " mode=\"lines\",\n", + " name=\"Piping concentration\",\n", + " )\n", + ")\n", + "\n", + "fig.update_layout(\n", + " legend_title=\"Components\",\n", + " font=dict(size=10),\n", + " # yaxis_type=\"log\",\n", + " # xaxis_type=\"log\",\n", + " xaxis_title=\"Time [day]\",\n", + " yaxis_title=\"Block value [Bq]\",\n", + ")\n", + "fig.update_xaxes(exponentformat=\"power\")\n", + "fig.update_yaxes(exponentformat=\"power\")\n", + "fig.write_html(\"baby.html\")\n", + "fig.show()" + ] + } + ], + "metadata": { + "kernelspec": { + "display_name": "baby-1l-run-3", + "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.12.9" + } + }, + "nbformat": 4, + "nbformat_minor": 2 +} diff --git a/analysis/tritium/multi_system_tritium_model.py b/analysis/tritium/multi_system_tritium_model.py new file mode 100644 index 0000000..6a2594a --- /dev/null +++ b/analysis/tritium/multi_system_tritium_model.py @@ -0,0 +1,162 @@ +######################################################################################### +## +## ARC fuel cycle +## +######################################################################################### + +# IMPORTS =============================================================================== + +from pathsim import Simulation, Connection +from pathsim.blocks import ODE, Source, Scope, Block +from pathsim.solvers import RKBS32, RKF21, RKCK54 +from pathsim.blocks.mixed import Pulse, Step, SquareWave + +from tritium_model import ( + k_top, + k_wall, + baby_model, + gas_switch_deltatime, + irradiations, +) + +# PARAMETERS ============================================================================ + +A_wall = baby_model.A_wall +A_top = baby_model.A_top +V_salt = baby_model.volume +TBR = baby_model.TBR +neutron_rate = baby_model.neutron_rate + +TBR = TBR.to("particle/neutron").magnitude + +A_top = A_top.to("m^2").magnitude +A_wall = A_wall.to("m^2").magnitude +V_salt = V_salt.to("m^3").magnitude + +k_top = k_top.to("m/s").magnitude + +k_top *= 1.1 + +k_wall = 4 * k_top +gas_switch_deltatime = gas_switch_deltatime.to("seconds").magnitude +neutron_rate = neutron_rate.to("neutrons/s").magnitude + +total_duration = 50 * 24 * 3600 +irradiation_time = ( + sum([irr[1] - irr[0] for irr in irradiations]).to("seconds").magnitude +) + +f = 0.05 # fraction of tritium that sticks to the piping (this is minimal compared to IC) +initial_piping_IV_concentration = 5.5e9 # initial concentration in the piping IV +initial_piping_OV_concentration = 2e9 # initial concentration in the piping OV +k_exchange = 10 * k_top +k_exchange_OV = 8 * k_top + + +# BLOCKS ================================================================================ + + +def salt_ode(x, u, t): + # unpack states + c_salt = x + # unpack inputs + neutron_rate_in, h2_in = u + dc_dt = ( + 1 + / V_salt + * (TBR * neutron_rate_in - c_salt * (k_top * A_top + h2_in * k_wall * A_wall)) + ) + return dc_dt + + +def piping_ode(x, u, t): + # unpack states + c_piping = x + # unpack inputs + c_salt_in, h2_in = u + dc_dt = f * c_salt_in * k_top * A_top - c_piping * h2_in * k_exchange + return dc_dt + + +def piping_OV_ode(x, u, t): + # unpack states + c_piping = x + # unpack inputs + h2_in = u + dc_dt = -c_piping * h2_in * k_exchange_OV + return dc_dt + + +salt = ODE(func=salt_ode) + +piping = ODE(func=piping_ode, initial_value=initial_piping_IV_concentration) + +piping_wall = ODE(func=piping_OV_ode, initial_value=initial_piping_OV_concentration) + +neutron_rate_block = Pulse( + amplitude=neutron_rate, + T=total_duration, + duty=irradiation_time / total_duration, +) +h2 = Step(amplitude=1, tau=gas_switch_deltatime) + +Sco = Scope(labels=["salt", "piping", "piping_wall", "neutron_rate", "h2"]) + +blocks = [neutron_rate_block, salt, piping, piping_wall, h2, Sco] + + +# CONNECTIONS =========================================================================== + +connections = [ + Connection(salt, piping[0], Sco[0]), + Connection(piping, Sco[1]), + Connection(piping_wall, Sco[2]), + Connection(neutron_rate_block, salt[0], Sco[3]), + Connection(h2, salt[1], piping[1], piping_wall[0], Sco[4]), +] + + +# SIMULATION ============================================================================ + +Sim = Simulation( + blocks, + connections, + log=True, + Solver=RKBS32, + tolerance_lte_rel=1e-4, + tolerance_lte_abs=1e-9, +) + + +if __name__ == "__main__": + + Sim.run(total_duration) + + # Sim.save("baby.mdl") + + # fig, ax = Sco.plot() + + # plt.show() + + import plotly.graph_objects as go + + # plot total plant inventory + time, data = Sco.read() + + fig = go.Figure() + for p, d in enumerate(data): + lb = Sco.labels[p] if p < len(Sco.labels) else f"port {p}" + fig.add_trace(go.Scatter(x=time, y=d, mode="lines", name=lb)) + + fig.update_layout( + legend_title="Components", + font=dict(size=10), + # yaxis_type="log", + # xaxis_type="log", + xaxis_title="Time [s]", + yaxis_title="Block value [#]", + ) + fig.update_xaxes(exponentformat="power") + fig.update_yaxes(exponentformat="power") + fig.write_html("baby.html") + # fig.show() diff --git a/environment.yml b/environment.yml index f98d80a..9228c12 100644 --- a/environment.yml +++ b/environment.yml @@ -23,3 +23,4 @@ dependencies: - openmc_data_downloader - matplotlib-label-lines - h-transport-materials + - pathsim