{ "cells": [ { "cell_type": "markdown", "metadata": {}, "source": [ "# Comparison between RS and WPM" ] }, { "cell_type": "code", "execution_count": 1, "metadata": { "ExecuteTime": { "end_time": "2020-05-01T22:03:09.074155Z", "start_time": "2020-05-01T22:03:08.024478Z" } }, "outputs": [], "source": [ "from diffractio import plt, np\n", "from diffractio import degrees, um\n", "from diffractio.scalar_sources_X import Scalar_source_X\n", "from diffractio.scalar_fields_XZ import Scalar_field_XZ" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**RS**" ] }, { "cell_type": "code", "execution_count": 2, "metadata": {}, "outputs": [], "source": [ "length = 150 * um\n", "wavelength = 1 * um\n", "x0 = np.linspace(-length / 2, length / 2, 1024)\n", "z0 = np.linspace(2 * um, 500 * um, 512)\n", "\n", "# source\n", "f1 = Scalar_source_X(x0, wavelength)\n", "f1.gauss_beam(A=1, x0=0 * um, z0=0 * um, w0=10 * um, theta=0 * degrees)" ] }, { "cell_type": "code", "execution_count": 3, "metadata": {}, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAn00lEQVR4nO2de7QnVXXnP7tvvwgyILQiARSMjQYVkDSoo1mCotM+RmatOA6OMaDM9EpGHRNNFEKWyzgrK2rWaMjEGHvEmBhHxEekR5kgEDCTTEQa5I1oB5GHaNsG0ES6m9t3zx9VRdetrsep16/qV/39rPVbt+rUObt23V/9vmfXrlOnzN0RQggxLVYM7YAQQojukbgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEWTm0A31gZhoCJIQIZYe7P2FoJ7pGkbsQYn/nu0M70AcSdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCASdyGEmCCjEnczO8TMPmdm3zSzO8zs+WZ2qJldYWbfjv8+fmg/hRBi7IxK3IELgb9292cAJwJ3AOcBV7n7euCqeF0IIUQJ5j6Od0mb2cHAjcBTPeWUmd0JnObuD5jZEcA17v70ClvjOCghxDxwvbtvGNqJrhlT5H4s8EPgz8zsG2b2MTM7EDjc3R+I63wfODyvsZltMrOtZrZ1Rv4KIcRoGZO4rwROBj7i7s8B/oVMCiaO6HOjcnff7O4bptgDCyFEXcYk7vcB97n7tfH654jE/gdxOob47/aB/BNCiLlhNOLu7t8H7jWzJJ/+EuB2YAtwdlx2NnDpAO4JIcRcsXJoBzK8FfiUma0G7gLeSNQBXWJm5wLfBV47oH9CCDEXjGa0TJdotIwQogYaLSOEEGI+kLgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEkbgLIcQEGdv0A0IIMXds3LjRd+zYkbvt+uuvv9zdN87YJYm7EEK0ZccPf8jWa/8+d5utOmDdjN0BJO5CCNEBDkuLQzuxDIm7EEK0xZdgcefQXixD4i6EEG2RuAshxARxl7gLIcT0UOQuhBDTw3VDVQghpofSMkIIMUWUlhFCiOkxwsh9dHPLmNmCmX3DzL4Urx9rZtea2TYz+4yZrR7aRyGEWEaSc8/7DMToxB14G3BHav39wIfc/WnAg8C5g3glhBCFxGmZvM9AjErczewo4JXAx+J1A14MfC6u8ufAvxvEOSGEKCJJy4xI3MeWc/9D4J3AQfH6YcBD7p5c29wHHJnX0Mw2AZv6dlAIIfZhhE+ojiZyN7NXAdvd/fom7d19s7tvcPcNHbsmhBDl+PjSMmOK3F8AvNrMXgGsBf4VcCFwiJmtjKP3o4D7B/RRCCFyGN9DTKOJ3N39fHc/yt2PAc4C/sbdXw9cDbwmrnY2cOlALgohRD4jzLmPRtxLeBfwdjPbRpSDv2hgf4QQYjlKy4Th7tcA18TLdwGnDumPEEKUMsKHmEYp7kIIMV+ML+cucRdCiLaMcCikxF0IIdqitIwQQkwRRe5CCDE9RpiWmYehkEIIMW5azgppZhvN7M549tvzcrY/2cyujmfMvTl+2LMURe5CCNGa5jl3M1sAPgy8lGj+rOvMbIu7356q9jvAJe7+ETM7HrgMOKbMrsRdCCHa0i4tcyqwLX6mBzO7GDgTSIu7E03JAnAw8L0qoxJ3IYRoS/lomXVmtjW1vtndN6fWjwTuTa3fBzw3Y+M9wFfM7K3AgcAZVS5J3IUQ+x2WWvYuDJZH7js6mK32dcAn3P2/m9nzgU+a2bPcfamogcRdCDEZrLpKT7R6QvV+4OjUet7st+cCGwHc/R/MbC2wDtheZFSjZYQQo8cCP4PRblbI64D18fuiVxPNirslU+ce4CUAZvbzRNOi/7DMqCJ3IcRgDCrIRNHtni4Mtbih6u6LZvYW4HJgAfi4u99mZu8Ftrr7FuAdwP80s98gyiSd4+6lGSWJuxCic4YU7UHSES0fYnL3y4iGN6bL3p1avp3ohUbBSNyFEEFMTbA7Px7NCimEGBNDifZcCHYoI5x+QOIuxAQZQuTGKNYzS9G4w+KuWe0tCIm7EHPAVKLrocS69/+fA+PKykjchRiKKUTX8yLWM4ngJe5CTJP9XaxnIdR19jHT78OBwmdFh2E04m5mRwN/ARxO9K/a7O4XmtmhwGeIZkC7G3ituz84lJ9i/2HexXrWUXWf0XSo7RCbvXyvSsuUsgi8w91vMLODgOvN7ArgHOAqd39fPM/xecC7BvRTzCHzLtQw26i6L6GepUjPXOgl7vm4+wPAA/HyT8zsDqLZ0s4ETour/TlwDRL3/Zp5F+pZpz76SHuE2ByLQJfZ6OxcUuQehpkdAzwHuBY4PBZ+gO8TpW3y2mwCNs3EQdEZsxbqsUTTYxDpsQh0n+I8E2EHiXsIZvY44PPAr7v7j832fgXu7maWO59CPD/y5thGJ7N4inDmWaTHGknPMoruU6CbCuys27VGN1SLMbNVRML+KXf/Qlz8AzM7wt0fMLMjKJniUnTDLIR6LBE0zIdAt42g+xDnrgV2DG0ao8i9GItC9IuAO9z9g6lNW4CzgffFfy8dwL25Yt7EeYix0mNMcYwliu2yzVD1m9hqhcS9lBcAbwBuMbMb47LfJhL1S8zsXOC7wGv7csBY/laW5OTI5niKTpqiXNDQ05rm0dUJPgVhDrE5RN543kW2ru26+8wrb3P10gqJezHu/ncUfzcvqWuvqejktQu11SZ6C03XDS3KY3uicN4i5rGLbJ26ddrXEeIhOpBOUM59NszbK6a6HifcxHbTfQzxMMo8Rspt685KjPvYdx1xb9tplNnuTRcc9ihynw1DpkKqTqCyDn4MN/e6tjnESIw+I+SuxXgoIR67f0V169gtKy/bb11c4j4bjDDB8NQy5OfbydmWt70OCxXbxzRGOdTWLCPpLn7wTUWnS6Frur8Q3/v0qU69orpDROd9BnwS9xlhVB+Ys/zLzgp9tm7etjb58yHHHYe0r7IxJrFuGm3OUoCbiG/IcXVVJ7RdHx1oUXnZOdb18MzWSNxngwGrSrYvsfyLzq6ny8lsKzt5kk4gLzLvMrVQ1iakbVX7rnztOpLrS+yaCF2VjS78mOX/JHR/dex1UbdNm7J2IW3r4rqh2j8GrC7Ylhby9HI6TZMV9RBBL/Ijj75P+DZt6/gRKt5dpUCaiF1dQa3rQ939d7GPpvvpo15R3Sb1y9q0aRfSvi1t0zJmthG4kCg2/Ji7vy+nzmuB9xDJzk3u/h/LbE5W3NcUbFvB8rRJWuiz9RJbkC/ieXbyfClbL2vbx4netsPpKpc7CxEsu+Iag4h3JdB9iHMTYW4aXbe9Cg2xUcdWE9qIu5ktAB8GXgrcB1xnZlvc/fZUnfXA+cAL3P1BM3tild1JivsK4MCc8qKbqOlIPXszNi3gIQ8pNb0Bl9c2dFsTm33eqBs6Oh7b9qZ1israXEWV1S2qX1YeYjNke8g+6tpru486OLDYPHI/Fdjm7ncBmNnFRLPh3p6q85+BDyfvsnD3ymlYaou7mR0I7HT3PXXbzooF4KDUenLzNPmbRO9poS9Lz5SNoimKDuv+qEK2d5UL7yv6DqnTZnvXtkPadJn+6EOom9QvaxPStk6dkH21sR1KrzdTAbw0577OzLam1jfHEx0mHAncm1q/D3huxsZxAGb290QS9x53/+sylyrF3cxWAGcBrwdOAXYBa8xsB/Bl4KPuvq3KzixZAA5OrSdinkTme1LrSVlC2duy8oQm5MdfZKOMpnb6zr12HXWPZR95derU6yq10SRFV2Uz1EYdW120CaV3Ye6C8rTMDnff0HIPK4H1wGnAUcDfmtmz3f2hsgZVXA1cSZTvudU96p/i19+dDrzfzP7K3f+yne/dsRJYFy+nI3aIhB2iDiCdmslG67Bv3r1I0IvKQugi+mmb+ukifxxap6sUxizEuah9WXmZrZC2Ie2b1q3jg6iH0+qG6v3A0an1o+KyNPcB17r7o8B3zOxbRGJ/XZHREHE/Iza4DHf/J6LpeT8fT9U7GlYB6bsNyQ9gKV5eSK0nwp4dOZP8zZZ1Qd0fVlc3SIvqt0kXzFJQ66a6xjDyoo2I9hkNj5WuXsQw61GJLUfLXAesN7NjiUT9LCA7EuaLwOuAPzOzdURpmrvKjFaKe1rYzWxVgdDvUzYkq4Aj4uV0fn0htZ6kZkjVS8S8jaDXOTm7vIHU5CbtLHK4feZ9+0wzDBHZ5vk59Ftn2ohkXd/r7CvU9sxEvoW4u/uimb0FuJxIpj7u7reZ2XuBre6+Jd72MjO7nUi+fsvdf1RmN/iGqpl9DHilmS0C3wNuBm529//R7JD6YzXLr3GcvQe6yF6xT6L4dN697EdddKJ0+QPsOiLsKhfbV6Q65BVRX3QtKPMgklV2+mxf1rZpuya0eYjJ3S8DLsuUvTu17MDb408QdUbL/CJwlLvvMbMjgROBE2q0nxmrgWNWwlL8z165GhZ3w56lvRH8wkqwFdFnReqR0qXUGKD0l5X3xS3NICxY0YFi2QxUr8999PXkXxO7od95iO02gtdEtMra1A1c6tavek6k6T67qNsF8z63zLXAYcB2d7+fKDd0WXmTYVh9ADz5abF4G+zeGf3YFlZFor6wsFyMkh9i8sN1X16e0FTMlzKDRldUzRzWE1k/+qCuYIbWD/nfl9laKvlll7Yr2Jbb2Rf8f0MDg9x6OX6H2svzJ9u2qW8h4tm0TrYspGPIq9O0XRPmXdw/CnzVzC4iEvqb3f3hftxqx+oD4MnPAlbA0iNx9LuG5XMSLKU+yfqe1HK6aseRo3cYQtTxLSiarKhT1kEU2S8qzxWRnP9NsBCGCFfG/6o2+2yvENts/ez+qsS1zJ+s0NdqW+FHuu0+20r22+Z4K//3ZLYTvj27rS9RTxucZ3H/S2Bz3Oa/ACeY2Vp3/7lePGvDWoPjV8HO3axYuxrWPh5WrISlRXj0kejvnp1R3aWl5d98jiAEZRyadgAh0XSJ7YUurm2rfgmwr59Vber+Msu+g/T6nuJ6ywTK88uhWGBCUnIh+yisXyKEQe29uk7tYyjoPEJFuGh/oZ1FaEdRliIt3VdqOeSqoSnzHrnf5+6/ny4ws6IpXIZl7cFw7L+Gg54UifriLljcGQn7Y6IefxNd5yqWAr7hojp55XtyBiJl63lOu3SdfcLVxeL19P6W2SjZ57J6SY4rbb9gOU+oi36NedsztlYUbF9oYjPPzxwhqhL4POHNE6o8Ea8S/CKRzd2e429VB9G0cwiqW7NTCOkQimyG7qMts7gHV4c64n6jmb3N3S9MCtx9Vw8+tWf1QXDSObBnFzzyIDz6z9E3HyK8bakj3NmyrJBXiXiRgBeJcp5wp7fnCXae3Tz7ie08e3l+Zm2V+Za1s+wXXfA3L8WWLGefWAsR9MzfFTmdxEKVL6myrNiWdQx5ohXSKYR2CFWdUt49qBDfCm1XCHXjTmBVte2lPcuvxH0J2E1r5j1yPxw4w8zeBdwA3ATc6O6f7cWzDCFTYj7GmoPhoCNg50OwKjOF2ELO81YrAv8NVZ1Dkyi7LKruRKQLxDSvLM/uPsJaIc5Fwlxko8pOV/vPdgyhnUJoh1DRGWTbr8hsXyiwU9QJlAltVx1AkPiv2lf8a12tLOwtW4iX01cACyvzfVhI2U7bXEgfR2qcbTKAYmFlzjF3IO5znXN399fCY6mYZwLPJprcpndxD5kSc7mzi9GQlAMPj4Q7Ee9E2EPFvIgika+KvIsEGMIj5PRyqRiX2AsWxAIhrSPIdcQ4RMh732eJ+LcR/kDR3+cKISP6CxkBdd8raLnCRb7oLyzsFddEnC0roPFw4gX2FfxlIppsi8V0RcpmMpQ3T/Cz+0vqJYPJEnFeclhYkaq3Iq6XLks6hj1ReVJvH5tx3aQDWehoCG/L6Qd6obbKxamYG+LPrAiZEjOFwc88AVaujc40Sx1mG2EvEvU60Xh2PS8iL9peJLrp5WDhrRDuEOEPEf2QCLxSfNfEZWv27nfl2vz92WKsRjntku8+b38Wl60oEftECIqi9WQio/RMdem/eW1C2yZ/FyIbK1I2l2Kx9rh94nJWt5Yy64+JYNZGqm72kNNtsnbT+1vKFOZp6GP929LeyPqx/aXLUumcsnq+tHeYcdp2uoPJ2uyM8lkhB6HOE6q3ED+VGn9uAc5299/rybc0lVNimtkmYBPAk4980gxcEq0pS3PN4v6IEB3hwJ6RTYJeJ4x9EdETqScQTWzzaeA2YBbiXkk8P/JmgA0n/rzz0x9GEXteWgaaR/ChKZm8ulV58WydLtIyITn2PtIydXLkdfLrXaVl6uTk2+Tja+biQ3P6Vbn4srx30UN7Zbn1bFqmap+t8vDZspybpFWjcOrc1O0CJ5raZEzUybn/E3BN/Ele+/Q7vXi1LyFTYu7FVkbf8s4fRcMfYbmYN7mpWhVJ5ol7XruQXHveepv0TXp7WT4/ZMRL6xuaNdNBY7ih2vZmKgE2Gor4skNqcUM1VwAbjqZpI+RA6ZDKOiNpmozhb8PIsjK10jLHufu3knV3/7aZndCPW/sQMiXmXnY9BD++N/r2ZzUUsm6KIa+szlDIZeIc0EGEjF8PGbqYXq7qcIqEOc9GWYfVxVDIOkJbNv6+KOKuaW+MQyFDx8VXjS1vK7RdjIdvMha+DXMduQMfNbOfIxLXm4G1wK1mdqC7/0sv3sUUTYlZ2GD3P8ONn4DHHQEr1+z7EFOe2DSlTodRVLdNqqds7Dvse4xlVwdFHUB2veoBJsgXvfRy0ZOmZSJcx1YTmzX8LHuIqfV48ADB36d9SZRd5G+on6W+Ju0CIuOhHmKq80RrG4aenjlLyGv2zCNOj9efTDQj5Enx3xvMDHd/ep+O5k2JWcjOh+GuK2H3EiTTD8De0RuJwBddolfRxQnRctqB0u1tph7IK8vzNVunar3OxB9Fwpz1peTyuu7UAGVzolSJcOg+qp4oDd1XWVqh7+kH9hHKFiJd1qbJ97HPNjLbKN7WFifsJz1Lgl6zZ2afBy5193vc/R7gHjO7nGga4LOBrX06WZudDt90WAF7frobsx+w4gCio02P7Up/oPQx87p0OTlYlro+hUYnVXbnadKwvHp5/peJQ+72mpN3dTqRVo19Vx170/22Od7K/212O5ntFG+vM3VRXv0umEdx3wi8Cfh0nPN+iCglswB8Bfigu9/Yl4NN2P0I3HNrtLxioZspf8uEbxZT6Y6Zppe2XXU6VbZmMd0v5J8HtTq3vE6rwyl/89qF2Mqtk7Wzb5Xc6HjIKX9D/GnKXEbu7r4T+BPgT+J3pa4DHil76/bQ7H4E7k5l5JOXdSzG364RPYyRvLAjjzoRxyzp4uUdMJsXeIyRph1R6HceYj/EVJnolLVv0q5udq+oTVH9OlnCovqzrtuEuRP3NPG7Uh/oyZfO2A3ck0ThADuXv2ZvD8AS2O4oS7NA9Wv2+uz18yh8/VzGka5ekdflq/ZC7TW13ZcPTemrr29yftXxJdR+284oxMYsO7Oqdk3wHmy2peUkK+NkN8sfZ4VI0FcR9oJsUn9DRKeLL7ULAdZLstvZamt/SrTtsOr+JobqlLpknodCzg3Zy4tk2o49xO9PjcuX2NvjJidCWtSz0Xwf0WXaftM6XQh0Xlle+7yy0Lah+y1q36R+k/9bVbuQ9k3sdbGfeaXLqHeIrGlaQ5oQOuutmf0S8DngFHcvHchS5yGm47OzMJrZae5+TaiNWfEosD1TlvxAkrxYOhWTFfiErKhbzra8+k1oG7n3Ldhd1QnxKaROm06mqO6sOo6qtiHtm9at44MIp80N1dBZb83sIOBtRK85raRO5H6JmX0S+ADRaJkPABuA59ewMRMWgR2p9fSkfMnyUrycpGMSqnrgvCg+RJzybFTR1E6bSLrqWJoIeJUwV9Vv0kE06TDa/N/qtC+zUXZudNEhdHXudd0uhLF3Ri0i99BZb/8b8H7gt0KM1hH358aG/x9wEPAp4AU12s+MPUD6zd3O3tQMxHNRs/xScIF9Sb6s7CVjkbA3EaG627vKeXcVHbcV6qG3h3ReXXUUdc6HOp1KX6mnEKEOFfMmwjxPN9krIvd1ZpZOoWyOJzpMCJn19mTgaHf/spl1Lu6PAo8ABxBF7t9xH9sMxhF7gJ/klKcFPp2SSUf12Tx7kcBn7aZpEnkWtQ3dVma3jVCEiFQTQWwrumXtZ91hNLERamfWHULZudn2inQs9zH6urooEfcd7r6hqV0zWwF8EDinTrs64n4dcClwCtFY9z81s19y939fZ4ezYAkomuwmnaJJ1qE63x4yRrbrH2vVti5HrrSJ+vuKdvvsAMZgv6v7GH3UK6pbVn/WQUtdO3Xt1aHlQ0xVs94eBDwLuMbMAJ4EbDGzV5fdVK0j7uemDD0AnGlmb6jRfmY4UPTm7iTXnl1OR/V5I2fK9lVG2x9YUf0+2/URLfaV2uhbcPtI4/RxTyFkP6G2i+p1UbesflmbqnZVbUNttKGFuJfOeuvuDxMF1ACY2TXAb3Y2WibPkLt/MrT9LHHK33mbFXjIPzHytnUt9Ald/xD6atvFFUGbfHRfN0rbRu0hbZqkbsZ4T6DOOdDnFWRVm6p2VW3r0OYhpqJZb83svcBWd9/SxO4kx7k70Q2CMhZZ/sUnDyDknQxF24pOjNAbEUOI8yz2P+uOocsOpKuIeKh7F03rFJWFinvbK4K650wfwVBb2jzElDfrrbu/u6DuaSE2JyvuVf/oopurZMrI2ZZsbzuXRFH7ri4vu8xBtu1MxtQZ1Knb5X2Urm6yNhXrWQh1m46kqP0sov+2tH2IqQ8mK+5D/qOrhLBqHH3dTqPopM3aqRO5ZG227YiqjmlWN46L2oytwwhtOyZxbtN5dFG3rH5dO3WZ9zcxzRVDTuLTJqKvatsk+khO4Dp+he4nsdn2KqFuR5TYadLpdNlZ9JV/blt3bIJc53iGiOq7QBOHzYixXSIltM3Tl9koIk8A21yi5u2/604p9IcS0lnUSSnlHUeZgNTtXPpMQRXV76tuXyLc5/9ipKNlemGy4t6kF03n4ZP1Ilt5J0nRPtN1u+h0uuq4mkYx6ZO47Y+lyIeyH0qdfYamucpshl6dlHUSVVcS89BJdFV/qJuqfUftEvcRkxXnsg6iTucxxOVaiAC26SSSH0rbY9tDfx1EYj+Puh1Ema0Qe3uo1zmk26V9qGrXRacy1M3vsvpN2sxS5CXuBZjZHwD/lmh4+j8Cb0ze9GRm5wPnEv3v/qu7Xz6Un/NEG9Htu2OA5T+woTuIqh97m7RW2vZYO4ds27odRJM2Ta5Uurq6KWvTBol7PlcA58eD+d8PnA+8y8yOJ3pa65nAzwJXmtlx7j62/+OkaCq2dQS2y84B2nUQXVzBFNFVx1BkK9RemwfsmnQsofc8umjXRWfUBbqhmoO7fyW1+jXgNfHymcDF7r4L+I6ZbSOaHvMfZuyiCKDpfY4mNBXkvB/+vHQMTVJJRbbq2Ky6b1Flo2nHUjcllW3XpG1TNBQyjDcBn4mXjyQS+4T74rJ9MLNNwKZ+XRNdM8sOAbrtFKDdVU4TX0LTCU1veoeOhKrz4FubjiHkhnjbq44uGPrZmjxmJu5mdiXRbGZZLnD3S+M6FxB1gJ+qaz+eH3lzbGdsV0iiQ2bdIcA4OoWmHUKZH1kSIa37/+qyU0jbaztqqu1VRx3Gliuembi7+xll283sHOBVwEvcPTn3q6bCFCKIIToEqC/GZSLc5hj6vEqAZp1CnVFOXd6n6GOsu0bLFBC/HPadwIvc/aepTVuA/2VmHyS6oboe+PoALor9kLpiOrbOANodQxNf6uyvy6uErM2EPh9ayrLfpmUq+GNgDXBFPBn919z9V+NpLy8hepfgIvBmjZQRY0WdQX3qdgjQbOhr33l3Re4FuPvTSrb9HvB7M3RHiJmwv3cG0NyfuvudRVQtcRdCNKKOoHWVjpj3zgCaXR3URZG7EGImDNERQPcT4HVxdTOrXPjYhuhJ3IXYz2k7A2dT+pgJdahOTQ8xCSHmlqGEE4bvCELQaBkhxOQZ6moA+n03QhFjzLnP4slcIYTIxQM/fdFVtJ2Ie94nBDPbaGZ3mtk2MzsvZ/vbzex2M7vZzK4ys6dU2ZS4CyFGz9CdQAhLBZ8qzGwB+DDwcuB44HXxjLhpvgFscPcTgM8BH6iyK3EXQkyGoTqBlpH7qcA2d7/L3XcDFxPNiLvXvvvVqaf3v0Y0FUspEnchxH5HH0JfIu7rzGxr6pOdvfZI4N7UeuHstzHnAv+nyh/dUBVCiJZU3FDd4e4butiPmf0ysAF4UVVdibsQQrSk5WiZoNlvzewM4AKiCRZ3VRlVWkYIITqg6Q1V4DpgvZkda2ariV4tuiVdwcyeA3wUeLW7bw8xKnEXQogOaCru7r4IvAW4HLgDuCSeEfe9ZvbquNofAI8DPmtmN5rZlgJzj2F734sxHfQmJiFEDa5vmxM/yMyfU7Dt/3ZgvwnKuQshREv263eoCiHElJG4CyHExBjj3DISdyGE6ACJuxBCTAxF7kIIMVHGJu6jGuduZu8wMzezdfG6mdkfxdNg3mxmJw/toxBCZBnjTJWjidzN7GjgZcA9qeKXA+vjz3OBj8R/hRBiVChyL+ZDwDtZ3tmdCfyFR3wNOMTMjhjEOyGEKKDtyzr6YBTibmZnAve7+02ZTcFTYZrZpmRKzZ7cFEKIQsYm7jNLy5jZlcCTcjZdAPw2UUqmMe6+Gdgc70vTDwghZsZ+/YSqu5+RV25mzwaOBW4yM4imu7zBzE4lcCpMIYQYGuXcM7j7Le7+RHc/xt2PIUq9nOzu3yea9vJX4lEzzwMedvcHhvRXCCGyjDHnPprRMgVcBrwC2Ab8FHjjsO4IIcS+6CGmAOLoPVl24M3DeSOEEGFI3IUQYoKMbRSHxF0IIVriwOLQTmSQuAshREv266GQQggxZZRzF0KIiaHRMkIIMVGUlhFCiImhyF0IISaKxF0IIabH5cC6gm07ZulIgsRdCCFa4u4bh/Yhy+AThwkhhOgeibsQQkwQibsQQkwQibsQQkwQibsQQkwQibsQQkwQibsQQkwQibsQQkwQibsQQkwQibsQQkyQ0Yi7mb3VzL5pZreZ2QdS5eeb2TYzu9PM/s2QPgohxLwwirllzOx04EzgRHffZWZPjMuPB84Cngn8LHClmR3n7mObgE0IIUbFWCL3XwPe5+67ANx9e1x+JnCxu+9y9+8A24BTB/JRCCHmhrGI+3HAL5rZtWb2VTM7JS4/Erg3Ve++uGwfzGyTmW01s609+yqEEKNnZmkZM7sSeFLOpgtiPw4FngecAlxiZk+tY9/dNwOb4315O2+FEGK+mZm4u/sZRdvM7NeAL7i7A183syWiie/vB45OVT0qLhNCCFHCWNIyXwROBzCz44DVRG8v2QKcZWZrzOxYYD3w9aGcFEKIeWEUo2WAjwMfN7Nbgd3A2XEUf5uZXQLcDiwCb9ZIGSGEqMYiDZ0WyrkLIWpwvbtvGNqJrhlLWkYIIUSHSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCjOVNTF2zA/ju0E50zDqi45oaUzyuKR4TTPe4nj60A30wSXF39ycM7UPXmNnWKb4tZorHNcVjgmkf19A+9IHSMkIIMUEk7kIIMUEk7vPD5qEd6IkpHtcUjwl0XHOFufvQPgghhOgYRe5CCDFBJO5CCDFBJO4jwcw+bmbbzezWVNmhZnaFmX07/vv4uNzM7I/MbJuZ3WxmJw/neTFmdrSZXW1mt5vZbWb2trh83o9rrZl93cxuio/rd+PyY83s2tj/z5jZ6rh8Tby+Ld5+zKAHUIKZLZjZN8zsS/H6FI7pbjO7xcxuTIY9zvs5GILEfTx8AtiYKTsPuMrd1wNXxesALwfWx59NwEdm5GNdFoF3uPvxwPOAN5vZ8cz/ce0CXuzuJwInARvN7HnA+4EPufvTgAeBc+P65wIPxuUfiuuNlbcBd6TWp3BMAKe7+0mpcfrzfg5W4+76jOQDHAPcmlq/EzgiXj4CuDNe/ijwurx6Y/4AlwIvndJxAT8D3AA8l+jpzZVx+fOBy+Ply4Hnx8sr43o2tO85x3IUkdC9GPgSYPN+TLF/dwPrMmWTOQeLPorcx83h7v5AvPx94PB4+Ujg3lS9++Ky0RJftj8HuJYJHFecvrgR2A5cAfwj8JC7L8ZV0r4/dlzx9oeBw2bqcBh/CLwTWIrXD2P+jwnAga+Y2fVmtikum/tzsIpJTj8wRdzdzWwux62a2eOAzwO/7u4/NrPHts3rcbn7HuAkMzsE+CvgGcN61A4zexWw3d2vN7PTBnana17o7veb2ROBK8zsm+mN83oOVqHIfdz8wMyOAIj/bo/L7weOTtU7Ki4bHWa2ikjYP+XuX4iL5/64Etz9IeBqopTFIWaWBExp3x87rnj7wcCPZutpJS8AXm1mdwMXE6VmLmS+jwkAd78//rudqCM+lQmdg0VI3MfNFuDsePlsopx1Uv4r8Z395wEPpy4xR4NFIfpFwB3u/sHUpnk/rifEETtmdgDRfYQ7iET+NXG17HElx/sa4G88TuiOBXc/392PcvdjgLOIfHw9c3xMAGZ2oJkdlCwDLwNuZc7PwSCGTvrrE32ATwMPAI8S5fnOJcphXgV8G7gSODSua8CHifK8twAbhva/4JheSJTvvBm4Mf68YgLHdQLwjfi4bgXeHZc/Ffg6sA34LLAmLl8br2+Ltz916GOoOL7TgC9N4Zhi/2+KP7cBF8Tlc30Ohnw0/YAQQkwQpWWEEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNyFEGKCSNzFpDCzA8zsq2a20IGt1Wb2t6nH74WYGyTuYmq8CfiCRxN7tcLddxM9xfgfWnslxIyRuIu5wcx+NX6bzo1m9h0zuzqn2uvZO08IZnaNmT0jXj7M4jddmdlnzeyPzezvzOy7ZvZCM/ukmX3LzC5K2ftibFOIuULiLuYGd/9Tdz8JOIVo/p30ZGTEr4B7qrvfnSp+GvCtePkEovlCAJ4N3OXuLyR6QcNFRHOZHw+80szWxPVujfcnxFwhcRfzyIVEsxD+70z5OuChZMXMngLc7+7JyydOAG42s7XAIUQvp4BocrOL3P0Bj148sQfYDY/N2747mVlQiHlB4i7mCjM7B3gK8Ls5mx8hmq0w4USimRsTfiFefyZwQ0r0TyR6QxRmdhTwPV8+o94aYGcX/gsxKyTuYm4ws18AfhP45ZQwP4a7PwgsxJE5RC+vXhu3XQ+cSZSWeTbRFLAJJ7C3E1jWIZjZYcAOd3+004MRomc0xEvME28BDgWujl/Vt9Xd/1OmzleI5pG/kkiod5rZTUSCfTvRixkOIZqDnLgjOCDuGGC50AOcDny5j4MRok80n7uYFGZ2MvAb7v4GM/s2cLK7/6SFvS8A57n7tyorCzEilJYRk8LdbyCK7A+OVlsJ+2rgixJ2MY8ochdCiAmiyF0IISaIxF0IISaIxF0IISaIxF0IISaIxF0IISaIxF0IISaIxF0IISbI/wcn3S+lFRTYIgAAAABJRU5ErkJggg==", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u_rs = Scalar_field_XZ(x=x0, z=z0, wavelength=wavelength)\n", "u_rs.incident_field(f1)\n", "u_rs.RS()\n", "u_rs.draw(kind=\"intensity\", logarithm=False, colorbar_kind=\"vertical\")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**WPM**" ] }, { "cell_type": "code", "execution_count": 4, "metadata": { "ExecuteTime": { "end_time": "2020-05-02T19:20:47.379068Z", "start_time": "2020-05-02T19:20:46.821378Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXcAAAEKCAYAAADpfBXhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAmAElEQVR4nO2de9QlV1nmf09/nU6YkEkkDTHmYoI0apAAsQmwwEWQ6LSXIbOWDAODGjQzvXTAQUUhERdLmeUawDVgHBHpIXhBxhAukh5kDElMdHQkpBOSkIuBNtwSg6ExgKPk0t3v/FH1dVefrvvlVJ3q57fWWd+pXXu/+63znfPsXW/teksRgTHGmHmxYWwHjDHG9I/F3RhjZojF3RhjZojF3RhjZojF3RhjZsjGsR0YAkleAmSMqcueiHj82E70jWfuxpgjnc+P7cAQWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGWNyNMWaGTErcJZ0g6QOS/kbSXZKeI+lxkq6W9Jn07zeN7acxxkydSYk7cCnwpxHxHcDTgLuAi4FrI2ILcG26bYwxpgRFTONZ0pKOB24BnhgZpyTdDZwXEfdLOhm4PiK+vcLWNA7KGLMK3BQRW8d2om+mNHM/E/gy8LuSPinpXZKOBU6KiPvTOl8CTsprLGm7pF2Sdi3JX2OMmSxTEveNwDnAOyLiGcA/sRCCSWf0ubPyiNgREVvnOAIbY0xTpiTu9wL3RsQN6fYHSMT+79NwDOnfB0byzxhjVobJiHtEfAn4oqT1ePoLgTuBncCFadmFwJUjuGeMMSvFxrEdWOBngPdK2gTcA/wEyQB0haSLgM8DLxnRP2OMWQkms1qmT7xaxhjTAK+WMcYYsxpY3I0xZoZY3I0xZoZY3I0xZoZY3I0xZoZY3I0xZoZY3I0xZoZY3I0xZoZY3I0xZoZMLf2AMcasHNu2bYs9e/bk7rvpppuuiohtS3bJ4m6MMV3Z8+Uvs+uGv8rdp6Mes3nJ7gAWd2OM6YGA/XvHduIQLO7GGNOV2A97Hxrbi0OwuBtjTFcs7sYYM0MiLO7GGDM/PHM3xpj5Eb6gaowx88NhGWOMmSMOyxhjzPyY4Mx9crllJK1J+qSkj6TbZ0q6QdJuSe+TtGlsH40x5hDWY+55r5GYnLgDrwbuymy/GXhbRDwJeBC4aBSvjDGmkDQsk/caiUmJu6RTgR8C3pVuC/he4ANpld8H/s0ozhljTBHrYZkJifvUYu6/AbwWOC7dPhH4akSsn9vcC5yS11DSdmD70A4aY8xhTPAO1cnM3CX9MPBARNzUpn1E7IiIrRGxtWfXjDGmnJheWGZKM/fnAi+S9IPAMcC/BC4FTpC0MZ29nwrcN6KPxhiTw/RuYprMzD0iLomIUyPiDOClwJ9FxMuB64AXp9UuBK4cyUVjjMlngjH3yYh7Ca8Dfl7SbpIY/GUj+2OMMYfisEw9IuJ64Pr0/T3AuWP6Y4wxpUzwJqZJirsxxqwW04u5W9yNMaYrE1wKaXE3xpiuOCxjjDFzxDN3Y4yZHw7LGGPMDJngk5hWYZ27McZMnG43MUnaJunuNLX5xTn7T5d0XZoO/bb0Tv5SPHM3xpiudAjLSFoD3g58H0lyxBsl7YyIOzPVfhm4IiLeIeks4KPAGWV2Le7GGNOVbqtlzgV2pzdsIuly4AIgK+5Bkm8L4Hjg76qMWtyNMUckSv9GH8bKZ+6bJe3KbO+IiB2Z7VOAL2a27wWetWDjV4CPSfoZ4Fjg/CqXLO7GmFmh6ioDUHpBdU8PqchfBvxeRPw3Sc8B3iPpuyJif1EDi7sxZiUYR7Rr0i0scx9wWmY7L7X5RcC2pKv4a0nHAJuBB4qMerWMMWZUVPM1BL0JYLeskDcCWySdKWkTScrznQt1vgC8EEDSd5I88+LLZUY9czfGDMJYM+1RZqwdVstExF5JrwKuAtaAd0fEHZLeCOyKiJ3Aa4D/IennSC4TvCIiSi8XWNyNMY0YQ7SHEOzej6PDTUwR8VGS5Y3Zsjdk3t9J8rS62ljcjTEW7K44/YAxZpksW+z6Fuyu/i8tRBMBex9eVm+1sLgbs2KssmCPJdaDf2YBTCu1jMXdmClwpAr2ssR6KTN4i7sxRwbLFOy+xGuMmXWTPpvYX+qAGUDh7UTjMBlxl3Qa8AfASSQf1Y6IuFTS44D3kSTJ+Rzwkoh4cCw/zZHLkSbWqyjUdWwO8n90WKaUvcBrIuJmSccBN0m6GngFcG1EvClNhXkx8LoR/TQz4UgS66kIdR2bfdlautBb3POJiPuB+9P3/yjpLpKEOhcA56XVfh+4Hou7WcBC3U8/yxbprgI9dPvaeOZeD0lnAM8AbgBOSoUf4EskYZu8NtuB7Utx0AzGql5YXIZQT1WkxxTosrZLE3awuNdB0mOBDwI/GxFflw7+CyIiJOXecpum0NyR2ugli6dpxyrfEDO12XRfoYUxBXoIcR6iXWd8QbUYSUeRCPt7I+JDafHfSzo5Iu6XdDIlWdBMd8bMvGeB7mZnqFlsUdu+BXYKbVrjmXsxSqbolwF3RcRbM7t2AhcCb0r/XjmCe5NlKmlQp3Jn4pQFelVmv32KeV99NK3fZaBrhcW9lOcCPwZ8StItadkvkYj6FZIuAj4PvGQc9/KZiri2YWq5PYZc0TE3cR5LNPPqDynIRTaGPs7GWNyLiYi/pPizfmFTe6ssuussKy/GHG5cmZowT0mUu9YdSoynMoD0hmPuy2EOTyEZ4os4pdu9x4wzz0Vkm9is275L310Fu+nn3MdA0AsB+zxzXw5Tn7lPIUY9xh2CQyx7W7YgdxWwxXpTEeIuAtxVfIcaJMvK+/wNhsV9OYh+/nHLGiC6+DrFGfXUYsh9z4wX67UV4rbC39fgMMZnVVRW1L6v70tZmz6wuC8JMa0DW0YoZOi49NCCPZXZcJsZbht/2vjS12DQ1p+6torqFdXtIxRTZqesTW9Y3JeDgKNG7LspfYZH2oZFphDaqBMm6EPw2gwETYV3iAGgr8+jTl91+2vavqhum/plbbq0a0v4gurwCNi0xL6aMrVZdtfVEH3Fc6cgul0Fd4yBpchOX8LfV92y+mVturSr074PHJZZEgKO7snWkCGVLgI+Zix6WSGFJkLbt8hOYeY+Vky8ad2mdprs73qm2sRWF7qKu6RtwKXAGvCuiHhTTp2XAL9Csqr+1oj492U2ZynuG4Bje7K17AuWbWYoy14xMbTwdbE/9Mx9ahdFi8qatC+rX1ZeZa/O/ir7bex1sd+WAPa2FHdJa8Dbge8D7gVulLQzIu7M1NkCXAI8NyIelPSEKruNxV3SscBDEbGvadtlsQYcN5Dtul+uqi/U2CLex+x71ff3Ef5oez2hqG0foY+yNl3aNalT1U8Xu3322RvRKeZ+LrA7Iu4BkHQ5SarzOzN1/iPw9vUHFUVEZY6tSnGXtAF4KfBy4JnAw8DRkvYAfwK8MyJ2NzuWYVkDjm/RbohZwdCx8SHjr1OMiS9TkPsW7qK6ZfWr9tURxKFi0n3OiEcT5b4oD8tslrQrs70jzWK7zinAFzPb9wLPWrDxZABJf0Uicb8SEX9a5lKdmft1wDUkpwS3RyTjU/r4uxcAb5b0xxHxhzVsLYWNwOaONpp8cbv8wJr+cIeeKQ4VG+5rwBkivtxUcNtcLCxrV7d903p1+23LMsIdbRgj33dQKu57ImJrxy42AluA84BTgb+Q9NSI+GpZgyrOj4hHFwsj4h9I0vN+ME3VOxmOAioDUg3oMybYVzimyFaTun2HCroMPE3aF9UtKy+zVdWuTvumttrY7ULfAr/sVX99CPaQPne8oHofcFpm+9S0LMu9wA2pFn9W0qdJxP7GIqOV4p4VdklHFQj9YWVjchRw8kC2m/xIugpG30LVZ+x2jHhvX9c72tjss88x6Hs220Yom/pQt4+6dgcdkLqJ+43AFklnkoj6S4HFlTAfBl4G/K6kzSRhmnvKjNa+oCrpXcAPSdoL/B1wG3BbRPz3ujaWxSYOHQaHZqiZWp9itooz05WPw1bQRWzGFMoqW1U2urQva9u2XV+0vaAaEXslvQq4iiSe/u6IuEPSG4FdEbEz3ff9ku4E9gG/GBFfKbPbZLXM9wCnRsQ+SacATwPObnMwQ7MJOH0kZdjQc7+auL2p0/oH16Dd/hp169gbQriK2pW1KdrXxFZR3T5sDFW3C13XuUfER4GPLpS9IfM+gJ9PX7VoIu43ACcCD0TEfSSnDx8tbzIOmx4DZzxpbC+mRR0BamW3xYLYusJZSxAr6pTZ2F/wKy9qU9RXXv2izyW3bk5Zbr0cf+vYy/Mlr91iWR2/6ohn2zqLZXn/rjp16vTVlVW/Q/WdwJ9LuoxE6G+LiK8N41Y3Nh0Dp3/nCB33+K1pKsZRczrSSQzL9jUQsybimXdci/XqCFwdAa5qc9j+Bd+q6i/2l61f2Vdme1HkG7Ut8WGx7WH7Svod6lihXMCrxH1xf53BoBUrLu5/COxI2/wn4GxJx0TEtw3iWRceAzylYfyh69S2bvOqeun+XO+LZslFNnPK1+pMf+pMoxZ9qbKR3W7yCyx6v+hDgdAsCnCReGXLs8JUJEpFfdSxXyS2pcJa0F9vx1AweDQ9hjZ9lLWpO0iUDmgs7KN4X1tWfeZ+b0T812yBpL5SuPTL0cfDlueX19nf8j8RFe3K7OZNb/Pq1ynb92j5/kU/s/sPm7LuzX+f7aOoThSU7y/4deW9LxDpA++Lfo37c8r2pQNjQZ9rZTZr2F4sWz/MPIHPE5880a0aNPIEskpoK9tUCHiTwaGvgaGvgafpwNaXwg8V+mxLE3G/RdKrI+LS9YKIeHgAn7pz9HFwxnnt21cJf5HAF8UmFu013S4S2UJh3VdQXiHcebbz7ObZXP+bZ7eOrbL2RXbyBpDFv4uCnHcGUdW2pM6GHBtrRW1qDA51BoY8waozyOQKfs6ZQ66A54hqnqCW+XaY7QqbVYNAYd2NBeVFZwQP0ZlVn7mfBJwv6XXAzcCtwC0R8f5BPFugTta0Axz1WDjt2Qe321z1K52B9zDbLptp1xHnOmKYtZtnM6/+ot1CYa0Q50phLmi/tP7XFYL6A0LdwaBof077DQvba1k7mX1Vg0CdM4MiG7l1Ss4GcsU/R1CrfFofBPfv58DGgbYLwr+28XB/D7RfGFjWsse2lilfO2h3bWPmOHoQ95WOuUfES+BAKOYpwFNJ8h8MLu51sqYdwsaj4ZtKlstUhVbWaTsTXywrC6EUzsRrzJKz7xdFLWuvqZAeJogV7csEuYkY1xHyrn02Ef+mIh45ZR2EP/t3w4KdtQUBrSP6G1gQ17VMnbXExlqOjfXls2tZu6kfyhN8Hdouuzw4T/APsbEhp2ztUHvZevv3w9qGQ+vt35eUZeuu95s9BjJi35WK9AOj0DgrZBqKuTl9LYs6WdMOsmEjHPv45H2b2HpRm7zyMuE+LO494Iy802y8gXjmta0j+mXCXVR/jYN/Yy9wdOpjxt+Nxxwu3ovt9q//6jfm18/2rfTvhkyd7CwfDl7t3k8ithsy+5RTh0ydtYU22bbr+yNn/+Lf1M6GTPv9qWBH2j4rfBvWBW4t+fi04aDIHXB1fXtdONPytUzdxY/gABkBXbSXrZy3UGD/wva6v9l9sR82pJ/t/g0V9TIif2CAyfiUHSQ29HVrcbeskIPQ5A7VT5HelZq+PgVcGBG/NpBvWSqzpknaDmwHOP3009uJOrRv1wdD9D3m8cDhg58xGdpETKdIAPsmdixNZu7PJ7kj9WyS3Ad/BNwBLEPcK0lTaO4A2HrO04NvPFhQseeQTF5Zk1h73fh49v0qhGWaxMjrxOf76rNOWGZxpl4nxDJAPL6Onaax+KI4fB0bjePwUbyv7qqhvAuwlfWaXHztMSwzsahMo5j7PwDXp6/1J4P88iBeHU6drGkH2fsQfOXu9r2VzXbrhmzKBH5xf9PwTVE8v0n4Jq9+nt06cfDs3ylfUM3G06G5EI8k4FmXm4p4Xp1BLqgOJOR5vhX512UpZR9MLCrTKCzz5Ij49Pp2RHxG0tnDuHUYdbKmHeTRf4L7CjNhFlMnhNH3MsiqOnUHhUNEu8FSyKL9TZYvZt/XWQpZtnKnzK+yFS8H6i78rRLqvH0lIltbsGvYqzNrzpYVzb4PaddwbXzZUsg8/6qWFrZZCtlkPfwQa+H7YKVn7sA7JX0bibjeBhwD3C7p2Ij4p0G8SynKmlbY4OGvwz3X9OhAjX9b6Wy/Q4gnr6zLDUyL/hSJ+2I/VYPAYXWyv6SFv4vvy9aTQ7kQZ99X3QyVZ7OqXontujccZctrz2Jr3hSVLa+627XqJqYmNwL1fXdrU9Eua9P4JqaeGOMhIWXUecyeIuEF6fbpJBkhn57+vVkSEfHtQzqalzWtkIe+Dp9ZUk6zJl+Sqrpl+3tIPVAr81Je2eJ2l9QDi36U1S3bV+JD09QAi/umln6gyZ2ldf0v7a+nu0zb9LHYpsx2mVgfso+FfRTva0tQ/BMdi1qP2ZP0QeDKiPhCRHwB+IKkq0jSAF8I7BrSycY8BNw1oP0BRn1oN5tYRsKwqvZNkoYVlefZr5M0LM9erTqLGRha2C1LHLZYtyphV11RhPIZbGXbEj/KBHKx36rPZqjEYVVpiUZLHMZqivs24CeBP0pj3l8lCcmsAR8D3hoRtwzlYBse+QZ8rjho05q+c7X3wRj5LNrGK5tcwKpTt8uANWS6X8j/jOoOaoV1K0S9yF4dX2oPuCXCe6DOop2WdeoIcZ06S0n5ywqKe0Q8BPw28Nvps1I3A98oezDr2Dz8EHx+QdzHfkhFnwPDEMcy9ufTF11XQDRpX3dgrTwrqmrfsm2bdkXlTR6GUVa/bmSwaZ991O3Kyol7lvRZqfcP5EtvPAJ8YbGw5FfQ101qpfqY078ft1efmYw9QPeZY1NxqttfHbt1bFXZaTsgdWk79AluVPQ/Bo3TD6wCj3Do7azr9PkQ4z4flN32wdB9PiS7TZu5PCS7id0h+h6KIcSmrUgONSA1sT20wK/yUsiVYS/wpQHsNvnR9i18TYS3Sd267fPK6rat22+T9kV1q/Z1HTyGOqtpav9Ioo9Bahkz9y591M16K+lHgA8Az4yI0oUsTW5iOmsxC6Ok8yLi+ro2lsWjwAM92WryY2srHH3MmNsKdx1bberU8adt32MONGXlZbbqtK3Tvmm9Jn13YexBaeyQSJcLqnWz3ko6Dng1yWNOK2kyc79C0nuAt5CslnkLsBV4TgMbS2EvsKdD+75jyW3EoOtMvU69IUS8qn5TH6rst+mjbZ2ien38D8vqV+0bK5TVtk0Zq3xtpcPMvW7W2/8CvBn4xTpGm4j7s1LD/xc4Dngv8NwG7ZfGPpL1mkX08YXs4wfTNDzTt+C3mSU3FeMu+4fue5kDRZPwWR+hqbbXPKraNqlTp6+2dofouwsVM/fNkrIhlB1posN16mS9PQc4LSL+RFLv4v4o8A2Sx08fA3w2YmoZjBP2Af+vZ5t9xdvr2Fum6PclVkPO0Pu2vez9besUlXUJVRXVLatfta/q+973mcNUr2eUiPueiNja1q6kDcBbgVc0addE3G8ErgSeSbLW/Xck/UhE/NsmHS6D/UCfyW6G/DJ1OaVeRuimqN4UZv1dBphlh5fq1Bk6fDTkANBW/Ps4ixgy7FSXjjcxVWW9PQ74LuB6SQDfDOyU9KKyi6pNxP2ijKH7gQsk/ViD9ksjgGU+uXuoi1tjx/KHXknTRhCHPmOYguj3MXDW7atrvaK6fdYva1OnbZ32fdBB3Euz3kbE10gm1ABIuh74hd5Wy+QZioj31G2/TIJkrftYDLmSoetMps0Mq8mPse8Lk32s7mkz0HQV8Tp1hrgmUKefItvLCgEt6/pBXxea69DlJqairLeS3gjsioidbezOcp37fpILBFOhy5dombH+KhtDz9CGnn0uUyyHuG7Rly9t/Wliq0u/Ze3LytvY6pMuNzHlZb2NiDcU1D2vjs1Zijv0d7dYH6dzdU7X+hgA+uxniEGgbbshBGWK1yDq9tfXQNanT03qNWlfVt7nGUBXut7ENASzFPcpftB5ZL+EXZMO1f3iLvbTZFDJ9lHkb9fY5yqcOTSpO4Uw1RD1hr4eU1a3id2y8rJ+m7LqT2JaKca+Y60OVYI+xExj8ctcd1BpOni0PUNoOvgop02dtmMMFEX1lzWojCnIU5/V98HUNGe24r4KM/cqhvgi5glh10Ekz8++Bq59Bfbr2Mv6UPeMIs/vvmO8y774OFTdriK87GMrK++DlU75u0oscxQd6guzjAFqA90/q8UvdZfPo+lAUbevOsdYZKvuIFEWtqrTrsnAMqUzj6b1mwrvqszaLe4zZMzTsa4DS18DSPaH0+XzWP+BtD2uoh9w0Q+vbj/rdtsONHXOQLqcfRQNEFVthj5TGat+G1tdsLgXIOnXgX9Nsjz9b4GfWH/Sk6RLgItIPrv/HBFXjeXnFBliYGkjrH0MEn0MEPvof2BYt7vIsgaGdao+kz4Gh8U2ddtNbVBp26YLFvd8rgYuSRfzvxm4BHidpLNI7tZ6CvAtwDWSnhwRU/scZ0UfA8ayB4i5DwxFdura2lDSvspGm7OOKYay1tsNJR6+oJpDRHwss/lx4MXp+wuAyyPiYeCzknaTpMf86yW7aBrS9oveVlzbDgyLQtDW7776X6Tt9Yy6q6KGHhig3RkHmT6bzs7bDCxd8VLIevwk8L70/SkkYr/OvWnZYUjaDmwf1jUzNKswKOQJRhu/1WP/i2QFrslnU/eCdt1Boah9XTtVA0uVjaqzjr4Efor31ixN3CVdQ5LNbJHXR8SVaZ3XkwyA721qP82PvCO1M7UzJDMwbcW1DXMdEGAag0LWXtfrFF3POpowtVjx0sQ9Is4v2y/pFcAPAy+MiPXvfVUqTGNa01RcuwhBU0EuEuO2Pg85IEC7VU5NrlP0GY4aYumyV8sUkD4c9rXA8yPinzO7dgL/U9JbSS6obgE+MYKLxhwRgwFMe0Co6qPtmUcfHLFhmQp+CzgauDpNRv/xiPipNO3lFSTPEtwLvNIrZcyq0ERYl32/QplAdhnE2vjRtL82q5qGvIEJPHMvJCKeVLLv14BfW6I7xiydZQ4E0EyEq4Rx2YNBmz6XMau2uBtjOrHKAwF0879tqGjoFRaeuRtjlsqRNBBA/jEsKxY+tSV6FndjDFBfnPq6UNn3QADLH8yy/fomJmPMSrPsQQDGHwjq4NUyxpgjgjEGARhmIKhiijH3oVcIGWNMKVHzNQR9zbbXxT3vVQdJ2yTdLWm3pItz9v+8pDsl3SbpWknfWmXT4m6MWQnGHATqsL/gVYWkNeDtwA8AZwEvSzPiZvkksDUizgY+ALylyq7F3RgzK8YYBDrO3M8FdkfEPRHxCHA5SUbcg/Yjrsvcvf9xklQspVjcjTFHJH2LfIm4b5a0K/NazF57CvDFzHZh9tuUi4D/XeWPL6gaY0xHKi6o7omIrX30I+lHga3A86vqWtyNMaYjHVfL1Mp+K+l84PUkCRYfrjLqsIwxxvRA2wuqwI3AFklnStpE8mjRndkKkp4BvBN4UUQ8UMeoxd0YY3qgrbhHxF7gVcBVwF3AFWlG3DdKelFa7deBxwLvl3SLpJ0F5g6gg8/FmA9+EpMxpgE3dY2JHyfFMwr2/Z8e7LfBMXdjjOnIEf0MVWOMmTMWd2OMmRlTzC1jcTfGmB6wuBtjzMzwzN0YY2bK1MR9UuvcJb1GUkjanG5L0m+maTBvk3TO2D4aY8wiU8xSOZmZu6TTgO8HvpAp/gFgS/p6FvCO9K8xxkwKz9yLeRvwWg4d7C4A/iASPg6cIOnkUbwzxpgCuj6sYwgmIe6SLgDui4hbF3bVToUpaft6Ss2B3DTGmEKmJu5LC8tIugb45pxdrwd+iSQk05qI2AHsSPty+gFjzNI4ou9QjYjz88olPRU4E7hVEiTpLm+WdC41U2EaY8zYOOa+QER8KiKeEBFnRMQZJKGXcyLiSyRpL388XTXzbOBrEXH/mP4aY8wiU4y5T2a1TAEfBX4Q2A38M/AT47pjjDGH45uYapDO3tffB/DK8bwxxph6WNyNMWaGTG0Vh8XdGGM6EsDesZ1YwOJujDEdOaKXQhpjzJxxzN0YY2bGFFfLjL7O3Rhj5sD+glcdJG2TdHeaAffinP1HS3pfuv8GSWdU2bS4G2NMR7rcxCRpDXg7SRbcs4CXSTprodpFwIMR8SSSJItvrrJrcTfGmB7ocIfqucDuiLgnIh4BLifJiJvlAuD30/cfAF6oNF9LEY65G2NMd64CNhfsO2YhW+2ONNHhOnnZbxefW3GgTkTslfQ14ERgT5FDFndjjOlIRGwb24dFHJYxxphxqZP99kAdSRuB44GvlBm1uBtjzLjcCGyRdKakTcBLSTLiZtkJXJi+fzHwZ2nurUIcljHGmBFJY+ivIonbrwHvjog7JL0R2BURO4HLgPdI2g38A8kAUIoqxH8l8ZOYjDENuCkito7tRN84LGOMMTPE4m6MMTPE4m6MMTPE4m6MMTPE4m6MMTPE4m6MMTNkMuIu6Wck/Y2kOyS9JVN+SZrm8m5J/2pMH40xZlWYxE1Mkl5AkvXsaRHxsKQnpOVnkSzWfwrwLcA1kp4cEVPLi2+MMZNiKjP3nwbeFBEPA0TEA2n5BcDlEfFwRHwW2E2SHtMYY0wJUxH3JwPfkz5h5M8lPTMtz0uFeUqeAUnbJe1aSK1pjDFHJEsLy0i6BvjmnF2vT/14HPBs4JnAFZKe2MR+mh95R9qX0w8YY45olibuEXF+0T5JPw18KM1y9glJ+0kS39dJhWmMMWaBqYRlPgy8AEDSk4FNJE8Y2Qm8NH047JnAFuATYzlpjDGrwiRWywDvBt4t6XbgEeDCdBZ/h6QrgDuBvcArvVLGGGOqccpfY8yRjlP+GmOMWQ0s7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Ms7sYYM0Om8iSmvtkDfH5sJ3pmM8lxzY05Htccjwnme1zfPrYDQzBLcY+Ix4/tQ99I2jXHp8XM8bjmeEww7+Ma24chcFjGGGNmiMXdGGNmiMV9ddgxtgMDMcfjmuMxgY9rpVBEjO2DMcaYnvHM3RhjZojF3RhjZojFfSJIerekByTdnil7nKSrJX0m/ftNabkk/aak3ZJuk3TOeJ4XI+k0SddJulPSHZJenZav+nEdI+kTkm5Nj+tX0/IzJd2Q+v8+SZvS8qPT7d3p/jNGPYASJK1J+qSkj6Tbczimz0n6lKRb1pc9rvp3sA4W9+nwe8C2hbKLgWsjYgtwbboN8APAlvS1HXjHknxsyl7gNRFxFvBs4JWSzmL1j+th4Hsj4mnA04Ftkp4NvBl4W0Q8CXgQuCitfxHwYFr+trTeVHk1cFdmew7HBPCCiHh6Zp3+qn8Hq4kIvybyAs4Abs9s3w2cnL4/Gbg7ff9O4GV59ab8Aq4Evm9OxwX8C+Bm4Fkkd29uTMufA1yVvr8KeE76fmNaT2P7nnMsp5II3fcCHwG06seU+vc5YPNC2Wy+g0Uvz9ynzUkRcX/6/kvASen7U4AvZurdm5ZNlvS0/RnADczguNLwxS3AA8DVwN8CX42IvWmVrO8Hjivd/zXgxKU6XI/fAF4L7E+3T2T1jwkggI9JuknS9rRs5b+DVcwy/cAciYiQtJLrViU9Fvgg8LMR8XVJB/at6nFFxD7g6ZJOAP4Y+I5xPeqGpB8GHoiImySdN7I7ffO8iLhP0hOAqyX9TXbnqn4Hq/DMfdr8vaSTAdK/D6Tl9wGnZeqdmpZNDklHkQj7eyPiQ2nxyh/XOhHxVeA6kpDFCZLWJ0xZ3w8cV7r/eOAry/W0kucCL5L0OeByktDMpaz2MQEQEfelfx8gGYjPZUbfwSIs7tNmJ3Bh+v5Ckpj1evmPp1f2nw18LXOKORmUTNEvA+6KiLdmdq36cT0+nbEj6TEk1xHuIhH5F6fVFo9r/XhfDPxZpAHdqRARl0TEqRFxBvBSEh9fzgofE4CkYyUdt/4e+H7gdlb8O1iLsYP+fiUv4I+A+4FHSeJ8F5HEMK8FPgNcAzwurSvg7SRx3k8BW8f2v+CYnkcS77wNuCV9/eAMjuts4JPpcd0OvCEtfyLwCWA38H7g6LT8mHR7d7r/iWMfQ8XxnQd8ZA7HlPp/a/q6A3h9Wr7S38E6L6cfMMaYGeKwjDHGzBCLuzHGzBCLuzHGzBCLuzHGzBCLuzHGzBCLuzHGzBCLuzHGzBCLu5kVkh4j6c8lrfVga5Okv8jcfm/MymBxN3PjJ4EPRZLYqxMR8QjJXYz/rrNXxiwZi7tZGST9VPo0nVskfVbSdTnVXs7BPCFIul7Sd6TvT1T6pCtJ75f0W5L+UtLnJT1P0nskfVrSZRl7H05tGrNSWNzNyhARvxMRTweeSZJ/J5uMjPQRcE+MiM9lip8EfDp9fzZJvhCApwL3RMTzSB7QcBlJLvOzgB+SdHRa7/a0P2NWCou7WUUuJclC+L8WyjcDX13fkPStwH0Rsf7wibOB2yQdA5xA8nAKSJKbXRYR90fy4Il9wCNwIG/7I+uZBY1ZFSzuZqWQ9ArgW4Ffzdn9DZJshes8jSRz4zrfnW4/Bbg5I/pPI3lCFJJOBf4uDs2odzTwUB/+G7MsLO5mZZD03cAvAD+aEeYDRMSDwFo6M4fk4dXHpG23ABeQhGWeSpICdp2zOTgIHDIgSDoR2BMRj/Z6MMYMjJd4mVXiVcDjgOvSR/Xtioj/sFDnYyR55K8hEeqHJN1KIth3kjyY4QSSHOSkA8Fj0oEBDhV6gBcAfzLEwRgzJM7nbmaFpHOAn4uIH5P0GeCciPjHDvY+BFwcEZ+urGzMhHBYxsyKiLiZZGZ/fLLZSdg3AR+2sJtVxDN3Y4yZIZ65G2PMDLG4G2PMDLG4G2PMDLG4G2PMDLG4G2PMDLG4G2PMDLG4G2PMDPn/CKSyZGJ1QPQAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u_wpm = Scalar_field_XZ(x=x0, z=z0, wavelength=wavelength)\n", "u_wpm.incident_field(f1)\n", "u_wpm.WPM(verbose=False)\n", "u_wpm.draw(\n", " kind=\"intensity\", logarithm=False, normalize=\"intensity\", colorbar_kind=\"vertical\"\n", ")" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "**Differences**" ] }, { "cell_type": "markdown", "metadata": {}, "source": [ "Here we show the relative difference between BPM and RS" ] }, { "cell_type": "code", "execution_count": 5, "metadata": { "ExecuteTime": { "end_time": "2020-05-02T19:20:47.822318Z", "start_time": "2020-05-02T19:20:47.381496Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAX0AAAEKCAYAAAD+XoUoAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAulUlEQVR4nO2dfdQlRX3nP9+ZYWZAlJcZRQR0IAy6EMEQRF018Q0zMVkne0LiZF1DEvZwkkjWbJKjEHc9xhN3g3uOxixGMyvsuugJEnxhomSRV3OSTZABeUd0RBQIigPDKC/DzDPz2z+qap6anu57u2/fl773/j7n9Hm6q6urqu/T/a1f/aq6SmaG4ziOMx8smXQBHMdxnPHhou84jjNHuOg7juPMES76juM4c4SLvuM4zhyxbNIFGAWSfEiS4zh12Wpmz510IcaFW/qO48w73510AcaJi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c0SnRl3SopMslfUPSPZJeJelwSVdL+lb8e9iky+k4jjOtdEr0gY8C/9fMXgKcAtwDnAdca2ZrgWvjseM4jjMAMuvGGuKSDgFuBY6zrFCS7gVeZ2YPSzoSuMHMXtwnrW7clOM408DNZnbapAsxLrpk6R8L/BD4X5K+LumTkp4FHGFmD8c43weOKLtY0jmSNkvaPKbyOo7jTB1dEv1lwKnAx83sp4AnKbhyYgug1Io3s41mdto81diO4zhN6ZLoPwg8aGY3xuPLCZXAD6Jbh/j3kQmVz3EcZ+rpjOib2feBByQlf/0bgbuBTcBZMews4IoJFM9xHGcmWDbpAhT4PeAzkpYD9wG/SaiYLpN0NvBd4FcnWD7HcZyppjOjd4aJj95xHKcBPnrHcRzHmU1c9B3HceYIF33HcZw5wkXfcRxnjnDRdxzHmSNc9B3HceYIF33HcZw5wkXfcRxnjnDRdxzHmSO6Ng2D4zjO1LFu3TrbunVr6bmbb775KjNbN+YiVeKi7ziO05KtP/whm2/8x9JzOuDA1WMuTk9c9B3HcVpjsGdh0oWohYu+4zhOW2wPLOyYdClq4R25juM4bUmiX7bVQNI6SfdK2iLpvJLzKyR9Np6/UdKa7NzJkv5J0l2S7pC0sldebuk7juO0xWxgS1/SUuBjwBmEFQRvkrTJzO7Oop0NbDOz4yVtAC4A3iZpGfBp4B1mdpukVcCuXvm5pe84jtOaVpb+6cAWM7vPzHYClwLrC3HWA5+K+5cDb5Qk4M3A7WZ2G4CZPWpmu3tl5qLvOI7TFosduWVbf44CHsiOH4xhpXHMbAHYDqwCTgBM0lWSbpH07n6ZuXvHcRynLb3dO6slbc6ON5rZxiHlvAx4DfBy4CngWkk3m9m1vS5wHMdxWtFz9M7WPssxPgQckx0fHcPK4jwY/fiHAI8SWgV/b2ZbASRdCZwKVIq+u3ccx3Hakiz9wXz6NwFrJR0raTmwAdhUiLMJOCvunwlcZ2GB86uAl0o6KFYGPwvcTQ86J/qSlkr6uqQvxeNj4xClLXHI0vJJl9FxHGcfWvj0o4/+XIKA3wNcZmZ3SfqApLfGaBcBqyRtAf4AOC9euw34MKHiuBW4xcy+3Cu/Lrp33kW48efE4wuAj5jZpZI+QRi69PFJFc5xHGd/2n2cZWZXAlcWwt6X7e8AfqXi2k8Thm3WolOWvqSjgV8APhmPBbyBMEQJwpClX5pI4RzHcapo594ZK12z9P8ceDfw7Hi8Cng8Nn+gfCgTAJLOAc4ZdQEdx3H2w6dhaI6kXwQeMbObB7nezDaa2Wl9eskdx3GGT8tpGMZJlyz9VwNvlfQWYCXBp/9R4FBJy6K1XzaUyXEcZ8JMzyybnbH0zex8MzvazNYQhixdZ2ZvB64nDFGCMGTpigkV0XEcp5wp8ul3RvR78B7gD+JQpVWEoUuO4zjdwd077TCzG4Ab4v59hAmJHMdxukmLWTbHTSdF33EcZ7qYHp++i77jOE5bpmjIpou+4zhOW9y94ziOM0+4pe84jjM/uHvHcRxnjjDvyHUcx5kjpsenPw0fZzmO43Sblh9nSVon6d64bsh5JedXxPVEtsT1RdbE8DWSnpZ0a9w+0S8vt/Qdx3Ha0mL0jqSlwMeAMwgzCd8kaZOZ5StgnQ1sM7PjJW0grDPytnju22b2srr5uaXvOI7TlnaW/unAFjO7z8x2ApcC6wtx1hPWE4Gwvsgb43ojjXHRdxzHaU3P5RJXS9qcbcV1P44CHsiOy9YN2Rsnzji8nTAXGcCxcYnZr0p6bb+SunvHcRynLb3dO1tHuM7Hw8ALzexRST8NfFHSSWb2o6oL3NJ3HMdpSzv3zkPAMdlx2bohe+NIWgYcAjxqZs+Y2aMAcQGqbwMn9MrMRd9xHKct7UT/JmCtpGMlLSesJ7KpEGcTYT0RCOuLXGdmJum5sSMYSccBa4H7emXm7h3HcZxhMODHWWa2IOlc4CpgKXCxmd0l6QPAZjPbRFhH5JK4rshjhIoB4GeAD0jaBewBftvMHuuVn4u+4zhOW1pOw2BmVwJXFsLel+3vAH6l5LrPAZ9rkpeLvuO0YKAxcwNiY8zLaYgZLDwz6VLUwkXfcQqMU8ibkJfLK4COYcB0TL3jou/MF10VdGcGcNF3nPEyL4LuVn4HMUI36hTQmSGbko6RdL2kuyXdJeldMfxwSVdL+lb8e9iky+pMDvXYZh3DBX/YDO25Se6dsq1jdEb0CT/PH5rZicArgXdKOhE4D7jWzNYC18ZjZ4aZZ2Evw8W+Ob2eoZE9T1Mi+p1x75jZw4RPijGzH0u6hzDfxHrgdTHap4AbgPdMoIidZFo799q8cOO0VCbRYp+m/+O46WzF7x257YhzRf8UcCNwRKwQAL4PHFFxzTlAcSKjkZEevn4v6LBEud/DPmjaxXSHLThNXtIuNTsTqUyjFv95F/rOinldXPQHR9LBhI8Nft/MfpTPHho/Oy59P8xsI7AxptH3HSoTO1Xsl11nhbDicTHtqryrClrnJagrFHVfqFFWSl0U9LqMQvDnReSnXsybMCUduZ0SfUkHEAT/M2b2+Rj8A0lHmtnDko4EHmmcbsM4VfsJY3/hropXlU6vigCqLcwmLYs6NBWfqvTrivqoRWCYYjqsd3jWBH6ahLzOc7l7GBlNkaXfGQMsLghwEXCPmX04O5VPNHQWcEWjdBuWo+wHKYYV01xaiGNZvGKHUeqUy8+l6/MN9hWdss68QTumjPL0+qVdVs60TaTjrKLcbdiTbYNgJdu00JX/YZGqZ67ONjamaPROlyz9VwPvAO6QdGsM+2Pgz4DLJJ0NfBf41boJ5g9okwdgCeGlz8W3KMTF9HKRWFoSngt9urbMzVMUekriluVfVo6ydIrUSXfcrps2FnZTkR1nXpOgK1Z5Z6zLUTFFln5nRN/M/oHqZ/SNTdNLCVU9bP3890uz4zyNJYX4ewrhZeKdX1OkGL/YSqj7stRxA/WqBMvK1+RFHWaHcKp069Ikr0FEvmviPikhnwbhnmgl5z79yVNmURcp+u9zwSlWDCm93AeYW/VllUBVi6DqA76UXr9KIk+rGJ5fv6QkrBhedr4t/fo4qqjz3tQV4VFVHKNgXGLVNeHuSkukNQa7p8TS79ozMDSKgl/H/7e0EDftL2OxMtjNon87xc/99AcAy+O2NMs/XZu25DJakqV1QLZf7CNI254sraL7J127rCS9Zdl+1TZK32kvUe3nR+/nI99TsfVLZ9R+97o+8kH6nSbl425yT5PsHxi339+i6JdtdZC0TtK9krZI2u8DVEkrJH02nr8xDmvPz79Q0hOS/qhfXjNt6SdLu8p1sadkP7loRBDuJNJJnJewKNopXqoUim6f3exv/eeCXubTJ0s77ZeVPb++6phCeFuqxLGXWDeJ3+uaftf1u3YYDOt3HKel1UVLehYtTWth6ceVrz4GnEFYFP0mSZvM7O4s2tnANjM7XtIG4ALgbdn5DwN/Vye/mRX9XPCr3Dy5pZ778JeyKNoQxJ94vCuml6zoVEkkod7Nvu6fZP1XuXqssF8UtrwyKe5TclxkEDdL3c7gOv0I46oQBqGNII5KuCYp0tMgxl2sxIC27p3TgS1mdh+ApEsJMxHkor8eeH/cvxy4UJLit0u/BHwHeLJOZjMr+smqLoptbsnDvqKbfoz0v1sez+9isWJYnqVXdNkkipZ/TrGCSOnkZS6rrMoqriYvaS6mZR3XxTijOF8Vr1/8ugziJmnLOESoa2LcWeGtwSh/S6u2WFZL2pwdb4wfkyaOAh7Ijh8EXlFIY2+cuLzidmCVpB2EaWnOAPq6dmCGRT/5qYuCadlfWPStLwd2xrAV8Xw6Xs6itZ4qgQUWRTu5gopCXxT43LJP1+V9Bv1aJ70e2CYfcdXpDK4TXpVP02GjZdQVlqYv8bAFaxKCPA2i27WKatT0ce9sNbPTRpT1+4GPmNkT+ewFvZhJ0c87inLLuczCX04Q8d3sL/YHsegCWgCeiXGJ4SvY10ef+/HzDtt89E9epryjt1hhlL00dQS313j/4rDQfvHK8hlE0Pt12tURiEGEbhTC0xXBnVZR7crvN2za+PSBh4BjsuOjY1hZnAclLQMOAR4ltAjOlPQh4FBgj6QdZnZhVWYzKfqwOHoliSzs33kLQcSTG2dnjH9wvGY3sIPFSmAp8Cz29eXn7p3kBkqth9QCyEfHlFnzZfQajpn7/qsEvdiq6Be/mGcxv0TZSztIv0LOIAI2KvHoqpjOqlh29fduigELg4v+TcBaSccSxH0D8O8KcTYRZiT4J+BM4DozM+C1KYKk9wNP9BJ8GED0JT0L2GFmQ5myYhQkoU/DFHNLPwlwsthXEoR9KaHqFMGif4Ig5MtYrARSpbFAqAgWWBT6fMhk3slbNVoH+gt2cUx/maD32s/TzMPqUGcEUK8XtolIDfPFnyZxnBXBGxbT9L/bD+vp0+99afDRnwtcRZCMi83sLkkfADab2SbCFDWXSNoCPEaoGAZCobLoEUFaEjN4O/BygiauALYCXwb+ysy2DFqAUbBUstUE8V3G/m6Thfg3Cf+B8e+TwNPx3AqClS6CWO7MtlTbpYol5ZNXMFUCnwt4k7/Fa6Hami9S7MiuOld2vixOv/i96KrQTbXgTDFdeB62w81tfe6nHCS78iXl547+evv0h0kdS/964BrgfOBOs1CfSToceD1wgaQvmNmnR1fMZqRhkmnLfe7pg4xdBMt+D/BjFt08zyWI9wKhBfA0QehTKyFVBknoi66aXKAX2PejqtzfX9bB20vse90rLPYRFMPK4lJxvipeL7rw0jr98UqtnGENAzam54vcOqL/JjPbVQw0s8cI0yB/Lk6J3BnEojgnl0v6knYXoYN2CfA4QZgPJlj7Segfj38hiPshBDdQsuaLI3RyMV+g/Ovb3ZSLftF3X2yV9BuX36vzt+pFd6F25oE63pahfbjYriN3rPQV/VzwJR1QUQHsFzZJRBDpNLoGFgX2YOBHBOFdRfgBnga2AU/FuCsIy3OtYLEzOE9jV/Y330/DOIuCXxwemsqY3ENV4/Gr3DKDfHE7LKGfxJxSo/7StgtMyVxdnSd/VsZq3MyS6CckfRL4BUkLwL8AtwO3m9n/GFXhBmUJwXJPom8sivsOggsHgtD/kPDCHQgcHv/mFcUugnsnCfwzhbCyj7PyoZlVY/DLxuL38r0Xz7ehjcC0EeBhC1tXKoNh3NfS/lHGQld+00HYw2TdWIN25I6bJqN3XgscbWa7JR0FnAKcPJpitWMJwaJPHbFLCKNxDovnf0Dw468AXsjieHwjiHqqHNLfZ9hX5BOp4zaf2KxM3HOBH/a8OHVe0n7P4jBntuwVr0rYhvWujEOwyso6TsHukih3RePSbzJJt+VMuXcybiR4RB4xs4cI40mvHEmpWpKLfnKvPJ/QPHma8AXDsQSx3kkYtbONUBE8wb6dt7Ao6KnlsIzFUUH5iB0ot9qbUvcDqLZfw6a4ZaI1jMqkbpzib9VG2IYtisXyD1vgJ1VRDYNBn5thkU+SOGlmVfT/CviqpIsIFcDtZrZ9NMVqx1LgOSx2vD5G+HTtRcCzCdb7Y4Ra63GC0KcPsJYRXDyHsO9QzOJwzLqiXvWRVZ0pD5qcS8d5uapewLwyq0qrLH6RdH2va4vrDTRhULFqK3J5OduI/KgroLr0uodRiPQoWwAT89n3Y0ZF/9PAxnjN7wInS1ppZj8xkpK1YCnBlWMEUX8RoWP3+8A9BD/+0yx2+D6PYMWvYP/x9rD/cEzo/UFUHtbrA6z0d0nhXL5fNd6/7FydqRuqwpbQWwD6uW2G6SIiptdUcAdt5rcRqV7X1rFAm/wmVb/HMIR7mEI9TDEulmsc7tCB0p1R0X/QzP5bHiBpxZDLMxQOIFj6K4AXAN8AthCE/kDCBBYHsTgMMxf4NEInCXeTj6jySqJYAeRTN+/JjvO4ZPv5F8Sw78NadMv0agGUWdrFl7LXC9+vCV3WaqiK04+qcrdNtyz9nH559cqnSuDG2UqpmjW1CW2FOi/3oMJcVv5hV+CjdAXt6UonRx+aiP6tkt5lZh9NAWb2zAjK1JrlBJ/9XcB9BHE/leDaSTecxsunztmy4ZZpyCUM1xJKJGEvzglUFNqy8fs5ZX7nfn78YvyUT5E8r0F9+GUvWllew7Bk67qaBkm7X/pleTVNv2gE1GXQkStVLcq61yUGneY7p24ZBqmE65ZhUGbV0j8CeJOk9wC3ALcBt5rZ34ykZAUkrQM+SnifPmlmf1YVdznBh38C8K9j2NPZ9hRhRE6+slWqANJ++shqTyFevy9ni5VElZsnDytSfPFViFu07MoqgUE+zOr3IvSqTIrp132pmnxAU0dABxHbJgLQVNyaiGlT91Re7jrX1K2Qc6p+x0Gvg3bPYL+86/7fh94fMIuib2a/CntdOicBLyVM6zly0a+5nNheDnk2/MdXwYPfhIfuD0M0nyZ02D5JEPw0e+ZCtuUWftmXs/lfSvbzv0Wxz8/1Cq+iKPwprJgW9LbOe1nddV7Gui97v5cqb9n0o2neTeM2EaEmgtNU3Pq50Io0dS0N4opqWqZe19UR5Lr9UnXzHddoImMGRT8RXTq3xG1c1FlObC9LDl3JYz/Ywf33w8OE4ZhPsGjhV31cNW6xL57rFdaLsgohUcfyLk4r0SudRN0mehPBqitWbfNuErdJ/0eT/Ju4zer8fwZNe5D4MFg/wqB9D01bj8PIszHG7H2cJekO4le4cbsDOMvMPjiisuX0XU5M0jnAOQAvPOwAtGr/RLowntdxnNnDgN2dnWx+X5q4tn4W+J8ET8kG4E7gLaMo1CCY2UYzO83MTjt8YReHroY1a8IXty8AVhOGcT6H8OHWQYQFUdIcPcvjlj68yr+yzbc0Rj/fzzcVNrJzFMKh/B/Q1N/Yy4pJ/RK9rs23fumMKr1+I4iGmXebcg4jzWK6/aib5iBpDxI/L1MTBrkGmpdtGHkOks9CxVYHSesk3Stpi6TzSs6vkPTZeP5GSWti+OmSbo3bbZL+bb+8mvj0HwNuiBuS1gL/ue71LamznNhefvRj+ItrgsAngz/Nmb+U4OZZRnDxLGNxMZR8FM8Ci+PPqzpylR0XOzmLD2lxlE7ukqkzYqPswW0yQqdX/H7XNUmjSVpN4jXJu2ncUZWhqVCNaiTRoPFhcMEc9Lq2HpJx+fDLGLTsNfsszwa2mdnxkjYAFwBvIxjfp8WFWI4EbpP0t2ZWWd80ce+cYGbfTMdm9i1JJze6u8Gps5zYXnYSLPt7CEvRrCSI/3MIFjzsO1Qtn4s+hYlF334S5ST+qvhbJf55nCrhp3BNjlXsV8WvilcVt981da5tEqdOXoPG7ZrQDmrZjuMaGL+oQzeEfdju92TpD0idPsv1wPvj/uXAhZJkZk9lcVZS4+dpNA2DpJ8giO7tMYM7JT3LzJ5skE5jqpYTq4q/E/gOQeRfAXwTuJXgl1pJGK9/MM0/zho2xc7gfL+YX79WBD3OV8XrF39Ylv+wWwdN0x00/UHyaJNXm+ug3TPa9vkehoh2Ucyb0KP8qyVtzo43mtnG7Lhvn2UeJ+rhdoItu1XSK4CLCZMPvKOXlQ81RD/WJmZmr4/HLyTMsPmy+PcWSZjZi/ul1QYzu5KaE7ztIsyZb4RpGFYDP0OYhuF7hFrr6Rh3JcH1k3z7ZdMw7FeW+Lco0PlInuIonl5xi+mWxSuL00/8y+JUhVVd3yu8TppN02mSXpv02+Y1jGthOELXJbEctnE0Cot8FBj7zsBbYOsol0s0sxuBkyT9K+BTkv7OzHZUxa+1XKKkzwFXmNn3zOx7wPckXUWYbvksYHPPFMbMbsIwzWUEIf8eoRI4AngxoTp8jLDI7zbCXDz5hGsHsG/H7qQmXGtq0depBCbl9qmTziBpDuu6nGEJQ9cFEEYjgqOytifpr69Di8E7dfosU5wHJS0jzAn5aB7BzO6R9ATwk/TQ5Dqivw74LeCvo0/9cYJhvBT4CvBhM7u1RjpjYzdB5POplVcRxP9pwq91GEH8jyR8sPUjwlj+NJ4/tRRg39E7ywrHxbVp206t3MQK77p/v0m8uvmO+voio3QXjEvExuXyGKcod21IfB9Lvx91+iw3EQzsfwLOBK4zM4vXPBBdPi8CXgLc3yuzOssl7gD+EvjLuBbuauBpM3u8yV2Nkz0E8c4XUXmM4Md/DuEL3QcJrYBDYvhhhHn20yIq6W/a30n4ire4QlZxOGfeIshbBil+/heqXUh1aSvOw3bLjKvjchj5NqULQtMla7cLvwd05zcZVPSr+iwlfQDYbGabgIuASyRtIcjZhnj5a4DzJKVZZH7XzLb2yk9mXfnJhsdKyX6S8uUSnyF05C4huHa2E/5ZK1ns3M1H+KS1b9OXvHWXS4TqlbQGXS6xGGdQxj3So22+VXTlye2K+LWlK7/noAz6f3gUbm7rcz9Osv9ace7XhpD+MGk8DcM0sIfgsklDKGFxFM5KgiNsF8Hqfz6hMtjOohNtBaFz90AWF0ZPHb0HMbyF0Ytz9veqCHp90FW3tTDM6XPHxbQLUT9mpcLoCpN8XqZk6p3ZFH1jsWM2uV52E8TxKRa/wH2U8I96FsG9k1bVeoLQuQtB9PMvdnPrPYXnc/ak2Tlz0a+a2yefujmVu0zcm7QK6lYInVp1yHGGzLgr02RUTgNNPs46sTirpaTXmdkNQy9VS4zghtlDEO18wqYlBOHfRRD75P/fRhDwAwm+/QUWF0ffzqIgJ/FPSynmFnqqEFIZil/zlk3o1mtyt6ox+zlNKoE6K22VxeuFVx7Th89BNfwWQcuO3LHSxNK/TNIlwIcI3o4PAacBrxpFwdqQfPFpvyiCC/Hv9hiWxuk/SbDwcyv+2QTx3Rm3ZwiVBiyO5snX0i123ubz8KT+hV4i30/8i+6itF9sNeT0sv7rtAx6CfsgAtLVisLFcHJ09ZlowsxZ+oQvxC4A/h9BCz8DvHoUhRoGyb2TrPt8ygOxWCmkDt0lBPFfzqKw/5jwAyW3zoGEyiDN1bPA4qIsSdSLW7HzNpG7bNIcP0UxryP2xX1KrqVwPqfKOhllZVDMc9gv/DSJ9yyI3TDIn4lp+v8lZtXS38XiMrMrge+YdXMG6dy/vpvFztzisoQQXDjJAt/OYqftwfHa5OMnxkuzbyYXT6oE0spbyaefRL5fBVA1y2axEzpffLxM1Ac9zsPy8Krz+QtZV7z7vcS7e1xbxaiEYdwi3EQoplEM61D8zae1834WRf8m4Arg5YSx+p+Q9Mtm9isjKVlL0j8gCboIgpxenIXsfBLtZTFuEvnl2bZAaD08mV1X3Jaxr99+d7wm5Z9XBMU+gGJLoJ8VXZzYrbhYRHHhjWKlkSqUorhXhVMSXjzf66HvJVhL+lzb7/pe6TalzovbBfGd5hZC8ffrpOXYkFm19M82s/Rp78PAeknvGEGZWmPZlkQY9p8VE4J1ngQ4rfKexP/p7DgfxZOP3U/pFv35ycefu2pSWRay6+qO4Yd6L3pVC6LYcijGHyS8eC7Pq0gv663sZRnGMNSql3CYoj0p8e26UPb6XabVku/HzIl+Jvh52CXDLc7wSP+A4vTFufDDYsWwi0VRLIp/6sDNrfo0xUNu0SeKE7blbp5iB23eAVvl/mlTCZSR33ueRi4kZYug5xUHFeehvHIpi1fEKvZ70cYXPMkP3ZrSNRdIv9+6y5XSKCrqZGROAzM5Th8WH7piB24ZufsHFi30fKw/LFYORTdNainkfQmJMuu9KOB5RVD1dW9VX0C/SqDsnqseznxoa/HaYgWRx+uVf1m8sgqlKm5Ov5eq6UvXRpjSvQ77Ra96RidVuVQxSYFr21Ib1W/pH2dNmCSkeedtTrFTNxfa9CFX+qhrJ/sKfS7uRddO7tZJ7px8pFCZ1d5vJE9eLtj/+l4fbuXhbSi2mPJ8yioJsvhllcUgH431e1mbClHT+HmZ2wpH1X2OQkyb/P/HUbm0tbS7aFHP5MdZ00j+jyg+aEWRh32FLbl1ilZ83ilLTL/MRw/7TwFRtOSrxDpPu+gq2cO+lYJqpJeHUXKuLVUVa1W/QK8WQ6+WSb+WRD9xaysWbSqJIsNoZdRlXK2RukxDxdKUvK+u68ys6CcrP4l03QetKGBJYHOhz4U7/xCrl1Vf5dLJ2V0Sj8JxlQ++OHogVQZVwl8mBKMclVLWUkj0cxv16zyu61bqdX+jsB6bpDkua3wYYjjM32pUz5zPE1XNNI/86ks+RLFsg33/UcXRPsV4Sexzoc/n1EkinyqO1AeQOmuLFm3eJ0AhveJ8PWUdRUuyTYUtH71UZyvOEVTMt99WpE6cKvo9lMV77ZVOcaubZp30h0mb37oJe1puw6bJfY/6t2lL1btVB0nrJN0raYuk80rOr5D02Xj+RklrYvgZkm6WdEf8+4Z+ec2spV+08Mte+DLhT+RW8p7Cfp5eMTyF5eF5ayG/RiVhZePjd7O/+BTvp9gCKLPq8/vsNeplSSFOP+FLLaph0VRc+rVgcvq1KvqlnzMpkRlVC6IOk25l9GOS/5NBh2xKWgp8DDiDsNTHTZI2FeY6OxvYZmbHS9pAmB3hbYQFAP+Nmf2LpJ8kzMl/VK/8Zlb0YV+hHcS9s6ckjMK5dD6PUybSxSGiVaNuerlvcsryqKKfe6Ru2lUvbZOyjJpebqQyBqkE8ryKTNraLDLJCqLIoBXGNLgj2og+cDqwxczuA5B0KbAeyEV/PfD+uH85cGFcv/zrWZy7gAMlrTCzZ6igE7+npP8u6RuSbpf0BUmHZufOj02aeyX9XNO0m76EZQ9mMaxMJMs6I8uanUX3i2XXF5vRRX98mXgP0ryt47ooS7usnHlLpgtN8Lbp13EF9WJS7qFh0JX/YZEuuaN60cO9s1rS5mw7p3DpUcAD2fGD7G+t741jZguEWWNWFeL8MnBLL8GH7lj6VwPnx2XDLgDOB94j6UTCsmAnAS8ArpF0gpk1qlTLHtIyAVXFfvG6ojVZdpz/rRL+PP1iPlUPbD9XRp17bZJenfSh/wuWi2jXrOEqylxugzANLYKmNDEwJsk4hb/Hb7J11CtnSTqJ4PJ5c7+4nbD0zewrsfYC+GfCavAQmjSXmtkzZvYdYAuhKdQ+T/a3XKr2y64rE+6ytGF/i6+OxVTHsqprTda12Nq8oP3S75pV1oRRvCTT2hpoShdbDqPACEM2y7YaPAQckx0fzeJCfvvFkbSMsLz3o/H4aOALwK+b2bf7ZdYJ0S/wW8Dfxf06zR4AJJ2Tmk8jLh9Q/2Ed1oPd72UZVDhG/fI1ecm72GwfV37zUAH0YtorB6PVs3oTsFbSsZKWE7wbmwpxNgFnxf0zgevMzKIr/MvAeWb2j3UyG5t7R9I1hCVpi7zXzK6Icd5LqBw/0zR9M9sIbIzpdPn5GCrTeqO9yl1H+LrcMmhLv473eabLbqVBO3KjW/tcwsibpcDFZnaXpA8Am81sE3ARcImkLcBjhIoB4FzgeOB9kt4Xw95sZo9U5Tc20TezN/U6L+k3gF8E3mhm6X9bp9njzBhtK4RZo8qN6FQz7srBaDfLppldCVxZCHtftr8D2G8aezP7U+BPm+TVCfeOpHXAu4G3mtlT2alNwIb4YcKxwFrga5Moo9MNprHpPyzm2f0zKob53EzaFVmXrozeuZAwXf3VkgD+2cx+OzZxLiOMV10A3tl05I4zP9R5gWdBNJt+i+CMnraW/jjphOib2fE9zn0Q+OAYi+PMMPNSMTjjx0XfcaaUJlb0OCsIt+67i1v6jjMnuBA7iWl5Flz0HcdxWpI+zpoGXPQdx3GGQBdH6pThou84jtMS9+k7juPMES76juM4c4a7dxzHceYEt/Qdx3HmDBd9x3GcOcEtfcdxnDlimkS/E7NsOo7jTDttZtmUtC6uA75F0nkl51dI+mw8f6OkNTF8laTrJT0h6cI6ebnoO47jDIFBRV/SUuBjwM8DJwK/FtcHzzkb2BYnp/wIYT1cgB3AfwH+qG45XfQdx3Fa0nKN3NOBLWZ2n5ntBC4lrA+esx74VNy/HHijJJnZk2b2DwTxr4WLvuM4Tkv6rJG7Oq3fHbdzCpfXWQt8bxwzWwC2A6sGKat35DqO4wyBHq6crWZ22vhK0hu39B3HcVqSRu+UbTWosxb43jiSlgGHAI8OUlYXfcdxnCHQQvRvAtZKOlbScmADYX3wnE3AWXH/TOA6MxtoCn937ziO47SkzTh9M1uQdC5wFbAUuDiuD/4BYLOZbQIuAi6RtAV4jFAxACDpfuA5wHJJvwS82czurspPA1YWnUbS7N2U4zij4ua2PveDpMqFvu8YQvrDpFPuHUl/KMkkrY7HkvQX8YOE2yWdOukyOo7jFLEeW9fojHtH0jHAm4HvZcE/D6yN2yuAj8e/juM4ncKnYWjOR4B3s2/luB74Pxb4Z+BQSUdOpHSO4zgVtBy9M1Y6IfqS1gMPmdlthVN1PlpIaZyTPn4YUTEdx3EqmRbRH5t7R9I1wPNLTr0X+GOCa2dgzGwjsDHm1UVXmuM4M0r6IncaGJvom9mbysIlvRQ4FrhNEoQPE26RdDr1PlpwHMeZOF206suYuHvHzO4ws+eZ2RozW0Nw4ZxqZt8nfJDw63EUzyuB7Wb28CTL6ziOU2SafPqdGb1TwZXAW4AtwFPAb062OI7jOPszTYuodE70o7Wf9g145+RK4ziOUw8XfcdxnDliWkaPuOg7juO0JC2iMg246DuO47TEh2w6juPMGe7TdxzHmRN89I7jOM6c4e4dx3GcOcEtfcdxnDnDRd9xHGd+uApYXXFu6zgL0g8XfcdxnJaY2bpJl6EuE59wzXEcxxkfLvqO4zhzhIu+4zjOHOGi7ziOM0e46DuO48wRLvqO4zhzhIu+4zjOHOGi7ziOM0e46DuO48wRLvqO4zhzRGdEX9LvSfqGpLskfSgLP1/SFkn3Svq5SZbRcRxn2unE3DuSXg+sB04xs2ckPS+GnwhsAE4CXgBcI+kEM5uWCe0cx3E6RVcs/d8B/szMngEws0di+HrgUjN7xsy+A2wBTp9QGR3Hcaaeroj+CcBrJd0o6auSXh7DjwIeyOI9GMP2Q9I5kjZL2jzisjqO40wtY3PvSLoGeH7JqffGchwOvBJ4OXCZpOOapG9mG4GNMS9rV1rHcZzZZGyib2Zvqjon6XeAz5uZAV+TtIewIMFDwDFZ1KNjmOM4jjMAXXHvfBF4PYCkE4DlhNVmNgEbJK2QdCywFvjapArpOI4z7XRi9A5wMXCxpDuBncBZ0eq/S9JlwN3AAvBOH7njOI4zOAraOlu4T99xnAbcbGanTboQ46Ir7h3HcRxnDLjoO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c4aLvOI4zR7joO47jzBEu+o7jOHOEi77jOM4c0ZWVs4bNVuC7ky7EkFlNuK9ZYxbvaxbvCWb3vl486QKMk5kUfTN77qTLMGwkbZ7F1X1m8b5m8Z5gtu9r0mUYJ+7ecRzHmSNc9B3HceYIF/3pYeOkCzAiZvG+ZvGewO9rJpCZTboMjuM4zphwS99xHGeOcNF3HMeZI1z0O4KkiyU9IunOLOxwSVdL+lb8e1gMl6S/kLRF0u2STp1cyauRdIyk6yXdLekuSe+K4dN+XyslfU3SbfG+/iSGHyvpxlj+z0paHsNXxOMt8fyaid5ADyQtlfR1SV+Kx7NwT/dLukPSrWl45rQ/g21w0e8O/xtYVwg7D7jWzNYC18ZjgJ8H1sbtHODjYypjUxaAPzSzE4FXAu+UdCLTf1/PAG8ws1OAlwHrJL0SuAD4iJkdD2wDzo7xzwa2xfCPxHhd5V3APdnxLNwTwOvN7GXZdwbT/gwOjpn51pENWAPcmR3fCxwZ948E7o37fwX8Wlm8Lm/AFcAZs3RfwEHALcArCF+rLovhrwKuivtXAa+K+8tiPE267CX3cjRBAN8AfAnQtN9TLN/9wOpC2Mw8g003t/S7zRFm9nDc/z5wRNw/Cnggi/dgDOsssfn/U8CNzMB9RTfIrcAjwNXAt4HHzWwhRsnLvve+4vntwKqxFrgefw68G9gTj1cx/fcEYMBXJN0s6ZwYNvXP4KDM5DQMs4iZmaSpHF8r6WDgc8Dvm9mPJO09N633ZWa7gZdJOhT4AvCSyZaoHZJ+EXjEzG6W9LoJF2fYvMbMHpL0POBqSd/IT07rMzgobul3mx9IOhIg/n0khj8EHJPFOzqGdQ5JBxAE/zNm9vkYPPX3lTCzx4HrCa6PQyUlQyov+977iucPAR4db0n78mrgrZLuBy4luHg+ynTfEwBm9lD8+wihgj6dGXoGm+Ki3202AWfF/bMIPvEU/utxpMErge1ZU7UzKJj0FwH3mNmHs1PTfl/PjRY+kg4k9FPcQxD/M2O04n2l+z0TuM6iw7grmNn5Zna0ma0BNhDK+Ham+J4AJD1L0rPTPvBm4E6m/BlsxaQ7FXwLG/DXwMPALoIf8WyCj/Ra4FvANcDhMa6AjxH8yHcAp026/BX39BqCP/V24Na4vWUG7utk4Ovxvu4E3hfDjwO+BmwB/gZYEcNXxuMt8fxxk76HPvf3OuBLs3BPsfy3xe0u4L0xfKqfwTabT8PgOI4zR7h7x3EcZ45w0Xccx5kjXPQdx3HmCBd9x3GcOcJF33EcZ45w0Xccx5kjXPQdx3HmCBd9Z6aQdKCkr0paOoS0lkv6+2waAseZelz0nVnjt4DPW5gQrRVmtpPw1ebbWpfKcTqCi74zNUj67bj60a2SviPp+pJob2dxHhUk3SDpJXF/leLKZJL+RtKFkv5B0nclvUbSJZK+KemiLL0vxjQdZyZw0XemBjP7hJm9DHg5YX6ifBI34lJ+x5nZ/Vnw8cA34/7JhPlUAF4K3GdmryEsnHERYS75E4FfkLQixrsz5uc4M4GLvjONfJQwq+PfFsJXA4+nA0kvAh4ys7QoyMnA7ZJWAocSFg2BMCncRWb2sIUFQXYDO2HvvPk700yNjjPtuOg7U4Wk3wBeBPxJyemnCbM/Jk4hzISZ+Ol4fBJwS1YZnEJY0QtJRwP/YvvORLgC2DGM8jvOpHHRd6YGST8N/BHw7zPB3ouZbQOWRksewqLlK+O1a4H1BPfOSwlT7SZOZrFy2KeikLQK2Gpmu4Z6M44zIXwomjNNnAscDlwfl1zcbGb/oRDnK4R5/K8hCPgOSbcRhPxuwoIZhxLmgCdWEAfGCgP2rQAAXg98eRQ34ziTwOfTd2YKSacC/8nM3iHpW8CpZvbjFul9HjjPzL7ZN7LjTAHu3nFmCjO7hdASOCQcthL85cAXXfCdWcItfcdxnDnCLX3HcZw5wkXfcRxnjnDRdxzHmSNc9B3HceYIF33HcZw5wkXfcRxnjnDRdxzHmSP+P6c7kxAHIit6AAAAAElFTkSuQmCC", "text/plain": [ "
" ] }, "metadata": { "needs_background": "light" }, "output_type": "display_data" } ], "source": [ "u_diff = Scalar_field_XZ(x=x0, z=z0, wavelength=wavelength)\n", "u_diff.u = np.sqrt(np.abs(np.abs(u_rs.u) ** 2 - np.abs(u_wpm.u) ** 2))\n", "u_diff.draw(kind=\"intensity\", logarithm=False, colorbar_kind=\"vertical\")" ] }, { "cell_type": "code", "execution_count": 6, "metadata": { "ExecuteTime": { "end_time": "2020-05-02T19:20:48.006772Z", "start_time": "2020-05-02T19:20:47.824894Z" } }, "outputs": [ { "data": { "image/png": "iVBORw0KGgoAAAANSUhEUgAAAZUAAAEGCAYAAACtqQjWAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAAsTAAALEwEAmpwYAAAkuUlEQVR4nO3dfVRUdf4H8PfAII86Diiog4owPiCimDNiVr9dYNX2qNPRNNlaspWiUjcXT5u1e45Zux213B7OoqflrBVtKp7cCnuQ1TB7IBWwsAQfAAF5UhAElIan8fv7w8PkBMoF7tyB4f06Z8483Zn7mdvBd5/7vfd7VUIIASIiIhm4OLoAIiJyHgwVIiKSDUOFiIhkw1AhIiLZMFSIiEg2akcXYC8qlQrBwcHQarWOLoWIaMAoKSnB5cuXe/15pw0VAHjppZcQGxvr6DKIiAYMg8HQp8879e6v9vZ2R5dARDSoMFSIiEg2Th0qbW1tji6BiGhQcepQYadCRKQshgoREcmGoUJERLJhqBARkWycOlQ4UE9EpCynDhV2KkREymKoEBGRbBgqREQkG4YKERHJxqlDhQP1RETKcupQYadCRKQshgoREcmGoUJERLJx6lDhmAoRkbKcOlTYqRARKYuhQkREsmGoEBGRbBgqREQkG6cOFQ7UExEpy6lDhZ0KEZGyGCpERCQbhgoREcmGoUJERLJx6lDhQD0RkbKcOlTYqRARKYuhQkREsmGoEBGRbJw6VDimQkSkLKcOFXYqRETKYqgQEZFsGCpERCQbhgoREcnGqUOFA/VERMpSNFTS09MxefJk6PV6bNmypdP7LS0tWLFiBfR6PSIjI1FSUmLz/oULF+Dj44Nt27ZJWh87FSIiZSkWKhaLBWvWrMGBAweQn5+PPXv2ID8/32aZnTt3QqvVorCwEImJidiwYYPN++vXr8dvf/tbyetkqBARKUuxUMnKyoJer0dwcDCGDBmC2NhYpKWl2SyTlpaGlStXAgCWLVuGjIwMCCEAAB999BEmTJiAsLAwyetsbW21fp6IiOxPsVCpqKjA2LFjrc8DAwNRUVFxy2XUajU0Gg1qa2tx7do1bN26Fc8//3yP1imEYLdCRKSgATFQv2nTJiQmJsLHx+e2yyUnJ8NgMMBgMEClUgG4MU5DRETKUCu1Ip1Oh7KyMuvz8vJy6HS6LpcJDAxEe3s7Ghoa4Ofnh+PHj2Pfvn145plnUF9fDxcXF3h4eGDt2rU2n09ISEBCQgKAG52OxWJBS0tLt2FERETyUCxUjEYjCgoKUFxcDJ1Oh9TUVOzevdtmGZPJhJSUFNx5553Yt28foqOjoVKp8PXXX1uX2bRpE3x8fDoFyi+xUyEiUp5ioaJWq5GUlIQFCxbAYrFg1apVCAsLw8aNG2EwGGAymRAfH4+4uDjo9Xr4+voiNTW11+tzcbmxZ4+hQkSkHJVw0sOjPDw80NLSgvz8fISGhjq6HCKiAcFgMCAnJ6fXnx8QA/W9wU6FiEh5ThsqHFMhIlKe04YKOxUiIuU5baiwUyEiUh5DhYiIZOO0odKx+6u5udnBlRARDR5OGyrsVIiIlOe0ocKBeiIi5TltqLBTISJSntOGCjsVIiLlOW2osFMhIlKe04YKOxUiIuU5bah0dCo8pJiISDlOHyrsVIiIlOO0ocLdX0REynPaUGGnQkSkPKcNFU7TQkSkPKcPFbPZ7OBKiIgGD6cPlaamJgdXQkQ0eDh9qPz0008OroSIaPBw+lBhp0JEpBynDxV2KkREyulxqDQ1NcFisdijFlm5uroCYKdCRKSkbkPl+vXr2L17NxYuXAh/f39MmTIFo0ePxtSpU/HnP/8ZhYWFStTZY+xUiIiU122oREVFoaioCJs3b8bFixdRVlaG6upqfPPNN5gzZw42bNiA9957T4lae4RjKkREylMJIcTtFmhra4Obm9ttv0TKMkozGAw4ceIEAMBisVhDhoiIbs1gMCAnJ6fXn+/2X1opYdHfAqWDp6cnAJ5VT0SkFLXUBXNycvDSSy+htLQU7e3tEEJApVLhhx9+sGd9feLl5QWz2YympiZ4eXk5uhwiIqcnOVQeeughvPLKKwgPDx8wu5K8vb1RW1vLwXoiIoVIDpWRI0fCZDLZsxbZdXQnHKwnIlKG5FB54YUX8OijjyImJgbu7u7W15cuXWqXwuTg7e0NgIcVExEpRXKovP322zhz5gza2tqsu79UKlW/DhV2KkREypIcKtnZ2Th79qw9a5FdR6fCUCEiUobkEfe5c+ciPz/fnrXIbtiwYQCAxsZGB1dCRDQ4SO5Ujh07hhkzZiA4OBju7u4D4pBihgoRkbIkh0p6enqfV5aeno5169bBYrHg0UcfxbPPPmvzfktLCx5++GGcOHECfn5+2Lt3L4KCgpCVlYWEhAQAgBACmzZtwpIlS7pdH0OFiEhZknd/LVq0CH/5y1+QmpqK/Px8uLq69mjOL4vFgjVr1uDAgQPIz8/Hnj17Ou1O27lzJ7RaLQoLC5GYmIgNGzYAAKZNm4acnBzk5uYiPT0djz/+ONrb27tdp0ajAQA0NDRIrpOIiHpPcqh8+eWXeOyxx+Dp6YnU1FRMmzYNn332meQVZWVlQa/XIzg4GEOGDEFsbCzS0tJslklLS8PKlSsBAMuWLUNGRgaEEPDy8oJafaOpam5uhkqlkrROdipERMqSHCq+vr749a9/jaeeegopKSnIzs6GXq+XvKKKigqMHTvW+jwwMBAVFRW3XEatVkOj0aC2thYAcPz4cYSFhSE8PBxvvvmmNWRuh6FCRKQsyaFy7tw5m+cTJ05UdJA+MjISeXl5yM7OxubNm7ucJDI5ORkGgwEGgwE1NTXc/UVEpDDJA/WPP/44ioqKoNPpMH36dDQ3N2PatGn46aefJE3WqNPpUFZWZn1eXl4OnU7X5TKBgYFob29HQ0MD/Pz8bJYJDQ2Fj48PTp06BYPBYPNeQkKCdUDfYDCwUyEiUpjkTuWLL77AhQsXsHfvXixatAh6vR5msxkRERGYMmVKt583Go0oKChAcXExWltbkZqa2mkuMZPJhJSUFADAvn37EB0dDZVKheLiYuvAfGlpKc6cOYOgoKBu18lQISJSluROpcO4ceMwbtw4LF682PratWvXul+RWo2kpCQsWLAAFosFq1atQlhYGDZu3AiDwQCTyYT4+HjExcVBr9fD19cXqampAIBvvvkGW7ZsgZubG1xcXLBjxw6MGDGi23Vy9xcRkbK6vfLjQGUwGLB//37odDqMGjUKVVVVji6JiKjfs/uVHweyjk6lvr7esYUQEQ0SknZ/nTlzBmlpadZDgHU6HUwmE0JDQ+1aXF95eXlhyJAhaG5uhtlstl5emIiI7KPbTmXr1q2IjY2FEAKzZ8/G7NmzIYTA7373O2zZskWJGntNpVJZjx7rON+FiIjsp9tOZefOncjLy4Obm5vN6+vXr0dYWFin+bv6G19fX1RVVaGurg6BgYGOLoeIyKl126m4uLigsrKy0+tVVVUD4lr17FSIiJTTbafy+uuvIyYmBhMnTrROoXLhwgUUFhYiKSnJ7gX2la+vLwCgrq7OwZUQETm/bkPl3nvvxblz55CVlWUzUG80GuHq6mr3Avuqo1NhqBAR2Z+ko79cXFwwZ84ce9diFx2dCnd/ERHZn1MfUgywUyEiUpJTH1IM/Bwqly9fdnAlRETOz+kPKfb39wcAXLp0ycGVEBE5P6c/pDggIAAAQ4WISAlOf0gxQ4WISDlOf0hxR6hUV1dDCCH5+vZERNRzfTqk+OLFixg1apTsRcnJ09MTQ4cOxdWrV1FfXw+tVuvokoiInFafBkXi4+PlqsOuuAuMiEgZfQqVTz/9VK467Kqjm7p48aKDKyEicm5Of/IjAOvsxGVlZQ6uhIjIuTn9yY8AMH78eAA3jlojIiL7cfqTHwFg3LhxAIDS0lIHV0JE5Nyc/uRH4OdQYadCRGRfTn/yI8BOhYhIKU5/8iNgO6bCEyCJiOzH6a+nAgAajQbDhg1DY2Mj6urqrDMXExGRvAbGoIgMuAuMiMj+Bk2o8LBiIiL7GzShwiPAiIjsr0ehcvjwYZv7gaQjVEpKShxbCBGRE+tRqDz99NM29wPJpEmTANyYcoaIiOyjV7u/hBBy12F3HfOUnT592sGVEBE5r0EzpqLX66FWq1FaWoqmpiZHl0NE5JQGTai4ublh4sSJEELg7Nmzji6HiMgpDZpQAYCpU6cC4C4wIiJ76VGo+Pj4AACGDh1ql2LsjeMqRET21aNQ+eqrr2zuB5qOTiU/P9/BlRAROadBtfuro1NhqBAR2YeioZKeno7JkydDr9d3edXIlpYWrFixAnq9HpGRkdYTFQ8dOoRZs2YhPDwcs2bN6vXJl6GhoXBzc8O5c+fQ2NjYl59CRERdUCxULBYL1qxZgwMHDiA/Px979uzp1DHs3LkTWq0WhYWFSExMxIYNGwAAI0aMwMcff4wff/wRKSkpiIuL61UN7u7uiIiIgBACJ06c6PNvIiIiW5JDpatdRkeOHJG8oqysLOj1egQHB2PIkCGIjY1FWlqazTJpaWlYuXIlAGDZsmXIyMiAEAIzZ87EmDFjAABhYWEwm81oaWmRvO6bzZ4921oPERHJS3KoPPDAA9i6dSuEEDCbzfjjH/+I5557TvKKKioqrFeOBIDAwEDrRb+6WkatVkOj0aC2ttZmmf/+97+444474O7u3mkdycnJMBgMMBgMqKmp6bKOjlDJzs6WXDsREUkjOVSOHz+OsrIyzJ07F0ajEWPGjEFmZqY9a+skLy8PGzZswL/+9a8u309ISEBOTg5ycnIwcuTILpcxGo0A2KkQEdmD5FBxc3ODp6cnzGYzmpubMWHCBLi4SB+S0el0KCsrsz4vLy+HTqe75TLt7e1oaGiwXqWxvLwcS5YswbvvvouQkBDJ6/2lyZMnY+jQoSgrK0NVVVWvv4eIiDqTnApGoxGenp7Izs7G119/jT179mD58uWSV2Q0GlFQUIDi4mK0trYiNTUVJpPJZhmTyYSUlBQAwL59+xAdHQ2VSoX6+nosXLgQW7ZswV133SV5nV25+dLIA/V8GyKi/kpyqOzcuRMvvvgi3NzcMHr0aKSlpXUKhdtRq9VISkrCggULEBoaigceeABhYWHYuHEj9u/fDwCIj49HbW0t9Ho9Xn31Vethx0lJSSgsLMSLL76IiIgIREREoLq6uoc/9WcxMTEAgM8//7zX30FERJ2pxECcx14Cg8GAnJycLt/LycmB0WjEhAkTcP78eYUrIyLqv273b6cUg+qM+g4zZ86EVqtFcXExQ4WISEaDMlRcXV0RFRUFgLvAiIjk1OtQqaqq6vUJiP3BggULAAAff/yxgyshInIevQ6VuLg4TJkyZUBerx64caSZSqXCwYMHcfXqVUeXQ0TkFHodKp9//jnOnz+PP/zhD3LWo5hRo0bhrrvuQmtrKz777DNHl0NE5BT6NKaiUqkQFhYmVy2KW7p0KQDggw8+cHAlRETOQXKoxMTEdPo/+oSEBNkLUlJHqOzfvx/19fWOLYaIyAlIDpXi4mJs3boVL7zwgvW1vhzL3B+MHz8eMTExaG5uRmpqqqPLISIa8CSHyvDhw5GRkYFLly5h8eLFaGhosGddilm1ahUA4K233nJwJUREA5/kUBFCQK1WY8eOHbj//vtx991392mqlP5iyZIl0Gg0yM7O5oW7iIj6SHKoPPHEE9bHjzzyCN555x3Mnz/fLkUpydPTE/Hx8QCAV1991cHVEBENbN3O/SWEgEqluu2XSFlGaT2Zv+bChQsIDg4GcGPs6OaLiRERDSZ2n/srKioK//znP3HhwgWb11tbW3H48GGsXLnSOl39QDVu3DgsX74cFosFL7/8sqPLISIasLrtVJqbm/HWW29h165dOH/+PLRaLcxmMywWCxYsWIDVq1dj5syZStUrWU/TNi8vD+Hh4VCr1SgoKMD48ePtWB0RUf/U106lR1Pft7W14fLly/Dy8oJGo+n1SpXQmw3z+9//Hrt27cLDDz884LsvIqLesHuoDB069LbjJY2Njb1euT31ZsMUFRVh6tSpaG1txVdffYV77rnHTtUREfVPdh9TuXr1KhobG295cyYhISF49tlnAQBPPvkk2traHFwREdHAMiivp3I7zz33HEJCQpCXl4dt27Y5uhwiogGFofILHh4e2LFjBwBg48aNyMrKcnBFREQDB0OlC/Pnz8dTTz2F9vZ2xMbGOs2UNERE9sZQuYWXX34Zd9xxB4qLixEXFweLxeLokoiI+j2Gyi24u7tj79690Gq1+Pjjj7Fu3Tr04OhrIqJBiaFyG3q9HmlpaRgyZAi2b9+OzZs3O7okIqJ+jaHSjXvuuQcpKSlQqVT461//ir/97W+OLomIqN9iqEgQGxuLt99+GyqVChs3bsQzzzyD69evO7osIqJ+h6Ei0cqVK/Hee+9BrVbjlVdewfLly9HU1OTosoiI+hWGSg88+OCDSE9Ph0ajwQcffACj0Ygff/zR0WUREfUbDJUeiomJwdGjRxEaGorTp0/DaDRi+/bt3B1GRASGSq+EhoYiOzsbjz32GFpaWrB27Vr83//9H06dOuXo0oiIHIqh0kve3t5ITk7G+++/j4CAAGRmZmLmzJlYv349Ll++7OjyiIgcgqHSR8uWLcOZM2ewevVqWCwWvPbaawgODsaLL76I+vp6R5dHRKQohooMhg8fju3bt+PEiRO49957cfXqVTz//PMYO3YsEhMTUVpa6ugSiYgUwVCR0cyZM3HgwAEcOXIEMTExuHbtGl5//XWEhIRg6dKl+OSTT9De3u7oMomI7IahYge/+tWv8Pnnn+O7777DQw89BAD48MMPsXjxYowdOxZPP/00jh07xiPGiMjp9Oga9QNJXy+JKaeqqir85z//wVtvvYWzZ89aX9fpdFiyZAnuu+8+3H333fDw8HBglUREClxOWE7p6emYPHky9Ho9tmzZ0un9lpYWrFixAnq9HpGRkSgpKQEA1NbWIioqCj4+Pli7dq2SJcti9OjReOaZZ3D69GlkZmZi3bp1GDt2LCoqKpCUlIR58+ZBq9ViwYIF2LZtG06ePMkuhogGJMU6FYvFgkmTJuHQoUMIDAyE0WjEnj17MHXqVOsyO3bswA8//IA333wTqamp+PDDD7F37140NTXh+++/x6lTp3Dq1CkkJSV1u77+1Kl0RQiBnJwcfPDBB0hPT0dubq7N+xqNBpGRkbjzzjsxd+5cREZGQqPROKZYIho0+vpvp1rGWm4rKysLer0ewcHBAG5M0piWlmYTKmlpadi0aROAG4fqrl27FkIIeHt74+6770ZhYaFS5dqdSqWC0WiE0WjE5s2bcenSJWRkZODQoUPIyMhAWVkZDh48iIMHD1qXnzhxImbMmGG9RUREQKfTQaVSOfjXEBHdoFioVFRUYOzYsdbngYGBOH78+C2XUavV0Gg0qK2txYgRIyStIzk5GcnJyQCAmpoamSpXRkBAAB588EE8+OCDAIDy8nIcPXrUejtx4gTOnTuHc+fO4f3337d+ztfXF6GhoZg8ebLNLSQkBG5ubo76OUQ0SCkWKkpISEhAQkICgBst3EAWGBiI5cuXY/ny5QBujDedPn0aubm5OHnypPVWV1eHzMxMZGZm2nze1dUVwcHB0Ov1CAoKQlBQECZMmGB9PGLECHY4RCQ7xUJFp9OhrKzM+ry8vBw6na7LZQIDA9He3o6Ghgb4+fkpVWK/5u7ujoiICERERFhfE0KgsrISZ86cwdmzZ21upaWlKCgoQEFBQZff5+XlZQ2YcePGQafTYcyYMdDpdNbHWq2WwUNEPaJYqBiNRhQUFKC4uBg6nQ6pqanYvXu3zTImkwkpKSm48847sW/fPkRHR/MftdtQqVTWEIiJibF5z2w2o7CwEMXFxSgpKel0u3LlCvLz85Gfn3/L7/f09MSYMWM6hc3o0aMREBBgvfn6+sLFhac8EZGCoaJWq5GUlIQFCxbAYrFg1apVCAsLw8aNG2EwGGAymRAfH4+4uDjo9Xr4+voiNTXV+vmgoCA0NjaitbUVH330EQ4ePGgzyE+2PD09ER4ejvDw8C7fr6+vR2lpKYqLi1FeXo7KykpUVFRYb5WVlWhsbERRURGKiopuuy5XV1f4+/vbBM2tbiNGjICrq6s9fjIR9QM8+ZFu6dq1a9aAuTlwLl26ZHO7cuWK5O90cXHBiBEjEBAQAH9/f4wcOfK29xqNht0qkYIGzCHFNPD4+PhYjya7ndbWVlRXV+PSpUu4ePFip9C5+VZbW4vq6mpUV1dLqsHNzQ0jR468bfjc/HjYsGEMISIHYqhQnw0ZMgSBgYEIDAzsdtm2tjbU1NTg0qVLqKmpQXV1tc39L1+7evUqKisrUVlZKbmW7rqfm+99fHwYQkQyYqiQotzc3KyD/1I0Nzd3GTa3um9qarLuppPCw8Oj207o5ntvb+++/Hwip8dQoX7Nw8MDY8eOtTlx9nbMZrPkAKqurobZbEZZWZnN4e634+npKWk3XMe9l5dXX34+0YDDUCGn4unpiXHjxmHcuHGSlm9qapK0G67j3mw2o7S0VPKF17y9vSV1QR03T0/Pvvx8IodjqNCg5u3tDW9vbwQFBXW7rBACTU1Nkrqgm3fHNTU1WWfcllLPzSHzyw7olzdvb2+OCVG/wlAhkkilUsHHxwc+Pj7WiVFvRwiBq1evStoN19Ed9TSEbh4TknLjIdpkbwwVIjtRqVQYNmwYhg0bhpCQkG6XF0KgsbHRGjBSbj0dE3Jzc8OIESO67YA6blqtlrMlUI8wVIj6CZVKBY1GA41GA71eL+kzHWNCt7rd3AXV1NTg2rVrqKqqQlVVlaTvd3V1hZ+fn+ROyM/PD2o1/1kZzPhfn2gA68mYEGB7iLaUW319fY9OVlWpVNBqtRg5ciRGjBgBPz8/6/0vH3c89/X1ZRA5Ef6XJBpEenqIdltbGy5fviy5E6qrq7Pezp49K7kujUZzy/C51WMPD4/ebgayI4YKEd2Sm5sbRo8ejdGjR0ta3mKxoLa2FjU1NaitrbXeLl++fMvHdXV1aGhoQENDQ7eTl97M29v7th1Qx2NfX1/4+vpCq9VCo9FwQlM7Y6gQkWw6Zqz29/eX/Jnr16+jvr7+tsHT1eOOI+UuXLggeV0d41YdISPlvuOxl5cXj5yTgKFCRA7l4uJi/cdbqo7DtaWE0JUrV1BXV4crV66goaEB9fX1qK+v73GdQ4YMkRxEN99rNBq4u7v3eH0DFUOFiAacmw/XnjBhguTPdVxRtiNkenLf3NxsnW27pzw8PDB8+HBoNBoMHz68x48H0kmuDBUiGjTUarV1rKWnzGYzrly50ilsbhdEHV1Rc3MzLl68iIsXL/aqbldX114FkkajwbBhw6DRaODm5tardfcUQ4WISAJPT0/rJbZ7QggBs9lsDZibd8FJffzTTz9ZA6y3PDw8rCHTETRd3fcVQ4WIyI5UKhW8vLzg5eXV40Dq0Nraaj1CrieB1NjYiMbGRjQ0NKC5udm6C+92Zs2a1asaOzBUiIj6uY6Lz40cObJXn+/olhoaGqwhc6v7L7/8sk+1MlSIiJzczd1Sd+ccGQyGPq2LM8UREZFsGCpERCQbhgoREcmGoUJERLJhqBARkWwYKkREJBuGChERyYahQkREsmGoEBGRbBgqREQkG4YKERHJhqFCRESyYagQEZFsGCpERCQbhgoREclG0VBJT0/H5MmTodfrsWXLlk7vt7S0YMWKFdDr9YiMjERJSYn1vc2bN0Ov12Py5Mn43//+p2DVREQklWKhYrFYsGbNGhw4cAD5+fnYs2cP8vPzbZbZuXMntFotCgsLkZiYiA0bNgAA8vPzkZqairy8PKSnp2P16tWwWCxKlU5ERBIpFipZWVnQ6/UIDg7GkCFDEBsbi7S0NJtl0tLSsHLlSgDAsmXLkJGRASEE0tLSEBsbC3d3d0yYMAF6vR5ZWVlKlU5ERBIpdjnhiooKjB071vo8MDAQx48fv+UyarUaGo0GtbW1qKiowJw5c2w+W1FR0WkdycnJSE5OBgCcOXOmz5fFdAY1NTW9vq61s+G2sMXt8TNui5+dOXOmT593qmvUJyQkICEhwdFl9CsGgwE5OTmOLqNf4Lawxe3xM26Lnw2Ya9TrdDqUlZVZn5eXl0On091ymfb2djQ0NMDPz0/SZ4mIyPEUCxWj0YiCggIUFxejtbUVqampMJlMNsuYTCakpKQAAPbt24fo6GioVCqYTCakpqaipaUFxcXFKCgowOzZs5UqnYiIJFJs95darUZSUhIWLFgAi8WCVatWISwsDBs3boTBYIDJZEJ8fDzi4uKg1+vh6+uL1NRUAEBYWBgeeOABTJ06FWq1Gtu3b4erq6tSpQ9o3B34M24LW9weP+O2+Flft4VKCCFkqoWIiAY5nlFPRESyYagQEZFsGCoD3KpVq+Dv749p06ZZX6urq8O8efMwceJEzJs3D1euXAEACCHw1FNPQa/XY/r06fjuu+8cVbZdlJWVISoqClOnTkVYWBjeeOMNAINzezQ3N2P27NmYMWMGwsLC8PzzzwMAiouLERkZCb1ejxUrVqC1tRXA7adIchYWiwUzZ87EokWLAAzebREUFITw8HBERERYDx+W82+EoTLAPfLII0hPT7d5bcuWLYiJiUFBQQFiYmKs86wdOHAABQUFKCgoQHJyMp588klHlGw3arUa//jHP5Cfn49jx45h+/btyM/PH5Tbw93dHYcPH8bJkyeRm5uL9PR0HDt2DBs2bEBiYiIKCwuh1Wqxc+dOALeeIsmZvPHGGwgNDbU+H8zb4osvvkBubq713BxZ/0YEDXjFxcUiLCzM+nzSpEmisrJSCCFEZWWlmDRpkhBCiISEBLF79+4ul3NGJpNJHDx4cNBvj6amJjFz5kxx7Ngx4efnJ9ra2oQQQnz77bdi/vz5Qggh5s+fL7799lshhBBtbW3Cz89PXL9+3WE1y62srExER0eLjIwMsXDhQnH9+vVBuy3Gjx8vampqbF6T82+EnYoTunTpEkaPHg0AGDVqFC5dugSg66lyupruxhmUlJTg+++/R2Rk5KDdHhaLBREREfD398e8efMQEhKC4cOHQ62+cSbBzb/3VlMkOYs//elPePnll+HicuOfvNra2kG7LVQqFebPn49Zs2ZZp7WS82/EqaZpoc5UKhVUKpWjy1DUtWvXcP/99+P111/HsGHDbN4bTNvD1dUVubm5qK+vx5IlS/o8p9NA9cknn8Df3x+zZs3CkSNHHF2Ow33zzTfQ6XSorq7GvHnzMGXKFJv3+/o3wk7FCQUEBKCqqgoAUFVVBX9/fwDSpsoZ6Nra2nD//ffjoYcewtKlSwEM7u0BAMOHD0dUVBSOHj2K+vp6tLe3A7D9vbeaIskZZGZmYv/+/QgKCkJsbCwOHz6MdevWDcptAcD6O/39/bFkyRJkZWXJ+jfCUHFCN093k5KSgvvuu8/6+rvvvgshBI4dOwaNRmNteZ2BEALx8fEIDQ3F+vXrra8Pxu1RU1OD+vp6AIDZbMahQ4cQGhqKqKgo7Nu3D0DnbdHVFEnOYPPmzSgvL0dJSQlSU1MRHR2NXbt2Dcpt0dTUhKtXr1ofHzx4ENOmTZP3b0TG8R9ygNjYWDFq1CihVquFTqcT//73v8Xly5dFdHS00Ov1IiYmRtTW1gohhLh+/bpYvXq1CA4OFtOmTRPZ2dkOrl5eX3/9tQAgwsPDxYwZM8SMGTPEp59+Oii3x8mTJ0VERIQIDw8XYWFh4oUXXhBCCFFUVCSMRqMICQkRy5YtE83NzUIIIcxms1i2bJkICQkRRqNRFBUVObJ8u/niiy/EwoULhRCDc1sUFRWJ6dOni+nTp4upU6eKv//970IIIevfCKdpISIi2XD3FxERyYahQkREsmGoEBGRbBgqREQkG4YKERHJhqFCZEfLli3D+fPnZfmu3/zmN9bZY4n6K4YKkZ3k5eXBYrEgODhYlu+Li4vDjh07ZPkuInthqBD1wptvvomIiAhERERgwoQJiIqK6rTMrl27rGcmA4CPj4/18b59+/DII48AuHH5gieffBJz5sxBcHAwjhw5glWrViE0NNS6DHDj7OY9e/bY7TcRyYGhQtQLTzzxBHJzc5GdnY3AwECbaWE6ZGZmYtasWZK+78qVKzh69Chee+01mEwmJCYmIi8vDz/++CNyc3MBAFqtFi0tLU41Yy45H4YKUR+sW7cO0dHRWLx4caf3qqqqMHLkSEnfs3jxYqhUKoSHhyMgIADh4eFwcXFBWFiYzZUH/f39UVlZKVf5RLLj1PdEvfTOO++gtLQUSUlJXb7v6emJ5uZm6/ObJyW8+XXgxpUaAcDFxcX6uON5x0y6HZ/z9PSUpX4ie2CnQtQLJ06cwLZt2/Dee+9ZL/z0S6GhoSgsLLQ+DwgIwOnTp3H9+nV8+OGHPV6nEAIXL15EUFBQb8smsjuGClEvJCUloa6uDlFRUYiIiMCjjz7aaZmFCxfaXBRqy5YtWLRoEebOndurKfZPnDiBOXPmWK9WSNQfcZZiIjsxm82IiopCZmYmXF1d+/x969atg8lkQkxMjAzVEdkHOxUiO/H09MQLL7wg23Xvp02bxkChfo+dChERyYadChERyYahQkREsmGoEBGRbBgqREQkG4YKERHJ5v8BD9etKGfBfEsAAAAASUVORK5CYII=", "text/plain": [ "
" ] }, "metadata": {}, "output_type": "display_data" } ], "source": [ "profile_error = u_diff.profile_longitudinal()\n", "plt.ylim(ymin=0)" ] } ], "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.10.12" }, "toc": { "base_numbering": 1, "nav_menu": { "height": "4px", "width": "254px" }, "number_sections": true, "sideBar": true, "skip_h1_title": false, "title_cell": "Table of Contents", "title_sidebar": "Contents", "toc_cell": false, "toc_position": {}, "toc_section_display": "block", "toc_window_display": false }, "toc-showtags": false, "varInspector": { "cols": { "lenName": 16, "lenType": 16, "lenVar": 40 }, "kernels_config": { "python": { "delete_cmd_postfix": "", "delete_cmd_prefix": "del ", "library": "var_list.py", "varRefreshCmd": "print(var_dic_list())" }, "r": { "delete_cmd_postfix": ") ", "delete_cmd_prefix": "rm(", "library": "var_list.r", "varRefreshCmd": "cat(var_dic_list()) " } }, "types_to_exclude": [ "module", "function", "builtin_function_or_method", "instance", "_Feature" ], "window_display": false } }, "nbformat": 4, "nbformat_minor": 4 }