{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# 5x3, 6x3, 6x4 assemblies\n", "\n", "In reality, rectangular LWR fuel assemblies are symmetric, having the same dimension both in x and y direction. Nevertheless, this example is just to illustrate through plotting the geometry how the origin (0,0) will be in the \"center\" of the assembly even when the number of rows and columns are not equal in the assembly." ] }, { "cell_type": "code", "execution_count": 1, "metadata": {}, "outputs": [], "source": [ "from feign.geometry import *\n", "from feign.blocks import *" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Materials\n", "\n", "The fuel pins are made of uranium-dioxide, in zirconium cladding. The pool around the assembly is filled with water. Outside the pool air is considered (in reality there is ofc concrete walls, but the direct gamma ray won't pass through that, so no need to model those).\n", "\n", "In front of the detector, three absorber sheets are placed made of lead, aluminium and copper." ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "uo2=Material('1')\n", "he=Material('2')\n", "zr=Material('3')\n", "h2o=Material('4')" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## Pins\n", "\n", "There are two types of pins in the assembly:\n", "\n", "- fuel pins\n", "- empty water channel" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [], "source": [ "fuel=Pin('1')\n", "fuel.add_region(uo2,0.52)\n", "fuel.add_region(he,0.53)\n", "fuel.add_region(zr,0.615)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 5x3" ] }, { "cell_type": "code", "execution_count": 4, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Warning: no pool in the problem, the surrounding of the Assembly is filled with coolant material\n", "No absorbers in the problem\n", "Warning: elines missing; only distance travelled in various materials will be computed\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAAD8CAYAAABUzEBbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xl0VOeZ5/Hvc6uk0ooktCKEQEIbOwaMWYyN19iJO07STsdOOnGc5DiJnZnTyemezkzO6fTp/qe7Z3pyZuIkPe7EWebE2brjtidxHGODzWLAyGwWQhsSSEJCC1rQXqq67/yhkqylJJXQvSDg+ZzDUdW9t169pXL9fJf3Pq8YY1BKqfmyrncHlFI3Bw0TpZQjNEyUUo7QMFFKOULDRCnlCA0TpZQjHAkTEXlBRFpFpGzcssUiskdEqkM/U6Z57ZOhbapF5Ekn+qOUuvac2jP5CfDQpGXfBN40xhQCb4aeTyAii4FvA3cAW4FvTxc6SqmFzZEwMcbsBzomLX4U+Gno8U+Bj4V56YeAPcaYDmNMJ7CHqaGklLoBeF1sO9MY0wxgjGkWkYww2ywFGsY9bwwtm0JEngaeBvDFxW3OzC90uLtKqVH1ZafajTHpc3mNm2ESCQmzLOz4fmPM88DzAMvXbTTffOUNN/ul1C3tmfz0C3N9jZtXc1pEZAlA6GdrmG0agWXjnucATS72SSnlEjfD5BVg9OrMk8DLYbb5I/CgiKSETrw+GFqmlLrBOHVp+BfAYaBYRBpF5IvAPwAPiEg18EDoOSKyRUR+CGCM6QD+HjgW+vd3oWVKqRuMI+dMjDFPTLPqvjDblgJfGvf8BeAFJ/qhlLp+dASsUsoRGiZKKUdomCilHKFhopRyhIaJUsoRGiZKKUdomCilHKFhopRyhIaJUsoRGiZKKUdomCilHKFhopRyhIaJUsoRGiZKKUdomCilHKFhopRyhIaJUsoRGiZKKUe4GiYiUiwiJ8f9uyIifzFpm90i0j1um79xs09KKXe4Om+OMaYS2AggIh7gIvBSmE0PGGMecbMvSil3XcvDnPuAc8aYOU/uo5Ra+K5lmDwO/GKaddtF5JSI/EFE1lzDPimlHHJNwkREooGPAr8Js/o4sNwYswH4LvAf07TxtIiUikhpb8dl9zqrlLoq12rP5GHguDGmZfIKY8wVY0xv6PGrQJSIpIXZ7nljzBZjzJaExanu91gpNSfXKkyeYJpDHBHJEhEJPd4a6pPueih1g3H1ag6AiMQxMj3ol8ct+wqAMeZfgMeAr4pIABgAHjfGGLf7pZRyluthYozpB1InLfuXcY+fA55zux9KKXfpCFillCM0TJRSjtAwUUo5QsNEKeUIDROllCM0TJRSjtAwUUo5QsNEKeUIDROllCM0TJRSjtAwUUo5QsNEKeUIDROllCM0TJRSjtAwUUo5QsNEKeUIDROllCM0TJRSjtAwUUo5wvUwEZHzIvJ+aB7h0jDrRUT+t4jUiMhpEdnkdp+UUs5zvaB0yD3GmPZp1j0MFIb+3QH8IPRTKXUDWQiHOY8CPzMjjgDJIrLkendKKTU31yJMDPC6iLwnIk+HWb8UaBj3vDG0bAKdHlSphe1ahMlOY8wmRg5nnhWRuyatlzCvmTIJl04PqtTC5nqYGGOaQj9bgZeArZM2aQSWjXueAzS53S+llLNcDRMRiReRxNHHwINA2aTNXgE+F7qqsw3oNsY0u9kvpZTz3L6akwm8FJqX3Au8aIx5bdJcw68CHwZqgH7gKZf7pJRygathYoypBTaEWT5+rmEDPOtmP5RS7lsIl4aVUjcBDROllCM0TJRSjtAwUUo5QsNEKeWIa3Wjn7qBGNum72I9gb5exOMhJi0DX8rcRx0H+vvob24k6PfjiYkhIWc5VlT0nNsZaG1mqLMDgOikZGIzswkNN1ALiIaJGhMcHKD12Du0HH6bluZm+vp6EbFIS11MekExWXfeS1Lhqlm/yP3NF7l0cC+t75+gtb2NYf8wMbExZGZkkLHpDjJ37J41nIwdpOP0cZoP7qWt/jxd3d0YA8lJi0jLyWXJzntIve12LI/+J7xQ6CehABjq6qDihec4e/IEVdXVDA0NjYWGbdukpb7L+rLTFN77ILkf+QQi4Y+Q208co/Lff86pkydoamqCUBvGGDweD7knT7Lh+FFKPvs0iSsKwrZhDw9T8+KPqHr3Hc6UldF95QqWZY31JTEhgdVnz1C0eStFn/syHl+MC38RNVcaJorAQD8VLzzHob1vcrGpCdu2p2zT1t7OW2/tY8g/hBUdzbIHPzplm67KM5T/+mfs3fsmfr+fkfGIH7Btm7q6OtpaW7GDNuu+8nXisibeIG6M4dyvfsyJfXs4ffr0WF+CweDYNld6enj33Xe5cuUKYlkUP/UMYnmc+FOoedATsIqWI/s5c/z4tEEyKmjbHD58mPr9b+Lv7pywztg251/+FQcP7GdoaGhKkIyybZue3l6OHDlM/asvTVnfU1fN+ePHJgTJdO1UVFRw7kQpXRWTb/dS14OGyS3OBIO0HDlAdXXVjF/ese2NobqqitZ3D01YfqWmgqaGevr6+yNqo6WlhbaaSoY6Jhbgu3ToLc6cORNRX2zb5mx5OZcO7pt1W+U+DZNbXM+FWi41N+EfHo5oe9u2OX/hAu0n3p2wvO29I1RXRRZIACLCuZoaLr9/fGxZcGiQzpoK2traIu5/R2cnl8/X4u/pjvg1yh0aJre44d4r9PT0hK1QNR2/38/wpD0Q/5Uu+vr6Im5j9HDH3/XB4VKgv4/+gcE59GQklPoG+gn09szpdcp5Gia3OLE8WGJNLW0302tEsDwT/9MRyzPnsR8ey0LGXdoVy8LjsaY93zIdy7L0BOwCoGFyi4tJyyAtLXVOX+CEhAR8yYsnLIvLyiY1NTXiQLEsi9S0NGIzs8aWeeMTiYuLw+uJPBhEhEWJiUQnJUf8GuUODZNbXFxWNum5K0hatCii7S3LorCgkMwduycsz9h6J0XFJRH/XhEhLy+P1HUfTJNkeb1kbLqD3NzciEJJRMhesoT0tRvxxMRG/LuVOzRMFEt23cf69RvGBoZNR0SIjY1lRV4ei9fdNmFdbOYSMlcWkL1kyaztWJZFcVER6Ru2TAmBjG13UbJqFV7v7EOgPB4Pa9auJXP73bNuq9ynYaJYvH4ThTvvYvPmzSPnH8LsFViWRVxcHLt376bw8c/jifZN2Sb/k59j+513kp6ePm2gWJZFXl4eG+7YTu5Dj05ZH5OaTsEjf8rdd99NdFRU2HYsy8Lr9bLrzjspePAR4pfmXsW7Vk5zbQSsiCwDfgZkATbwvDHmf03aZjfwMlAXWvRbY8zfudUnFZ6IxcpPfg4rOprFi1OpqqygsbGR4UAAESE+Pp6CggLy8/MpevwpkgpXhW3Hl5LKmqe/jifKR33tOaoqK+jo7MQYg2VZZKSnU7J6NTlFJRR97ivTHppkbb97LLxqamqora1lcHDkKo/P52PF8uUUFhWR98AjLLnrftf+LmpuZK5nziNueGRWviXGmOOhCvXvAR8zxpSP22Y38JfGmEfm0vbydRvNN195w9H+qhEDLc20HNlPR9kJhgeHEEuISUohc+du0sIcloRjgkG6Ksq4dGgfvU0N2LaN5fGQlF9E1s7dJOYVTHtvz3j+nm7aSg/T+u4h/L09gBAdH0/65m2k375jyklg5Zxn8tPfM8ZsmctrXAuTKb9I5GXgOWPMnnHLdqNhsmDZwQBiWRF98adjjMHYwXnf3WvskXtz9BLwtXE1YXJNzpmIyArgNuBomNXbReSUiPxBRNbM0IZOD3qNWR7vvIIERsekzP9oWiyPBskC53qYiEgC8O/AXxhjrkxafRxYbozZAHwX+I/p2tHpQZVa2FwtQSAiUYwEyc+NMb+dvH58uBhjXhWR74tImjGmffK26towxqanrob2E8cYvtKFeLzEZy8j/fbtRCelRNxOcHCA9lOlXKmpIDg0hDcujpTVG0hZvR7LGxVxOwMtzbSVvsNAeysYgy81nfTN24nPzrmat6dc5ObVHAF+BJw1xvzPabbJAlqMMUZEtjKyp6THMNdJV1U5F/7fb7jUUE91VRV9/f1YIqSlp1P01h9JL15N3sc/TVRC4rRt2MEADa+9TEvpYepqa2m8eJHA8DDRPh8r8vJYtiyXpbsfIGvnvTMOTBtou0Tdb1+kNXRVqKurC4BFixZRfGAv6bkryPvEZzRUFhA390x2Ap8F3heRk6Fl/w3IhbFZ/R4DvioiAWAAeNxcqzPCaoL2E+9S/uv/y8GDB+jr65tw929XdzfV1dXkl5fTd6mJNV/+BtGLkqa0YQcCVP3sXyg7tJ+TJ09ijJnQTnt7O6dOnGBbWytDl9tZ/tE/CxsofU2NlP/ou7xz6CAtLS0AY8P9u7q7qW9oIC01lZ3trax+6hkSl690+s+hroJrYWKMOQgz34xqjHkOeM6tPqjI9F2sp/K3L45USBsamnLT32gg1J0/z3AgQFRsHGue/S9TgqD+1d9SdnA/J06eCFuKwLZthvx+Dhw4AAhxWUvJuOPOCdsEBweo+Mn3eWvfXrpCY1TGG33efvkye/fuxfJ42fj1bxGVENntAMo9OgJW0bz/DY6Xlo6UWpxhO9u2aWxs5OK5anrqaiasC/T3cem9I5w8dXLWmia2bXPs3aM07d+DmbRt+6lSakKHNfYMO6nGGHp7e6k4W07bscOzvkflPg2TW9xwbw9t5ae51NIS0Z3Dtm1TUV7OpUMTq5u1nzhKbW1txHcfDw4N0dTQwJWairFlxhguHdpHTU1NxJXWas+d49KR/WPjUNT1o2Fyi+ttOE9zS8ucapG0tbfTW187YVlX5VkaGxoirrRm2zb15+u4Uls9tizQ10NPe3tEpR9HDfn9dHd2MKRjj647DZNbnO0fYjhMJfkZX2PbBCeVebT9gwwHAnP63X6/n8DgB8ER9PsZDgbmFGwiwpB/mODQ3Cq0KedpmNziPL4Yon2+OX2BPR4PnqiJY0U8vhiioiIfPwLgi4nBGxs/rg0f0V7vnILNGIMvOlrrmSwAGia3uITcPJZkZc3pC5yZkcGilcUTliWvXs/y3NxZa5mMsiyLFXl5E+5A9sYlkJiRSWJCQsR9iYmJISU9/aqmL1XO0jC5xXnj4klfu5Gl2ZHN32uJUFxSQtakSmtpG7aQl58fcZjExsaSlbOMxLwPZvUTEZbsvJfCwqKI2rEsi5UrV5K5/S4kwt+r3KOfgGLJXfezafNmYmNiZgwUy7LIz88nu3gV8ctWTFjniYkl645dbN6yZdYg8Hg8bN+2jaX3Pjzl9y1ev4mVJSWkp6XN2I5lWSQnJ1Oyeg3pm7bN/iaV63R6UEVcZjbFj3+eYNDmnXcO0dk5Mv3E6KHPaPW1woICbtu+g+LPfjls6Cx78BGGOtoQ4Pjx4wQCgQlXdyzLIiYmhm3btlH8oT8h7batU9rwRPsoeepZ7ECAo4ff4eLFi8AHA+csa6R6fVZmJtt27KTk81/FGxc/pR117WmYKABSVq1n/Zf+E0lZ2bTV143cm9PXh1gWGRmZ5OfnkbFuE7mP/CneaU52iuWh4PEvEJe1lKU5y2hsqKehoYFhvx+fz0feypVk5yxj6X0Pk755+7R9iU3PZO0zf0lizm9oqyqnuqqKrq4ujDEkJSdTWFhIekExKz76Z8SkZbj1J1FzdM2KIzlJiyO5q6+pkcunSxnu7kK8XuKW5JC28fY57QHYgQCd5ae4UltFcHAQb1w8KavWsWhl8ZzOb/i7O2k7fpSh9laMgZjUNNJu24pvcdrVvDUVoaspjqR7JmqK+Oyced+Na3m9pK7fTOr6zfNqJzophaX3PDSvNtS1oSdglVKO0D0TFZYdDBAcHEA8Xjw+31WVbzTGYA/7sf1DeHyxWHMc1DbWjm0THBwAwBMTo+UbFygNEzXGGENvfR0t77xFZ+UZhoaHsUSITUgk4447ydiyg6jE2W/1D/qH6Dh9nOZDe+lrbyMQDBLtjWJRTi5ZO+8huWQtEsEUoIOX22g9sp/W40cZGhwCINoXTfqGLWRuv4vYjCXzfcvKQRomCgB7eJhzv/4pF44f42z5mbG7iI0xxPh85J06SdFbeyh87DMTpvScrL+5kYqf/IDa6mqqq6u40tMzdhl5cUoKRSdKWVGyipLPPzvt/MDGGJrfep1ze35HZUUF9Q0NDIfuBfJ6veQcP86qIwdYcdf9LHvoo/Mueq2coWGiMHaQmhd/yIl9b3D69Okpd/4ODg1RWVlJXV0dfr+ftZ+xWLx245R2BlovceaH3x0pbNTVNdbO6BXDyx0dHD16lIuNjdh+P2ue+Sui4qcOnW9+64+UvfwbDhw8SDAYnDDUPxAIcOHCBRoaGth+5QrYNrkf+YSTfw51lTTSFe0njlFx9J2wQTLKtm0GBwd566191Pz7z8fOYYx37jc/49CB/XR2ds7YTuPFixw/eoSGP7w0Zf1ASzM1r/+eAwcPEggEwt4zZIwhGAxy+MgRzr31Or31dVO2UdfetZjq4iERqRSRGhH5Zpj1PhH5VWj90dAcO+oaMcbQfOBNzrz/fkS1SIaGhqirq6X9VOmE5X1NDbRdqKO1rW3WmwZt26a2ro7W908Q6O+bsK7lyH7Olp8hGJy92JEdDHLmzBkuvfPWrNsq97kaJiLiAb4HPAysBp4QkdWTNvsi0GmMKQC+A/yjm31SE/U3N9LaUE9Pb29E29u2TU1NDa1HDkxY3vruISrOno349xpjqKut4/Lp4x+0HRim7dR7NDY2RnQXswGaL12ivfx9AgORF1RS7nB7z2QrUGOMqTXG+IFfAo9O2uZR4Kehx/8G3CdzKa6h5mWoo53Ozo6I7/YF6O/vZ6i7c8KywbZLY0PeI2HbNh0dlxlsbxlbNtzbQ39fL8EIq7WN6unrxX+la06vUc5zO0yWAg3jnjeGloXdxhgTALqBKcUpdHpQdxh75IoNc7ytYkrVeNvMWIx6ujYmFJQ2BoPMqVATMFJ42r7xbgu52bgdJuH+q5j8qUeyjU4P6pLoRUksWpQ0pyCIjo4matJ9Or6UxSTMoaiRZVksWrQIX/LisWXe+ARiY3xz6MlIICXExhGVOP3EYOracDtMGoFl457nAE3TbSMiXiAJ6HC5XyokITePrKVL8fki+xKP1jTJ2LJjwvL0LdspLi6O+HDJGMPKggJSN3xw744n2kdqyVqyMjMj7n9aaiqpKwt13pwFwO0wOQYUikieiEQDjwOvTNrmFeDJ0OPHgL06q9+1I5ZF5va7KSkpibi6WcHKlaRvmVhCIDGvgKxluSQmJs56mGJZFkuzs0kvWjVl/uLMHbtZvWZNxH1ZvWYtWTvvmXVb5T5XwyR0DuRrwB+Bs8CvjTFnROTvROSjoc1+BKSKSA3wDWDK5WPlrow77mTVbZtZsWLFjF9ir9fLrjt3sfyBj0wZVi9ikffxT7Nr113ExsZOGyiWZZGclMTWbdvJ/fDUwWYJuXnkbbuTzZs2z1ppbf369eTffgdJRaum3U5dO66PgDXGvAq8OmnZ34x7PAh80u1+qOl5on2UfOFrIELq6ZNUVlbS09MzVtUMYElWFmvXraPwwUdYsuv+sO0syi9kzZ9/Ca/Xw/unT1Pf0EAwGEREMMYQHR1Nfn4+q9espeSpZ8IWNhIRVnz8CUzQJiEhgbKy92m/fHksWGzbZnFKCmvWrGXl1m0UfOrzOpx+gdDiSGqMPTzM5dPv0XxwLx1NjfT19WN5LFKSkklbNXI4kZCbN2s7gx3ttB4euUHvcsdlhoeH8UX7SE1PJ2vbXaTfvoOohJlPmBpj6K4amTmw/VwVV670YIBFiQmkrVg5csPgqnVaSNolV1McScNETWGMYehyG8O9PYjXS8zitKuqs2oP+xlovYQ97MfjiyU2M+uqygf4r3Th7xqZxDw6KXnCFSDlDq20phwhIsSkZcy7vqoVFU380tx59yd6UTLRi8LfYawWDt1HVEo5QvdM1BSB/j66KsoY7utFPB5i0zPnXAgaRgpT99bXYvv9eGJiWFRQQswcC0HbgWG6K8sZ6uqA0GFOcslarKjoObWj3KdhosYMdrRz8c1XaSs7yYUL9fT19mBZHjKyMsnOWUbWjt1kbr971ippnWdOcXHfa7TVn6ehvp7hQACfz8fy5ctJLyhm6X0Pk7h85YxtBAcHaHp7Dy3H3qGxsYHOjpEwSUpOITc3l4xNW8m+5yGdM2cB0TBRAPRdrOfsj79H6bFjNDSM3E41enK+uqYan8/HupoaVlefpfCzT2N5p9ZzNcZwcc/vqHz9d5QeO0b3lStjl4VFhPLyctLSjrH1XBUlj/152Em4YOSGv7M/+i5lx45SWVU1oa6JiHDy5AlWlpWxsaKMVV/6z1MGvqnrQ8NE4e/upOLH32fvm2/S3d09pa6JMYaBgQFKS48xPOzHExdHwaeemtJO67sHOfvaK7z99ttj9UhGQ2D0Z3t7O2+88QbGwIZFSVMmQDfBIJU//QHvvr2PuvPnw/bFtm2qqqvo6+9DLA9rnv0rPNFzu6dHOU9PwCouHdzLyRPHwwbJeLZtc+rUKRqPlzLQ2jxxXWCYhj2/59ChQzMWNjLG4Pf7OXz4HRpee3nK+q6KMurKy8IGyeS+XLx4keoz79Px/okI3qVym4bJLS7oH6L1+FEuXLgQUaU1YwxVlRW0HN4/YXln+WkaGurHCj/P1kZXVxetF+roa2qcsO7Sob1UnD0bUV9s26aqqopLB/dGXEdFuUfD5BbXe6GWpubmiAsSGWNovHiRjvJTE5ZfPlVK7blzEYXAqHM1NXSOaycw0E9X/YWxidMj0dPTQ0dzE/4uvdH8etMwucUF+vsZ6J9bycNAIEBwaGjisr4+hgYHI27DGEP/wADDvT1jy4ID/QwO++dUHMkSYXBwUMs2LgAaJrc4K8qLxzu38/AigjXpNVaUFyuCibXG83q9E06cijcKj+WZ0yGLATxez1XPFqico2Fyi4vLWkpmRsacvsDJycnETZpNLyE3n8zMzIiLI1mWRXZ2NvE5y8eWRSUkkrBoEb7oyAekeTwekpOS8SVr9b3rTcPkFudbnEZ6QTFpqZF9GS3Lori4hKw7JxYkSr99B0WFRRGHktfrZdmyXFJWrx9bJpZF5rZd5OXnRxRKIkLusmVkbL5D90wWAA0TxZK77mfjbbfhmeUwRURYnJLCshV5JBevnbDOl7yYjNXrKCosnDUILMtiw/oNZG7bNeVwKWPLDoqLi4mJiZn13InP52P1mrVk3LFrxu3UtaGD1hRJBSUUPfgIgUCAg6EpOSdflbEsi9TUVHbt2kXxk18JO6Q+7xOfHik5YNvUnDsHTKxiPxoyGzdsZM3OXSy956EpbUQlLqL4iS8wPBxg//636evrC9uXmJgY7t69m6LHPkNMavq8/wZq/jRMFABL7/8I3rh4kpOTqauro/bcOQaHhhARUpKTKVm9miW5yyn89JeIzVwStg1vbByrv/x1opNTyK8op6qygkstLQSDwZHDmpwcCouKWHLbVlY8+qlp7/FJLl7DuqeeIS4hgcaGeqoqK+nt7QVjiIuPp6ioiJxlueR/4tMsHneYpK4vV4ojich/B/4E8APngKeMMVNmSRKR80APEAQCkRZj0eJI7gkODtB+qpT2944S6A/dNZyRReb23STmFUR82XagtZmWw/vpqavBDgzjifaRvGodGVt3RnwvjQkG6aoso+XIAfydHRgM0YuSybhjFymr1085RFLOWTCV1kTkQUaqzAdE5B8BjDF/HWa788AWY0z7XNrXMFHKXVcTJq6cgDXGvB6qTA9whJH5cpRSN7FrcTXnC8AfpllngNdF5D0ReXqmRnR6UKUWtqs+6BSRN4CsMKu+ZYx5ObTNt4AA8PNpmtlpjGkSkQxgj4hUGGP2h9vQGPM88DyMHOZcbb/VzAKDA1w+eYz240cJ9PchlkVsxhIyd9wdOmcS2f9/BlqaaTn8NlfO12CGh7GifaSsXk/67TsiLghtBwN0VZTReuQA/q4ODCPTmX5wzkTHliwkVx0mxpjwk6eEiMiTwCPAfdPN0GeMaQr9bBWRl4CtQNgwUe4yxnDp0D4a3nyVutpaamtrx67mJCcns6r0CFm5Kyj8zBeJy8yetp1Afx81v/wxlyrPUllxlkstLdjBIN6oKHJyDlD01h6W3HY7Kz72+IwnULtrKjn3m5/RWF9PVWUFvb29GCA+Pp7CY0dZlruc/E88weLVG1z4a6ir4crpcBF5CPhr4G5jTNg7sEQkHrCMMT2hxw8Cf+dGf9TsLr7xe8pf/Y+w40xaWlpoa2sjdXEZ/t4rrHn6L4jLWjqljcBAP+X/5zuUHjowZZyJf3iYmpoaampq2NDUxHBPN0WfCz9epavyDGU/e579b79FX3//hL74/X5KS0spKyvj7v4+Vn/qSdI2zOk8oXKJW9fWngN8jBy6ABwxxnxFRLKBHxpjPgxkAi+F1nuBF40xr7nUHzWD7uqzVL3+O/bv3z9tYSPbtrnc0cH+t9/G64tl419+e0oQ1P32Rd47fIhztbVhh9WPhsKpUyfxejzEL1tBzv0fmbDNcM8VKn/xAvv27WVgYGDadvr7+9n35pt4LA8JOct14NoC4EqYGGMKplneBHw49LgW0H3UBaD5wJucPHFixgppMPIl7ujspOF8HSsqy0gZd4gx1NVBa/n7VFdXz3p/jm3bnDp9itwVK8je/aEJhzutpe9QWVHB4ODgrO0M+f2cLT/DsqMHws5brK4tvTfnFjfU0U5bTSXtlyO7QmbbNpWVFVw6uG/C8rZj71BVXRXxoLZAIEBDQz2d5afHlhnbpuXIAWrr6iKu+lbf0EDre0ewI6jwptylYXKL6790kZbW1jkVJOrq6qJ/Ug3Y3vpaWlpaIq60Zts2zU1N9DVeGFs23NtD75Ur+P3+iPsSDAbp6u5mqEuHC1xvGia3OHs4QDAQmH3DcYwx2JNeYw8HsGc5TJosEAgQ9H9Qsc0EhgnawTkFmwDBQFD3TBYADZNbnDcujti4uLm9xuvF45s4tYQ3Ph5fTEzEbYgIcbGxRCUkji3zxMaKBVAUAAAJIUlEQVQRExU9p0JNtjHExMTgjZ3be1DO0zC5xSWuWEn2kiV4IqyQJiLkLF06ZXxH6oYt5K9cGXGlNYCVBQUsXvNBO97YOJJzV5CSEvmkWomJiSxekk10hAPhlHs0TG5xVlQ0GZvuYPny5RFXNysqLiFz+10TlqesXs+yZblERVDxbHQgXMbyfOKWTLxtK2vnPaxatSqivliWRVFREVl33junQyPlDg0TRdad97Lxtk0kJSXN+CW2LIsNGzaQs+l2YifVgLW8USx74CPs3LlzxoptIkJ0dDQ7tu9g2UMfnbI+uWQtK1avJW/Filn7snTpUgrXrCN1/aYI3qVymxaEUEQnpVDy1DMYDKXHjtHY2IgxZmyOYIAYn4+169azesed5P/pZ8K2k7H1TlZd6cKyLEpLS+nu7p4w17AxhvS0NG7fupWSx/58ytSgAOLxUPzkVzGBAHFx4eca9ng8FKwsYMPtt1PyhWexoiIvQK3co2GiAIhfmsu6Z/6KxXv/QOvpE9Q31NPb24tlWWRkZpGdk0PWzt1kbr8bscLveYgIOQ/8CfFLc8nKe53WC7U0NDQQGB7GFxNDbm4u6QXFLL3vYRKXr5y2L1EJiaz+yjdIXvU6Re++w8XGBjo7Rm70S0pKJjc3l4xNd5B9z4fwxsW79BdRc+VKcSS3aXEkdwX6++iqPMNwbw+W10tMeiaL8ouQOZxcBehraqS3vg572I/HF8OigmJiFqfNqQ07EKC76gxDnR1gDNFJKSSXrNG9EZddTXEk3TNRU3jj4km7beu824nPziE+e351sSyvd8KwfbVw6QlYpZQjNEyUUo7QMFFKOULDRCnlCA0TpZQjNEyUUo7QMFFKOcK1MBGRvxWRiyJyMvTvw9Ns95CIVIpIjYh8063+KKXc5fagte8YY/7HdCtFxAN8D3gAaASOicgrxphyl/ullHLY9T7M2QrUGGNqjTF+4JfAo9e5T0qpq+B2mHxNRE6LyAsiEq7izVKgYdzzxtCyKXR6UKUWtnmFiYi8ISJlYf49CvwAWAlsBJqBfw7XRJhl083+97wxZosxZkvC4tT5dFsp5YJ5nTOZbYrQUSLyr8DvwqxqBJaNe54DNM2nT0qp68PNqznjS3F9HCgLs9kxoFBE8kQkGngceMWtPiml3OPm1Zx/EpGNjBy2nAe+DDB+ilBjTEBEvgb8EfAALxhjzrjYJ6WUS1wLE2PMZ6dZPjZFaOj5q8CrbvVDKXVtXO9Lw0qpm4SGiVLKERomSilHaJgopRyhYaKUcoSGiVLKERomSilHaJgopRyhYaKUcoSGiVLKERomSilHaJgopRyhYaKUcoSGiVLKERomSilHaJgopRyhYaKUcoSGiVLKEa6UbRSRXwHFoafJQJcxZmOY7c4DPUAQCBhjtrjRH6WU+1wJE2PMp0Yfi8g/A90zbH6PMabdjX4opa4dV+caFhEB/gy4183fo5S6/tw+Z7ILaDHGVE+z3gCvi8h7IvK0y31RSrnoqvdMROQNICvMqm8ZY14OPX4C+MUMzew0xjSJSAawR0QqjDH7p/l9TwNPAyzOzrnabiulXHLVYTLb1KAi4gU+AWyeoY2m0M9WEXkJ2AqEDRNjzPPA8wDL120MOx+xUur6cfMw536gwhjTGG6liMSLSOLoY+BBwk8hqpS6AbgZJo8z6RBHRLJFZHT2vkzgoIicAt4Ffm+Mec3F/iilXOTm9KCfD7NsbGpQY0wtsMGt36+UurZ0BKxSyhEaJkopR2iYKKUcoWGilHKEholSyhEaJkopR2iYKKUcoWGilHKEholSyhEaJkopR2iYKKUcoWGilHKEholSyhEaJkopR2iYKKUcoWGilHKEholSyhEaJkopR2iYKKUcMa8wEZFPisgZEbFFZMukdf9VRGpEpFJEPjTN6/NE5KiIVIvIr0Qkej79UUpdP/PdMyljZG6cCXPdiMhqRqrTrwEeAr4vIp4wr/9H4DvGmEKgE/jiPPujlLpO5hUmxpizxpjKMKseBX5pjBkyxtQBNYxMsDUmNA/xvcC/hRb9FPjYfPqjlLp+3JrqYilwZNzzxtCy8VKBLmNMYIZtxoyfHhQYeiY//WacsCsNaL/enXDBzfq+4OZ9b8VzfcGsYRLhnMJTXhZm2eQpPSPZ5oMV46YHFZFSY8yW6ba9Uen7uvHcrO9NRErn+ppZw2S2OYWn0QgsG/c8B2iatE07kCwi3tDeSbhtlFI3CLcuDb8CPC4iPhHJAwoZmQJ0jDHGAPuAx0KLngSm29NRSi1w8700/HERaQS2A78XkT8CGGPOAL8GyoHXgGeNMcHQa14VkexQE38NfENEahg5h/KjCH/18/Pp9wKm7+vGc7O+tzm/LxnZQVBKqfnREbBKKUdomCilHHHDhMl8h+7fKETkb0XkooicDP378PXu03yIyEOhz6VGRL55vfvjFBE5LyLvhz6jOV9GXUhE5AURaRWRsnHLFovIntCtLntEJGW2dm6YMGH+Q/dvJN8xxmwM/Xv1enfmaoU+h+8BDwOrgSdCn9fN4p7QZ3SjjzP5CSPfnfG+CbwZutXlzdDzGd0wYTKfofvqutkK1Bhjao0xfuCXjHxeagExxuwHOiYtfpSRW1wgwltdbpgwmcFSoGHc8xmH5d8gviYip0O7n7PuXi5gN+NnM8oAr4vIe6FbPW42mcaYZoDQz4zZXuDWvTlXxcWh+wvKTO8T+AHw94y8h78H/hn4wrXrnaNuuM9mDnYaY5pEJAPYIyIVof/D37IWVJi4OHR/QYn0fYrIvwK/c7k7brrhPptIGWOaQj9bReQlRg7pbqYwaRGRJcaYZhFZArTO9oKb4TBn1qH7N5LQBzfq44yceL5RHQMKQ0Wwohk5Uf7Kde7TvIlIvIgkjj4GHuTG/pzCeYWRW1wgwltdFtSeyUxE5OPAd4F0RobunzTGfMgYc0ZERofuBxg3dP8G9U8ispGRw4HzwJevb3eunjEmICJfA/4IeIAXQrda3OgygZdGSvLgBV40xrx2fbt09UTkF8BuIC10e8y3gX8Afi0iXwTqgU/O2o4Op1dKOeFmOMxRSi0AGiZKKUdomCilHKFhopRyhIaJUsoRGiZKKUdomCilHPH/AZA/m+m1p1+FAAAAAElFTkSuQmCC\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "assy=Assembly(5,3)\n", "assy.set_pitch(1.62)\n", "assy.set_source(uo2)\n", "assy.set_coolant(h2o)\n", "assy.set_pins(fuel)\n", "\n", "fuelmap= [['1', '1', '1'], \n", " ['1', '1', '1'], \n", " ['1', '1', '1'],\n", " ['1', '1', '1'],\n", " ['1', '1', '1']]\n", "\n", "assy.set_fuelmap(fuelmap)\n", "\n", "det=Detector('D')\n", "det.set_location(Point(174.726, 174.726))\n", "\n", "ex5=Experiment()\n", "ex5.set_assembly(assy)\n", "ex5.set_detectors(det)\n", "ex5.set_materials(uo2,he,zr,h2o)\n", "uo2.set_color(\"#0f0e0e\")\n", "he.set_color(\"#474545\")\n", "zr.set_color(\"#635f5f\")\n", "h2o.set_color(\"#96deeb\")\n", "ex5.Plot(xl=[-10,10],yl=[-10,10],detectorSize=0.4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 6x3" ] }, { "cell_type": "code", "execution_count": 5, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Warning: no pool in the problem, the surrounding of the Assembly is filled with coolant material\n", "No absorbers in the problem\n", "Warning: elines missing; only distance travelled in various materials will be computed\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAAD8CAYAAABUzEBbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Wl0XNd55vv/e6owkxgIFgCCBEESIwnOojiTEmXLGpJYUWJbst22YjuWJdv3trs7veLuvivJSr6kh3RuEnecJTuKhziW7eTKVmJZ1kBSJCHOAyiMBECAAIiRmOcazr4fUIBQQAGoIqtAUHx/a2Gh6pxTu3aphIdn2OfdYoxBKaXulHW3O6CU+nDQMFFKRYSGiVIqIjRMlFIRoWGilIoIDROlVEREJExE5GUR6RSR8mnLVojIWyJS6/+dNsdrn/NvUysiz0WiP0qpxRepPZPvAY/PWPZN4B1jTAHwjv95ABFZAfwxsAfYDfzxXKGjlFraIhImxpgTQM+MxU8B3/c//j7w20Fe+hjwljGmxxjTC7zF7FBSSt0DnFFsO9MY0wZgjGkTkYwg26wGmqc9b/Evm0VEngeeB4hLTHwgc0NBhLurlJrUVF52yxjjCuc10QyTUEiQZUHH9xtjXgJeAsjdst1887W3o9kvpe5rX93guhHua6J5NadDRFYB+H93BtmmBciZ9nwN0BrFPimloiSaYfIaMHl15jngF0G2+TXwMRFJ8594/Zh/mVLqHhOpS8M/Bk4DRSLSIiJfAv4ceFREaoFH/c8RkV0i8l0AY0wP8GfAef/Pn/qXKaXuMRE5Z2KM+fQcqz4SZNsLwO9Pe/4y8HIk+qGUunt0BKxSKiI0TJRSEaFhopSKCA0TpVREaJgopSJCw0QpFREaJkqpiNAwUUpFhIaJUioiNEyUUhGhYaKUiggNE6VURGiYKKUiQsNEKRURGiZKqYjQMFFKRYSGiVIqIjRMlFIREdUwEZEiEbky7WdARL4xY5uHRaR/2jZ/FM0+KaWiI6rz5hhjaoDtACLiAG4CrwbZ9KQx5jej2RelVHQt5mHOR4B6Y0zYk/sopZa+xQyTZ4Efz7Fun4iUicivRKRkEfuklIqQRQkTEYkFPg78LMjqS0CuMWYb8DfAz+do43kRuSAiF4Z6uqPXWaXUbVmsPZMngEvGmI6ZK4wxA8aYIf/j14EYEVkZZLuXjDG7jDG7lq1Ij36PlVJhWaww+TRzHOKISJaIiP/xbn+fdNdDqXtMVK/mAIhIIhPTg35l2rIXAIwxfwd8AnhRRLzAKPCsMcZEu19KqciKepgYY0aA9BnL/m7a428B34p2P5RS0aUjYJVSEaFhopSKCA0TpVREaJgopSJCw0QpFREaJkqpiNAwUUpFhIaJUioiNEyUUhGhYaKUiggNE6VURGiYKKUiQsNEKRURGiZKqYjQMFFKRYSGiVIqIjRMlFIRoWGilIoIDROlVEREPUxEpFFE3vfPI3whyHoRkb8WkToRuSoiO6PdJ6VU5EW9oLTfEWPMrTnWPQEU+H/2AN/2/1ZK3UOWwmHOU8APzIQzQKqIrLrbnVJKhWcxwsQAb4rIRRF5Psj61UDztOct/mUBdHpQpZa2xQiTA8aYnUwcznxNRA7PWC9BXjNrEi6dHlSppS3qYWKMafX/7gReBXbP2KQFyJn2fA3QGu1+qfn53OO4+3vxDA1gbN9ttWGMwTsyzHhfD96x0dvui+314h7oxz3Qh+313HY7KrqiegJWRJIAyxgz6H/8MeBPZ2z2GvB1EXmFiROv/caYtmj2SwVnjM3g9Tra3ztGb201424PiJCUlETG7gNk7D5AXOqKBdvxjo3SfeU87aXHGeq5hdfrIzYmhuTs1WQdeIS0kq1YzpgF2xnpaKXz9Am6rl5kbGwMYwwJ8Qmkl2wjc/9DJGXnLNiGWjzRvpqTCbzqn5fcCfyTMeaNGXMNvw48CdQBI8AXotwnFYTPPU79j/+BxssXqKyqpLOzc2pdTEwM666WUXzyKAVPP8vKnXNfbBtuuUH1975NXe01auvqGBkZmVqXkpJCUdll1hdtpPgLXyUuLfjhqjGGm2//kutHf01VZQXNLS3Ytg2AZVlkXzjPxvPvsf7gEXJ/43cRaylcR1BRDRNjzHVgW5Dl0+caNsDXotkPNT/j81H7j9/h0rG3qaiomPrDneR2u6mrq+PGjRuMucfZ4rBYue3BWe2MtLdS8fff4tixowwMDMxqp6+vj/Pnz9PS3Izt9bD5xT8gZnnyrHZuvv1L3v/Xf6G0tBTbtpk+j73P56OlpYXWmzfZPTAAxpD7W5/E/w+Wuos00hVdl85Qc+500CCZZNs24+PjnHj3XepffQXv6EjAemMM9T/7AaWnTtLf3z9vO23t7Vw8c5qmX706a/1Ieyv17/yK0tJSfD5fQJBMfy+fbXP27FmunzzK0I3rt/GpVaRpmNznjDG0nzpKRXn5nAEwndvtpqGhge4r5wOWj9xsoutGA123bgUNgOls26ahsZGuijI8w0MB6zpOH6eqsjKkvti2TXl5Oe2lxxbcVkWfhsl9buRmE503WxgaHg5pe9u2qauro+PsyYDlnedKqa6pCeu9r1+/Ts/Vix+07fFw6+plWm7eXDCQJnV0dNBdU4l3JLT+q+jRMLnPjff20NvTgxXGSczR0VHcA/0By8ZuddLf1xdyCNi2TW9PD2O3PjjR6xkeZGRkJKS9kukGh4dm9UctPg2T+5wxJuQAmPm6GQtmjzQM9739bYR7KtUYA7fxGVRkaZjc52JTUklJTQ0rUOLi4ohNWhbYzop0kpcvD7kNy7JISUkNuDzsTFpGQnwchHFlxhjDssQkYpJnXxVSi0vD5D63bO06MrNXEx8fH9L2lmWxYcMGXA8eCFiesWs/RcXFWGEEQV5eHunbHph67oiNY+WmraxatSrkvROXy0V6fiExSaEHmYoODZP7nIhF1v6H2FhcHNJ5E4fDQX5+Pq4H9gYsX7Yuj4w1a0lJSVlwzIdlWazOzsZVtInY5NSAdZn7H2bTpk0hDUSzLIuSTSVk7X94wW1V9GmYKDJ2H6T4gQfJy8ubM1BEBKfTyaFDh1j32G8Rs2z5rPUbfvczHDx0mMTExDnbsSyLFWlp7Nm7j7W/8Tuz1ietyWX9/od4cPfuecPNsix2bN/Bhj37SS4oDuPTqmhZrOJIagmzYmIp/sLXELFYmb6Squoq+vv7sSwLYwwiQnZ2Npu3bCH/8Y+Ttf9I0HaW5+ZR8vnncTgcVFZW0NjYiNfrRUSwbZuEhATy8vIoLtlM8e99lbgVK2e1ISKse+pTGGOTlJhERUU5XV1dCBO3khtjWLlyJSUlJeTt3s+GT34OEf03cSmQ2zmTf7flbtluvvna23e7Gx86ttdLb0UZ7aVH6W5uYmRkBMthkZKcjGvzDjL3P0xS9poF23H399Jx9hSd50rp7evF4/ESHxdLmiuDrH0PsXLnHpyJSfO2YYxh8Hot7aXHuHWtisHhIYyB5UlJrMwvJHP/w6QUbtQgiZKvbnBdNMbsCuc1GiYqqPG+HrxDg4jTSVzqChzxCWG3YXu9jHd34fO4ccTFE7/SdVt//J7hQdz9/WBsYpNTg97PoyLrdsJED3NUUHGpK0IqNzAfy+kkIfPOK3DGJC3XqzX3AN1HVEpFhO6ZqFk8QwP0Vr6PZ2gQy+kk3pVJauEmxOEIuQ1jDMPNjQw1NeBzj+OITyAlv5iEjKyw+uJzj9NXVc54XzcYQ2xKGmmbtuKIC21cjFo8GiZqymhXBy1v/Rtdle/TeKOR4aEhLMtBRlYmq9fkkLn3EKsOfxTLMff/NsYYuq9epPXYr+lsbqK5qQmv10NcXBy569bjWr+BNR/9TZLzCufti3d0hJvvvE7X5XM0NTXR292NAVJSU8nNzSVj2y5Wf/Q3Zl2iVnePhokCYKi5karv/S3nzpyhtW2iaqZt2whwveE6sTEX2Fxfx6b6axQ99wJWTOysNowxNP/qVWre/hXnz59naGhoarllWVTX1LAiLY3djQ1s/J1P43pwf9C+uAf7qfruX3P13Flqa2vx2XZApbWysjI2lJez41olm7787+es2KYWl4aJYry3m6rvfZujb7/NwOBgwF27BjC2zdj4OJcuXcLtduNMSKTgs78/q52O945T9eYvOXnyJD5fYBHqyTZ7enp45523QSA2NY2Ugo2B2/m81Hzv25w5foym5uZZdxBPPq+vr2dkeBjL4WTz//WHOGLjIvBfQt2JqJ2AFZEcETkmIlUiUiEi/z7INg+LSL9/6tArIvJH0eqPmlvbqaNcuXxpVpDMNFmM6GbZRUbaAycQsD0eWo6+wXvvvTcrSKYzTBRYOl1aStMbv5h1g2Ff5fs0VJYHDZKZfWlrb+daxft0l82adVbdBdG8muMF/pMxZiOwl4k5czYF2e6kMWa7/2dm5XoVZT73OF2Xz9PU1BRSHRFjDNXV1XScfjdgeU/FFZqabuD1ekN63/6BAbqaGhlpbQ5Y3l56lOrq6pArrdXW1tJ+6thtlVFQkRW1MDHGtBljLvkfDwJVBJmpT91dg431tLW1hlyQyBhDa1sbvVXvByzvuXqRhuvXwypsdL2+nt7KD9rxjgzT19JEX19fyG0MDQ3R29mOu68n5Neo6FiUcSYisg7YAZwNsnqfiJSJyK9EpGSeNnR60CjwjY4wOhreBFlerxefezxw2cgwY+Pjc7xiNmMMI6OjeIYHPmhjbJRxjyesSvMiwtjo2KwC12rxRT1MRGQZ8C/AN4wxAzNWXwJyjTHbgL8Bfj5XOzo9aHRYMTE4nc6wqptZloXlDDx3b8XE4ghjHAqA0+kMGC9ixcTgsBxhH7I4nE6smIUn9VLRFdUwEZEYJoLkR8aY/2/memPMgDFmyP/4dSBGRGbfSqqiJnHVGjIyMrDD+ANOS00lKSvwiHX5ujyysrJCriVrWRZr1qxhWc66qWUxSctZnppKXFzoV2YcDgepqSl6eXgJiObVHAH+HqgyxvzvObbJ8m+HiOz290ePYRZRXFo6GQXFuFaGluGWZVFUvJHMA4FlCFwP7qewYP6BaNPFxMSwOieH1I1bppaJZZG599C8dVWmExHW5eaS8cC+kKYbVdEVzT2TA8DngEemXfp9UkRemJweFPgEUC4iZcBfA88aPS2/6FYdfpQdO3fidM4/7MiyLFamp7N63TpSiwIvzMUmp5JRspWiwsIFg8CyLLZt386qfQ/NGk3remAfRYVFJCYmznvuRESIj49nU0kJmXsOLvAJ1WKI2qA1Y8wpFig0boz5FvCtaPVBhSY5r5DCJ57C5/Vy8uRJPF5vwFUZ/64jLpeLA4cOUfzci4g1+/zI+qc/zWhXB7YxXLt2baoo0qTJkNm5cyclBx5i1UOPzmojZtlyCj/7JdweDydPvMvg0NCsK0SWZZGYmMjhww9R+KnPBy2ypBafjoBVAGQ//BjOhCRSUlNpbGykvr6esbExRIQVaWkUb9xEZu46Cj7zJeJXZgRtwxGfwMbnv0Fc2gryqsq5VlNDW1sbPp8PZ0wMuWvXsmFDHtkP7iP3Nz8RNJAAUvKL2fqlr5OUnMzN5iauVVczODiIAZYlJVFYXMzqNTls+N3PkloYbOiSuhu0OJIK4HOP01t+ha4Lp/EMDyEOB4mZq8jc9xBJOetCvmw71t1F59mTDFyvxfYXR0rbtBXXrn3ELAutuJGxbfrrquk8c5Lx3ltgJobgZ+w+SGpRSVh3MavwaKU1pVRE3E6YaHEkpVREaJgopSJCT8CqAN6RYbounaH70jm8o8OI5SDelUnmvodIKdgY0uRYAMM3m+g4/S6DjdcxXg9WbCypxZvJ2HOI+HRXSG3YXg+9FWV0nDmBu3/ifp2Y5clk7D7Iii07tOzAEqNhooCJe2Xa3n2LpqNvcL3hOtevX8c9Pg4ipKamUnzxHNlrcyn47Jfnne7CMzRI7Y++S3tdDdXV1XR2dk5czXE6WbP6FIWnjrFq6042/O6/m3cIfF9NBXU/+yHNTU1cu1bD8PAwAIkJCRRcOE/uulzWP/UM6VsfmLMNtbg0TBQAzW/8nMo3/pXS0lLsaZXNALq6uuju7iatooLx4WE2f/n/Jmn12llteIaHqPi7/8350pM0NDYCTN1n4/V6ud7QwPWGBja3tOAZGqDo97426x4fgJ6KMip+9F3effddRkdHA/rS7/Fw6fIlKirKOTw0RInHM2uqUnV3aJgo+moqqH3nDU6dOjVnYSPbtunp7eXdd48TEx/P9v/8J7NGrzb8y4+4eOY9Gm/cCHqz3mQolJeXExMTw7K168n52McDtnEP9FH7sx9w7NgxxkZHCXat0fZXfjt+7BgxMbEsz90w59gXtXj0BKyi7eTbXL58ed4KaTCxl9Hf30/LjUb6ZtQzGe/tpqO6nPr6+gVrmti2zdWrV+k4V4rt8QSs6zr3HlVVVYyPjwcNkul9cXs8VFaU03HmxLzvpxaHhsl9bqy7i676Wnp6QisuZNs21dXVtJceC1jeea6U2trakAe1+Xw+mm7coLeybGqZsX10nDtFY2NjyFXfmpqb6bx0DtvjDul9VfRomNznRttb6ezsDLl0AEB/fz+jXR0By4abG+ns6Ai50ppt27S3tzN8s2lqmWdoiKHBQdzu0IPBtm36B/oZ79VKa3ebhsl9zvZ68fl8YRUkMsZgz6j1as+4OTAUPq8Xe1pwGK8Hn+0Lr9Ia4PP6ZvVHLT4Nk/ucMzGRhMTEsF4T43TinDGjXsyyZcTHhz7LnoiQmJgYMAm5IyGR+JjYsILNNoaEhAScYX4GFXkaJve55evyyV6VHXLJRRFhzZo1rNiyI2B5+rYHycsvCOtwKS8vj7SS7VPPnQmJpK5bz4oVoU+YnpycTFr2amJT0kJ+jYoODZP7nBUTQ8aufaxftz7k6mZFRcVk7j0csDx10xbWrM0hNnb2TH/B2khLSyNjfT6JWdkB67IOPMLGjRtD6otlWRQWFrHq4CNhHRqp6NAwUWQdOMLWHTtYkZY27x+xZVns2LGDNbv2zBrXYTmcrP3Yb3Hw4MF5K7aJCHFxcezff4Ccx5+atT61aBPrSrYuWLrRsizW5uRQsGULKzbvmHM7tXh00JoiNjmFjV/8GhjDpYsXaLpxA9sYjDET/+IbQ2JiIlu2baN43yHWPf3poO24du1n4+AAlmVx8cIFenp7Z+0xZLhc7HpwN8XPfI7l6/NntSGWg6LnXsD4vCQlJVFdVYV7ciyKMSCC0+mksKCALbt2T4yi1cr0S0LUw0REHgf+CnAA3zXG/PmM9XHAD4AHmCgm/YwxpjHa/VKBErNWs/lr/5n042/SWXaB5uZmhoaGcDgcuDIyyVqTQ9bBI2Q8eGDem/2yjzxOYnYOmXlv0nm9jtbWm3g8HuLi4sjJycFVtInsI4+zbE3unG04E5PY+Pw3SNt0lIKzp2htaaGnpxsMpKalsnr1GjJ272fV4UdxxidE4z+Hug1RLY4kIg7gGvAo0AKcBz5tjKmcts1Xga3GmBdE5FngaWPMM/O1q8WRoss3Pkb/tSo8w4OIw0GCK4tluRvCPi8x2tnOUHMjttuNIz6e5A2FxKakhtWG8fnor6/B3dsDGGJT0kjOLw56T4+KnNspjhTtb2Q3UGeMuQ4gIq8ATwGV07Z5CvgT/+N/Br4lIqJV6u8eR1z8rKs1tyMhI4uEjKw7akMcDq3zeo+I9gnY1cD0malbmD3f8NQ2xhgv0A/MmlFJpwdVammL9p5JsP3imXscoWyDMeYl4CWYOMy5866puRhjGG1vnTrMiU/PIDY5Jex2fGOjjHS0+g9zEkhctea2Dk/Gem4FHOboHcJLU7TDpAXImfZ8DdA6xzYtIuIEUgC90eIu8LnHuXXxLO3vHaer7SbDw8NYYrEiLQ1X0SayDh4hecPCs/aNdrbT8d5xOssu0HWrG4/bTXxCPK6VLjJ27ydz72Fik+c/d2KMTW/FVdpLj9LVcJ3+gX4wkJySjCtnHVkHj7Bi8w6tUL+ERDtMzgMFIrIeuAk8C3xmxjavAc8Bp5mY4e+oni9ZfO6Bfqpf/hY1ZVeouVbD6MjI1FUbYwyuC+fZ+v5lCj7yOGsee2rOk7Hd5Ze59tMfcLWsjOaWlqnLy8YYnE4n66+WsflcKcXPvRgwz/B0ttdD/U++T93ZUsrLy+nt7Z0ac2LbNikpKZRUvE/Bg3vI/+zva/nGJSKqYWKM8YrI14FfM3Fp+GVjTIWI/ClwwRjzGhPzEf9QROqY2CN5Npp9UrP5xkapfvlvOH3sKE3NzR/csDetvklHRwdHjx7F7XZjxcax+pEnZrXTX1tF1T/9A+8cfWeiHsmMfxPcbje1dXV0dHRgbJvNL/xHEjJWBWxjjOH6z35I2TtvcvnK5am+TK+10tfXx+kzpxkY6AcRCj//Qsi1aVX0RP36mjHmdeD1Gcv+aNrjMeCT0e6HmlvH2VNUXroYGCRB+Hw+SktLSU5OxrVrf8B5FGNsGn7xE06ePBE0SCbZts3A4CCnT58mee3PKfq9FwPWD92op+H8mYAgmaudispK0tJWkLn3MKnFm8P81CrSNM7vc8b20XHmBDXXroVUQsC2bWpra+k6XxqwfKD+Gu3NTQwNDy941+9kLZNbtVWM9wZemWsvPU5VVWXIfamsrKD91NEFt1XRp2Fynxu60UD7zZshFySybZvGxka6Lp4NWN514TQ1NTUh1zQREerq6ui+emlqmc89Ts+1Sjo6OuZ5ZaDunh66G+rwDA6E/BoVHRom9zn34ACDgwNBr8/PZXx8HM/ocGA7/X1T01GEYvJwx933wYU77/AQo2PjYfRkIpSGRkfxDGmY3G0aJvc5sQQRa97izbNeIzLrjl6xrLACCSbu/J1+aVcsC8uSsIojAVgioCdg7zr9Bu5z8eku0tNXhPUHnJSYSNyMYkQJGVmkpaWFfP+OZVmkp6cT7/pguL0zaTkJCYkhF2qatDxpmRZHWgI0TO5zCVmrceXksnz58pC2tyyLgsJCMvYFFkfK2HOQ4uLikN9XRFi/fj3pW3d+0LbTiWvHg+SsyQkplESE7OxsXCXb9O7hJUDD5D4nIqw69BG2btm6YHUzAeLj41m3bt2saTkTs1aTsSGfrMysBduxLIv8/Hwytu7EmRBYuzVz72E2btqEM4S9E8uyKCkpIfPAwwtuq6JPw0SRvm0XhfsPsmPHjonzGEH2CizLIiEhgYcffpiCTz2HI2528egNn/w8+w8eIH3FijkDxbIs1q5dy449+8h54ulZ6xNcmeQ/8RSHDh8mJiYmaDuWZeF0ODhw4AB5H3li3tooavFoUQiFWBZ5z3wByxlLeno61VVV3GxtnRp1mpiYSF5eHnl5+RQ+89ycA8Ti012UfPkbOGLiaKyv5dq1a/T19U28hwiulSsp2riRtYXFFD334qy9kklZBx8Bf/X6a9eu0djYOHXpOiYmhty1ayksKmL9I4+z+iNPRuG/iLodUS2OFC1aHCl6hm820XH6BD0VZXi8HiwR4pYnk7XvIVbu3IMzMWnBNmyvZ+omvaG2VgzgcFgk524g88ARUvKLQxr+Pt7XQ9f5UjrPn8Y9OoIB4uITWLlzNxm7DxKf7rrzD6yCup3iSBomKihjDLZ7HHE476iqmbF92B4PVkzsbd8/Y4zBeD0YYyba0Ur0UbcUK62pe5SIBD0vEnY7lgNH3J2VCRARJGbhKTTU3aUnYJVSEaF7JiqAMYbh5kZuXTmPe6APy+EkMTsH1849AVN5LsTnHqfn/csM1NfgGxvDmZhE2qYtpBZtDqug0Vh3F10XTzPW1QkY4tJduHbumVW6QN19GiZqykB9DY3/+s90NN2grvYaw8PDiGWRkZFB3oY8Mkq2su6pZ+Y9CWtsHy1v/5KOs6dobGikpaV5aqqL9Xl5rF6Tw+ojj5Gx59C85z7Gurto/PkrdNVf41pNDX19fRggJTmZgsK3yNiQz7rffobEzOw521CLS8NEAdD9/kWqXvk+padO0T8wAMZM3a/T09NDdXU1eTXVDLe2sOkr3yBm2ey9FOPzUfuPL1FReoLLly7hs+2Au4g7OjuJjb3Ivo52irq7WPvk7wQNlJGOViq/81ecee89WtvaECYmKJ/sS0NjIxmuSxxob2XjF782Z8U2tbg0TBQjbS3U/PSHvPNO8Appk4FQX1+P1+PBkZBAyYt/MCsImt74OeWn3uXixYtBSxHYts3Y2BgnTp6cGEeSuQrXrv0B2/jGx6j+h7/l3ePH6enpmbiSM239ZN+6bt3i6LGjiMNi23/4f4hJCu12ABU9egJW0XbibS5dujRvhTTw1zK5cYO2umsM3agPWOcdHaH9/HtcunRpwZomPp+Pc2fPcvP4mxgTuG132UXqaqrp6elZsNLawMAA1VVVdF04E8KnVNEWlTARkf8pItUiclVEXhWRoKXIRaRRRN4XkSsiciEafVHz8wwP0VVxlba2tpDuHDbGUFVVRXvp8YDlty6fo+F6Q8h3H4+OjdHe3MxA/bWAtttLj1JXWxtypbXr9fV0nDmBCbEok4qeaO2ZvAVsNsZsZWJ60P8yz7ZHjDHbwx0goyJjuLmR1va2sAaCdXV2MthYF7Csv6aCpuamkCut2bbNjcaGgDDxDg8x0NXFUBhFlsbGx+nrvsV4z62QX6OiIyphYox50z87H8AZJubLUUuQb3wMj9sdVj0Tn23j83hmtzNj2ULG3W680yq2+dzjeHzesIJNRBh3e/C5w6vQpiJvMc6ZfBH41RzrDPCmiFwUkefna0SnB40OR1w8MbHhDVF3WBaOmJjZ7cxYtpC42FicCR9cZnbExhHjcIYVbMYY4mJjdO6cJeC2w0RE3haR8iA/T03b5r8BXuBHczRzwBizE3gC+JqIHJ5jO4wxLxljdhljdi1bMWsqYnWbknLWkZ21Kqw/YFdGBsvX5QcsSykqYW3O2gVrmUyyLIvcdetJzvtghkBn0jKSXS6WJS18M+Gk+Lg4UleuJG7FypBfo6LjtsPEGPNRY8zmID+/ABCR54DfBD471wx9xphW/+9O4FVg9+32R92emKRluEq2smrVqpD2TiwRNm4MRN8gAAAN/UlEQVTcSNaMgkQrd+xm/Yb1Ie/hJMTHk5WTExAmIkLWgUfILygIKZQsy2JDXh6Zew7rJFxLQLSu5jwO/CHwcWPMyBzbJInI8snHwMeA8mj0R81v1eFH2blzJ3FxcfOGgWVZ5Obmkl1QxLLcvIB1zoREVj14gJ07dy4YBA6Hgz1797L6yGOIBG6bvu0B8ouKWTFPgaXJviQnJ7Nx4yZcu/aG8ClVtEVr0Nq3gDjgLf//nGeMMS+ISDbwXWPMk0Am8Kp/vRP4J2PMG1Hqj5pH4qrVFH3q89i2TWlpKf39/QEjYCf/qPPz8tm5bx+Fn/tK0NDJefwpxro7EZGgI2AtyyI2NpZ9e/dS9NEnWblzdgg44uIp/sJXMV4vp0+/R2tra8AI2Mn3dblcHDhwgKLnXtABa0tEVMLEGJM/x/JW4En/4+vAtmi8vwpf+padbE1cRrIrk87mJmqv1TA8PIw1eW9OXh6ukm2s+/in5rw3RxwOCv7dl0nMWk1OTg43Gm/QPP3enA0bJu7NeeRxMnYfnHMvKDEzm5IX/xNJ2WvoqquZqthmjCElOZnComJcG/JY/9SzJGTqDX9LhRZHUgGMMQy33ODWlfN4+vsQp5PE7DUTdw0HuR9nLj73OL3lV+ivq8HnHsOZkETapq2kFm4K+67hW5fOMOq/azh+hYuVD+whYdoUGSrytDiSumMiwrKcdXd885wjNo6VO/ewcueeO2onPt3Fmkd/647aUItDT4ErpSJCw0QFZYzBNz6G7Q1vVOusdnw+vGOjd3TvjDEGn3scn3v+GxHV3aWHOSrARHX6d+mpKMPr9YFA/PJkMvcexvXA3jCq05fRfuooQ+1t2Bgc1kR1+qwDR0gp2BhydfrOc6founAG99joxGjX+HhWbt9Nxt5DWp1+idEwUQDYXi8N//KPNF08S1VlJa1tbVPz5iQkJJBXdoW8/Ncp/NRzrNi0dc52RjvaqP7e39JYXzdRIa2/H5g4F7MyPZ3iyxdYW7SRoudenLMM5MTdw8e4/sYvqKmp4caNGwHz5qy9eJGi0mNT8+ZotfqlQcNEYWyb+ldepuzo21y+cnmiING0w4nR0VEqKyu5Xl+Px+1my+e/EnQirrFbnZR/5//l5PHjdM+oR2KMoevWLbpLS2lubsbnHqfkxT8IOhFX+6l3KH/1J5w8eRKfzxfQjsfjoaGhgRs3brB/aAjj85Hz2Mcj/F9E3Q49Z6LoLrtAzelSLl+5jG3bQc9L2LbN6Ogox48fp/an38c3PjZrm+s/+yGnT52aFSQz22lqauLy2TM0vf7qrPWjXe3U/eo1Tp44gcfjCV6xzRi8Ph+lpaXUH32DoebG8D+0ijgNk/ucMYa2k+/w/vtXF6xFYoCxsTEaGxvpvnoxYN1I2006Gupo7+hYsB3btqmrq6Pr/Ut4RwPvtug4c5Kqykq8/kOshdqpKK+g4713F9xWRZ+GyX1utP0mXc03GBwcDGl727apvXaNjtMnApZ3njtFTXV1yO9rjKGhoYHuq5c+aNvrpevyeZpbmkOu+tba1kpXRRnesdGQ31tFh4bJfW6su4vu7p6wTmIOj4zg7u8NWDba2U5vb2/Il25t26a7u5uxrvapZd7hQUZHR6ZO/IZqcHhoVn/U4tMwuc8Z22CMTTjXQ4wxsw5ljG0T7ggQ27Yx04LD2Da2bcK+OmMbA1oD9q7TMLnPxS5PJjk5OawgiIuLIyYhcLxJbEpqWEWNLMsieflyYlNXTC1zJi0jIT68imnGGJYlJIR135CKDg2T+9yy3PVkZq8mNja0icEty2LdunW4Hgi858a1ax+FRUUhV1ozxpCfn0/61p1Tyxyxcawo3ERmZmbI/U9PTyd9fX5YU5eq6NAwuc+J5SBz72GKQgwCy7IoKCjA9eCBgOXJeYVk5axl2bJlCx6mWJZFVlYWKws2EpcWWIIz68DDbNy4KeS+bNq4iayDjyy4rYo+DRNF5t5DbNrxAGtzcub9I3Y4HBw4cIDcI48Rm5wSsE7EYv1vP8uhg4eIn6di22SFtH379rH2yadnrV+Wm8f63fvYsX3HgpXWSkpK2PDAblIKN4X4SVU06QhYNVHd7ItfB4EVZVe4VlPD8MjI1B+zMQaXy8W2rdvI/+jjZB95PGg7KfnFbPzMF3A4HVwtK6OlpQXbTJxQNcbgdDpZv24dm7dupfi5F0nImF2TRETY8InPYWybpGVJlJeX09vbO9UX27ZJSUmhpKSEgl17yP/MF7X+6xKhxZHUFJ97nFuXz9Feeoxbba0MDQ3hsCzSUtNwFZeQdfAIyesLFmxntKudjtLjdJZdoOtWN16Pm7j4eFwrXWTuOUjGnkOz9mxmMsamr/J92k4dpauhnoGBAYwxJCcn48pdT9aBI6zYsh2xQi+0pEJ3O8WRohYmIvInwJeBLv+i/2qMeT3Ido8DfwU4mKgP++cLta1hEl3GGEY7WvEMDWE5HcSlu4hdPv8ffzC+8TFG2luxPW4ccfEkrlqD5Qx/Z3i85xbjfb1gDLEpqcSvzAi7DRWepVhp7S+NMf9rrpUi4gD+D/Ao0AKcF5HXjDGVUe6XmoeIkJi1+o7bccTFszx3wx23E7dC58W5F9ztg83dQJ0x5roxxg28Ajy1wGuUUktQtMPk6yJyVUReFpG0IOtXA83Tnrf4l82i04MqtbTdUZgsMEXot4E8YDvQBvxFsCaCLJtr9j+dHlSpJeyOzpkYYz4aynYi8h3g34KsagFypj1fA7TeSZ+UUndH1A5zRGT67EhPE3zqz/NAgYisF5FY4FngtWj1SSkVPdG8mvM/RGQ7E4ctjcBXAKZPEWqM8YrI14FfM3Fp+GVjTEUU+6SUipKohYkx5nNzLJ+aItT//HVg1vgTpdS95W5fGlZKfUhomCilIkLDRCkVERomSqmI0DBRSkWEholSKiI0TJRSEaFhopSKCA0TpVREaJgopSJCw0QpFREaJkqpiNAwUUpFhIaJUioiNEyUUhGhYaKUiggNE6VURGiYKKUiIiplG0XkJ0CR/2kq0GeM2R5ku0ZgEPAB3nCnI1RKLR1RCRNjzDOTj0XkL4D+eTY/Yoy5FY1+KKUWT1TnGhYRAT4FPBLN91FK3X3RPmdyCOgwxtTOsd4Ab4rIRRF5Psp9UUpF0W3vmYjI20BWkFX/zRjzC//jTwM/nqeZA8aYVhHJAN4SkWpjzIk53u954HmAFdlrbrfbSqkoue0wWWhqUBFxAr8DPDBPG63+350i8iqwGwgaJsaYl4CXAHK3bA86H7FS6u6J5mHOR4FqY0xLsJUikiQiyycfAx8j+BSiSql7QDTD5FlmHOKISLaITM7elwmcEpEy4BzwS2PMG1Hsj1IqiqI5PejvBVk2NTWoMeY6sC1a76+UWlw6AlYpFREaJkqpiNAwUUpFhIaJUioiNEyUUhGhYaKUiggNE6VURGiYKKUiQsNEKRURGiZKqYjQMFFKRYSGiVIqIjRMlFIRoWGilIoIDROlVERomCilIkLDRCkVERomSqmI0DBRSkXEHYWJiHxSRCpExBaRXTPW/RcRqRORGhF5bI7XrxeRsyJSKyI/EZHYO+mPUuruudM9k3Im5sYJmOtGRDYxUZ2+BHgc+FsRcQR5/X8H/tIYUwD0Al+6w/4ope6SOwoTY0yVMaYmyKqngFeMMePGmAagjokJtqb45yF+BPhn/6LvA799J/1RSt090ZrqYjVwZtrzFv+y6dKBPmOMd55tpkyfHhQY/+oG14dxwq6VwK273Yko+LB+LvjwfraicF+wYJiEOKfwrJcFWTZzSs9QtvlgxbTpQUXkgjFm11zb3qv0c917PqyfTUQuhPuaBcNkoTmF59AC5Ex7vgZonbHNLSBVRJz+vZNg2yil7hHRujT8GvCsiMSJyHqggIkpQKcYYwxwDPiEf9FzwFx7OkqpJe5OLw0/LSItwD7glyLyawBjTAXwU6ASeAP4mjHG53/N6yKS7W/iD4H/KCJ1TJxD+fsQ3/qlO+n3Eqaf697zYf1sYX8umdhBUEqpO6MjYJVSEaFhopSKiHsmTO506P69QkT+RERuisgV/8+Td7tPd0JEHvd/L3Ui8s273Z9IEZFGEXnf/x2FfRl1KRGRl0WkU0TKpy1bISJv+W91eUtE0hZq554JE+586P695C+NMdv9P6/f7c7cLv/38H+AJ4BNwKf939eHxRH/d3SvjzP5HhN/O9N9E3jHf6vLO/7n87pnwuROhu6ru2Y3UGeMuW6McQOvMPF9qSXEGHMC6Jmx+CkmbnGBEG91uWfCZB6rgeZpz+cdln+P+LqIXPXvfi64e7mEfRi/m0kGeFNELvpv9fiwyTTGtAH4f2cs9IJo3ZtzW6I4dH9Jme9zAt8G/oyJz/BnwF8AX1y83kXUPffdhOGAMaZVRDKAt0Sk2v8v/H1rSYVJFIfuLymhfk4R+Q7wb1HuTjTdc99NqIwxrf7fnSLyKhOHdB+mMOkQkVXGmDYRWQV0LvSCD8NhzoJD9+8l/i9u0tNMnHi+V50HCvxFsGKZOFH+2l3u0x0TkSQRWT75GPgY9/b3FMxrTNziAiHe6rKk9kzmIyJPA38DuJgYun/FGPOYMaZCRCaH7nuZNnT/HvU/RGQ7E4cDjcBX7m53bp8xxisiXwd+DTiAl/23WtzrMoFXJ0ry4AT+yRjzxt3t0u0TkR8DDwMr/bfH/DHw58BPReRLQBPwyQXb0eH0SqlI+DAc5iillgANE6VURGiYKKUiQsNEKRURGiZKqYjQMFFKRYSGiVIqIv5/9SYqmlkpMCMAAAAASUVORK5CYII=\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "assy=Assembly(6,3)\n", "assy.set_pitch(1.62)\n", "assy.set_source(uo2)\n", "assy.set_coolant(h2o)\n", "assy.set_pins(fuel)\n", "\n", "fuelmap= [['1', '1', '1'], \n", " ['1', '1', '1'], \n", " ['1', '1', '1'],\n", " ['1', '1', '1'],\n", " ['1', '1', '1'],\n", " ['1', '1', '1']]\n", "\n", "assy.set_fuelmap(fuelmap)\n", "\n", "det=Detector('D')\n", "det.set_location(Point(174.726, 174.726))\n", "\n", "ex5=Experiment()\n", "ex5.set_assembly(assy)\n", "ex5.set_detectors(det)\n", "ex5.set_materials(uo2,he,zr,h2o)\n", "uo2.set_color(\"#0f0e0e\")\n", "he.set_color(\"#474545\")\n", "zr.set_color(\"#635f5f\")\n", "h2o.set_color(\"#96deeb\")\n", "ex5.Plot(xl=[-10,10],yl=[-10,10],detectorSize=0.4)" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "## 6x4" ] }, { "cell_type": "code", "execution_count": 6, "metadata": {}, "outputs": [ { "name": "stdout", "output_type": "stream", "text": [ "Warning: no pool in the problem, the surrounding of the Assembly is filled with coolant material\n", "No absorbers in the problem\n", "Warning: elines missing; only distance travelled in various materials will be computed\n" ] }, { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAARMAAAD8CAYAAABUzEBbAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDMuMC4zLCBodHRwOi8vbWF0cGxvdGxpYi5vcmcvnQurowAAIABJREFUeJzt3Xd8XNd95/3PuVPQeyUBEADR2ZtYRFKiqm3ZsWKvbMlVcnckbXaTzSvx7j6vJK/knzz7PHmym9hxVk4Up9ixE2/kKLYsWxIr2DtI9EKAKEQjehlMuef5YwAQgxkAM8AMC/h7v154AXPvnYMzA86Xt5z7O0prjRBCrJRxrzsghFgdJEyEEGEhYSKECAsJEyFEWEiYCCHCQsJECBEWYQkTpdSbSqlepdT1OctSlVLvKaUap7+nLPDcl6e3aVRKvRyO/ggh7r5w7Zl8H/jwvGXfAj7QWpcAH0w/9qGUSgX+ANgD7Ab+YKHQEULc38ISJlrr48DAvMXPA383/fPfAb8e4KkfAt7TWg9orQeB9/APJSHEA8AawbaztNa3ALTWt5RSmQG2yQHa5zzumF7mRyn1deDrAFGxsTuz1peEubtCiBk3r1/t11pnhPKcSIZJMFSAZQHH92ut3wDeAMjfvE1/6+33I9kvIR5qr67PaAv1OZG8mtOjlFoDMP29N8A2HUDenMe5QFcE+ySEiJBIhsnbwMzVmZeBfwuwzS+BZ5VSKdMnXp+dXiaEeMCE69LwPwGngTKlVIdS6ivAnwDPKKUagWemH6OU2qWU+msArfUA8MfA+emvP5peJoR4wITlnInW+jMLrHoqwLYXgK/Oefwm8GY4+iGEuHdkBKwQIiwkTIQQYSFhIoQICwkTIURYSJgIIcJCwkQIERYSJkKIsJAwEUKEhYSJECIsJEyEEGEhYSKECAsJEyFEWEiYCCHCQsJECBEWEiZCiLCQMBFChIWEiRAiLCRMhBBhEdEwUUqVKaWuzPkaUUr953nbHFJKDc/Z5vcj2SchRGREdN4crXU9sA1AKWUBOoG3Amx6Qmv9sUj2RQgRWXfzMOcpoFlrHfLkPkKI+9/dDJOXgH9aYN0+pdRVpdQvlFIb72KfhBBhclfCRCllBz4O/EuA1ZeAfK31VuAvgJ8u0MbXlVIXlFIXxgZuR66zQohluVt7Jh8BLmmte+av0FqPaK3Hpn9+B7AppdIDbPeG1nqX1npXfGpa5HsshAjJ3QqTz7DAIY5SKlsppaZ/3j3dJ9n1EOIBE9GrOQBKqVi804N+Y86ybwJorf8KeAH4DaWUG5gEXtJa60j3SwgRXhEPE631BJA2b9lfzfn528C3I90PIURkyQhYIURYSJgIIcJCwkQIERYSJkKIsJAwEUKEhYSJECIsJEyEEGEhYSKECAsJEyFEWEiYCCHCQsJECBEWEiZCiLCQMBFChIWEiRAiLCRMhBBhIWEihAgLCRMhRFhImAghwkLCRAgRFhEPE6VUq1Lq2vQ8whcCrFdKqT9XSjUppaqUUjsi3SchRPhFvKD0tCe01v0LrPsIUDL9tQf47vR3IcQD5H44zHke+HvtdQZIVkqtudedEkKE5m6EiQZ+pZS6qJT6eoD1OUD7nMcd08t8yPSgQtzf7kaY7Nda78B7OPOaUuqxeetVgOf4TcIl04MKcX+LeJhorbumv/cCbwG7523SAeTNeZwLdEW6Xw8z0+3GOTKEc2QI0+1edjsexyRTQwO4xsdY7iSM2jRxjY0wNTSAxzm17L6YLifO4UFcoyNo07PsdsTyRfQErFIqDjC01qPTPz8L/NG8zd4GXldK/QjviddhrfWtSPbrYTVxq5Oe00fpv3aFyclJAGJiYkjfvI2sfYeIXeN3dOnHdLsZqq2iu/IIw503cbrcWK1W4lNSyXr0cdK378YaE7tkO86RIXrPnaT33EnGx0YxPR6io6JILiole/8TJBaVotTi/9dprRltbaLn5FEGGmqYcrkARWxcLFm7HiVzzwGiUmQv9m6J9NWcLOCt6XnJrcAPtdbvzptr+B3gOaAJmAC+FOE+PXS0Nml/921ajr1HbU0NnZ2dmKYJgGEY5Fy8QMW5U6x//GnyPvz8gh9i5/AQdX/7HVrra6mrq2NoaGh2XWxMDEVVVygpLaXsC98gIX/9gv25XXWRhp/8gPr6OlpbW3E6nbPrMjMyKL98gYIt2yn53FexREUHbMN0uWj+8fdpvXSe2ppqenp7Z/eObFYr+VeuUlZ5mJLnP03mI/tDfs9E6CIaJlrrFmBrgOVz5xrWwGuR7MfDrv0XP+Xaz3/K6TNnME3T55DE4/HQ3t5OZ2cn+0ZHQWvWPfdJvzZc42PUvPFnnDp+zCeMZoxPTHD9+nXaWlvxuN1s+tp/Ii5nnV87A9VXqP7R9zly+DBOp9OvnZ7eXvr6++nv7UV73JR96XUMq+8/U216aPrh33D5yHtUVVX5teF0uWhqbqLtZhtOpxNlGGTs3Bfy+yZCcz9cGhYRNN7RRsux9zl95gwejyfguQ2tNR6Ph9NnztBy7H3GOtr8tun45b9x+dy5gEEywzRNRkdHOXH8OM3//Pd+v8sz5aDlX3/I0SNHcDgci7ZTW1dH7bkz9J6r9Ft/++pF6s+eDBgkc1+T0+nk2NGjNL/9L7jGRwNuJ8JHwmSV6z51lOrq6gU/dHOZpkl1dTU9p476LHc7JumtukRzS/OS7ZhaMzA4SM/NNsbaWnzW3a66xI0bN5iaWvpEq2ma1FRX03PqGFrf+Z1aa25VHg76Nbncbpqbm+m/eHbJbcXKSJisYp4pB/3VVXR1dQV1tUVrTVdXF/3VVXimHLPLB69fpvXGjaCv2Gitqaur9dur6D17gsbGxqBCAGDS4aC7s52xtht3lvXcor+9jZGRkaDaME2T5uZmes767+GI8JIwWcVcoyOMT0yE/LzxiQlco3c+rI7bfQwODgYdAgDDIyM4+vt8ljmHB5kIsT+DtweYGhqYfTw1eJvBwUEMI/h/ug6HA9f46LIvX4vgSJisYlprzGV8gDTa54OnTb3MD6Lvc5bThNba94l6eX3xPkfCJJIkTFYxW0Ii8bExIT8vLiYGW0Li7OOolDSSkpND2huIj4/Hnpzqs8wen0B0dOBLvQtJSknGnph8p42kZJKSk9Eh7CXZ7XZssXFLjlsRKyPv7ipmjY4hrXQDWZmZQT8nKzOTtNINWKPvhFDq5m2sX7/wuJH5DMOgvLyczEce9VmesXs/ReuLgg4lu93OmpxcEgqLZpfFrs0jY20usXFxQfeloKCATLk0HHESJqtc1v5DbNiwMagPsGEYVGzYSNajh3yW2+ISSK/YzLq8vCXbUUoRHx9Pdm4eicVlPuvSt+9mfdF6rNalhzfNBFLWngMow+LTfvb+Q1RUVAT9mkpLS8mYF2wi/CRMVrmEgmIK9x1g546di374DMNg546drN93gITCYr/1eR9+nkf27CU9LW3BdgzDICYmhscee5z1n/yc32GFNSaWwuc+wWOPPYbNZmN6ZHTAdgoKCtiwYxdZ+x73W5+xcy9lOx+htLhkwb4opbBarRw8cICCpz+KPSk54HYifO5WcSRxjyilKPzkZ0FrEhITqL5+nd6+vtkPstaazIwMNm7axPq9Byj85GcDfsijklOp+OpvoiwW6quraWpqZGJyEsMw0FpjsVgoyC9g46aNlL30JRKLSgP2J3P3ATa6Xdhsdqqrr9PZ2YnWGoV3jEpiYiLl5eUUb95K2SuvBhxOb1htlL38GxiGhdT0dGpraxgaGprtC8DatWvZtGkTxc9+jDWPPR2+N1QsSD2Il8vyN2/T33r7/XvdjQeK1pqRpjq6Tx7ldlM9o2NjACTEx5NWXEb2/kMkFpcvuLcww+2YpP/SWbpPHWWwtxfHlAOrxUpKaipZj+wjc3dwN9dNdHfSc+oovVWXGBkZxePxEBMTQ2pOLmsOPkXqpm0YVtuibZgeN0M13psOb7e3Mj4+gWExSIxPIGPTNrIffZy43Pxg3yIxx6vrMy5qrXeF8hwJk4eQa2wU54j3Jj17YjK2+ISQ29DaZOp2P27HJBabnai09CU//IF4phxMDd5Gu91Y4+KXfZevc3gQ19goyrBgT04J6s5lsbDlhIkc5jyEbPEJywqQuZQyiE4P/irRQixR0cRmL136YCn2pBTsSSkrbkcsn5yAFUKEheyZPGQ8Uw6Gaq/NDlGPSk4luWLzgnVDFuLo72W4oQa3YxLDZic+r4D4/PVLnnOZS5sehhvrmOy5hel2YYtPILl8M/bEpJD64hofZajmGs6xEQyLlej0DJLKNmJY5J/33STv9kPCPTFOx/s/p+/KBdpaWxkaGgQgOTmF/IICMrbtIvfpj2KNXXww2OiNJm87zY20trUy5XBgtVrJzcsjK28dax5/lvTtuxcNFe3xcOvkYXpOH6ervZ3u7luYHg+xcXEUFhSSUb6R3Gc/Rkzm4pMUOG730fHez+ivvcaNGzcYGxvDUIqMzCxvf/YcYM3jzyzrXI4InYTJQ8A5PEjtX/85l8+dpbnZW0ZgbqW1q1VXKaq+zvbGWiq++psLnnvov3yOup/8I2fPnOH2wABKKUzTRClFY1MT8XFx7GprpbS9lfyPfzpgoJhuF43/8D1qTp+gqqoKp9OJnr7fxjAMampqWHPpInuaGyh/+RskFPiPeQEY77xJ7d9+h/Nnz9LR2elte/o1tdy4waVLF9nU1MiG5nrvJWZ7VDjeSrEICZNVznS5qPvb73Dy6JGAhY1mgqWxsZGJiQmUxcKm134Pw+b7v/lISwN1//KPvP/B+z4BAMx+Hx0b49jx43hME3tSCmsPPevXnxv/54dcPX6Yq1evBuwLwK3ubt5//z1QsPm13yU6LcNnO+fIEHV/+5cc/uADhoeHA7YzNTXF5StXcExNYYmKpvSL3wzpEEyELmInYJVSeUqpI0qpWqVUtVLqPwXY5pBSanh66tArSqnfj1R/HlYD1y/TXFO9aIU08H4AOzs7aa6pZuD6Zb/1Hb/8d06fOT0bJIu1c/r0aTqPv+9TEwVgsq+HzsvnAwbJ/DbGxse5ePEit46/57e+++RRrl69EjBI5rdTW1tL5/WrTHTeXHA7ER6RvJrjBv6L1roC2It3zpwNAbY7obXeNv01v3K9WKFblYepr68PutJafX09tyoP+yyf6O6ir7WZwcHBoG7/93g8tLbe4HbVJZ/lvWdPUF9fF1QbpmnS0dFBX9Ul3I7JO8tdLnovnqa1tTWo16S1pq62lu7Tx5bcVqxMxMJEa31La31p+udRoJYAM/WJyHGODjPc3RV0VTKAkZERhru7cI4Ozy4bqr9Oc0vLIs/yZZomN1paGKi66LN8sOba7PD5YHV0dDLa0jD7eKz9Bj09PXg8wc2No7XmVnc3Q/XVQf9OsTx3ZZyJUqoA2A4EKsS5Tyl1VSn1C6XUxkXakOlBQ+SZnGDS4QjpXIFSCsfUFJ7JOxXR3GNjOCYnQwqBKacT96RvVTXPlANXCJN+aa1xOCZxT97ZM3FPTjAxPhHSa3K73XiczpBqoIjQRTxMlFLxwP8B/rPWev5/kZeAfK31VuAvgJ8u1I5MDxo6w2bHZrWFFAJaa6wWK4bNfqedqKigygbMZbFY/E7iKqslpAJLABbrvL7Y7FhtoV3qVUqhLBZUiL9bhCai765SyoY3SH6gtf7X+eu11iNa67Hpn98BbEqp9Ej26WFiT0wmMSkJu92+9MYzz7Hbvc+ZU90sPjefnNyla5nMMAyD7KxsEvKLfJbHrc0jNTV1gWf501qzJiubuLW5s8tis3PIysgIKSCTk5OJy1p8zIpYuUhezVHA3wC1Wuv/b4Ftsqe3Qym1e7o/cgwTJspiIXP3fgoKCoIuJFRQUEDm7v0oy52CREmlFazNyws6lLTWlJSWkLnbdya97EcPUV5eHnQopaamklFc4nMPkD0xiYzy4KvHGYZBWVk52fufDGp7sXyR3DPZD3wBeHLOpd/nlFLfnJkeFHgBuK6Uugr8OfCSfhBvY76PZe7eT0VFBdHR0YueZ1BKER0dTUXFBr/pNJVhIXv/IbZv375kEBiGQUlxMVkVm4hK9d3JTCwuY23hejIzMpZsx2KxsGPHTtY89ozfuuyDT7Nt2zZsSxx6KaVISUlhXWEhyRWbF91WrFzEBq1prSuBRc+Saa2/DXw7Un0Q3kOdspe+hNvt4dixo4yPj/tdUjUMg7i4OB577HHKXnolYFWy7EefoOJmKx6Ph4sXL3or389pZyaoitavZ+ejB1j/whf82lDKoOyL38Q9Ps6pk5Xc6u4G8DlkMQwDq8XC/gMHKH32oySX+o8mSMhfT8nH/gNuj4fjx4/jcrn8XpNSioz0dA4cPEjZy9/0m2JUhJ+8ww+B5PJNbP7SbxCbkEBXezv1dbU+xZHKyr2HMUUvfJ7EorKAbSjDoHg6aNasWUNzU/PsXL4Wi4Xs7GzKSsvIKt9A0adfXrCeSHRaBhtf/R2iklPoudFCXX0dt2/fRmtNVFQURUVFFBQUkvfUR/xq0c615sCTWGNiSExM4ubNNpoaG2evXKUkJ1NeUUHWugJKPvvlJe/xEeEhxZEeItr0MNxQS8/ZEzgHvXcN21NSydpzkKTSCp/CzYtxjY/Sf/EsA9cv45lyYFhtJBQWk7X3saBrnGitGe+8Sc/pY0zc6kR7PFhj40jfsYe0rTuDvpfGdDkZuH6FvguncY2PogyDmIxssh59nPh1od3FLO6QSmtCiLBYTpjIhXchRFhImAghwkJOwD5EZs4v9J49gXPEe++NPTGJzD0HvdXgbcGNI5ka6Kf3XCWDNdcwnVMoq5X4dYVkPXqIuJx1QZ2n0NpkpKmenlPHmOy9hTY9WGJiSd+2m/Sde7DFBVej1j05Qf/lc/RfPIN7YhxlGESlZZD96KGQzgOJlZMweUgMVF+l5V9/SFtbK02NjYxPeO+biYuNpfjCOfLzC1j/yc+SunHrgm2Ybhetb/2IzsvnaGxopKOjHZfbjWEYZGZmUn7uNNlFxZR8/mvYExYuvTjR3UXjD75HV1srdXV1DA0OYk5fzSk8e5ai4nfIPfgUOU89t2Awaa3prjzMzQ/eoaWlhZaWFqYcDlCKpKQkyi9dIGfdOoo/+xXiZbqLu0LC5CHQf/UCNT/+O44dOYJjaspnTMbQ8DCXLl2iprqaxycn2PDiy6Rv9T/vpj0eGv7+f1N14hhXr14B8Gmns7OTzs5O1jU24BwdYdOrv4MtPtGvnYnuTqrf+J+cOH6c2wMDPm243W5qamuoq69j7+3buCfGyf+1TwUMlM73f07Nz39K5clKPB6PTzv9/f2cOnWSpOrrTE2Ms/Er/5H4vIJlvXcieBImq9zU4G0af/IDjhw+jMPhCLiNaZpMTE5y5PBhbBYrCesK/eav6Tr2K2pOV3L16pWAdURmrgq2t7ejTp0kJj2T8i+/7ruN6aHhH97g+LFjDAwOBmxnpvLbqVOnsNvtJBWVkTJvb2mkuYHGX/2M4yeOL1iKwDRNhoaGOHrkCFabne2/+0d+Nx6K8JITsKtc79lKamtrmZqaWnLbqakpamtr6T1b6bPc9LjpOeOt2bpUQSLTNGlra6O3oQ5Hf6/PuuGGWjpaWxcMkvntXL58mVsn/IcA3DrxPleuXl2yponWmpHRUdpabzBQfWXRbcXKSZisYqbbTc+FU7S2tQZ1l63Wmta2VnounMKcU3dkqO46He3tuFyuoH93fX0dPWdO+CzrPnmE+vq6oCqkAQwNDdHb2sJkz63ZZVNDA/Q31tHX1xdUG6Zp0tDQQM/JI0H3XSyPhMkq5hoZYnR0FHcIBYncbjejo6O4pqcPBZjo6uBW1+I1ZOfSWtPX18d4R5vP8omeWwwNDS3wLH+GUvR0dzPR0zW7bLLnFr39/SHVRRkZGWGirzeksgUidBImq5jpduPxeEKutObxeHz2TEyXE0+IVco8punTBnjPmQQbSAAabz1Zc84ekXa78YQQjuANN216QMIkoiRMVjFrbCzR9qiQK61F26Owxt65Uc8an0hsTAxGCKEUHR2NNc53Qi9rVAy2EE6CKqWIiY31mRjMGhtHbGxcSK/JarVijYqSSmsRJu/uKmaLTyQ1L5/kZP+SAgtJTk4mNS/f57Ju6oYtrF9fRLAfX8MwKCoqIn3bIz7LUzdvY13euqAPUZRS5ObkklhUOrssbl0B2WvWBB1KSily1q4lZcOWIHsvlkvCZJXLPvBk0NXNDMOgvLyc7AO+Vcmi0zPJKCohPS0tqEMmm9XKunX5pG7a7rM8c/cBSkpLguq3YRjkr1tHxvZHfO4gNixWMnfvZ31hYVCvSSlFeXk5WfseD+r3iuWTMFnlkjdsZv3GzRQuUbpxpmTj+o2bSd7gX5Us70MfZ8/efUtWbLNMFzbKffo5v3EdUanp5O3ez86dO5fsS1JSEtt27GDNwaf81mftfYxN27aTnpa2ZDtbt24ld8cjxGbLLCuRJoPWVjnDYqXslVfRHg9xceeob2jwXuKdOeegFDabjbLSUjbt2k3ZK69iWPz/WcSvK6Tis1/GMBTnz5+np6fHZ73WmuTkZHbt2kXpsx8la+9jAfuT/2ufwj0xjt1u5+rVq4yPj8+G08z3dXl57Ni1i4ovvUZUsn8Balt8AhVffg1Mk8sXL9La5p2QS2s920ZMdDSbN2+hfN/+gFXfRPhFvJ6JUurDwP8CLMBfa63/ZN76KODvgZ14i0m/qLVuXaxNqWcSOs+Ug1snPqD33Ek6OtoZGhwEIDklhdzcPDJ372fNwaewREUv2s54VwddR35BX2017R3tTE5OYrPZyMnNJbNgPWsef5aUJeqtam3Sd+E03ZWH6Wm/OTupVlx8PHm5uWRu3cnaQx/yG4U7n3NkiK6jv6LvygXa29sZHR3BYhikZ2SyJm8d2fsPkbnngNzstwz3XXEkpZQFaACeATqA88BntNY1c7Z5Fdiitf6mUuol4BNa6xcXa1fCZPlMj5uRpnqcQ94wsSenkFhcFnBvZDHOkWFGmuvxOCYxbHbi8vKJzVobUhtaa8Zu3sDR243pcWONjSeptAJrdExI7XicUww31OAaG0VZLESnZ5JQUIRSchS/XMsJk0gf5uwGmrTWLQBKqR8BzwM1c7Z5HvjD6Z9/AnxbKaWkSn1kGBYryWULTpwYNHtiEunbd6+oDaUUCfnrSchfv6J2LPYov5O94u6LdHTnAO1zHnfgP9/w7DZaazcwDPjt38r0oELc3yK9ZxLotP/8PY5gtkFr/QbwBngPc1betYeX43YfzqHpgtLJqUSnZYTchulxM9HVMXuYE5O9NuTDEwDX6AiT/T1otxtrXAKxa9aGfHiitWayp8t7mGNYiE7LCDhdh4isSIdJB5A353Eu0LXANh1KKSuQBAxEuF8PHW2aDFy/QvfJw/S13mBkxDvtc2JiIhkFhWTvf5LUTduWHCXqGh2h58wJes5X0tfbi2PSgd1uIy0tjcwtO8ja/0RQ505GW5vorjxMX10NA4MDsydg07PXkP3oITJ27l3yZLDHOUX/5XN0nzzK7a5ORsfHMJRBakoK6SXe8TJJxYGn7hDhF+kwOQ+UKKUKgU7gJeCz87Z5G3gZOI13hr/Dcr4kvEyXk6Z/epPGc2eorr7O0NDQ7PgM0zRJTk5mY/V1Sh7ZS/Fnv7xg+cbxzpvUff+7XKu6SmtrKy6Xy2fMSc6lS2y5dJ6yFz7nN/p1htaazvd/TtP771B19So9vb2zbWitiY6OpvT6Ncq3nqB8gUvDAK6xEer+9i+pv3KJ+vp6xsfHfV5TRsY5tly/QvGhZ1n30U/Iydi7IKJhorV2K6VeB36J99Lwm1rraqXUHwEXtNZv452P+B+UUk1490heimSfHjZamzT9+PtcPvyeTz2SubVABgcHOXXqJKOjI2AxKPncV/0+fI7+Xmre/A5HDn/A8PBwwBv22tvbZy/zbomOIbl8k9823Sc+oObnP+XEieO4A9QjmZiYoKqqitv9/WBqNr76O34TenmcU9S9+R1OH/mAtra2gK+pt7d3tiCUYbOR96GPh/CuieWI+KA1rfU7wDvzlv3+nJ8dwKci3Y+H1UhzA60Xzi1Z2Mg0Ta5du0ZqaipZex8jqbjcZ337u//GuTNnFgwS8O5ZTE1Ncfz4MRKSk9k2r6Cza3yUtg/eofJkZcAgmduXzq4url26QPrJo+Q8/ZzP+r4Lp6m7fMknSALxmCanT58mKTGRzN37lxy3IlZG9v1Wue6TR6iuqQ7q1n/TNKmprqZ7XiEh58gwfXXVdN3qCqqdyclJutrbGW6o9Vned/EMTU1NS1ZIm+lLU1MTPecqMT13Sg5obdJ96ih1QRZZ0lrT0NhI77nKJbcVKyNhsoq5JycYaKoPuioZQG9fHwNN9bgnxmeXDV6/TMuNlqDbME2Thvo6es+f8lnef/EMN27cCLqmicvloqurk9GWxtll4x036evqXLCebaC+tLW10X/pXND9F8sjYbKKucZGGZ90hFwcaXzSgWt8dHbZ1PAQIyMjIRU2Gp+Y8KnWBuAaHw+qFu1cI8MjuEZH7rQxOsLIyEhIr8npdOKanEDr0Ao8idBImKxiyjBCKmg0w2Ion3Mdy2lHAcy7zLycScQNQ/lcrlaGEXKRI6Vm2pBJzCNJwmQVsyUkEh8fF/KeSVxcHLaEO8WRYjKySEvPCKmoUXJKCjEZ2T7Lo1JSiY+PD7ovWmvS0tKImjOoLiotndTU1JD2kmJjYrAnJi0rzETwJExWMYs9iozN28nJyQnqg6SUIicnh4zN230KEqVs3ErBEvVQ5qsoryBz70GfZVn7HqekpCToduLj4sjKW0fcnBn5YjKyycgvDLp6nGEYFJeULFgSQYSPhMkql/XoITZu3ITFsvRt+BaLhQ0bN5K175Dv8qhoMrY/Qllp6ZJBYChFRkYGGfkFxOWs81mXumkb+QUFxMbELBluhmGwafNmsvc/6bftmgNPsXnz5qBCKSoqisLC9Su+KVEsTcJklYvNzqHomec4eOAAVqs14AfQMAysVisHDxyg+JmPErvGvypZ3jMfY8vuvRQuUi7RMAxSUlI4cPBNJqgQAAAWwElEQVQgxS++4hcChs1OyYuv8PihJ4iNjQ3YjlIKwzDYsmULZXv3k7Frr982KZu2UvLoY+zc4a3YFiiYDMMgJjqaQ088QckLn/cb+CbCTyqtPQTWPvFhlMVCXHw89fX1tLW14XQ6AbDb7eTn51NWVkbh0x9lzePPBGzDEh1Dxdd+E4s9itz6Gupqa+nt65utEp+UlERpaSkFRcWUv/wbxGRmB2wnqaSCTS9/g6goO83NzT6TqFssFtauWUN5RQV52x+h6FNfDFjYSCmDok99EcNmIzUtlbraWjo7O2cHwsXExFC0voii4mJKXvg8qZu2rfg9FEuLeKW1SJDiSMszNXib3nOV9F08w9TkJABRMTFk7NxL5u4DQY0Q1dpkpNk7Q97wjSbcbg+GUsRmZrPm4JOkbtwW1Jy+7skJ+i+fo+fUMRwjQ2itsVqtpGzYQva+x4nNWRfUeZ6JW530nDnG7WuXcblcGEoRFZ9A1t7HSN+5B1tcwtJvjPBz31VaixQJk5XRWqPd3omtlNW27Ksc2jQxnVMYNjsqiHMyCzHdbrTHjWGPWn5ftMZ0TqEsFgyrTFC+UvdjpTVxH1JKoRa4MzikdgwDyzJqmMxnWK1gXdk/RaXUkiULRGTJCVghRFjInslDZrKvh/5LZ3D09wLeCbbSd+wlJiMr6Da06WGovobBmirc42NYoqJIWF9C2padIe0duMZG6b90lvHOm5huF/bEZNK27iQ+vyjowx2tNeMdbfRfOY9zeBDDsBC7Jpf0nXuxJyYF3RexchImD4nJvh5af/oj+pobaWioZ2h4GIDkpCRKSt8js6iEgl9/0W/U6nx950/R/sE7dLW309LSzJTDgdVmIzc3l4KCQrIe2Ufus7+2aLV79+QEbf/+L/Rdv0pTcxO9PT14TJO42FhKSkvJyltH/kf/A0mlGxbty+iNJlrf/md6brbR2NjA+NgYyjBIT0+npPgXZFRsouD5F7HFy0nYu0HC5CEw3tVOzd98m9MnK+menjxr5sT7wMAALTdukJ11mX3dXWz4yuvErc0L2E77r96m7t1/5/Tp00xNTfkMae/v76eqqoqtba1M3Oqk9Ivf9J4Lmcc1PkbtG/+TS2dO0djovRt4pp2BgQE6OjtJSEhgf28vGz79hQUHmw3UVFH7wzc5ebKSoaEhv9fU0NBAYU01Y13tbPzGb2FPlJqwkSZhssq5Jyeo+/53OXrkMENDQ8y/ejfzuKenh6OHP8CwWNj6W/+X3yCv/qvnafjlzzh27FjAeiSmaWKaJpcvX0ZrTVRqGoW//hm/39X4j9/jfOVxWhYoRaC1ZnRkhMOHP8AwFNsysoifM5weYLL3Fg0//j4ffPA+DofD7zXNtHujtRWX2401OoZNr/+ulG6MMHl3V7n+y+doqKtlaGho8UprWjM0PExDXS39l31rf2it6Tr8LmfOnl2ysJFpmly5coWeS+d9yhiAt4ZsV0MdzS0ti/ZF4y0bcOniRW4de89v/a0Th7l8+XLAIJnfl46ODjqbGhhpbli032LlIhImSqn/RylVp5SqUkq9pZQKuI+plGpVSl1TSl1RSl2IRF8eZjNVyZqbmoKutNbc1ET3qaM+tT/G2lroab/J+Pj4wk+ep7m5mf4LZ3yW9Zw6Sn1dXZB91/T09tJXdx3n6PDscvfkBH3XLtHZ2blokMwwTZO62lq/6nEi/CK1Z/IesElrvQXv9KD/dZFtn9Babwt1gIxYmnN4iOH+fiamR7sGY2JykuH+fpzDdwobjbQ00Nbaemey8yWYpsnN9psM1Vf7LB9paaCntyeoEJjR2dnFWOudKm/jnTfp6ekN+vkAff39jLYFXylOLE9EwkRr/avp2fkAzuCdL0fcZaZzCqfLGXI9E6fLiem8UxHNMzmB0+n0nxltEW63G4/Tt6qa6XIHVf91Lqdzyqcdz5TDW60thEAyTdM7yjaEGigidHfjnMmXgV8ssE4Dv1JKXVRKfX2xRmR60NAZ9ijsNntIewJaa+w2O8aceiaWmFjsdntIdcqsVqtPTRQAw2YNqhTCXHZ7lE87lqhooqKiIISANAwDw2oNuUKbCM2y312l1PtKqesBvp6fs81/B9zADxZoZr/WegfwEeA1pdSCFWy01m9orXdprXfFp8qUBcGwJyWTlJ5ObEzwQ95jY2JITk/3mV4zcX0p+QUFQX+ADcNg3bp1fhOkJ64vJSsrK6Q9pZyctcQX3JnYPC5nHVlZmUE/HyAjPX3Fk6OLpS07TLTWT2utNwX4+jcApdTLwMeAzy00Q5/Wumv6ey/wFiAVbMJIKYPsRw9RVFwcVCEhwzAoKi4m69FDPpdR4/PXk5W3jvi4uKB/d1FREenzapFkPXqIivKKIPuuyMrMJKN8E/aEOyNZrTGxZGzeEXT1OMMwKK+oIHv/E0H3XSxPpK7mfBj4PeDjWuuJBbaJU0olzPwMPAtcj0R/Hmbp23dTWl5BcnLyooFiGAbJycmUllf4DRRTSrH2yY+wZ8+eJQ9TDMNg+7btZG/f7Xf7f1zOOrJLyylaX7RoXxTeOis7du4KWF9lzcGn2L59OzHR0YsGimEY5OXmklNcSmJR6aL9FisXqUFr3waigPem/9hntNbfVEqtBf5aa/0ckAW8Nb3eCvxQa/1uhPrz0LLGxFL+yquYHk/AEbAzH8aszEz27T9A+SuvBqxKlr51F2UfuoWpdcARsIZhYBgGW7duY/OBx1j3ay/4taGUovTzX8M9MYbVaqFh3gjYmW0SExLYf+AA5S98zm/AGkBMZjZlL72Cx2MGHAE7E1SFhYXs2ruPspd/Qwas3QURCROtdfECy7uA56Z/bgG2RuL3C19xa3PZ9M3fJn5NDn3NDTQ0NDA8/QFMSk6mtLSUjKISCp5/adEb/nKf/TWiUtJIz8r23pvT3IRjagqb1UpObh6FhQVk7nqU3Gc/tuC9OdbYODZ847eJW7uOsuorNDc109PTjWmaxM69N+djL5BUsvAhUUrFFrZ89XUSM7Pou9lGfUM942NjGNP35hSXlJBRsZmC5z8tBZLuEimO9JBx9PfSd/EMjtveWf6i0zLI2LmX6PTgT2pq08NwQy0D1VdxT4xhiYq+c9fwvCs4i3GNj9J/6RzjnTfRbje2hETStu0ift36kO4anui86b1reGgQZbUSm73We9dwgtw1vFxSaU0IERbLCRM5kBRChIWEyUNIa41nemTpSvZMtWnidkyiQxzVOp/pduFxTK5oLmCtNZ4pB6bLtaK+iOWTEgQPkamBfnrOVtJ/6SxTjunq9NExpO/YQ9aeA0Slpi/ZhtYmI031dE9Xp/eYprc6fUYWaw48SeqmbRhB1Jd1T07Qf+ks3aePMTUyjNZgsVpILd9M1v5DxAVdnb6D7lPHGKi+gtvlAqWIjk8gc89BMnbtlZOvd5GEyUNAa82tY+/R8t7PaJiZN2f6f3C7zUb+5UuUVh5m/TMfY83jzyz4IXaNj9Lwd3/Fzbpa6upq6Zs3b07J1csUFpdQ/vI3iV2z8O1YQ/XVNP74+zQ1NtHU3MTEzLw5hsGa8+epuHB6dt6chabNMD1uWt/6EW3nT83OmzNz3090dDRFV69QfPhdSl74PGmbty/7vRPBkzB5CHQdeZfqt3/CicrK2SJGM5wuF80tLdxobeXA+Dja9JDz5Ef82nA7Jqn93p9z7vhRbrS2+pU0GB4e5tKlS9xoacHjcrLpG79FTOYav3aGG2u5/vf/m6NHjjDpcPi04zFNOjs76erqYktfH6bTSekXv+43EZfWJi3/8g9Uvf9LLl666J26Y87hmsPhoLa2lubmZqamHGz6/FdJ2ySBEmkSJqvcxK1Omt97hxOVlbjd7oDbzARMZWUl0dHRpFRs8ZsitONX/87Vc2cCBsncdgaHhqg8cYLopBQ2/cdv+ezlmC4njT/+PseOHmVicjLg+ZqZYKiqqvL2ZcMWMnfv99lm8PoVGk8d5+Kli4v2xeFwcOzYMaKjY0gqKpMpQiNMTsCucj2nj1JdfT2oW/89Hg81NTX0nD7qu3zKQd+VCzQ0NCxZZMk0Tfr6+ui72cZ4502fdQPXr9DW2rpgkMxv5/r163SfPOy37a3Kw1y7di2ogk9TU1PcuNHiVz1OhJ+EySrmcU7Rd+1y0FXJtNZ0dnbSd+2yTw2RweqrtC6yRxJIbV0tvWdO+CzrOX2MxsbGoNsZHx/3VnjraJtdNtnXTV/bjdkh9EsxTZOmxkZ6Th8Puu9ieSRMVjHX6AhjY+Mh1zMZHx/HNToyu2yyr4fb/X1Bh4DWmqHBQSb7un2WTw0OMDY2FnRfDKW4ffs2U9OjdQGmbvczMDAQ1F3QMyYmJ3GODq/oMrhYmoTJKqZNE3MZHyCPqdHmncOi5bSjAeaFT6gfZg2YpvapkKZNM+SKaVrPtCFhEkkSJquYLT6BuJjokPdM4mKiscUnzi6LSk4hMTExpL2BuNhYbPPmqrHFxXurpIUgMSkRW8KdvtgSEklMTAzpNUXZ7dhiYuXO4QiTd3cVs8bEklpcRkZGRtDPyczIILXY98pHysZtrC8MvlKZYRiUlpWT+cijPsvTd+6hsLAw6FCy2WysXZtDwvqS2WVxuevIWJtDdHRw05AahsG6/HzSd0jdrUiTMFnlsvc/wcYNG4OutLZhw0a/qmT2xCQyyjeSs3ZtUO3ExsayNi+PpFLfEgIZO/dSXFwcVB1YwzAoLi4ma/cBn3IGM9XjysvKg+qLUorS0lIydx9YcluxMhImq1xiUSmFu3azZcuWJSutbdmyhcJduwNWJcv7yK/zyJ69JCUlLdiOUoqoqCgOHnyMwudf9BtsZotLIP/p5ziw/wDWRQLFMAxy1q5l885dZO0/5Lc+Y9c+ynfspCC/YNHXZDEM9u3bR/7jTxOVInWDI00Gra1yShkUvfgltKmJj0+guvo6Q0NDsx9C0zRJTk5m48ZNlOzeS9GLrwQ8txCdlsGGL7+GYRhcv1bFjRs3cLlcKKXQWqOUIicnh61bt1L6wudILt8UsD/ZB57CMzWFPcpOVVUVPT09swPbtNbEREdTUlpKxdZtlH3ptYADzSz2KMq/9BoAKVcv01BXx9j4uM9ryszIYPOWrRQ/8Qy5z3wsLO+lWJzUM3lIaNNksPoq3ScP09d6g5ER7yx5iYlJZBQUkr3/SVI2bl1yOgjX6Ai9ZyvpOVdJX28Pkw4HNpuN9LQ0MrfuJPvRJ4jJ8h9GP99oWzPdlUfoq73OwNAgHo+HuLg4MrLXkr3/EOk79mCJWvy8iOly0n/5HN0nj9Df2cno+BgWwyAlJYWMknKyDzxJYlFZ8G+SmHVfFUdSSv0h8DVgZpDAf9NavxNguw8D/wuw4K0P+ydLtS1hsjKO2304hwcBsCelEJ0W/AnaGabHzcStTjyOSQybnZisNVijg59SY4ZrbARHfy+m24MtLo6Y7LUhX3XRWjPZewv32CjKsBCVlo49MeCMtCJIywmTSB/m/JnW+v9daKVSygJ8B3gG6ADOK6Xe1lrXRLhfD7XotIxlBchchsUasNhzqGzxiT6XoZdDKUVs1lpviXJxz9zrE7C7gSatdYvW2gn8CHh+iecIIe5DkQ6T15VSVUqpN5VSKQHW5wDtcx53TC/zI9ODCnF/W1GYLDFF6HeBImAbcAv400BNBFi20Ox/Mj2oEPexFZ0z0Vo/Hcx2SqnvAT8LsKoDyJvzOBfoWkmfhBD3RsQOc5RSc68PfoLAU3+eB0qUUoVKKTvwEvB2pPokhIicSF7N+R9KqW14D1tagW8AzJ0iVGvtVkq9DvwS76XhN7XW1RHskxAiQiIWJlrrLyywfHaK0OnH7wB+40+EEA+We31pWAixSkiYCCHCQsJECBEWEiZCiLCQMBFChIWEiRAiLCRMhBBhIWEihAgLCRMhRFhImAghwkLCRAgRFhImQoiwkDARQoSFhIkQIiwkTIQQYSFhIoQICwkTIURYSJgIIcIiImUblVI/BmYmeU0GhrTW2wJs1wqMAh7AHep0hEKI+0dEwkRr/eLMz0qpPwWGF9n8Ca11fyT6IYS4eyI617BSSgGfBp6M5O8RQtx7kT5nchDo0Vo3LrBeA79SSl1USn09wn0RQkTQsvdMlFLvA9kBVv13rfW/Tf/8GeCfFmlmv9a6SymVCbynlKrTWh9f4Pd9Hfg6QOra3OV2WwgRIcsOk6WmBlVKWYFPAjsXaaNr+nuvUuotYDcQMEy01m8AbwDkb94WcD5iIcS9E8nDnKeBOq11R6CVSqk4pVTCzM/AswSeQlQI8QCIZJi8xLxDHKXUWqXUzOx9WUClUuoqcA74udb63Qj2RwgRQZGcHvSVAMtmpwbVWrcAWyP1+4UQd5eMgBVChIWEiRAiLCRMhBBhIWEihAgLCRMhRFhImAghwkLCRAgRFhImQoiwkDARQoSFhIkQIiwkTIQQYSFhIoQICwkTIURYSJgIIcJCwkQIERYSJkKIsJAwEUKEhYSJECIsJEyEEGGxojBRSn1KKVWtlDKVUrvmrfuvSqkmpVS9UupDCzy/UCl1VinVqJT6sVLKvpL+CCHunZXumVzHOzeOz1w3SqkNeKvTbwQ+DPylUsoS4Pn/N/BnWusSYBD4ygr7I4S4R1YUJlrrWq11fYBVzwM/0lpPaa1vAE14J9iaNT0P8ZPAT6YX/R3w6yvpjxDi3onUVBc5wJk5jzuml82VBgxprd2LbDNr7vSgwNSr6zNW44Rd6UD/ve5EBKzW1wWr97WVhfqEJcMkyDmF/Z4WYNn8KT2D2ebOijnTgyqlLmitdy207YNKXteDZ7W+NqXUhVCfs2SYLDWn8AI6gLw5j3OBrnnb9APJSinr9N5JoG2EEA+ISF0afht4SSkVpZQqBErwTgE6S2utgSPAC9OLXgYW2tMRQtznVnpp+BNKqQ5gH/BzpdQvAbTW1cA/AzXAu8BrWmvP9HPeUUqtnW7i94DfVko14T2H8jdB/uo3VtLv+5i8rgfPan1tIb8u5d1BEEKIlZERsEKIsJAwEUKExQMTJisduv+gUEr9oVKqUyl1ZfrruXvdp5VQSn14+u/SpJT61r3uT7gopVqVUtem/0YhX0a9nyil3lRK9Sqlrs9ZlqqUem/6Vpf3lFIpS7XzwIQJKx+6/yD5M631tumvd+51Z5Zr+u/wHeAjwAbgM9N/r9Xiiem/0YM+zuT7eD87c30L+GD6VpcPph8v6oEJk5UM3Rf3zG6gSWvdorV2Aj/C+/cS9xGt9XFgYN7i5/He4gJB3urywITJInKA9jmPFx2W/4B4XSlVNb37ueTu5X1sNf5tZmjgV0qpi9O3eqw2WVrrWwDT3zOXekKk7s1ZlggO3b+vLPY6ge8Cf4z3Nfwx8KfAl+9e78LqgfvbhGC/1rpLKZUJvKeUqpv+H/6hdV+FSQSH7t9Xgn2dSqnvAT+LcHci6YH72wRLa901/b1XKfUW3kO61RQmPUqpNVrrW0qpNUDvUk9YDYc5Sw7df5BM/+FmfALviecH1XmgZLoIlh3vifK373GfVkwpFaeUSpj5GXiWB/vvFMjbeG9xgSBvdbmv9kwWo5T6BPAXQAbeoftXtNYf0lpXK6Vmhu67mTN0/wH1P5RS2/AeDrQC37i33Vk+rbVbKfU68EvAArw5favFgy4LeMtbkgcr8EOt9bv3tkvLp5T6J+AQkD59e8wfAH8C/LNS6ivATeBTS7Yjw+mFEOGwGg5zhBD3AQkTIURYSJgIIcJCwkQIERYSJkKIsJAwEUKEhYSJECIs/n9P3Id0SGtE2AAAAABJRU5ErkJggg==\n", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "assy=Assembly(6,4)\n", "assy.set_pitch(1.62)\n", "assy.set_source(uo2)\n", "assy.set_coolant(h2o)\n", "assy.set_pins(fuel)\n", "\n", "fuelmap= [['1', '1', '1', '1'], \n", " ['1', '1', '1', '1'], \n", " ['1', '1', '1', '1'],\n", " ['1', '1', '1', '1'],\n", " ['1', '1', '1', '1'],\n", " ['1', '1', '1', '1']]\n", "\n", "assy.set_fuelmap(fuelmap)\n", "\n", "det=Detector('D')\n", "det.set_location(Point(174.726, 174.726))\n", "\n", "ex5=Experiment()\n", "ex5.set_assembly(assy)\n", "ex5.set_detectors(det)\n", "ex5.set_materials(uo2,he,zr,h2o)\n", "uo2.set_color(\"#0f0e0e\")\n", "he.set_color(\"#474545\")\n", "zr.set_color(\"#635f5f\")\n", "h2o.set_color(\"#96deeb\")\n", "ex5.Plot(xl=[-10,10],yl=[-10,10],detectorSize=0.4)" ] } ], "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.7.3" } }, "nbformat": 4, "nbformat_minor": 2 }