ogs-jupyter-lab-vpf1.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "id": "7cea3d8c",
   "metadata": {
    "tags": []
   },
   "source": [
    "**ogs-jupyter-lab: DECOVALEX 2023: Task G: Benchmarking**"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "0b601bd1",
   "metadata": {
    "tags": []
   },
   "source": [
    "<!--- ![tu-dresden-blue.png](attachment:tu-dresden-blue.png) --->\n",
    "<img src=\"ogs-jupyter-lab+decovalex2023.png\" alt=\"drawing\" width=\"300\"/>"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6064292a",
   "metadata": {},
   "source": [
    "**Description**\n",
    "\n",
    "..."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "7fdafdb4",
   "metadata": {},
   "source": [
    "**Governing equations** (template for TeX integration)\n",
    "\n",
    "The governing equations ... are given by\n",
    "\n",
    "$$\\phi \\rho_L\\alpha_p\\dfrac{\\partial p}{\\partial t}-\\nabla\\left(\\rho_L\\frac{\\mathbf{k}}{\\mu}(\\nabla p -\\rho_L \\mathbf{g})\\right)-\\rho_L\\left(\\phi \\alpha_T^L +3 (1-\\phi) \\alpha_T^s\\right) \\dfrac{\\partial T}{\\partial t}=Q_H$$\n",
    "for the fluid mass balance, and ...\n",
    "\n",
    "where \n",
    "* $p$ is the water pressure,\n",
    "* ..."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "6293ebdb",
   "metadata": {},
   "source": [
    "**Boundary conditions**\n",
    "\n",
    "..."
   ]
  },
  {
   "cell_type": "markdown",
   "id": "aa2e71bd",
   "metadata": {},
   "source": [
    "**Input data:** Inspecting vtu-file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "5907758b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "inspecting vtu-file:filename\n"
     ]
    },
    {
     "data": {
      "text/html": [
       "<table><tr><th>Header</th><th>Data Arrays</th></tr><tr><td>\n",
       "<table>\n",
       "<tr><th>UnstructuredGrid</th><th>Information</th></tr>\n",
       "<tr><td>N Cells</td><td>2500</td></tr>\n",
       "<tr><td>N Points</td><td>3636</td></tr>\n",
       "<tr><td>X Bounds</td><td>0.000e+00, 1.000e+00</td></tr>\n",
       "<tr><td>Y Bounds</td><td>0.000e+00, 5.000e-02</td></tr>\n",
       "<tr><td>Z Bounds</td><td>0.000e+00, 5.000e-02</td></tr>\n",
       "<tr><td>N Arrays</td><td>1</td></tr>\n",
       "</table>\n",
       "\n",
       "</td><td>\n",
       "<table>\n",
       "<tr><th>Name</th><th>Field</th><th>Type</th><th>N Comp</th><th>Min</th><th>Max</th></tr>\n",
       "<tr><td>ispatch</td><td>Fields</td><td>int32</td><td>1</td><td>0.000e+00</td><td>0.000e+00</td></tr>\n",
       "</table>\n",
       "\n",
       "</td></tr> </table>"
      ],
      "text/plain": [
       "UnstructuredGrid (0x7fe4ce15f0a0)\n",
       "  N Cells:\t2500\n",
       "  N Points:\t3636\n",
       "  X Bounds:\t0.000e+00, 1.000e+00\n",
       "  Y Bounds:\t0.000e+00, 5.000e-02\n",
       "  Z Bounds:\t0.000e+00, 5.000e-02\n",
       "  N Arrays:\t1"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import pyvista as pv\n",
    "filename = \"bar.vtu\"\n",
    "mesh = pv.read(filename)\n",
    "print(\"inspecting vtu-file:\" 'filename')\n",
    "mesh"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 26,
   "id": "df1db57b",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAABI8AAABZCAYAAAC33QQjAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAbuElEQVR4nO3deZxU1Z338e+PplsEWWRRGRYbZREMotgCalxwidA9ji+X5wmYVZOHh0TGaDKZIBo1KtJJJqOSiVE0xsdMEpKJy5CAKGoUHTcao6CA2GILuCIYUJCl6fP8UdXVdW9Xdd1be3V93q9Xvaxf3XPuPXWk6p769bnnmnNOAAAAAAAAQCJdCt0AAAAAAAAAFC+SRwAAAAAAAEiK5BEAAAAAAACSInkEAAAAAACApEgeAQAAAAAAICmSRwAAAAAAAEiK5BEAAACKmpndY2YfmtmrSbabmc03s0YzW2Vm4/PdRgAAOjOSRwAAACh290qa0sH2qZJGRB8zJP0yD20CAKBskDwCAABAUXPOLZe0rYMi50m6z0U8L6mPmQ3MT+sAAOj8ugYpZGZTJN0mqULS3c65et92i26vlbRL0tedcy9FtzVJ+kTSfknNzrmaVMfr37+/q66uDv4uAAAoYStXrvzIOTeg0O0AStggSZvi4s3R197zFzSzGYrMTlKPHj2OP+qoo/LSQAAAikG6486UySMzq5D0C0lnK3IiXmFmi5xza+KKxU8VnqjIVOGJcdsnO+c+Ctqo6upqNTQ0BC0e2vbP9qn3gZVJ4zD1P93TrOb9LerTvUqStKd5v3bvbVHv7sn35z/e6+9/olGH9YzFb2/dqcP79Uha/7O9+9Wli3RA1wpJ0qZtu9Sne6V6dquMxUP6dg/8Hva3OG3ctkvD+rcd8/3tu3VY726B38PGrbs0tF/bMbft3Ku+PaoC13/375/pH/ocGIs/2b0v9n6C1N++a58OqOyibpWRPmlpcdq7vyUWB9nH1k/3qN9BB8TifftbVFmRfHKev/7HO/fq4Lj3vL/FqaKLhXoP8f9unHOK5GUD1v9sn3p16+qp09E+nHPasbs5a5+FUqy/t7lFzS0t6l4V+SpsaXH6ZE9zSb2HbNSP/3eze99+OScdWJX8s5OLNpRy/Wwws7fzekCg80l0snOJCjrnFkhaIEk1NTUul2NOAACKTbrjziAzjyZIanTObYgeaKEiU4Pjk0exqcKSnjezPmY20DnX7q89+VY9e7Ekqam+Trv2NmvMtY/EYv/2MPFb82plZvrcdZH93TbtWJ137CCNumapJOn4ww/W/d86KVa+dR+5jk+uf0Lv/P2zpO9h5dvbdOEvn4vFcx5crd+9sDEWz75/lRau2BS4TxYsf1M3L1knSWqcO1W3Pf6Gfv5EoyRp1fVfUK9ulbHyj155qkYe2jMW9+rWVbOnjtacB1fH3sO/fGGk/u3R9YHf808vOkbf/9OqWHx17WjNXbLWU/4vq97VrN/9LRb/ZOk63f7km5KkW744Tlf+4ZVY+evPHaPr/7zGU7+jPrht2rH6zsKXY+V/9E9H67pFrwWu//1zRumnj7yujqTqg1SyUT++zU31ddrf4nTknCVp14+PU7nirBG64qyRRdcHTR/t1On/9mTa9ePjVFq/b1rLX3nWSF1+5nANuyrY/4PWY8Yf7615taHrx7d5w821OuNnT6pp66606r95c63uenqD6h9el1b9xrlTdfczbwWu37qP+D741ymj9JOlHX/+/PVv+PMa3fM/b3naBKAobJY0JC4eLOndArUFAIBOJ0jyKNE04IkByrROFXaSHjUzJ+nO6F972omfQjx06NBAjQ+rdXaBJP3uhY26eGLbcb75/1bo7q+dEHhfw65aouXfnxyLv7PwZQ3o2TZr5YqzRqTcx7yH13ric3/+TODjS+1/eAb5IXr84X0D108188VffvjVD3vilU0fa/JRh8Tivc0tnu07djd7EkeSPImjIOITR5I8iaNWu/bs98TnHzcoljyKTxxJ8iSOkpn/+Bux5/GJI0mexJEk7di9T73iZlDt29+iY65/NBanShwlEibp0Rnq/3HFJl1x1siM9vHpnuaM6r+08WNdcPuzade//clGT5IibP2JNz+uDz/ZE4tveWy9bnks3GfFf8wwiaNE9Y8ImDxMVj9o8jFZff/3TTr7CJM4avXBjt2h6wDIi0WSZkX/yDlR0vZi+CMmAACdRZDkUZBpwB2VOdk5966ZHSJpmZmtiy566C3sm0IcoF0ZmfPgao067KBY/NjaD3XxXc+H2sepP/2rJ774rhdizzu6XKrVnU9t8MSr39ke6vi5lihxtNP3I7wj8YkjSepaYTrN12dhpZO4OKCy7dKzfCc+enWr1DNvtF2xOSILP3jzUX/pFadoyq1Pp13/ji8fr5n/uTLt+oMPPlCbP47MoHt3++6M+6B1hmC64hNH6UgnSREvPnGE/Jly9GFa+tr7hW4GAElm9ntJp0vqb2abJV0nqVKSnHN3SFqiyNqbjYqsv3lJYVoKAEDnFORua0GmASct45xr/e+Hkh5U5DK4otB6+VarZ9/cWqCWpO/QXgd44hvOO7rD8ql+hL8xd2qHdY9O8SN82glDkm6bcuvTejvgJS6ZOCpu/ajxNy5rNzsolfhLUapnLw6duPj119tmsFXPXqwv/+qFDkoXh7/98GxP3Jo4Cuqcow/1xK2Jo3S1Jo6AXEp02Vn8Wm0kjoDi4Zyb7pwb6JyrdM4Nds79yjl3RzRxpOhd1i5zzh3pnBvrnGMhIwAAsihI8miFpBFmNszMqiRNU2RqcLxFkr5qEZMUnSpsZj3MrKckmVkPSV+Q9GoW219Wtn+2r91rH+zwzkjY6btEK16QJEg6M2Pita6XJEUu0crHrJkzfDOc1r3/Sez5tp17U9bv51vYO9M2X3Lviozq59qgPgfqe2d7LwM77sZlgetfcNygdq898toHGbcLyDXWKAIAAADSkzJ55JxrljRL0iOS1kr6o3PuNTObaWYzo8WWSNqgyFThuyR9O/r6oZKeMbNXJL0oabFzbmmW30PZGPejR1OW+fHS4IvHpvLvj76eUSIl00RUUE+s+zBU+cWXf94Tbw2QYGr160tOyOgHaFN9nX77Tf+SYbl14fjBnvidv3+mny0Lt15OvAf+9k6mTQIy9uKcM3OSDAqScAYAAADKTZA1j+ScW6JIgij+tTvinjtJlyWot0HSuAzb2Gmc+bMn83q8R197XzN+k/7lQ/Ojd00rJmHvkpWo/uaPw1069+zsM3RS/ROSpEt+HW5W0fqbpmrrzj06cV6kfqazmtJx/0ub835MINv8d1s7pFe3QjYHAAAAKCtBLltDlry5ZWdej5dJ4iiZJ753Wuz59AnJ1zdKZOkVp2Q8U2BPc/LL8hJpqq/TbdOOjcXVsxfr8z8Ovmj3ZZOPjCWO0jHymodjiSMAbdbfNNXzfTCwtzcZNOPUI0LtL36tMUm67twx6TcOAAAAgEfZJo+WXH5KVvfX0ULT6Uo1S6VnN+/EsSDJnEevPDX2/PWbpoRqT1N9nY4Y0HaHur6+tYISmXfBWE8c9jbX/mTTqGvCXfVYPXtxqAWzn7/qTE/8i7++Gep4+VBZkejmhkBx8X8nNtXXeT7PVV29p5/nfJ+9ObWjQx2vogufCwAAACBXyjZ5dOEvw996+5PdbQtW+38YVVZ0aZfMSeU33/DeeG79TeESUKuvP8cT//7FTUlKRjTV12nkoW13Ijuga0XKY1xTl/wHXKrESlN9naZPGBqLp9z6tCbe/HiHdZZe4U3qhb3Mq6P2JuJPTk2a13H7JOk/v5HfNYv89u13BT0+ILVPBvk/S5UV+T+9XHJydey5c9KvvlYTqr7/PbDANgAAABBhkeWKiktNTY1raMjOHVZbkw9N9XX6ZPc+jb0+9aLTreXj6zfOnarhHSwAffLwfvrtNyd5kh2p1ufxbw+7nk/Y8m/Nq9Wwq5akLhi3fyn9dXoWfOX4nFw6l28De3fTe9vDzZgC0Pkk+s7OFjNb6ZwLl+0CkLFsjjkBACgF6Y47y2rmUc9ulYHLtrR4k2qpZnv8T+NWz4+Kgw7omjLp4t8eNkmTqvzJw/tp3Y1tl6aFSRy18vdDGJ0hcSSJxBEAAAAAoKyFu86qxL0fIgnQxbd+xuhrw62107dHlT7d0xyqTrb99puTdNfyDWnXL8SdwQCgWFR17aK9zS2e1+rGDtTi1e8VqEUAAABAYZTVzKObFq/J27E2bgt3O/hs+dLEtjWGqmcv1twla0PVfzzubmoA0JkEWdMo/kYA/sQRAAAAUK7KKnn0l1WF/WvxPx4zMOfHCDO7KpEzf/ZUlloCALl12sgBWd8nCSMAAACgvbK6bK3Q8pG8enzdhzk/BgDkm38B/0wX9E+m0JcbAwAAAMWorGYeAQAKZ9bk4bHnTfV1Wb1bGQAAAIDcIXnUgf0d3Gns6trRnti3vnZRmnv+50qinQA6n6b6Ov3LOaOSbp8+YYgnvv1L40Pt/75LJ3ji684dE6o+AAAAgORIHnXgyDnJb23vX4g6gzva583VD75aEu0EUHz8M4X88fqbprYrH8a8C47xxLVjc79GXCok2wEAAIAIkkdJtBRhlmXdjVNClff/JR8AWqW681iQ5I9zbd+TVV0LezppbmnR/OnHxeLtn+3TEyHvHul/zxvmte+DnXtZEwkAAADlx+IH/8WipqbGNTQ0ZGVf8YurplpYdd2NU3TUD5dm5bgAgM7Ffx7J5ppNZrbSOVeTtR0CCCSbY04AAEpBuuNOZh7FmfGblYVuAgCgiHWvqih0EwAAAIC8K5vk0UsbP05ZZvn6LXloCQCgGB3Wq5v69qjqsMzkUYfkqTUA/Mxsipm9bmaNZjY7wfbTzWy7mb0cfVxbiHYCANAZlU3yaM4DqwvdBABAAaVa1+n9Hbvz2RwAIZhZhaRfSJoqaYyk6WaW6LaKTzvnjo0+bshrIwEA6MTKJnm07v1PCt0EAECaUt3tLVu27dyb9X0CyIoJkhqdcxucc3slLZR0XoHbBABA2Sib5BEAoHDylfwB0GkNkrQpLt4cfc3vRDN7xcweNrOj89M0AAA6P5JHAICUvnX6kXruqjNiMckfAHlmCV7z3zL4JUmHO+fGSfq5pIcS7shshpk1mFnDli2sdwkAQBAkjwCgDG24uTZU+R9MOUoDex+Yo9ZIt3xxnCd+7Lunhap/36UTPPF15yZaCgVACdssaUhcPFjSu/EFnHM7nHOfRp8vkVRpZv39O3LOLXDO1TjnagYMGJDLNgMA0GmQPAKAEtSne2XseWWF9w/yJ1Qf3K68f5ZQly7eOheOH+yJ/cmcVFItRp3K+cd5jz/8kINC1c+Fs8ccWugmAGizQtIIMxtmZlWSpklaFF/AzA4zM4s+n6DIOHdr3lsKAEAnRPIIAArMn+x5cc6ZnnjDzbWeZMzYQb01a/LwWLxvv9PN54+NxSuaPta1/9jxzJs/rNjoic8a7b0FvfNfDCJpf0vbiw+8tFmX3rsiFj+x7gNP2T+/4pkQ0M6yNR/ou394ORY/++ZHnu1Pvv5hh/WfWr9F1zzUdhfNZWs+0D+N+4dYvGD5Bv2v4wcnqhpz01/WeGL/bKXRA3u1q/Pmlk873CeA3HDONUuaJekRSWsl/dE595qZzTSzmdFiF0l61cxekTRf0jTnEn2bAQCAsKwYz6k1NTWuoaEhK/uqnr04K/sBAJS3pvo6zzklm2s+mdlK51xN1nYIIJBsjjkBACgF6Y47mXkEAEBA/XpUFboJAAAAQN6RPAIAlIVXrvuCJ040c6hviuTQpCP6ZbVNAAAAQCkgeQQAKHpN9XWeZI8/DqL3gZWpCwEAAABoh+QRACDnspH8yYdtO/cWugkAAABA0SF5BABlyH8JVyqlkvwBAAAAkH0kjwCgyP3XzBP1xtypsfjLk4aGTtz4y3MJFwAAAICgSB4BQJaNG9LHE6eatbPmhnPalY93QnVfVVYk/7qurDBf+YNDt7nQ/vLPn/fEq68PNzPqvksneOLrzh2TcZsAAAAARJA8AlD27v/WSZ547Q1TPHGq5I8/2fPfl50c6vjdq7qGKp8L8y4Y64n9yZxU/H0QdmbU5wb19sQ9u+V/ZtSEYX098TV1o/PeBgAAAKAYmXOu0G1op6amxjU0NGRlX9WzF2dlPwBK11VTj9K8h9elXf/yM0folU1/11Prt6RVf9bk4Vr1znYtT7P+/z31CK3avF3PbdiaVn1JOmZwb63avD3t+mMG9tKa93akXf/wft319tZdadfv16NKWwu8mHVTfZ1O+ckT2rTts1icLWa20jlXk7UdAggkm2NOAABKQbrjzsL/uRsAciyTxJEkzX/8jYzq/8dfGzOqf+fyDRnVl5RR4khSRokjSRkljiQVPHHUqjVxBAAAAJQTLlsDACCgsb7L6wAAAIByQPIIAFAUunbxLvy9+PK2dZdSrTslSX+aeaInvrrWu2ZRkHWZ+vao6rCNQ/t273A7AAAA0BmRPAKAMrRwxqR2r339pGpP/H9OGRZ7HiR589IPz/bE8cmcIPWbW7xr8NXNfybFu/C66I7nPPHcJWtD1QcAAACQGMkjAMixIImTMPW/PGlo6MWa/eUnHdGvXZl7n23yxHc9/VaoY4y/cZkn9idzSsG2IllbCQAAACgmJI8AlJyHLjs5VPlMbyOfafIHAAAAAEoZySMAHcr2rJmw9f91yqh25Y8d0idw/XLUrZKv9v4Hdbx2EQAAAIDgAv3CMLMpZva6mTWa2ewE283M5ke3rzKz8UHrAuXk3//3OC278tRY7E+K3PP1mnZ1sjFrxh/f+ZXjQ+2jI7d8cZwnfuy7p4Wqn+n7KwZnjT7UE8+7YGyo+m/Nq/XEYfugssK70HQx3BHshOqDQ5X/m2+9pA031yYpmdh9l07wxJdNHh6qPgAAAIDkUiaPzKxC0i8kTZU0RtJ0MxvjKzZV0ojoY4akX4aoC2TNrV881hM31dfptmnHemL/LBh/eX/c0ayZsPUvGD9YIw7tGfTtBHLh+MGeeP7040LVzzR5c/5x3uMPP+SgUPVzwZ/MSeXVH53jif3JnFTu/po36Td9wtBQ9c0sdaEce/DbJ3nil689O0nJxPz/bv5r5klJSiZ2sO8uZ1265L9Pvnf2SE/8hwSLigMAAADlyJxzHRcwO1HS9c65c6LxVZLknJsXV+ZOSU86534fjV+XdLqk6lR1E6mpqXENDQ3pvaM41bMXZ7wPAEB5aqqv85xHmurrdNlvX9Li1e/F4mwxs5XOufZTDwHkVLbGnAAAlIp0x51BLlsbJGlTXLw5+lqQMkHqSpLMbIaZNZhZw5YtWwI0CwAAAOUik2UUAABAZoIkjxJdO+CfrpSsTJC6kRedW+Ccq3HO1QwYMCBAs1IrxbVTAACZ63lAVx3er3ssrqro4okleeLxQ/uo/oKxuvKsyKVrv4peivjinDMlta1HVn/hWB07pI9+HvLyUACZyWQZBQAAkLmuAcpsljQkLh4s6d2AZaoC1M0pEkgAgDC+c9aI2PNDenXznEd6dqvUQ5edXIhmAeVugqRG59wGSTKzhZLOk7Qmrsx5ku5zkTUZnjezPmY20Dn3Xv6bCwBA5xIkebRC0ggzGybpHUnTJF3sK7NI0qzoiXyipO3OuffMbEuAuu2sXLnyIzN7O8T7SKW/pI+yuL9yRB9mjj7MHH2YHfRj5rLdh4dncV9AZ5RoKYSJAcoMkuRJHpnZDEVmJknSHjN7NbtNhQ/nnPygn3OPPs49+jg/RqVTKWXyyDnXbGazJD0iqULSPc6518xsZnT7HZKWSKqV1Chpl6RLOqob4JjZuW4tyswaWIg0M/Rh5ujDzNGH2UE/Zo4+BPIuk2UUvC84t0DSAonPcj7Qx/lBP+cefZx79HF+mFlad4oIMvNIzrkliiSI4l+7I+65k3RZ0LoAAABACJksowAAADIUZMFsAAAAoJBiyyiYWZUiSyEs8pVZJOmr0buuTVJ0GYV8NxQAgM4o0MyjTmBBoRvQCdCHmaMPM0cfZgf9mDn6EMijTJZRSIHPcu7Rx/lBP+cefZx79HF+pNXPFrniDAAAAAAAAGiPy9YAAAAAAACQFMkjAAAAAAAAJNVpkkdmNsXMXjezRjObnWC7mdn86PZVZja+EO0sdgH68UvR/ltlZs+a2bhCtLOYperDuHInmNl+M7son+0rBUH60MxON7OXzew1M3sq320sdgE+y73N7M9m9kq0D4OsDVJWzOweM/vQzF5Nsp3zClAiGCfmHmPI3GOMmR+MQ3OPcWru5WQc65wr+YciCye+KekISVWSXpE0xlemVtLDkkzSJEkvFLrdxfYI2I8nSTo4+nwq/Ri+D+PKPaHI4p4XFbrdxfQI+O+wj6Q1koZG40MK3e5iegTswzmSfhx9PkDSNklVhW57MT0knSppvKRXk2znvMKDRwk8GCcWTR8zhsxxH8eVY4yZw35mHJqXPmacmnk/Z30c21lmHk2Q1Oic2+Cc2ytpoaTzfGXOk3Sfi3heUh8zG5jvhha5lP3onHvWOfdxNHxe0uA8t7HYBfm3KEn/LOl+SR/ms3ElIkgfXizpAefcRklyztGPXkH60EnqaWYm6SBFTsrN+W1mcXPOLVekX5LhvAKUBsaJuccYMvcYY+YH49DcY5yaB7kYx3aW5NEgSZvi4s3R18KWKXdh++gbimQr0SZlH5rZIEnnS7ojj+0qJUH+HY6UdLCZPWlmK83sq3lrXWkI0of/IWm0pHclrZb0HedcS36a12lwXgFKA+PE3GMMmXuMMfODcWjuMU4tDqHPe11z2pz8sQSvuTTKlLvAfWRmkxU58X8+py0qPUH68FZJP3DO7Y8k0+ETpA+7Sjpe0pmSDpT0nJk975xbn+vGlYggfXiOpJclnSHpSEnLzOxp59yOXDeuE+G8ApQGxom5xxgy9xhj5gfj0NxjnFocQp/3OkvyaLOkIXHxYEWylGHLlLtAfWRmx0i6W9JU59zWPLWtVATpwxpJC6Mn9f6Sas2s2Tn3UH6aWPSCfp4/cs7tlLTTzJZLGieJk3ZEkD68RFK9i1z03Ghmb0k6StKL+Wlip8B5BSgNjBNzjzFk7jHGzA/GobnHOLU4hD7vdZbL1lZIGmFmw8ysStI0SYt8ZRZJ+mp0VfFJkrY7597Ld0OLXMp+NLOhkh6Q9BWy6wml7EPn3DDnXLVzrlrSnyR9m5O6R5DP839LOsXMuppZd0kTJa3NczuLWZA+3KjIX8xkZodKGiVpQ15bWfo4rwClgXFi7jGGzD3GmPnBODT3GKcWh9DnvU4x88g512xmsyQ9osjq7fc4514zs5nR7XcocseBWkmNknYpks1EnID9eK2kfpJuj/5Vo9k5V1OoNhebgH2IDgTpQ+fcWjNbKmmVpBZJdzvnEt6GshwF/Hd4o6R7zWy1ItNWf+Cc+6hgjS5CZvZ7SadL6m9mmyVdJ6lS4rwClBLGibnHGDL3GGPmB+PQ3GOcmh+5GMdaZCYYAAAAAAAA0F5nuWwNAAAAAAAAOUDyCAAAAAAAAEmRPAIAAAAAAEBSJI8AAAAAAACQFMkjAAAAAAAAJEXyCAAAAAAAAEmRPAIAAAAAAEBS/x8UxRZBkeSNAwAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 1440x72 with 2 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "import matplotlib.pyplot as plt\n",
    "import numpy as np\n",
    "import vtk\n",
    "from vtk.util.numpy_support import vtk_to_numpy\n",
    "import matplotlib.tri as tri\n",
    "\n",
    "#file\n",
    "reader = vtk.vtkXMLUnstructuredGridReader()\n",
    "reader.SetFileName(filename)\n",
    "reader.Update()  # Needed because of GetScalarRange\n",
    "data = reader.GetOutput()\n",
    "temperature = data.GetPointData().GetArray(\"ispatch\")\n",
    "#points\n",
    "points = data.GetPoints()\n",
    "npts = points.GetNumberOfPoints()\n",
    "x = vtk_to_numpy(points.GetData())\n",
    "triang = tri.Triangulation(x[:,0], x[:,1])\n",
    "#cells\n",
    "triangles=  vtk_to_numpy(data.GetCells().GetData())\n",
    "ntri = triangles.size//4  # number of cells\n",
    "tri = np.take(triangles,[n for n in range(triangles.size) if n%4 != 0]).reshape(ntri,3)\n",
    "#plots\n",
    "fig, ax = plt.subplots(ncols=2, figsize=(20,1))\n",
    "#plt.subplots_adjust(left=0.1, bottom=0.1, right=0.9, top=0.9, wspace=0.4, hspace=0.6)\n",
    "plt.subplots_adjust(wspace=0.5)\n",
    "#mesh\n",
    "ax[0].triplot(x[:,0], x[:,1], tri)\n",
    "#ax[1].tricontourf(x[:,0], x[:,1], tri, phasefield, 12);\n",
    "#contour_right = ax[1].tricontourf(triang, temperature)\n",
    "#fig.colorbar(contour_right,ax=ax[1],label='$T$ / [K]')\n",
    "plt.savefig(\"figures.png\")"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "8683ff4f",
   "metadata": {},
   "source": [
    "**Running ogs using [ogs6py](https://github.com/joergbuchwald/ogs6py)**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 41,
   "id": "bfd2a8b9",
   "metadata": {},
   "outputs": [],
   "source": [
    "#import ogs6py\n",
    "#ogs6py.__version__"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "id": "409ac11f",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['E',\n",
       " 'nu',\n",
       " 'k',\n",
       " 'gc',\n",
       " 'ls',\n",
       " 'rho_sr',\n",
       " 'displacement0',\n",
       " 'phasefield_ic',\n",
       " 'phasefield_bc',\n",
       " 'dirichlet0',\n",
       " 'Dirichlet_spatial',\n",
       " 'dirichlet_right_time',\n",
       " 'dirichlet_right']"
      ]
     },
     "execution_count": 64,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "from ogs6py import ogs\n",
    "PROJECT_FILE = \"AT1_iso_tensile.prj\"\n",
    "model = ogs.OGS(INPUT_FILE=PROJECT_FILE,PROJECT_FILE=PROJECT_FILE)\n",
    "#print(\"PROCESS: \")\n",
    "#model.tree.find(\"./processes/process/type\").text\n",
    "#print(\"PARAMETERS: \")\n",
    "[i.text for i in model.tree.findall(\"./parameters/parameter/name\")]\n",
    "#print(\"PROCESS: \")\n",
    "#[i.text for i in model.tree.findall(\"./processvars/parameter/name\")]\n",
    "#[i.text for i in model.tree.findall(\"./nonlinsolvers/parameter/name\")]\n",
    "#[i.text for i in model.tree.findall(\"./linsolvers/parameter/name\")]\n",
    "#[i.text for i in model.tree.findall(\"./parameters/parameter/name\")]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 65,
   "id": "85143f21-f2ef-43ab-a2a3-9412f1ca55fa",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "===============\n",
      ">>> run ogs <<<\n",
      ">>> ogs is running in background <<<\n",
      "OGS finished with project file AT1_iso_tensile.prj.\n",
      "Execution took 27.632313013076782 s\n"
     ]
    }
   ],
   "source": [
    "import multiprocessing\n",
    "from ogs6py import ogs\n",
    "#run ogs\n",
    "PATH_OGS = \"/home/ok/ogs/build/release/bin/\"\n",
    "PATH_SINGULARITY = \"/usr/bin/\"\n",
    "PATH_CONTAINER = \"/home/ok/ogs/container/gcc-10-openmpi-4.0.5-system-cmake-f0a156cb-4a3437aac8bc.sif\"\n",
    "PROJECT_FILE = \"AT1_iso_tensile.prj\"\n",
    "print(\"===============\")\n",
    "print(\">>> run ogs <<<\")\n",
    "model = ogs.OGS(INPUT_FILE=PROJECT_FILE,PROJECT_FILE=PROJECT_FILE)\n",
    "model.geo.add_geom(filename=\"bar.gml\")\n",
    "model.mesh.add_mesh(filename=\"bar.vtu\")\n",
    "model.mesh.add_mesh(filename=\"bar_left.vtu\")\n",
    "model.mesh.add_mesh(filename=\"bar_right.vtu\")\n",
    "#model.run_model?\n",
    "#model.run_model(path=PATH_OGS,logfile=\"console.log\")\n",
    "#model.run_model(path=PATH_SINGULARITY,container_path=PATH_CONTAINER,logfile=\"console.log\")\n",
    "# listing model as process for parallel processing\n",
    "ogs_process = multiprocessing.Process(\n",
    "    target=model.run_model,\n",
    "    kwargs={\"path\":PATH_SINGULARITY, \"container_path\":PATH_CONTAINER,\"logfile\":\"example.log\"})\n",
    "# start of ogs in background\n",
    "ogs_process.start()\n",
    "print(\">>> ogs is running in background <<<\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "c48523e9-9ae0-4bbb-baba-aad09a976cf4",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\u001b[32minfo:\u001b[m This is OpenGeoSys-6 version 6.4.1-200-gc33c20aba.\n",
      "\u001b[32minfo:\u001b[m OGS started on 2021-11-16 00:59:16+0100.\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: E\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: nu\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: k\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: gc\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: ls\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: rho_sr\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: displacement0\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: phasefield_ic\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: phasefield_bc\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: dirichlet0\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: Dirichlet_spatial\n",
      "[0] \u001b[32minfo:\u001b[m CurveScaledParameter: dirichlet_right_time\n",
      "[0] \u001b[32minfo:\u001b[m ConstantParameter: dirichlet_right\n",
      "[0] \u001b[32minfo:\u001b[m No source terms for process variable 'phasefield' found.\n",
      "[0] \u001b[32minfo:\u001b[m No source terms for process variable 'displacement' found.\n",
      "[0] \u001b[33m\u001b[1mwarning:\u001b[m Secondary variable type specification is deprecated and is ignored. All secondary variable types are 'static'.\n",
      "[0] \u001b[33m\u001b[1mwarning:\u001b[m Secondary variable type specification is deprecated and is ignored. All secondary variable types are 'static'.\n",
      "[0] \u001b[32minfo:\u001b[m The equations of the coupled processes will be solved by the staggered scheme.\n",
      "[0] \u001b[32minfo:\u001b[m Initialize processes.\n",
      "[0] \u001b[32minfo:\u001b[m Solve processes.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 0 took 0.0621726 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #1 and time 0.1 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m Calculate non-equilibrium initial residuum.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.1488 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00871912 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 303 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.427557 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=1.5738e+01, |dx|/|x|=9.8517e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.585541 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.118436 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000828393 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 4 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0159251 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.0087e-13, |x|=1.5738e+01, |dx|/|x|=1.2764e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.135351 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.120344 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000820027 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.011203 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.5738e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.132518 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -6.01492e-20\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.856445 s in time step #1  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.7612000000010184e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03599640599999998 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04419282099999999 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.051612856999999845 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.049909031999999964 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.206691 s in time step #1  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.128649 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000978125 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0133358 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.5738e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.143131 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -6.01492e-20\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.146126 s in time step #1  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=4.5000e-01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.4964000000071707e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.031081702999999905 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.043420400999999886 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035690922999999763 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.05224714399999986 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.173348 s in time step #1  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.000253125 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #1 took 1.38833 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 1 took 0.092069 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #2 and time 0.2 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.104048 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00100366 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 303 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.390684 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=3.1476e+01, |dx|/|x|=4.9259e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.495929 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.113097 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000747725 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 20 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0513308 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.4405e-12, |x|=3.1476e+01, |dx|/|x|=1.4108e-13\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.16538 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.109534 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000767087 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0112058 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=3.1476e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.121666 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -1.20298e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.785938 s in time step #2  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.9086999999705654e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.029596267000000065 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.041678102000000106 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.032314787000000234 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.038268903999999715 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.150062 s in time step #2  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0987623 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000910991 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0213097 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=3.1476e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.121122 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -1.20298e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.12442 s in time step #2  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=9.0000e-01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.7526000000221984e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03516180499999999 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036079285999999655 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.05832658100000021 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.056907786000000016 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.199183 s in time step #2  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.0010125 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #2 took 1.271 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 2 took 0.0890288 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #3 and time 0.3 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0989069 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000993042 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 298 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.377236 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=4.7213e+01, |dx|/|x|=3.2839e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.477304 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.106244 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000798444 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 17 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0305558 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.7266e-12, |x|=4.7213e+01, |dx|/|x|=3.6571e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.13777 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.114289 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000753666 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0113191 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=4.7213e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.126513 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -1.80447e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.744549 s in time step #3  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.0959999999691235e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.02993981400000001 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03363847899999994 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03805871200000022 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03854531799999972 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.150152 s in time step #3  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.106846 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00103468 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0158727 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=4.7213e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.123892 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -1.80447e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.127417 s in time step #3  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.3500e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.863899999994146e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03707012800000076 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03554856399999995 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.042201109000000514 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.05124199299999965 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.178154 s in time step #3  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.00227813 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #3 took 1.20604 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 3 took 0.0626093 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #4 and time 0.4 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.105192 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000971091 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 300 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.377893 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=6.2951e+01, |dx|/|x|=2.4629e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.484244 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.10846 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000812096 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 16 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0308438 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.2574e-12, |x|=6.2951e+01, |dx|/|x|=3.5859e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.140263 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.117557 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00324897 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0154317 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3024e-14, |x|=6.2951e+01, |dx|/|x|=2.0689e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.136587 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -2.40597e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.76476 s in time step #4  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 5.117200000004374e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04825873899999955 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.042078348999999626 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04570064400000007 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.041370275000000234 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.186731 s in time step #4  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.109603 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00175183 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0168842 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=6.2951e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.128397 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -2.40597e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.131999 s in time step #4  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.8000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.93260000000717e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04496691300000055 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.05133136400000016 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.060472294000000204 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034980031999999994 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.201341 s in time step #4  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.00405 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #4 took 1.29116 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 4 took 0.0743965 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #5 and time 0.5 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.115817 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000992558 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 298 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.400734 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=7.8689e+01, |dx|/|x|=1.9703e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.517714 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.108298 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000832903 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 11 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0313598 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.1507e-12, |x|=7.8689e+01, |dx|/|x|=1.4623e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.140769 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.128123 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000798076 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0241362 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3363e-14, |x|=7.8689e+01, |dx|/|x|=1.6982e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.153215 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -3.00746e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.815008 s in time step #5  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.7830999999700623e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0352782859999996 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04469962199999955 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.054499412999999386 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03761084499999967 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.185396 s in time step #5  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.112193 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.0013532 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0144155 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=7.8689e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.128134 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -3.00746e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.13469 s in time step #5  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=2.2500e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00013230500000016576 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04940320499999995 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.042711254999999504 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04335519199999993 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03753405399999998 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.183731 s in time step #5  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.00632812 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #5 took 1.32523 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 5 took 0.0565403 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #6 and time 0.6 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.119365 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00159765 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 300 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.433089 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=9.4427e+01, |dx|/|x|=1.6420e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.554255 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.108764 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000766186 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 13 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0341535 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5592e-12, |x|=9.4427e+01, |dx|/|x|=1.6512e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.143842 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.1204 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000821957 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0126555 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.2857e-14, |x|=9.4427e+01, |dx|/|x|=1.3616e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.134056 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -3.60895e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.8351 s in time step #6  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.8183999999598086e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04032171200000079 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033380671999998945 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036060642999999004 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.039733433999998624 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.159973 s in time step #6  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.10978 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00121831 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0159625 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=9.4427e+01, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.127103 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -3.60895e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.130922 s in time step #6  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=2.7000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.086799999962864e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.042681095000000724 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03869844299999947 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033952890000000124 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.040343527999999296 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.165907 s in time step #6  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.0091125 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #6 took 1.29774 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 6 took 0.0578221 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #7 and time 0.7 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.112138 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00169041 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 296 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.447535 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=1.1016e+02, |dx|/|x|=1.4074e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.561612 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.12205 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00107263 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 11 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0250214 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3221e-12, |x|=1.1016e+02, |dx|/|x|=1.2001e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.148313 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.120396 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000771061 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0124042 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.2793e-14, |x|=1.1016e+02, |dx|/|x|=1.1613e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.133775 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -4.21044e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.846591 s in time step #7  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.472500000116895e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.02961705400000092 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.032372738999999484 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04411965999999978 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03386196699999999 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.148567 s in time step #7  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.112614 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.0021198 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0220938 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.1016e+02, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.137052 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -4.21044e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.142394 s in time step #7  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=3.1500e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.632599999965237e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03955556999999921 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04114966000000031 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04101439899999981 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03609813699999975 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.168492 s in time step #7  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.0124031 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #7 took 1.3125 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 7 took 0.0547083 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #8 and time 0.8 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.115209 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00162135 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 299 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.59924 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=1.2590e+02, |dx|/|x|=1.2315e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.716252 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.119379 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000803514 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 14 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0322045 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3320e-12, |x|=1.2590e+02, |dx|/|x|=1.0580e-14\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.152567 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.112379 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000986353 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0140203 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4104e-14, |x|=1.2590e+02, |dx|/|x|=1.1202e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #3 took 0.127586 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -4.81193e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.999318 s in time step #8  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.4101999998981114e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04136958000000135 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03349003299999964 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035411746999999494 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03435480999999996 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.153417 s in time step #8  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.109279 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00187397 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 0 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0273621 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.2590e+02, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.138727 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -4.81193e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.143072 s in time step #8  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=3.6000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 7.229700000088712e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03521785899999941 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.038994731000000726 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03796789999999994 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03939244500000072 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.161677 s in time step #8  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.0162 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #8 took 1.46359 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 8 took 0.063698 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #9 and time 0.9 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.109872 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.0010554 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 298 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.592739 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=1.4164e+02, |dx|/|x|=1.0946e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.703867 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.119843 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000916598 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 8 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0280563 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=7.5439e-13, |x|=1.4164e+02, |dx|/|x|=5.3261e-15\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.14906 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -5.41342e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.856508 s in time step #9  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 5.952000000064572e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04325937000000124 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0435224630000004 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0358954839999992 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03256576299999914 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.164644 s in time step #9  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.118112 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000983543 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0175526 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3348e-14, |x|=1.4164e+02, |dx|/|x|=9.4239e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.136854 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -5.41342e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.140484 s in time step #9  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=8.8818e-16, |x|=4.0500e+00, |dx|/|x|=2.1930e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.641500000028941e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.043357549999999634 s.\n",
      "  0 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03674489900000033 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0343178489999989 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03752306100000169 s.\n",
      "  1 SNES Function norm 0.000000000000e+00 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.161308 s in time step #9  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=0.0000e+00, |x|=1.0000e+00, |dx|/|x|=0.0000e+00\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 0.0205031 Surface energy: 2.65325e-33 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #9 took 1.32866 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 9 took 0.064621 s.\n",
      "[0] \u001b[32minfo:\u001b[m === Time stepping at step #10 and time 1 with step size 0.1\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.108089 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00612235 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 303 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.585852 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5504e+01, |x|=1.5738e+02, |dx|/|x|=9.8517e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.700497 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.118615 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000840602 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 12 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.045731 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4447e-12, |x|=1.5738e+02, |dx|/|x|=9.1798e-15\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.165396 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: -6.01492e-19\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.869784 s in time step #10  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 5.536399999961361e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04163715999999873 s.\n",
      "  0 SNES Function norm 6.716165200470e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.049117755999999346 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03602839200000041 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04359306300000121 s.\n",
      "  1 SNES Function norm 1.069517467844e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.193897 s in time step #10  coupling iteration #0\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.102459 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00120034 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 273 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.532068 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.6961e-01, |x|=1.5765e+02, |dx|/|x|=3.6132e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.635906 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.120149 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.0013276 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0159863 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3669e-14, |x|=1.5765e+02, |dx|/|x|=8.6707e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.137651 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.000790744\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.776976 s in time step #10  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4780e-02, |x|=4.5000e+00, |dx|/|x|=3.2846e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.643199999991964e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.037877673000000556 s.\n",
      "  0 SNES Function norm 3.381192434862e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034998733000000115 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.032754799000001 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04647594299999902 s.\n",
      "  1 SNES Function norm 3.653110572006e-06 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0443613850000002 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03356572500000077 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.043803211000000175 s.\n",
      "  2 SNES Function norm 1.072486884155e-06 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036974803000001444 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034546600999998844 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03272463400000092 s.\n",
      "  3 SNES Function norm 9.175467397608e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.424079 s in time step #10  coupling iteration #1\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.0054e-02, |x|=1.0000e+00, |dx|/|x|=4.0054e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.111873 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00140303 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 276 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.518927 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.0047e+00, |x|=1.5810e+02, |dx|/|x|=6.3547e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.632374 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.129547 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000988897 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0156024 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4997e-14, |x|=1.5810e+02, |dx|/|x|=9.4854e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.146312 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00088077\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.782242 s in time step #10  coupling iteration #2\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.6003e-02, |x|=4.5000e+00, |dx|/|x|=5.7785e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.081900000052485e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03377088799999939 s.\n",
      "  0 SNES Function norm 4.861844557462e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04434351400000125 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.048433790999999005 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03387656299999975 s.\n",
      "  1 SNES Function norm 1.373820206958e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03939523900000097 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03691925799999751 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03496954299999899 s.\n",
      "  2 SNES Function norm 1.274109059977e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.299756 s in time step #10  coupling iteration #2\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.6379e-02, |x|=1.0000e+00, |dx|/|x|=4.6379e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.110973 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00150812 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 277 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.533665 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.6961e+00, |x|=1.5881e+02, |dx|/|x|=1.0680e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.646362 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.120156 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00099737 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.019497 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3777e-14, |x|=1.5881e+02, |dx|/|x|=8.6751e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.140844 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00106288\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.791002 s in time step #10  coupling iteration #3\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.4562e-02, |x|=4.5000e+00, |dx|/|x|=9.9026e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.613799999868661e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03724708000000021 s.\n",
      "  0 SNES Function norm 6.670609256774e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0334813849999982 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0337235540000016 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035554773000001205 s.\n",
      "  1 SNES Function norm 1.469227100668e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.043029872999998275 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03352749399999766 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0360774159999977 s.\n",
      "  2 SNES Function norm 4.802779938447e-06 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03718399000000261 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.044952053999999464 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03515473600000263 s.\n",
      "  3 SNES Function norm 9.445085585250e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.42765 s in time step #10  coupling iteration #3\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.6594e-02, |x|=1.0000e+00, |dx|/|x|=5.6594e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.112668 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00167377 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 280 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.542485 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.8291e+00, |x|=1.5993e+02, |dx|/|x|=1.7690e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.657032 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.117518 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000916678 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0138422 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5297e-14, |x|=1.5993e+02, |dx|/|x|=9.5649e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.132409 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00141312\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.792848 s in time step #10  coupling iteration #4\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=7.5873e-02, |x|=4.5000e+00, |dx|/|x|=1.6861e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.8502999999346912e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.049202477999997996 s.\n",
      "  0 SNES Function norm 1.062084643379e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04075983699999952 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03495436999999768 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03527365800000126 s.\n",
      "  1 SNES Function norm 2.226851568766e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034473284000000604 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03652481499999993 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.037312789999997875 s.\n",
      "  2 SNES Function norm 1.214737270257e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.299553 s in time step #10  coupling iteration #4\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=9.4383e-02, |x|=1.0000e+00, |dx|/|x|=9.4383e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0996377 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00150928 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 281 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.542888 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.9432e+00, |x|=1.6185e+02, |dx|/|x|=3.0541e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.64423 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.137363 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000970608 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0146655 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5142e-14, |x|=1.6185e+02, |dx|/|x|=9.3553e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.153151 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00205739\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.800878 s in time step #10  coupling iteration #5\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.3528e-01, |x|=4.5000e+00, |dx|/|x|=3.0063e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.5362000001603064e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04185888199999965 s.\n",
      "  0 SNES Function norm 1.741749978545e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04543599299999812 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03877931999999973 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036280427000001225 s.\n",
      "  1 SNES Function norm 4.590918481631e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03528239199999916 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0329142460000007 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03552358099999964 s.\n",
      "  2 SNES Function norm 7.123391485989e-06 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03977782300000143 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04350298199999969 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03318988200000206 s.\n",
      "  3 SNES Function norm 1.004060532190e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.426946 s in time step #10  coupling iteration #5\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5475e-01, |x|=1.0000e+00, |dx|/|x|=1.5475e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.107949 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00137253 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 288 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.542126 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=8.8680e+00, |x|=1.6540e+02, |dx|/|x|=5.3616e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.651647 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.134231 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000848249 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.015594 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5237e-14, |x|=1.6540e+02, |dx|/|x|=9.2119e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.150844 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00328218\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.807609 s in time step #10  coupling iteration #6\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.4801e-01, |x|=4.5000e+00, |dx|/|x|=5.5114e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.265200000119762e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04294709700000254 s.\n",
      "  0 SNES Function norm 3.093312026126e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04084205300000221 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04110325399999937 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0388607659999991 s.\n",
      "  1 SNES Function norm 6.871417971535e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035606720000000536 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03291571299999774 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035311105999998205 s.\n",
      "  2 SNES Function norm 5.548207329722e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.294702 s in time step #10  coupling iteration #6\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.1129e-01, |x|=1.0000e+00, |dx|/|x|=2.1129e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.102189 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00121708 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 301 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.575034 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4886e+01, |x|=1.7173e+02, |dx|/|x|=8.6682e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.67861 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.122063 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00118194 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0167546 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.9798e-14, |x|=1.7173e+02, |dx|/|x|=1.1528e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.140201 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00549594\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.823373 s in time step #10  coupling iteration #7\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.3375e-01, |x|=4.5000e+00, |dx|/|x|=9.6389e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.537299999933907e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03477341299999992 s.\n",
      "  0 SNES Function norm 5.111423952799e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0392484080000024 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036467766000001234 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035969894000000835 s.\n",
      "  1 SNES Function norm 1.066092159707e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03308899299999979 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034421173000001914 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03733613200000008 s.\n",
      "  2 SNES Function norm 1.621278701286e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04231733199999965 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03468828199999763 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036224170999997085 s.\n",
      "  3 SNES Function norm 1.440731753431e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.396992 s in time step #10  coupling iteration #7\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=3.1116e-01, |x|=1.0000e+00, |dx|/|x|=3.1116e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.105437 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00129971 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 310 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.578726 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.8574e+01, |x|=1.7993e+02, |dx|/|x|=1.0323e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.685689 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.124595 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00114249 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0174709 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=2.3188e-14, |x|=1.7993e+02, |dx|/|x|=1.2887e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.143417 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.00836121\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.832811 s in time step #10  coupling iteration #8\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.0866e-01, |x|=4.5000e+00, |dx|/|x|=1.3526e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 5.604100000056178e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03533022300000255 s.\n",
      "  0 SNES Function norm 6.351235346874e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04746567100000121 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03346419500000053 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.041415321000002336 s.\n",
      "  1 SNES Function norm 2.278328192572e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03587765199999993 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03615852400000108 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0333642679999997 s.\n",
      "  2 SNES Function norm 5.187405436210e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03425133300000027 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.040260213999999905 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04906292899999798 s.\n",
      "  3 SNES Function norm 2.959168761059e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.419479 s in time step #10  coupling iteration #8\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.0966e-01, |x|=1.0000e+00, |dx|/|x|=5.0966e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0991618 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00150028 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 316 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.588113 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5217e+01, |x|=1.8567e+02, |dx|/|x|=8.1958e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.689004 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.116926 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.000933392 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 2 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0256013 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.4832e-14, |x|=1.8567e+02, |dx|/|x|=2.9532e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.143675 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0102152\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.838709 s in time step #10  coupling iteration #9\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.3794e-01, |x|=4.5000e+00, |dx|/|x|=1.4176e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 6.174099999967098e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035656022999997816 s.\n",
      "  0 SNES Function norm 5.647617984353e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03387298099999825 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03985082499999848 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04193467600000034 s.\n",
      "  1 SNES Function norm 2.019742428902e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03326628599999992 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03296799499999992 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.036000135000001876 s.\n",
      "  2 SNES Function norm 2.728850930470e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0344920860000002 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.042937925000000376 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0364770039999982 s.\n",
      "  3 SNES Function norm 2.150587532717e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 3 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.399601 s in time step #10  coupling iteration #9\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.6366e-01, |x|=1.0000e+00, |dx|/|x|=6.6366e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.106525 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00119118 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 335 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.625723 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.1885e+01, |x|=1.8864e+02, |dx|/|x|=6.3003e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.73361 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.120842 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00354079 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 2 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0208243 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.5352e-14, |x|=1.8864e+02, |dx|/|x|=2.4042e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.145411 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0109003\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.884586 s in time step #10  coupling iteration #10\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.6223e-01, |x|=4.5000e+00, |dx|/|x|=1.4716e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.4951999999653935e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04053638799999959 s.\n",
      "  0 SNES Function norm 5.122143931150e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034813942999999625 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033472198999998426 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.045856889999999595 s.\n",
      "  1 SNES Function norm 1.829678023264e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03393854300000143 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04369161700000035 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03698821400000085 s.\n",
      "  2 SNES Function norm 1.012002939742e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.289562 s in time step #10  coupling iteration #10\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=7.0474e-01, |x|=1.0000e+00, |dx|/|x|=7.0474e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.11721 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00139933 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 360 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.679296 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=9.9469e+00, |x|=1.9032e+02, |dx|/|x|=5.2265e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.798124 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.114628 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00145253 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.017879 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=3.2199e-14, |x|=1.9032e+02, |dx|/|x|=1.6919e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.134165 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0111392\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.935612 s in time step #10  coupling iteration #11\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=7.1823e-01, |x|=4.5000e+00, |dx|/|x|=1.5961e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 2.702100000107066e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04268573199999892 s.\n",
      "  0 SNES Function norm 4.875094759968e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034995283000000654 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035407087999999476 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03574232700000124 s.\n",
      "  1 SNES Function norm 1.032400444127e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033340372999997925 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04585725100000104 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03940144700000303 s.\n",
      "  2 SNES Function norm 6.829950449544e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.287618 s in time step #10  coupling iteration #11\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.0715e-01, |x|=1.0000e+00, |dx|/|x|=6.0715e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.111164 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00158942 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 401 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.74535 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.8195e+00, |x|=1.9102e+02, |dx|/|x|=3.0465e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.858308 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.125811 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00399542 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 2 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0212432 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.9224e-14, |x|=1.9102e+02, |dx|/|x|=3.1003e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.151283 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0112205\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 1.01354 s in time step #10  coupling iteration #12\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.6554e-01, |x|=4.5000e+00, |dx|/|x|=1.4790e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.978599999854282e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.045100658000002625 s.\n",
      "  0 SNES Function norm 4.042491443363e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033852677999998804 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03358423999999971 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03674211600000277 s.\n",
      "  1 SNES Function norm 8.572565694816e-05 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0331664219999972 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03317765499999936 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04892505599999808 s.\n",
      "  2 SNES Function norm 3.359416704661e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.284046 s in time step #10  coupling iteration #12\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.1530e-01, |x|=1.0000e+00, |dx|/|x|=4.1530e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.101057 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00210486 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 400 iterations (relative convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.770149 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=8.8325e+00, |x|=1.9166e+02, |dx|/|x|=4.6084e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.87353 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.117105 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00142516 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0146126 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=6.1869e-14, |x|=1.9166e+02, |dx|/|x|=3.2280e-16\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.133339 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0112392\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 1.0103 s in time step #10  coupling iteration #13\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.1607e+00, |x|=4.5000e+00, |dx|/|x|=2.5793e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 4.409299999963423e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03720593699999952 s.\n",
      "  0 SNES Function norm 5.347348474189e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03424898700000156 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04750657700000005 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03499568100000161 s.\n",
      "  1 SNES Function norm 1.209803641258e-04 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03864321000000004 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033783989999999875 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.0328336070000006 s.\n",
      "  2 SNES Function norm 3.590400715002e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 2 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.278619 s in time step #10  coupling iteration #13\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.1972e-01, |x|=1.0000e+00, |dx|/|x|=4.1972e-01\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.103604 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00555007 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 394 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.740326 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=7.9205e-01, |x|=1.9166e+02, |dx|/|x|=4.1327e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.84981 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.121981 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00128159 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0165007 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5315e-14, |x|=1.9166e+02, |dx|/|x|=7.9909e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.139967 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0112187\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.993415 s in time step #10  coupling iteration #14\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.4856e-01, |x|=4.5000e+00, |dx|/|x|=3.3013e-02\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 5.8651999999881355e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03498630500000033 s.\n",
      "  0 SNES Function norm 4.740432411912e-06 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.037819906999999375 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.03743865299999882 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035523156 s.\n",
      "  1 SNES Function norm 7.956800461649e-19 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.155333 s in time step #10  coupling iteration #14\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.0141e-03, |x|=1.0000e+00, |dx|/|x|=5.0141e-03\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.124852 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00205379 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 371 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.695704 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=5.1991e-03, |x|=1.9166e+02, |dx|/|x|=2.7127e-05\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #1 took 0.822823 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.126049 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 0.00130536 s.\n",
      "\n",
      "================================================\n",
      "Linear solver cg with bjacobi preconditioner\n",
      "converged in 1 iterations (absolute convergence criterion fulfilled).\n",
      "================================================\n",
      "[0] \u001b[32minfo:\u001b[m [time] Linear solver took 0.0156027 s.\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=1.5630e-14, |x|=1.9166e+02, |dx|/|x|=8.1555e-17\n",
      "[0] \u001b[32minfo:\u001b[m [time] Iteration #2 took 0.143165 s.\n",
      "[0] \u001b[32minfo:\u001b[m Integral of crack: 0.0112189\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #0 took 0.969834 s in time step #10  coupling iteration #15\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #0 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=9.2502e-04, |x|=4.5000e+00, |dx|/|x|=2.0556e-04\n",
      "[0] \u001b[32minfo:\u001b[m [time] Applying Dirichlet BCs took 3.0286000001211733e-05 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.04042054199999967 s.\n",
      "  0 SNES Function norm 1.924088341236e-08 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.034373963000000174 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.033754085000001766 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Assembly took 0.035764857000000205 s.\n",
      "  1 SNES Function norm 2.284427915539e-18 \n",
      "[0] \u001b[32minfo:\u001b[m PETSsSNES convergence reason 2.\n",
      "[0] \u001b[32minfo:\u001b[m PETScSNES used 1 iterations.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Solving process #1 took 0.15402 s in time step #10  coupling iteration #15\n",
      "[0] \u001b[32minfo:\u001b[m ------- Checking convergence criterion for coupled solution of process #1 -------\n",
      "[0] \u001b[32minfo:\u001b[m Convergence criterion: |dx|=4.3986e-07, |x|=1.0000e+00, |dx|/|x|=4.3986e-07\n",
      "[0] \u001b[32minfo:\u001b[m Elastic energy: 9.8768e-06 Surface energy: 0.00343721 Pressure work: 0 \n",
      "[0] \u001b[32minfo:\u001b[m [time] Time step #10 took 18.9636 s.\n",
      "[0] \u001b[32minfo:\u001b[m [time] Output of timestep 10 took 0.0602107 s.\n",
      "[0] \u001b[32minfo:\u001b[m The whole computation of the time stepping took 10 steps, in which\n",
      "\t the accepted steps are 10, and the rejected steps are 0.\n",
      "\n",
      "[0] \u001b[32minfo:\u001b[m [time] Execution took 31.8034 s.\n",
      "[0] \u001b[32minfo:\u001b[m OGS terminated on 2021-11-16 00:59:50+0100.\n"
     ]
    }
   ],
   "source": [
    "#for running ogs in the jupyter container\n",
    "#!ogs AT1_iso_tensile.prj"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "442ab8bb",
   "metadata": {
    "tags": []
   },
   "source": [
    "**Extracting results using [VTUInterface](https://github.com/joergbuchwald/VTUinterface)**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "ddae910b-4fc6-4384-bebf-723113924613",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stderr",
     "output_type": "stream",
     "text": [
      "ERROR:root:Cannot find UnstructuredGrid element in file.\n",
      "ERROR:root:Algorithm vtkXMLUnstructuredGridReader(0x33d6910) returned failure for request: vtkInformation (0x33d5c60)\n"
     ]
    },
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "=====================\n",
      ">>> print results <<<\n",
      "./AT1_iso_tension.pvd\n"
     ]
    },
    {
     "ename": "AttributeError",
     "evalue": "'NoneType' object has no attribute 'GetData'",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mAttributeError\u001b[0m                            Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-6-fa713380b881>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m\u001b[0m\n\u001b[1;32m      9\u001b[0m \u001b[0mtime\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0;34m[\u001b[0m\u001b[0;36m0.0\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     10\u001b[0m \u001b[0;32mfor\u001b[0m \u001b[0mt\u001b[0m \u001b[0;32min\u001b[0m \u001b[0mtime\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 11\u001b[0;31m     \u001b[0mfunction_xaxis_t\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mpvdfile\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mread_point_set_data\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0;34m'phasefield'\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mpointsetarray\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mxaxis\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     12\u001b[0m     \u001b[0mplt\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mplot\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mr_x\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mfunction_xaxis_t\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mlabel\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m't='\u001b[0m\u001b[0;34m+\u001b[0m\u001b[0mstr\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mt\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     13\u001b[0m \u001b[0;31m# plot formatting\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/.local/lib/python3.8/site-packages/vtuIO.py\u001b[0m in \u001b[0;36mread_point_set_data\u001b[0;34m(self, timestep, fieldname, pointsetarray, interpolation_method)\u001b[0m\n\u001b[1;32m    506\u001b[0m                 \u001b[0mfilename\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mvtufilenames\u001b[0m\u001b[0;34m[\u001b[0m\u001b[0mi\u001b[0m\u001b[0;34m]\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    507\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0mfilename\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0;32mNone\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 508\u001b[0;31m             vtu = VTUIO(os.path.join(self.folder,filename),\n\u001b[0m\u001b[1;32m    509\u001b[0m                     \u001b[0mnneighbors\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnneighbors\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mdim\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdim\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    510\u001b[0m                     \u001b[0mone_d_axis\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mone_d_axis\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m~/.local/lib/python3.8/site-packages/vtuIO.py\u001b[0m in \u001b[0;36m__init__\u001b[0;34m(self, filename, nneighbors, dim, one_d_axis, two_d_planenormal)\u001b[0m\n\u001b[1;32m     49\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpdata\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mGetPointData\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     50\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcdata\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mGetCellData\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m---> 51\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mpoints\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mvtk_to_numpy\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0moutput\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mGetPoints\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mGetData\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     52\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mdim\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mdim\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     53\u001b[0m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mnneighbors\u001b[0m \u001b[0;34m=\u001b[0m \u001b[0mnneighbors\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mAttributeError\u001b[0m: 'NoneType' object has no attribute 'GetData'"
     ]
    }
   ],
   "source": [
    "# read and process (point interpolation) vtu- and pvd-files \n",
    "import vtuIO\n",
    "import matplotlib.tri as tri\n",
    "print(\"=====================\")\n",
    "print(\">>> print results <<<\")\n",
    "pvdfile=vtuIO.PVDIO(\"./AT1_iso_tension.pvd\", dim=3)\n",
    "xaxis =  [(i,0,0) for i in np.linspace(start=0.0, stop=1.0, num=100)]\n",
    "r_x = np.array(xaxis)[:,0]\n",
    "time = [0.0]\n",
    "for t in time:\n",
    "    function_xaxis_t = pvdfile.read_point_set_data(t, 'phasefield', pointsetarray=xaxis)\n",
    "    plt.plot(r_x, function_xaxis_t, label='t='+str(t))\n",
    "# plot formatting\n",
    "titlestring=\"TH process: FEBEX type repository\"\n",
    "plt.title(titlestring)\n",
    "#plt.xscale('log')\n",
    "plt.xlabel('x')\n",
    "plt.ylabel('Temperature')\n",
    "plt.legend()\n",
    "plt.grid()\n",
    "plt.savefig(\"profile.png\")\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "bd15424c-f926-4b11-8f74-b1d591b3571b",
   "metadata": {},
   "source": [
    "**Vertical cross-sections**"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "d211f74c-179c-4838-a34b-07a285b8477e",
   "metadata": {
    "scrolled": true,
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAApMAAAJDCAYAAABXKJ/SAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjIsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8li6FKAAAgAElEQVR4nOzdfbTlVX3n+ffHKotChicBW0KBRQLGFhpLoRHHNkFQqFYbNOIMzhDoFWchjvbYnbY1jOlGRXoF84DtijFNlIiabmRpUBofCFErJhkFQQQhYFMQ0BIMwSqeEnko/M4fZ9/kcLm36t5zz/N5v9Y6656zf3vv3y6S7HzO3r/f76SqkCRJknrxtFEPQJIkSZPLMClJkqSeGSYlSZLUM8OkJEmSemaYlCRJUs8Mk5IkSeqZYVLS0CVZm+SaJDckuTnJe1v5e5L8MMl32utVrfyVSa5L8t3297jR/gskzbokq5Jcn+SK9vnTXXPXnUm+01X37CSbk3wvyYld5Ue2eW1zkg8lSSvfpfW3OcnVSdZ3tTkjyW3tdUZX+cGt7m2t7Zph/HcAiM+ZlDRsbcLcraoeTvJ04C+AtwMbgYer6rfm1X8h8DdVdXeSw4Erq+qAoQ9ckpokvwocBexRVa+Zd+y3gQeq6n1Jng/8d+Bo4GeAPwWeW1VPJLmGztz3TeCLwIeq6ktJ/m/giKo6K8mpwOuq6n9P8kzg2nbeAq4DjqyqbUkuBf64qi5J8vvADVX1kSH8p3BlUtLwVcfD7ePT22vRb7ZVdX1V3d0+3gysTbLLgIcpSQtKsg54NfDRBY4F+N/oBEiAk4FLqurRqvprYDNwdJL96QTRb1RnZe8TwGu72lzc3n8GOL71eyJwVVVtraptwFXAxnbsuFaX1naur4EzTEoaibZF9B3gXjqT49Xt0NuS3JjkoiR7L9D09cD1VfXo0AYrSU/2QeCdwE8XOPYyOjspt7XPBwA/6Dq+pZUd0N7PL39Sm6raDjwA7LODvvYB7m915/c1cKuHdaKlWPO0XWvXVbuPehjS2Htw+9/eV1X7LafNy45dW9u2LjTv9d/N3338ZuCRrqILq+rC7jpV9QSwIclewGVt+/ojwLl0VinPBX4b+JW5NkkOA84HThjsv6A3zmHS0ix3Dhun+SvJa4B7q+q6JMcu0MUb+cdVSYAsUKd2UN5Lmx31NXBjFSZ3XbU7/+u+bxj1MKSx9+Uf/d5dy22zbetP+ewX9h3EcJ7ieQfd80hVHbWUulV1f5JNwMbuayWT/AFwRdfndcBlwOlVdXufh9wXzmHS0ix3Dhuz+eulwEntBsG1wB5JPlVVpyVZDfwScGRX/S3AgV2f1wF3t/J1C5R3t9nS+twT2NrKj53XZhNwH7BXktVtdbK7r4Fzm1vS0CXZr61IkmRX4BXAre0aojmvA25qdfYCvgCcXVV/OezxStKcqjq7qtZV1XrgVOCrVXVaO/wK4Naq6t6+vhw4td2hfTBwKHBNVd0DPJTkmHbN4+nA57vazN2pfUo7RwFXAick2btdBnQCnRsSC/haq0trO9fXwI3VyqSkmbE/cHGSVXS+1F5aVVck+WSSDXS2Z+4E3tzqvw04BPiPSf5jKzuhqu4d8rglaUdO5clb3FTVze1O678CtgNvbZf5ALwF+DiwK/Cl9gL4GPDJJJvprEie2vramuRc4Fut3vuqamt7/y7gkiTvB65vfQyFYVLS0FXVjcALFyj/5UXqvx94/6DHJUnLUVWb6Gwzz33+14vUOw84b4Hya4HDFyh/BFjwmpmqugi4aIHyO+g8fmjo3OaWJElSzwyTkiRJ6plhUpIkST0zTEqSJKlnhklJkiT1zDApSZKknhkmJUmS1DPDpCRJknpmmJQkSVLPDJOSJEnqmWFSkiRJPTNMSpIkqWeGSUmSJPXMMClJkqSeGSYlSZLUs9WjHkC3WruGx563bsn119y6ZYCjkSRJ0s6MVZhcruUEz0Ey1EqSpFk10WFyXAwj1BpYJUnSODJMTojlBlbDpyRJGgbD5JRytVSSJA2DYVI982YpSZJkmNRQLCV4GjglSZo8hkmNjX5tzRtKJUkaHsOkpo43K0mSNDyGSc28UT6v1CArSZp0hklphHoOsj/q7zgkSeqVv80tSZKknhkmJUmS1DPDpCRJknpmmJQkSVLP+hImk1yU5N4kN3WVvSfJD5N8p71e1Y9zSVI/OX9J0sr0a2Xy48DGBcovqKoN7fXFPp1Lkvrp4zh/SVLP+hImq+rrwNZ+9CVJw+T8JUkrM+hrJt+W5Ma2jbT3gM8lSf3k/CVJSzDIMPkR4OeADcA9wG8vVCnJmUmuTXLtY4//3QCHI0lLtqT5C+bNYT/9ybDGJ0ljY2Bhsqr+pqqeqKqfAn8AHL1IvQur6qiqOmrN03cb1HAkacmWOn+1uv84hz1t1+ENUpLGxMDCZJL9uz6+DrhpsbqSNE6cvyRp6fry29xJ/jtwLLBvki3AOcCxSTYABdwJvLkf55KkfnL+kqSV6UuYrKo3LlD8sX70LUmD5PwlSSvjL+BIkiSpZ4ZJSZIk9cwwKUmSpJ4ZJiVJktQzw6QkSZJ6ZpiUJElSzwyTkiRJ6plhUpIkST0zTEqSJKlnhklJkiT1zDApSZKknhkmJUmS1DPDpCRJknpmmJQkSVLPDJOSJEnqmWFSkiRJPTNMSpIkqWeGSUmSJPXMMClJkqSeGSYlSZLUM8OkpKFLsjbJNUluSHJzkvfOO/6OJJVk33nlByV5OMk7hjtiSXqyJKuSXJ/kiq6yf5Pke21e+0ArW5/kJ0m+016/31X/yCTfTbI5yYeSpJXvkuTTrfzqJOu72pyR5Lb2OqOr/OBW97bWds0w/juAYVLSaDwKHFdVLwA2ABuTHAOQ5EDglcD3F2h3AfCloY1Skhb3duCWuQ9JXg6cDBxRVYcBv9VV9/aq2tBeZ3WVfwQ4Ezi0vTa28jcB26rqEDrz3vntHM8EzgFeDBwNnJNk79bmfOCCqjoU2Nb6GArDpKShq46H28ent1e1zxcA7+z6DECS1wJ3ADcPa5yStJAk64BXAx/tKn4L8BtV9ShAVd27kz72B/aoqm9UVQGfAF7bDp8MXNzefwY4vq1anghcVVVbq2obcBWdL+MBjmt1aW3n+hq41cM6kSR1S7IKuA44BPhwVV2d5CTgh1V1Q9vtmau7G/AuOiuWY7vFXWvX8Njz1o16GNL4+9GoB7BiH6TzpXf3rrLnAi9Lch7wCPCOqvpWO3ZwkuuBB4Ffr6o/Bw4AtnS139LKaH9/AFBV25M8AOzTXT6vzT7A/VW1fYG+Bs4wKWkQ9k1ybdfnC6vqwu4KVfUEsCHJXsBlSY4A3g2csEB/76WzffNwd8iUpAHY4fyV5DXAvVV1XZJju+qtBvYGjgH+OXBpkp8F7gEOqqofJzkS+FySw4CFJrO5HZnFji23fCgMk9KM2PbEM/jMgy8a0tm+cF9VHbWUmlV1f5JNdLZ1DgbmViXXAd9OcjSd64NOaRe07wX8NMkjVfW7Axm+pLEyZvPXS4GTkrwKWAvskeRTdFYD/7htWV+T5KfAvlX1t3SuE6cF0NvprGJuoTPPzVkH3N3ebwEOBLYkWQ3sCWxt5cfOa7MJuA/YK8nqtjrZ3dfAec2kpKFLsl9bkSTJrsArgOur6llVtb6q1tOZNF9UVT+qqpd1lX8Q+M8GSUmjUFVnV9W6Nh+dCny1qk4DPkfnukWSPBdYA9zX5rtVrfxn6dxoc0dV3QM8lOSYds3j6cDn22kuB+bu1D6lnaOAK4ETkuzdbrw5AbiyHftaq0trO9fXwLkyKWkU9gcubhPs04BLq+qKnbSRpHF2EXBRkpuAx4AzqqqS/ALwviTbgSeAs6pqa2vzFuDjwK50nlQx97SKjwGfTLKZzorkqQBVtTXJucDctZjv6+rrXcAlSd4PXN/6GArDpKShq6obgRfupM76RcrfM4AhSdKyVdUmOtvMVNVjwGkL1Pks8NlF2l8LHL5A+SPAGxZpcxGd4Dq//A46jwsaOre5JUmS1DPDpCRJknpmmJQkSVLPDJOSJEnqmWFSkiRJPTNMSpIkqWeGSUmSJPXMMClJkqSeGSYlSZLUM8OkJEmSemaYlCRJUs8Mk5IkSeqZYVKSJEk9M0xKkiSpZ4ZJSZIk9cwwKUmSpJ4ZJiVJktQzw6QkSZJ6ZpiUJElSzwyTkiRJ6plhUpIkST0zTEqSJKlnhklJkiT1zDApSZKkno1VmHxibXjg53bhgZ/bZdRDkSRJ0hKsHvUAFtPPQLnn7Y/2rS9JkiT9o7ENk/20WDA1ZEqSJK3MTITJxSwUMg2YkiRJSzfTYXIhrmJKkiQtnWFyiQyZkiRJT2WYXCG3yiVJ0iwzTA5Ad8A0WEqSpGlmmByw+SuXhktJkjRNDJND5qqlJEmaJobJEXLVUpIkTTrD5Bhx1VKSJE0aw+SYMlhKkqRJ8LRRD0A798DP7dLX3yqXJEnqF1cmJ4irlZIkadwYJieUwVKSJI2DvmxzJ7koyb1Jbuoqe2aSq5Lc1v7u3Y9z6anmtsHdCpeWz/lLklamX9dMfhzYOK/s14CvVNWhwFfaZw2YoVJato/j/CVJPetLmKyqrwNb5xWfDFzc3l8MvLYf59LSuFopLY3zlyStzCDv5v4nVXUPQPv7rAGeSztgqJSWzflLkpZo5I8GSnJmkmuTXLv9J3836uFMNUOl1H/dc9hjjzuHSZo9gwyTf5Nkf4D2996FKlXVhVV1VFUdtXrX3QY4HM0xVEo7taT5C548h615unOYpNkzyDB5OXBGe38G8PkBnks9MFRKi3L+kqQl6stzJpP8d+BYYN8kW4BzgN8ALk3yJuD7wBv6cS71n8+s1Cxz/pKklelLmKyqNy5y6Ph+9K/hmQuWhkrNCucvSVoZfwFHC3K1UpIkLcXI7+bW+PPaSkmStBjDpJbMUClJkuYzTGrZDJWSJGmO10yqZ15XKUmSXJlUX7haKUnSbDJMqq8MlZIkzRbDpAbCUClJ0mwwTGqgDJWSJE03w6SGwlApSdJ0MkxqqAyVAkiyNsk1SW5IcnOS97byc5PcmOQ7Sf4kyc90tTkiyTda/e8mWTu6f4GkWZdkVZLrk1zRVfZvknyvzVMf6Co/O8nmduzErvIj23y2OcmHkqSV75Lk06386iTru9qckeS29jqjq/zgVve21nbNoP8bzDFMaiTmQqXBcmY9ChxXVS8ANgAbkxwD/GZVHVFVG4ArgP8EkGQ18CngrKo6DDgWeHwkI5ekjrcDt8x9SPJy4GTgiDZP/VYrfz5wKnAYsBH4vSSrWrOPAGcCh7bXxlb+JmBbVR0CXACc3/p6JnAO8GLgaOCcJHu3NucDF1TVocC21sdQGCY1cobK2VMdD7ePT2+vqqoHu6rtBlR7fwJwY1Xd0Nr/uKqeGNqAJalLknXAq4GPdhW/BfiNqnoUoKrubeUnA5dU1aNV9dfAZuDoJPsDe1TVN6qqgE8Ar+1qc3F7/xng+LZqeSJwVVVtraptwFV0vowHOK7VpbWd62vgfGi5xsZcoPQB6LOhfTO/DjgE+HBVXd3KzwNOBx4AXt6qPxeoJFcC+9GZmD/w1F5H64m18YuRtBSbRj2AFfsg8E5g966y5wIva3PYI8A7qupbwAHAN7vqbWllj7f388tpf38AUFXbkzwA7NNdPq/NPsD9VbV9gb4GzpVJjR1XKqfCvkmu7XqdOb9CVT3RtrPX0fmWfngrf3dVHQj8EfC2Vn018C+A/7P9fV2S44fyL5E0a3Y4fyV5DXBvVV03r91qYG/gGOA/AJe2FcMscI7aQTk9tNlRXwPnyqTGliuV/fXg9rVc9aPnDelsX7ivqo5aSs2quj/JJjrXCt3Udei/AV+gc33QFuDPquo+gCRfBF4EfKWfo5Y0nsZs/nopcFKSVwFrgT2SfIrOPPXHbcv6miQ/BfZt5Qd2tV8H3N3K1y1QTlebLe2a8T2Bra382HltNgH3AXslWd1WJ7v7GjhXJjX2XKmcPkn2S7JXe78r8Arg1iSHdlU7Cbi1vb8SOCLJM9rE+ovAXw1zzJIEUFVnV9W6qlpP58aar1bVacDn6Fy3SJLnAmvohLzLgVPbHdoH07nR5pqqugd4KMkxbQXzdODz7TSXA3N3ap/SzlF05sITkuzdbrw5AbiyHftaq0trO9fXwLkyqYnhSuVU2R+4uF03+TTg0qq6Islnk/w88FPgLuAsgKraluR3gG/R2br5YlV9YURjl6SFXARclOQm4DHgjBbybk5yKZ0vwNuBt3bdQPgW4OPArsCX2gvgY8Ank2ymsyJ5KkBVbU1yLp25EOB9VbW1vX8XcEmS9wPXtz6GwjCpiWOonHxVdSPwwgXKX7+DNp+i83ggSRoLVbWJdjtRVT0GnLZIvfOA8xYovxY4fIHyR4A3LNLXRXSC6/zyO+g8Lmjo3ObWxHL7W5Kk0TNMauIZKiVJGh3DpKaGgVKSpOHzmklNle5A6TWVkiQNnmFSU8tgKUnS4BkmNRMMlpIkDYZhUjPHYClJUv8YJjXTDJaSJK2Md3NLjY8YkiRp+VyZlOZxtVKSpKVzZVLaAVcrJUnaMcOktASGSkmSFmaYlJbBUClJ0pN5zaTUA6+rlCSpwzAprdDOVioNm5KkaWaYlAZsKdviBk5J0qQyTEpjwNVNSdKkGqsw+cQaeOg5WVab3e+qAY1GGh9PCZubRjIMSZKeYqzCZC8WCp8GTEmSpOGY+DC5EAOmJEnScExlmFxId8A0WEqSJPXHzITJbvNXLg2XkiRJvZnJMDmfq5aSJEm98ecU53noOVn2HeWSJEmzypXJRbhaKUmStHOuTC6Bq5WSJEkLM0wug6FSkiTpyQyTPTBUSpIkdRgmV8BQKUmSZp1hsg8MlZIkaVYZJvvIQClJkmaNYbLPXKWUJEmzxDA5IAZKSZI0CwyTA+QqpSRJmnaGySEwVEqSpGllmBwiA6UkSZo2hskhc5VSkiRNE8PkiBgoJUnSNDBMjpCBUpIkTTrD5Ii57S1JkiaZYXJMGCglSdIkMkyOEVcpJUnSpDFMjiEDpSRJmhSGyTFloJQkSZPAMDnGDJSSJGncGSbHnIFSkiSNM8PkBDBQSpKkcWWYnBAGSkmSNI4MkxPEQClJksbN6kGfIMmdwEPAE8D2qjpq0OecZg89J+x+V416GNJMcP6SpJ0beJhsXl5V9w3pXFPPQCkNlfOXJO2A29wTyi1vSZI0DoYRJgv4kyTXJTlzCOebGQZKaeCcvyRpJ4axzf3Sqro7ybOAq5LcWlVfnzvYJugzAVbvufcQhjNd3PKWBmqH8xc8eQ5bs5tzmKTZM/CVyaq6u/29F7gMOHre8Qur6qiqOmrVbrsNejhTyRVKaTB2Nn+1Y/8wh63e1TlM0uwZaJhMsluS3efeAycANw3ynLPKQCn1l/OXJC3NoLe5/wlwWZK5c/23qvrygM85s9zylvrK+UuSlmCgYbKq7gBeMMhz6MkMlFJ/OH9J0tL4aKAp5Ja3JEkalrEKk7WmePSgxxZ9aekMlBpnSdYmuSbJDUluTvLeVv6bSW5NcmOSy5Ls1cqfnuTiJN9NckuSs0f7L5A065KsSnJ9kiva5/ck+WGS77TXq1r5+iQ/6Sr//a4+jmzz2uYkH0q7ribJLkk+3cqvTrK+q80ZSW5rrzO6yg9udW9rbdcM67/FWIXJnTFgLo+BUmPsUeC4qnoBsAHYmOQY4Crg8Ko6AvifwFxofAOwS1X9M+BI4M3dk6skjcDbgVvmlV1QVRva64td5bd3lZ/VVf4ROo8WO7S9NrbyNwHbquoQ4ALgfIAkzwTOAV5M5+kS5ySZeybZ+e38hwLbWh9DMVFhciEGS2nyVMfD7ePT26uq6k+qansr/yawbq4JsFuS1cCuwGPAg8McsyTNSbIOeDXw0RX0sT+wR1V9o6oK+ATw2nb4ZODi9v4zwPFt1fJE4Kqq2lpV2+h8Ad/Yjh3X6tLazvU1cMP6be6h6A6Uu3x/aKu7Y80bcjSukqwCrgMOAT5cVVfPq/IrwKfb+8/QmVzvAZ4B/Luq2jqssQ6LuwnSxPgg8E5g93nlb0tyOnAt8O9b4AM4OMn1dL4E/3pV/TlwALClq+2WVkb7+wOAqtqe5AFgn+7yeW32Ae7v+jLe3dfATVWY7Gaw/EcGSo3Avkmu7fp8YVVd2F2hqp4ANrTrIi9LcnhV3QSQ5N3AduCPWvWjgSeAnwH2Bv48yZ+2O67HxhNrDITSFNjh/JXkNcC9VXVdkmO76n0EOJfOTsq5wG/T+VJ8D3BQVf04yZHA55IcBiw0Wcz9P+vFji23fCimNkx2mwuWsxwqDZR67LHV3Lllv2Gd7r6qOmopFavq/iSb6FwrdFO7oPw1wPFt6wfg/wC+XFWPA/cm+UvgKGCswqSkwRiz+eulwEntBpu1wB5JPlVVp81VSPIHwBUAVfUonevEaQH0duC5dFYP13X1uw64u73fAhwIbGmX9+wJbG3lx85rswm4D9gryeq2Otnd18BN/DWTy+F1ldJ4SLJf153auwKvAG5NshF4F3BSVf19V5PvA8elYzfgGODWYY9bkqrq7KpaV1XrgVOBr1bVae0ayDmvo/1iVpvvVrX3P0vnRps7quoe4KEkx7RrHk8HPt/aXw7M3al9SjtHAVcCJyTZu914cwJwZTv2tVaX1naur4GbiZXJbrO8SunqpMbI/sDFbYJ9GnBpVV2RZDOwC3BVe0LGN9udjx8G/pDO5BzgD6vqxtEMXZIW9IEkG+hsL98JvLmV/wLwviTb6Vyuc1bXNd9vAT5O58bCL7UXwMeAT7Y5cSud0EpVbU1yLvCtVu99XX29C7gkyfuB61sfQzFzYXLOowc9ZqCURqQFwRcuUH7IIvUfpvN4IEkaG1W1ic42M1X1y4vU+Szw2UWOXQscvkD5Iywy51XVRcBFC5TfQef68qGbqW3u+Wb1kULeICBJkvplpsPkHAOlJElSbwyTzSyuUhooJUnSShkm5zFQSpIkLZ1hcgGzFiglSZJ6ZZhcxCwFSlcnJUlSrwyTOzBL11EaKCVJUi8Mk0swK4FSkiRpuQyTSzQLgdLVSUmStFyGyWWYhUApSZK0HIbJZZr2QOnqpCRJWg7DpJ7CQClJkpbKMNmDaV+dlCRJWirDZI+mPVC6OilJkpbCMLkC0x4oJUmSdsYwuULTHChdnZQkSTtjmNQOGSglSdKOGCb7YJpXJ8FAKUmSFmeY7JNpD5SSJEkLWT3qAXRbs2Y769f97YLH7tyy35BHs3yPHvQYu3x/zaiHMRAPPSfsfleNehiSJGnMjFWY3JHukDkJwXIaGSglSdJ8ExMmu41zsJzm1UlJkqT5Jv6aycW2xUdpmq+f9GYcSZLUbeLDJHQC5TiGymlloJQkSXOmIkzOGadQOc2rk5IkSXOmKkzOGZdAOc1cnZQkSTClYRLGI1C6OilJkqbd1IZJGI9AKUmSNM2mOkzC6APlNK9OutUtSZKmPkzC6AOlJEnStJqJMAmjDZSuTkqSpGk1M2ESXKGUJEnqt5kKk2CgHARXJyVJml0zFyZHZZq3uiVJ0uyayTDp6mT/uTopSdJsmskwCaMJlK5OSpKkaTNWYXKP1Y/wymffyiuffeuoh6IeuDopSdLsWT3qASxmfqC86kfP6/s51q/7W+7csl/f+92RRw96jF2+v2ao55QkSRqUsVqZ3JFBrVZ6/WR/uTopSdJsmZgwCUzNFrjXTkqSpGkxUWFyTr9DpauT/eXqpCRJs2Miw+ScaVillCRJmmQTHSahf4Fy2KuT077V7eqkJEnjKckzl/Daa6n9je3d3MvxymffOpC7vSVJkqbQ3e21o5WfVcBBS+lsKsKkxtNDzwm731WjHoYkSXqyW6rqhTuqkOT6pXY28dvcc/qx3e1WtyRJmgEv6VMdYMpWJt3uliRJ2rGqemTufZK9gUOBtV3Hv95dZ2emZmWyX3xMUH95I44kSeMpyf8FfB24Enhv+/ue5fYzVmFy71V/zyl7fJtT9vh2z31M2uOC3OqWJEkj8nbgnwN3VdXLgRcCy15VG9tt7u5A+ZkHXzTCkUiaJEmeuYRqP62q+wc+GEkab49U1SNJSLJLVd2a5OeX28nYhslup+zx7WUFypVeO7l+3d9y55b9em6vJ/Oubg1ZXx95IUlTbEt7nuTngKuSbKMzfy7LRIRJWH6glDSz+vrIC0maVlX1uvb2PUm+BuwJfHm5/YzVNZM7s5JrKcfZLFw36Y04GqK+PvJCkqZNkrVJ/m2S303y5iSrq+rPquryqlp2KJmoMLkck3YjjqS+2XexA0n+FTz5sRiSNIMuBo4Cvgv8S+C3V9LZxIXJYa1O+oggaWJ9Jcn6+YVJfgX44NBHI0nj5/lVdVpV/VfgFOBlK+ls4sKkJpdb3RqSf0fnQvJD5wqSnN3Kf3Fko5Kk8fH43Juq2r7SziYyTC51ddKtbmn2VNUXgbOALyU5PMkHgdcAv1BVW0Y7OkkaCy9I8mCSh5I8BBwx9z7Jg8vtbGLu5p52jx70GLt8f82ohyFNhar6SpJ/DWwC/j/geK+TlKSOqlrVz/4MkxoqnzmpQWvfsovOcyZ3AY4H7k0SoKpqj1GOT5JGLcnlOzpeVSctp7+Bb3Mn2Zjke0k2J/m1fvU7rY8JkmZBeyzFNUluSHJzkve28t9McmuSG5Nc1h6mO9fm7DaPfC/JiYv1XVW7V9Ue7e+aqtqt6/OyguSg5i9Jky/JqiTXJ7mifX5Pkh8m+U57vaqr7oLzV5Ijk3y3HftQ+9JLkl2SfLqVX919U2GSM5Lc1l5ndJUf3Ore1truaLvzJcCBwJ8Dv0Xnbu7u17IMNEwmWQV8mM5t588H3pjk+YM8Zz95R/dgeCOOgEeB46rqBcAGYGOSY4CrgMOr6gjgfwJnA7R541TgMGAj8HttfhmYSZ+/JA3c24Fb5pVdUFUb2uuLsNP56yPAmcCh7bWxlb8J2FZVhwAXAOe3vp4JnAO8GDgaOCfJ3q3N+e38hwLbWh+LeTad+fVw4L8ArwTua8+a/LPl/ocY9Mrk0cDmqrqjPQTzEuDkAZ9T0pirjofbx6e3V1XVn0TEQ2UAACAASURBVHTdWfhNYF17fzJwSVU9WlV/DWymM788RZKdblsspQ7OX5IWkWQd8Grgo0uovuD8lWR/YI+q+kZVFfAJ4LVdbS5u7z8DHN9WLU8ErqqqrVW1jc4X8I3t2HGtLq3tXF9PUVVPVNWXq+oM4Jg2pk1J/s1S/xt0G/Q1kwcAP+j6vIVOmtYCZukmHK+dVPtmfh1wCPDhqrp6XpVfAT7d3h9AJ1zO2dLKFvJPk9y4o1PT+cmwnVn2/FVraiZ+0UoSHwTeCew+r/xtSU4HrgX+fQt8i81fj7f388uha/6pqu1JHgD2YeF56YB27P6uL+M7miOBzlY6nUD8RmA98CHgj3fUZjGDDpML7Wc+KUEkOZPOEi8/c0D/d61e+exbuepHz+t7v9KkyWMZ5peVfZNc2/X5wqq6sLtCVT0BbGjXRV6W5PCqugkgybuB7cAfteo7nUu6LOX/4J9YQp0lnbN7Dlu1z15PaSBp5cZp/kryGuDeqrouybFd9T4CnEtnnjiXzrWHv8Lic8mO5pjltlnOHEmSi+lscX8JeO/c3NurQYfJLXQu8JyzDri7u0L7H9CFAIcfsabvS1UGSWkk7quqo5ZSsaruT7KJzrVCN7ULyl9D53E+c3PCTueSrv7u6nnUT7akc3bPYbscvM7ldmny7Wz+eilwUrvBZi2wR5JPVdVpcxWS/AFwRfu42FyyhX+8lKe7vLvNliSr6eymbG3lx85rswm4D9ir/cb2dnYwRza/DPwd8Fzg/2n3/UAnlC77qReDvmbyW8Ch7Q6jNXQuQN3h7eiSpl+S/ebu1E6yK/AK4NYkG4F3ASdV1d93NbkcOLXd4XgwnQvVrxnwMJ2/JD1FVZ1dVeuqaj2deeGrVXVauwZyzuuAudW+BeevqroHeCjJMe2ax9OBz3e1mbtT+5R2jgKuBE5Isne78eYE4Mp27GutLq3tXF8L/Rue1p5wMff0i7nXsp96AQNemWz7/G+j849fBVxUVTcP8pyaHF43OdP2By5u100+Dbi0qq5IspnOsyGvat+Uv1lVZ1XVzUkuBf6Kzvb3W9s2+VMkeUlrt6L/5XL+krRMH0iygc728p3AmwF2Mn+9Bfg4sCudLecvtfKPAZ9sc+JWOqGVqtqa5Fw6X3YB3ldVW9v7dwGXJHk/cH3rYygG/tDydmv8Fwd9HkmTo6puBF64QPkhO2hzHnDeEro/A/hwkv8JfBn4clX9qMdxOn9JWlRVbaKzzUxV/fIO6i04f1XVtXSuXZxf/gjwhkX6ugi4aIHyO1jkKRfzJfl2Vb1opXXm+As4kqZKVZ0FkOR5dJ4R+fEke9LZAvoy8JeLrWpK0ozo11MvgAkOk595cElhWWPOrW4NSlXdCtwKXNCuy3w5nW/6vwMs6eYgSZpS/XrqBTDBYVKSlqqqfkJnu9ota0kzr49PvQCG8NvckiRJml4TGSbd4p4u/la3+inJb416DJI0SyYyTC6VDyyXZtJxox6AJI2zJC9J15PKV2riwqSrkpIkSStyBnBdkkuS/Oskz15JZxN1A45BUtISvCDJXwPfpfMLFHN/b62qx0c6MkkaA/1+hNrErEwuN0j2Y4v7zi37rbiP5Xr0oMeGfs5x4HWT6qMb6fx27u8CP6bzc2N/CNyX5KYdNZSkWVJVt1bVBVW1kc4lQn9B5xFqVy+nn4lYmXRFUtJyVNXdwN3An8yVteuDFv2FHUmaZSt5hNpYh0lDpKQefHihwvZb3bcNeSySNPXGKkxue+IZfQmQ3sU9mfw1HPVDVX101GOQpFkyMddMSpIkqX+S/MskVyf5XpJLk7ykl36mLkxO+qrkLt9fM+ohSJKk2fB7wK8CxwAXAr+Z5I3L7WTqwmS/jOJObnlXtyRJQ/Q3VfWXVbWtqv4UOBF493I7maowOemrkpIkSUN0Z5L3J5nbFn0ceGi5nUxNmDRISpIkLUsBvwT8IMlfAJuBTUkOXU4nY3U3tyRJkoajqt4IkGQtcDjwgvb6aJKfraoDl9LPVITJfq9Ker3kaPmIIEmShqeqHgGuba9lm/htbre3JUmSRmeiw6RBUpIkabQmcpt7kCHSLW5JkqSlm7iVSVcjZ4PPm5QkaTJMzMqkIVKSJGn8jGWYHFVwHPUWtz+lKEmSJs1YbXM/uH2tK5CSJEkTZKzC5CiNelVST+V1k5IkjT/DpCRJknpmmMRVSUmSpF7NfJg0SEqSJPVu5sOkxpvXTUqSNN5mOky6KilJkrQyMxsmDZKSJEkrN7Nhctz4wPLFudUtSdL4mskw6aqkJElSf8xcmDRISpIk9c9MhUmDpCRJUn/NVJgcV14vKUmSJtXMhElXJSVJkvpvJsKkQVKSJGkwVo96AIM0CSHSLW5JkjTJpnZlchKCpCRJ0qSbyjBpkJQkSRqOqdrmnrQQ6Ra3JEmadFOzMjlpQVKSJGkaTPzKpCFSkiRpdCYyTE5DgHSLW5IkTYOJCpPTECIlSZKmyViFycceW21glCRJmiBTcwPOJHGLW5IkTQvDpCRJknpmmJQkSVLPDJND5hZ3bx56TkY9BEmStADDpCRJknpmmBwiVyUlSdK0MUxKkiSpZ4bJIXFVUpIkTSPDpCRJknpmmBwCVyX7wzu6JUkaP4ZJSUOXZG2Sa5LckOTmJO9t5W9on3+a5Kiu+q9Mcl2S77a/x41u9JIESVYluT7JFfPK35GkkuzbPq9P8pMk32mv3++qe2Sb1zYn+VCStPJdkny6lV+dZH1XmzOS3NZeZ3SVH9zq3tbaDm0lyzA5YK5KSgt6FDiuql4AbAA2JjkGuAn4JeDr8+rfB/yrqvpnwBnAJ4c5WElawNuBW7oLkhwIvBL4/ry6t1fVhvY6q6v8I8CZwKHttbGVvwnYVlWHABcA57f+nwmcA7wYOBo4J8nerc35wAVVdSiwrfUxFIZJSUNXHQ+3j09vr6qqW6rqewvUv76q7m4fbwbWJtllSMOVpCdJsg54NfDReYcuAN4J1BL62B/Yo6q+UVUFfAJ4bTt8MnBxe/8Z4Pi2ankicFVVba2qbcBVdL6MBziu1aW1netr4FYP60SzyFVJaXFJVgHXAYcAH66qq5fY9PXA9VX16MAG16M1a7azft3fjnoY0ti7a9QDWLkP0gmNu88VJDkJ+GFV3dB2q7sdnOR64EHg16vqz4EDgC1ddba0MtrfHwBU1fYkDwD7dJfPa7MPcH9VbV+gr4EzTEozYtVjsPtdO/2y3C/7Jrm26/OFVXVhd4WqegLYkGQv4LIkh1fVTTvqNMlhdLZyTuj7iCWNrXGav5K8Bri3qq5Lcmwrewbwbhaem+4BDqqqHyc5Evhcm8sWuqt07h+52LHllg+FYXJAXJUcjIeek2FOKOrdfVV11M6rQVXdn2QTnWuFFg2TbVvpMuD0qrq9L6OUpKfa2fz1UuCkJK8C1gJ70LmO+2BgblVyHfDtJEdX1Y/oXCdOC6C3A8+ls3q4rqvfdcDc5TxbgAOBLUlWA3sCW1v5sfPabKJzXfleSVa31cnuvgbOayYlDV2S/dqKJEl2BV4B3LqD+nsBXwDOrqq/HM4oJempqursqlpXVeuBU4GvVtXrq+pZVbW+lW8BXlRVP2rz3SqAJD9L50abO6rqHuChJMe0ax5PBz7fTnM5nZsNAU5p5yjgSuCEJHu3G29OAK5sx77W6tLazvU1cIbJAXBVUtqp/YGvJbkR+BadC8qvSPK6JFuAlwBfSHJlq/82OtdW/seux2s8azRDl6Rl+QXgxiQ30LlB5qyq2tqOvYXOTTybgduBL7XyjwH7JNkM/CrwawCt3bl05s1vAe/r6utdwK+2Nvu0PobCbe4+M0hKO1dVNwIvXKD8Mjpb2fPL3w+8fwhDk6Qlq6pNdLaZ55ev73r/WeCzi7S/Fjh8gfJHgDcs0uYi4KIFyu+g87igoXNlso8MkpIkadYMLEwmeU+SH3ZtSb1qUOcaBwZJaXrM2vwlSSsx6G3uC6rqtwZ8Ds0Q7+TWEDl/SdISuM3dB65KSpKkWTXoMPm2JDcmuajrtyOnikFSmlpTP39JUj+sKEwm+dMkNy3wOpnOj5f/HLCBztPff3uRPs5Mcm2Sa5946O9WMpyhM0hKk6sf81fr5x/msMfv//shjV6SxseKrpmsqlcspV6SPwCuWKSPC4ELAXY5eN3EXBBnkJQmWz/mr9bPP8xhu//8sydmDpOkfhnk3dz7d318HTv4mTRJGifOX5K0dIO8m/sDSTbQ+aHxO4E3D/BcQ+WqpDT1pnb+kqR+G1iYrKpfHlTfo2SQlKbftM5fkjQIPhpoGQySkiRJT2aYXCKDpCRJ0lMZJpfAIClJkrQww+ROGCQlSZIWZ5jcAYOkJEnSjg3y0UATyxApSZK0NK5MzmOQlCRJWjrDZBeDpCRJ0vIYJhuDpCRJ0vLN/DWThkhJkqTezWyYNERKkiSt3MyFSUOkJElS/8xMmDRESpIk9d/Uh0lD5PTY/a4a9RAkSdI8UxsmDZGSJEmDN3Vh0hApSZI0PBMfJg2PkiRJozP2YdKwKEmSNL7GKkzmsRgeJUmSJog/pyhJkqSeGSYlSZLUM8OkJEmSemaYlCRJUs8Mk5IkSeqZYVKSJEk9M0xqIvi73JIkjSfDpCRJknpmmJQkSVLPDJMae25xS5I0vgyTkiRJ6plhUmPNVUlJksabYVKSJEk9M0xKkiSpZ4ZJjS23uCVJGn+GSUmSJPXMMClJkqSeGSY1ltziliRpMhgmJUmS1DPDpCRJknpmmNTYcYtbkqTJYZiUJElSzwyTkiRJ6plhUmPFLW5JkiaLYVKSJEk9M0xKkiSpZ4ZJjQ23uCVJmjyGSUmSJPXMMClJkqSeGSY1FtziliRpMhkmJQ1dkrVJrklyQ5Kbk7y3lT8zyVVJbmt/957X7qAkDyd5x2hGLkkdSVYluT7JFfPK35GkkuzbVXZ2ks1JvpfkxK7yI5N8tx37UJK08l2SfLqVX51kfVebM9oceVuSM7rKD251b2tt1wzy39/NMClpFB4FjquqFwAbgI1JjgF+DfhKVR0KfKV97nYB8KWhjlSSFvZ24JbugiQHAq8Evt9V9nzgVOAwYCPwe0lWtcMfAc4EDm2vja38TcC2qjqEzrx3fuvrmcA5wIuBo4Fzur50nw9c0ObPba2PoTBMShq66ni4fXx6exVwMnBxK78YeO1cmySvBe4Abh7iUCXpKZKsA14NfHTeoQuAd9KZz+acDFxSVY9W1V8Dm4Gjk+wP7FFV36iqAj7BP8553XPhZ4Dj26rlicBVVbW1qrYBV9H5Mh7guFYX5s2fg7Z6WCeSpG7tm/l1wCHAh6vq6iT/pKruAaiqe5I8q9XdDXgXnW/8Y7vFvcfqR3jls28d9TCksfdnox7Ayn2QTmjcfa4gyUnAD6vqhrZbPecA4Jtdn7e0ssfb+/nlc21+AFBV25M8AOzTXT6vzT7A/VW1fYG+Bs4wqZHz5pvhWPVIseftjw7rdPsmubbr84VVdWF3hap6AtiQZC/gsiSH76C/99LZvnl43iQtaQaM0/yV5DXAvVV1XZJjW9kzgHcDJyzQ30KTVu2gvJc2O+pr4AyTGimD5NS6r6qOWkrFqro/ySY61wr9TZL926rk/sC9rdqLgVOSfADYC/hpkkeq6ncHMXhJM21n89dLgZOSvApYC+wBfBI4GJhblVwHfDvJ0XRWCQ/sar8OuLuVr1ugnK42W5KsBvYEtrbyY+e12QTcB+yVZHVbnezua+C8ZlLS0CXZr61IkmRX4BXArcDlwNzdiWcAnweoqpdV1fqqWk9ne+k/GyQljUJVnV1V69p8dCrw1ap6fVU9q2ue2gK8qKp+RGdeO7XdoX0wnRttrmmX9DyU5Jh2zePptDmPJ8+Fp7RzFHAlcEKSvduNNycAV7ZjX2t1oWv+HAZXJjUyrkrOtP2Bi9t1k08DLq2qK5J8A7g0yZvo3A35hlEOUpJWqqpuTnIp8FfAduCt7TIfgLcAHwd2pfOkirmnVXwM+GSSzXRWJE9tfW1Nci7wrVbvfVW1tb1/F3BJkvcD17c+hsIwKWnoqupG4IULlP8YOH4nbd8zoGFJ0rJU1SY628zzy9fP+3wecN4C9a4FnnK9eFU9wiJfpqvqIuCiBcrvoPO4oKFzm1sj4aqkJEnTwTCpoTNISpI0PQyTkiRJ6plhUkPlqqQkSdPFMClJkqSeGSY1NK5KSpI0fQyTkiRJ6plhUkPhqqQkSdPJMClJkqSeGSY1cK5KSpI0vVYUJpO8IcnNSX6a5Kh5x85OsjnJ95KcuLJhSlL/OYdJ0sqt9Le5bwJ+Cfiv3YVJnk/nR8kPA34G+NMkz+36YXNJGgfOYZK0QitamayqW6rqewscOhm4pKoeraq/BjYzoh8f12i5xa1x5hwmSSs3qGsmDwB+0PV5SyuTpEngHCZJS7TTbe4kfwo8e4FD766qzy/WbIGyBZeokpwJnAmwes+9dzYcTRBXJTUOhjmH7bn/rj2NUZIm2U7DZFW9ood+twAHdn1eB9y9SP8XAhcCrD3gQNOHpL4a5hx2wGF7OYdJmjmD2ua+HDg1yS5JDgYOBa4Z0Lk0hlyV1IRzDpOkJVrpo4Fel2QL8BLgC0muBKiqm4FLgb8Cvgy81bsgJY0b5zBJWrkVPRqoqi4DLlvk2HnAeSvpX5PJVUlNCucwSVo5fwFHkiRJPTNMqq9clZQkabYYJiVJktQzw6T6xlVJSZJmj2FSkiRJPTNMqi9clZQkaTYZJrViBklJkmaXYVKSJEk9M0xqRVyVlCRpthkm1TODpCRJMkyqJwZJSZIEhklJkiStgGFSy+aqpCRJmmOY1LIYJCVJUjfDpJbMIClJkuYzTGpJDJKSJGkhhkntlEFSkiQtxjCpHTJISpKkHTFMalEGSUmStDOGSS3IIClJkpbCMKmnMEhKkqSlMkzqSQySkiRpOQyT+gcGSUmStFyGSQEGSUmS1BvDpAySkiSpZ4bJGWeQlCRJK2GYnGEGSUmStFKGyRllkJQkSf1gmJxBBklJktQvhskZY5CUJEn9ZJicIQZJSZLUb4bJGWGQlCRJg7B61APQYBkiJUnSILkyOcUMkpIkadBcmZxChkhJkjQsrkxOGYOkJEkaJlcmp4QhUpIkjYJhcsIZIiVJ0ii5zT3BDJKSJGnUXJmcQIZISZI0LlyZnCC731UGSU2FJAcm+VqSW5LcnOTtrfwFSb6R5LtJ/keSPbraHNGO3dyOrx3dv0DSrEuyKsn1Sa5on89NcmOS7yT5kyQ/08rXJ/lJK/9Okt/v6uPINp9tTvKhJGnluyT5dCu/Osn6rjZnJLmtvc7oKj+41b2ttV0zrP8WhskJYYjUlNkO/Puq+qfAMcBbkzwf+Cjwa1X1z4DLgP8AkGQ18CngrKo6DDgWeHwUA5ek5u3ALV2ff7OqjqiqDcAVwH/qOnZ7VW1or7O6yj8CnAkc2l4bW/mbgG1VdQhwAXA+QJJnAucALwaOBs5Jsndrcz5wQVUdCmxrfQyFYXLMuRqpaVRV91TVt9v7h+hMyAcAPw98vVW7Cnh9e38CcGNV3dDa/LiqnhjuqCWpI8k64NV0vgADUFUPdlXZDdjh//NOsj+wR1V9o6oK+ATw2nb4ZODi9v4zwPFt1fJE4Kqq2lpV2+jMkxvbseNaXVrbub4Gzmsmx5QBUrOibd+8ELgauAk4Cfg88AbgwFbtuUAluRLYD7ikqj4w9MHuxN6r/p5T9vj2qIchjb1fH/UAVu6DwDuB3bsLk5wHnA48ALy869DBSa4HHgR+var+nM4X6C1ddba0MtrfHwBU1fYkDwD7dJfPa7MPcH9VbV+gr4EzTI4hg6QGIY88xppbt+y8Yn/sm+Tars8XVtWFTxlT8r8AnwX+bVU9mORXgA8l+U/A5cBjrepq4F8A/xz4e+ArSa6rqq8M9F8haSyM0/yV5DXAvVV1XZJjuxtW1buBdyc5G3gbnS3pe4CDqurHSY4EPpfkMCALnHsuACx2bLnlQ2GYHCOGSE2R+6rqqB1VSPJ0OkHyj6rqjwGq6lY6W9okeS6dbSTofMv+s6q6rx37IvAiwDApqd92Nn+9FDgpyauAtcAeST5VVad11flvwBeAc6rqUeBRgBZAb6ez27IFWNfVZh1wd3u/hc7OzJZ2zfiewNZWfuy8NpuA+4C9kqxuq5PdfQ2c10yOAa+L1Kxp1/d8DLilqn6nq/xZ7e/T6OyEzd31eCVwRJJntIn1F4G/Gu6oJQmq6uyqWldV64FTga9W1WlJDu2qdhJwK0CS/ZKsau9/ls6NNndU1T3AQ0mOaXPi6XQu8YHOzszcndqntHMUnbnwhCR7txtvTgCubMe+1urS2s71NXCuTI6QAVIz7KXALwPfTfKdVvb/AocmeWv7/MfAHwJU1bYkvwN8i87WzRer6gtDHrMk7chvJPl54KfAXcDcXdu/ALwvyXbgCTpPpdjajr0F+DiwK/Cl9oLOl+1PJtlMZ0XyVICq2prkXDpzIcD7uvp6F3BJkvcD17c+hsIwOWQGSAmq6i9Y+BofgP+ySJtP0Xk8kCSNharaRGebmap6/SJ1Pkvnkp6Fjl0LHL5A+SN0bkJcqM1FwEULlN9B53FBQ2eYHDDDoyRJmmb/f3t3HmtXWa9x/PtYxKFXFMWJQUGts6b3QnAgKqKROsQpYtCo/GGsGkj0xuRekcThqonzlCiKIjijgqhxoOJsjBNIBSoSKk6VXqsWRaMChZ9/rPeEnXLanq599ni+n2TnrP3u9a79rp6TX5+91/AaJkfAAClJklYKw+QyMUBKkqSVyDA5BAOkJEla6QyTPRgiJUmSOobJJTJASpIk3ZxhchEGR0mSpKVZ0WHS0ChJkjScFRMmDY6SJEnLb27CpGFRkiRp/KYqTK66zlAoSZI0S24x6QFIkiRpdhkmJUmS1JthUpIkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUpIkSb0NFSaTHJdkU5Ibkxwx0H5okn8m2dge7x9+qJK0vKxhkjS8YadTvBR4JvCBRV77ZVWtHXL7kjRK1jBJGtJQYbKqLgNIsjyjkaQxsoZJ0vBGec7kYUkuSvKdJI8a4ftI0ihYwyRpCfb4zWSSrwN3W+SlU6rqC7vothW4R1X9OcnhwOeTPKiqrllk++uB9QD7rt5/6SOXpCUYZw078KBVyzVsSZoZewyTVfX4vd1oVV0LXNuWL0zyS+C+wAWLrHsacBrA6jsfUnv7XpK0O+OsYQ9+6L7WMEkrzkgOcye5c5JVbflewBrgylG8lyQtN2uYJC3dsLcGekaSLcAjgC8n2dBeejRwcZKfAWcDL6mq7cMNVZKWlzVMkoY37NXc5wLnLtJ+DnDOMNuWpFGzhknS8JwBR5IkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUpIkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUpIkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUpIkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUpIkSb0ZJiVJktSbYVKSJEm9GSYlSZLUm2FSkiRJvRkmJUmS1JthUtLYJTkkybeSXJZkU5KXtfa1SX6YZGOSC5Ic2dpvmeQjSS5pfU6e7B5IWumSrEpyUZIvteevT3Jxq19fS3LgwLonJ9mc5PIkxw60H97q2uYk70mS1n6rJJ9u7T9KcuhAnxOSXNEeJwy0H9bWvaL13Xcc/w5gmJQ0GTuAV1TVA4CHAycmeSDwFuB1VbUWeHV7DnAccKuqeghwOPDiweIqSRPwMuCygedvraqHtvr1JboaRqttxwMPAtYB70uyqvU5FVgPrGmPda39hcDVVXUf4J3Am9u27gi8BngYcCTwmiT7tz5vBt5ZVWuAq9s2xsIwKWnsqmprVf20Lf+NriAfBBSwX1vt9sBVC12A1Un2AW4DXAdcM9ZBS1KT5GDgycCHFtqqarAmraarWwBPA86qqmur6lfAZuDIJHcH9quqH1RVAR8Fnj7Q5yNt+Wzgce1by2OB86tqe1VdDZwPrGuvHdPWpfVd2NbI7TOuN5KkxbRvGP8T+BHwcmBDkrfRfdh9ZFvtbLriuhW4LfDfVbV97IOVpM67gP8BbjfYmOSNwAuAvwKPbc0HAT8cWG1La7u+Le/cvtDndwBVtSPJX4E7Dbbv1OdOwF+qasci2xq5qQqT//jTlr9fePorLp/0OMbgAOBPkx7EGLifo3PPve1wzY4/bjjv/993wCgGs4hbJ7lg4PlpVXXazisl+Q/gHODlVXVNkjfQBcVzkjwbOB14PN3hnBuAA4H9ge8l+XpVXTnyPdkLmy65/u/3v8dWa9j8cD9HZ69q2DTVryRPAbZV1YVJjh7sWFWnAKe087pPojsknUXeo3bTTo8+u9vWyE1VmAQur6ojJj2IUUtygfs5P2ZlP6tq3Z7XGp8kt6QLkp+oqs+15hPozkMC+Cw3HUJ6LnBeVV0PbEvyfeAIYKrCJNawueJ+To8pq19HAU9N8iTg1sB+ST5eVc8bWOeTwJfpwuQW4JCB1w6mO4VnS1veuZ2BPlva6T23B7a39qN36vNtug8Dd0iyT/t2cnBbI+c5k5LGrp3fczpwWVW9Y+Clq4DHtOVjgCva8m+BY9JZTXfRzi/GNV5JWlBVJ1fVwVV1KN2FNd+squclWTOw2lO5qUZ9ETi+XaF9GN2FNj+uqq3A35I8vNXEFwBfGOizcKX2s9p7FLABeEKS/duFN08ANrTXvtXWpfVd2NbITds3k5JWhqOA5wOXJNnY2l4FvAh4d/sk/i+6qxwB3gucAVxKdzjnjKq6eLxDlqTdelOS+wE3Ar8BXgJQVZuSfAb4Od2dLE6sqhtan5cCdbZukwAAApdJREFUZ9JdWPjV9oDuw/bHkmym+0by+Lat7UleD/ykrfd/A+eP/y9wVjtd6KK2jbFIF2anQ5L1i51XNW/cz/myUvZTe7ZS/hbcz/myUvZTozNVYVKSJEmzxXMmJUmS1NvUhMkk69o0Q5uTvHLS4xmVJL9uUydt3OnWAzMtyYeTbEty6UDbHZOc36Z2On/gLv0zaxf7+dokv2+/043tCj+tINav2Wb9sn5pOFMRJtu0Qu8Fngg8EHhOm35oXj22qtZO+60Y9tKZ3DQN1IJXAt9oUzt9oz2fdWdy8/2Ebgqrte3xlTGPSRNk/ZoLZ2L9sn6pt6kIk3Q3JN5cVVdW1XXAWXSzXWhGVNV36a44GzQ4HdRYp3YalV3sp1Y269eMs35Jw5mWMLmr6YHmUQFfS3JhkvV7XHu23bXdR4v28y4THs8onZTk4nYYaeYPh2mvWL/mk/VLWqJpCZMTnQZozI6qqv+iOyR2YpJHT3pAGtqpwL2BtXRzR799ssPRmFm/NMusXxratITJXU01NHeq6qr2cxtwLt0hsnn1hyR3B2g/t014PCNRVX+oqhuq6kbgg8z371Q3Z/2aT9YvaYmmJUz+BFiT5LAk+9Ld6f2LEx7TskuyOsntFpbppkG6dPe9ZtrgdFBjndppnBb+w2mewXz/TnVz1q/5ZP2SlmgqplOsqh1JTqKbc3IV8OGq2jThYY3CXYFzuyk42Qf4ZFWdN9khLY8kn6KbfP6AJFvoJrd/E/CZJC+km1v5uMmNcHnsYj+PTrKW7tDmr4EXT2yAGjvr1+yzflm/NBxnwJEkSVJv03KYW5IkSTPIMClJkqTeDJOSJEnqzTApSZKk3gyTkiRJ6s0wKUmSpN4Mk5IkSerNMClJkqTe/g3aWHIEsd3UiAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 720x720 with 4 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "m_plot=vtuIO.VTUIO(\"th_decovalex_ts_30_t_8.000000.vtu\", dim=2)\n",
    "triang=tri.Triangulation(m_plot.points[:,0],m_plot.points[:,1])\n",
    "u_plot = m_plot.get_point_field(\"T\")\n",
    "p_plot = m_plot.get_point_field(\"p\")\n",
    "fig, ax = plt.subplots(ncols=2, figsize=(10,10))\n",
    "contour_left = ax[0].tricontourf(triang, u_plot)\n",
    "contour_right = ax[1].tricontourf(triang, p_plot)\n",
    "fig.colorbar(contour_left,ax=ax[0],label='$T$ / [K]')\n",
    "fig.colorbar(contour_right,ax=ax[1],label='$p$ / [MPa]')\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "a823f8e3",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Wed Nov  3 21:37:20 2021\n"
     ]
    }
   ],
   "source": [
    "import time\n",
    "print(time.ctime())"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "54a4390e-6829-4e80-a9d1-27438146702d",
   "metadata": {
    "tags": []
   },
   "source": [
    "**Results have been shown using [matplotlib](https://matplotlib.org/)**"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "3a47f1f1-7d5f-4950-9284-26b1fc1ed88f",
   "metadata": {
    "tags": []
   },
   "source": [
    "**OGS links**\n",
    "- description: https://www.opengeosys.org/docs/benchmarks/hydro-thermal/decovalex-th/\n",
    "- project file: https://gitlab.opengeosys.org/ogs/ogs/-/blob/master/Tests/Data/Parabolic/HT/StaggeredCoupling/ADecovalexTHMCBasedHTExample/th_decovalex.prj"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "f3d0f9af",
   "metadata": {},
   "source": [
    "**Credits**\n",
    "- Jörg Buchwald for [ogs6py](https://github.com/joergbuchwald/ogs6py) and [VTUInterface](https://github.com/joergbuchwald/VTUinterface)\n",
    "- Wenqing and Jörg for this benchmark set up\n",
    "- Pierre Navaro for providing code https://perso.univ-rennes1.fr/pierre.navaro/read-vtk-with-python.html"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 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.8.10"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}