#30: keep temporary nc-files to check results

nansencenter · Jul 11, 2019 · cb4fcd8 · cb4fcd8
1 parent c83278f
commit cb4fcd8
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 10 deletions.
diff --git a/geospaas_wind/managers.py b/geospaas_wind/managers.py
@@ -45,8 +45,9 @@ def process(self, uri, force=False, *args, **kwargs):
 
         # Direct reprojection fails - gdal can't read the bands if we do w.reproject...
         # Workaround: Export wind to temporary file
-        fd, tmp_filename = tempfile.mkstemp(suffix='.nc')
-        os.close(fd) # Just in case - see https://www.logilab.org/blogentry/17873
+        tmp_filename = os.path.join(settings.PRODUCTS_ROOT,'TMP_WIND_'+os.path.basename(uri))
+        #fd, tmp_filename = tempfile.mkstemp(suffix='.nc')
+        #os.close(fd) # Just in case - see https://www.logilab.org/blogentry/17873
         w.export(tmp_filename)
 
         # Read temporary file
@@ -85,7 +86,7 @@ def process(self, uri, force=False, *args, **kwargs):
         ww.export2thredds(thredds_fn, mask_name='swathmask', metadata=metadata)
         wds, cr = super(WindManager, self).get_or_create(wind_uri)
 
-        os.unlink(tmp_filename)
+        #os.unlink(tmp_filename)
 
         return wds, cr
 
diff --git a/notebooks/Sentinel-1 wind.ipynb b/notebooks/Sentinel-1 wind.ipynb
@@ -32,21 +32,21 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 4,
+   "execution_count": 22,
    "metadata": {},
    "outputs": [],
    "source": [
     "#fn = 'http://nbstds.met.no/thredds/dodsC/NBS/S1B/2019/06/30/IW/S1B_IW_GRDM_1SDV_20190630T051108_20190630T051141_016921_01FD7F_21CF.nc'\n",
     "#fn = '/home/vagrant/project/S1A_EW_GRDM_1SDH_20190704T070553_20190704T070654_027964_032853_C95D.nc'\n",
     "#fn = '/vagrant/lib/project/products/WIND_WIND_UNPROJECTED_S1A_IW_GRDM_1SDV_20190704T052723_20190704T052755_027963_03284A_AF7D.nc'\n",
     "fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165618_20190709T165650_028043_032AC4_6010.nc'\n",
-    "fn = '/home/vagrant/project/products/\n",
-    "fn = '/home/vagrant/project/products/"
+    "fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190709T165553_20190709T165625_028043_032AC4_4501.nc'\n",
+    "#fn = '/home/vagrant/project/products/WIND_S1A_IW_GRDM_1SDV_20190704T052723_20190704T052755_027963_03284A_AF7D.nc'"
    ]
   },
   {
    "cell_type": "code",
-   "execution_count": 5,
+   "execution_count": 23,
    "metadata": {},
    "outputs": [
     {
@@ -63,12 +63,12 @@
   },
   {
    "cell_type": "code",
-   "execution_count": 17,
+   "execution_count": 24,
    "metadata": {},
    "outputs": [
     {
      "data": {
-      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARYAAAD8CAYAAACy5YsfAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAANcklEQVR4nO3df6jdd33H8efLtLYyHbbrbQlJXDLIH0tlVrlkQsdw6tZMxfSfQgRH/ijknw4qG0gyYcM/Cm5/iH/1j6CygD9CQKWhyLYQFRkM0xtttWmMvVrXXhKa60TU/VHX+t4f9ys73tz2njTvk/u93fMBl+/3fO73nPO+pH3yPT84J1WFJHV63UYPIOm1x7BIamdYJLUzLJLaGRZJ7QyLpHYzC0uSfUkuJFlMcnhW9yNpfDKL97Ek2QL8APhzYAl4DPhQVT3VfmeSRmdWZyx7gcWq+lFV/Qo4Duyf0X1JGpkbZnS724DnJi4vAX/8cgffdttttXPnzhmNImlWzp49+5Oqmlu9PquwZI2133rMleQQcAjgLW95CwsLCzMaRdKsJPnPtdZn9VBoCdgxcXk7cHHygKo6WlXzVTU/N3dF8CRtYrMKy2PA7iS7krweOACcnNF9SRqZmTwUqqoXk/w18K/AFuCzVXVuFvclaXxm9RwLVfVV4Kuzun1J4+U7byW1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFkntDIukdoZFUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGpnWCS1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqt25Yknw2yeUkT06s3ZrkVJKnh+0tE787kmQxyYUk98xqcEnjNc0Zyz8D+1atHQZOV9Vu4PRwmSR7gAPAncN1Hk6ypW1aSZvCumGpqm8CP121vB84NuwfA+6dWD9eVS9U1TPAIrC3aVZJm8SrfY7ljqq6BDBsbx/WtwHPTRy3NKxJ+n+k+8nbrLFWax6YHEqykGRheXm5eQxJG+nVhuX5JFsBhu3lYX0J2DFx3Hbg4lo3UFVHq2q+qubn5uZe5RiSxujVhuUkcHDYPwg8MrF+IMlNSXYBu4Ez1zaipM3mhvUOSPJF4F3AbUmWgH8APgGcSHI/8CxwH0BVnUtyAngKeBF4oKpemtHskkZq3bBU1Yde5lfveZnjHwIeupahJG1uvvNWUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGpnWCS1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFkntDIukdoZFUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVK7dcOSZEeSryc5n+RckgeH9VuTnEry9LC9ZeI6R5IsJrmQ5J5Z/gGSxmeaM5YXgb+tqj8E3gk8kGQPcBg4XVW7gdPDZYbfHQDuBPYBDyfZMovhJY3TumGpqktV9e1h/xfAeWAbsB84Nhx2DLh32N8PHK+qF6rqGWAR2Ns9uKTxuqrnWJLsBN4OfAu4o6ouwUp8gNuHw7YBz01cbWlYW31bh5IsJFlYXl6++skljdbUYUnyRuBLwEeq6uevdOgaa3XFQtXRqpqvqvm5ublpx5C0CUwVliQ3shKVz1fVl4fl55NsHX6/Fbg8rC8BOyauvh242DOupM1gmleFAnwGOF9Vn5z41Ung4LB/EHhkYv1AkpuS7AJ2A2f6RpY0djdMcczdwF8B30vy+LD2d8AngBNJ7geeBe4DqKpzSU4AT7HyitIDVfVS++SSRmvdsFTVv7P28yYA73mZ6zwEPHQNc0naxHznraR2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFkntDIukdoZFUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGpnWCS1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFknt1g1LkpuTnEnyRJJzST4+rN+a5FSSp4ftLRPXOZJkMcmFJPfM8g+QND7TnLG8ALy7qt4G3AXsS/JO4DBwuqp2A6eHyyTZAxwA7gT2AQ8n2TKL4SWN07phqRW/HC7eOPwUsB84NqwfA+4d9vcDx6vqhap6BlgE9rZOLWnUpnqOJcmWJI8Dl4FTVfUt4I6qugQwbG8fDt8GPDdx9aVhbfVtHkqykGRheXn5Wv4GSSMzVViq6qWqugvYDuxN8tZXODxr3cQat3m0quaran5ubm66aSVtClf1qlBV/Qz4BivPnTyfZCvAsL08HLYE7Ji42nbg4jVPKmnTmOZVobkkbx723wC8F/g+cBI4OBx2EHhk2D8JHEhyU5JdwG7gTPfgksbrhimO2QocG17ZeR1woqoeTfIfwIkk9wPPAvcBVNW5JCeAp4AXgQeq6qXZjC9pjFJ1xdMf1938/HwtLCxs9BiSrlKSs1U1v3rdd95KamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGpnWCS1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFkntDIukdoZFUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1G7qsCTZkuQ7SR4dLt+a5FSSp4ftLRPHHkmymORCkntmMbik8bqaM5YHgfMTlw8Dp6tqN3B6uEySPcAB4E5gH/Bwki0940raDKYKS5LtwPuBT08s7weODfvHgHsn1o9X1QtV9QywCOztGVfSZjDtGcungI8Cv55Yu6OqLgEM29uH9W3AcxPHLQ1rvyXJoSQLSRaWl5evenBJ47VuWJJ8ALhcVWenvM2ssVZXLFQdrar5qpqfm5ub8qYlbQY3THHM3cAHk7wPuBn43SSfA55PsrWqLiXZClwejl8CdkxcfztwsXNoSeO27hlLVR2pqu1VtZOVJ2W/VlUfBk4CB4fDDgKPDPsngQNJbkqyC9gNnGmfXNJoTXPG8nI+AZxIcj/wLHAfQFWdS3ICeAp4EXigql665kklbRqpuuLpj+tufn6+FhYWNnoMSVcpydmqml+97jtvJbUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGpnWCS1MyyS2hkWSe0Mi6R2hkVSO8MiqZ1hkdTOsEhqZ1gktTMsktoZFkntDIukdoZFUjvDIqmdYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVI7wyKpnWGR1M6wSGo3VViS/DjJ95I8nmRhWLs1yakkTw/bWyaOP5JkMcmFJPfManhJ43Q1Zyx/VlV3VdX8cPkwcLqqdgOnh8sk2QMcAO4E9gEPJ9nSOLOkkbuWh0L7gWPD/jHg3on141X1QlU9AywCe6/hfiRtMtOGpYB/S3I2yaFh7Y6qugQwbG8f1rcBz01cd2lY+y1JDiVZSLKwvLz86qaXNEo3THnc3VV1McntwKkk33+FY7PGWl2xUHUUOAowPz9/xe8lbV5TnbFU1cVhexn4CisPbZ5PshVg2F4eDl8CdkxcfTtwsWtgSeO3bliS/E6SN/1mH/gL4EngJHBwOOwg8MiwfxI4kOSmJLuA3cCZ7sEljdc0D4XuAL6S5DfHf6Gq/iXJY8CJJPcDzwL3AVTVuSQngKeAF4EHquqlmUwvaZTWDUtV/Qh42xrr/wW852Wu8xDw0DVPJ2lT8p23ktoZFkntDIukdoZFUjvDIqmdYZHULlUb/276JMvAfwM/2ehZ1nEbztjBGXuMYcbfr6q51YujCAtAkoWJj2QYJWfs4Yw9xjyjD4UktTMsktqNKSxHN3qAKThjD2fsMdoZR/Mci6TXjjGdsUh6jdjwsCTZN3ya/2KSwxs4x2eTXE7y5MTaqL6JIMmOJF9Pcj7JuSQPjm3OJDcnOZPkiWHGj49txon73ZLkO0keHfGMm/MbMqpqw36ALcAPgT8AXg88AezZoFn+FHgH8OTE2j8Bh4f9w8A/Dvt7hllvAnYNf8OW6zDjVuAdw/6bgB8Ms4xmTlY+mvSNw/6NwLeAd45pxolZ/wb4AvDoGP+9h/v+MXDbqrXRzbn6Z6PPWPYCi1X1o6r6FXCclU/5v+6q6pvAT1ctj+qbCKrqUlV9e9j/BXCelQ8qH82cteKXw8Ubh58a04wASbYD7wc+PbE8qhlfwejn3OiwTPWJ/hvomr6JYJaS7ATezsoZwajmHB5iPM7K5yCfqqrRzQh8Cvgo8OuJtbHNCDP4hozrYdpP6Z+VqT7Rf4Q2dO4kbwS+BHykqn4+fGzomoeusTbzOWvlo0jvSvJmVj7W9K2vcPh1nzHJB4DLVXU2ybumucoaa9fr37v9GzKuh40+Yxn7J/qP7psIktzISlQ+X1VfHuucAFX1M+AbrHwj5phmvBv4YJIfs/Lw+91JPjeyGYHN+w0ZGx2Wx4DdSXYleT0rX816coNnmjSqbyLIyqnJZ4DzVfXJMc6ZZG44UyHJG4D3At8f04xVdaSqtlfVTlb+m/taVX14TDPCJv+GjI14xnjVM9zvY+XVjR8CH9vAOb4IXAL+h5Xy3w/8HivfS/30sL114viPDTNfAP7yOs34J6yc2n4XeHz4ed+Y5gT+CPjOMOOTwN8P66OZcdW87+L/XhUa1YysvFr6xPBz7jf/f4xtzrV+fOetpHYb/VBI0muQYZHUzrBIamdYJLUzLJLaGRZJ7QyLpHaGRVK7/wXKnfjw51YqggAAAABJRU5ErkJggg==\n",
+      "image/png": "iVBORw0KGgoAAAANSUhEUgAAARcAAAD8CAYAAABdJ+AhAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjAsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy+17YcXAAANfElEQVR4nO3df6jdd33H8efLtFaZDtv1toQkXTLIH0tlVrlkQsdw6tZMxfSfQgRH/ijknw6UDSSZsOEfBbc/xL/6R1BZwB8hoNJQZFuIigyG6Y1tZ9MYe7WuvSQ014mo+6Ou9b0/7ld2TG56T5L7zv2e+nzA5XzP53zPOe9L2me/50e/SVUhSevtdRs9gKTXJuMiqYVxkdTCuEhqYVwktTAuklq0xSXJniTnkiwmOdj1PJLGKR3fc0myCfg+8OfAEvA48KGqembdn0zSKHUduewGFqvqh1X1S+AosLfpuSSN0E1Nj7sFeGHi+hLwx1fa+fbbb6/t27c3jSKpy+nTp39cVXOr3dYVl6yy9huvv5IcAA4A3HXXXSwsLDSNIqlLkv+60m1dL4uWgG0T17cC5yd3qKrDVTVfVfNzc6uGT9IM64rL48DOJDuSvB7YBxxvei5JI9TysqiqXk7y18C/ApuAz1XVmY7nkjROXe+5UFVfA77W9fiSxs1v6EpqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIarFmXJJ8LsnFJE9PrN2W5ESSZ4fLWyduO5RkMcm5JPd1DS5p3KY5cvlnYM8laweBk1W1Ezg5XCfJLmAfcPdwn0eSbFq3aSXNjDXjUlXfAn5yyfJe4MiwfQS4f2L9aFW9VFXPAYvA7nWaVdIMudb3XO6sqgsAw+Udw/oW4IWJ/ZaGtcskOZBkIcnC8vLyNY4haazW+w3drLJWq+1YVYerar6q5ufm5tZ5DEkb7Vrj8mKSzQDD5cVhfQnYNrHfVuD8tY8naVZda1yOA/uH7f3AoxPr+5LckmQHsBM4dX0jSppFN621Q5IvAe8Cbk+yBPwD8EngWJIHgeeBBwCq6kySY8AzwMvAQ1X1StPskkZszbhU1YeucNN7rrD/w8DD1zOUpNnnN3QltTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIarFmXJJsS/KNJGeTnEnykWH9tiQnkjw7XN46cZ9DSRaTnEtyX+cvIGmcpjlyeRn426r6Q+CdwENJdgEHgZNVtRM4OVxnuG0fcDewB3gkyaaO4SWN15pxqaoLVfWdYfvnwFlgC7AXODLsdgS4f9jeCxytqpeq6jlgEdi93oNLGreres8lyXbg7cC3gTur6gKsBAi4Y9htC/DCxN2WhrVLH+tAkoUkC8vLy1c/uaRRmzouSd4EfBn4aFX97NV2XWWtLluoOlxV81U1Pzc3N+0YkmbEVHFJcjMrYflCVX1lWH4xyebh9s3AxWF9Cdg2cfetwPn1GVfSrJjm06IAnwXOVtWnJm46DuwftvcDj06s70tyS5IdwE7g1PqNLGkW3DTFPvcCfwV8N8mTw9rfAZ8EjiV5EHgeeACgqs4kOQY8w8onTQ9V1SvrPrmkUVszLlX176z+PgrAe65wn4eBh69jLkkzzm/oSmphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdRizbgkeUOSU0meSnImySeG9duSnEjy7HB568R9DiVZTHIuyX2dv4CkcZrmyOUl4N1V9TbgHmBPkncCB4GTVbUTODlcJ8kuYB9wN7AHeCTJpo7hJY3XmnGpFb8Yrt48/BSwFzgyrB8B7h+29wJHq+qlqnoOWAR2r+vUkkZvqvdckmxK8iRwEThRVd8G7qyqCwDD5R3D7luAFybuvjSsSfotMlVcquqVqroH2ArsTvLWV9k9qz3EZTslB5IsJFlYXl6eblpJM+OqPi2qqp8C32TlvZQXk2wGGC4vDrstAdsm7rYVOL/KYx2uqvmqmp+bm7uG0SWN2TSfFs0lecuw/UbgvcD3gOPA/mG3/cCjw/ZxYF+SW5LsAHYCp9Z7cEnjdtMU+2wGjgyf+LwOOFZVjyX5D+BYkgeB54EHAKrqTJJjwDPAy8BDVfVKz/iSxipVl70dcsPNz8/XwsLCRo8h6SolOV1V86vd5jd0JbUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1mDouSTYleSLJY8P125KcSPLscHnrxL6HkiwmOZfkvo7BJY3b1Ry5fAQ4O3H9IHCyqnYCJ4frJNkF7APuBvYAjyTZtD7jSpoVU8UlyVbg/cBnJpb3AkeG7SPA/RPrR6vqpap6DlgEdq/PuJJmxbRHLp8GPgb8amLtzqq6ADBc3jGsbwFemNhvaVj7DUkOJFlIsrC8vHzVg0satzXjkuQDwMWqOj3lY2aVtbpsoepwVc1X1fzc3NyUDy1pVtw0xT73Ah9M8j7gDcDvJvk88GKSzVV1Iclm4OKw/xKwbeL+W4Hz6zm0pPFb88ilqg5V1daq2s7KG7Vfr6oPA8eB/cNu+4FHh+3jwL4ktyTZAewETq375JJGbZojlyv5JHAsyYPA88ADAFV1Jskx4BngZeChqnrluieVNFNSddnbITfc/Px8LSwsbPQYkq5SktNVNb/abX5DV1IL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWxkVSC+MiqYVxkdTCuEhqYVwktTAukloYF0ktjIukFsZFUgvjIqmFcZHUwrhIamFcJLUwLpJaGBdJLYyLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWhgXSS2Mi6QWU8UlyY+SfDfJk0kWhrXbkpxI8uxweevE/oeSLCY5l+S+ruEljdfVHLn8WVXdU1Xzw/WDwMmq2gmcHK6TZBewD7gb2AM8kmTTOs4saQZcz8uivcCRYfsIcP/E+tGqeqmqngMWgd3X8TySZtC0cSng35KcTnJgWLuzqi4ADJd3DOtbgBcm7rs0rEn6LXLTlPvdW1Xnk9wBnEjyvVfZN6us1WU7rUTqAMBdd9015RiSZsVURy5VdX64vAh8lZWXOS8m2QwwXF4cdl8Ctk3cfStwfpXHPFxV81U1Pzc3d+2/gaRRWjMuSX4nyZt/vQ38BfA0cBzYP+y2H3h02D4O7EtyS5IdwE7g1HoPLmncpnlZdCfw1SS/3v+LVfUvSR4HjiV5EHgeeACgqs4kOQY8A7wMPFRVr7RML2m01oxLVf0QeNsq6/8NvOcK93kYePi6p5M0s/yGrqQWxkVSC+MiqYVxkdTCuEhqkarLvjx744dIloH/AX680bOs4Xaccb3MwpyzMCNs7Jy/X1Wrfgt2FHEBSLIw8X9cj5Izrp9ZmHMWZoTxzunLIkktjIukFmOKy+GNHmAKzrh+ZmHOWZgRRjrnaN5zkfTaMqYjF0mvIRselyR7hhN5LyY5uMGzfC7JxSRPT6yN6kTkSbYl+UaSs0nOJPnI2OZM8oYkp5I8Ncz4ibHNOPG8m5I8keSxEc84myfIr6oN+wE2AT8A/gB4PfAUsGsD5/lT4B3A0xNr/wQcHLYPAv84bO8a5r0F2DH8HptuwIybgXcM228Gvj/MMpo5WTkb4ZuG7ZuBbwPvHNOME7P+DfBF4LEx/nkPz/0j4PZL1kY356U/G33kshtYrKofVtUvgaOsnOB7Q1TVt4CfXLI8qhORV9WFqvrOsP1z4Cwr5ygezZy14hfD1ZuHnxrTjABJtgLvBz4zsTyqGV/F6Ofc6LjMwsm8R3si8iTbgbezcmQwqjmHlxtPsnL60xNVNboZgU8DHwN+NbE2thlhRk+QP+0JurtMdTLvkdrQ2ZO8Cfgy8NGq+tlwpsBVd11lrX3OWjn74D1J3sLKmQzf+iq73/AZk3wAuFhVp5O8a5q7rLJ2o/681/0E+TfCRh+5THUy7w12XSci75DkZlbC8oWq+spY5wSoqp8C32TlL8gb04z3Ah9M8iNWXo6/O8nnRzYj0HOC/Btho+PyOLAzyY4kr2flb2o8vsEzXWpUJyLPyiHKZ4GzVfWpMc6ZZG44YiHJG4H3At8b04xVdaiqtlbVdlb+uft6VX14TDPCjJ8gfyPeRb7kXe/3sfKJxw+Aj2/wLF8CLgD/y8p/AR4Efo+Vv6722eHyton9Pz7MfQ74yxs045+wcpj7n8CTw8/7xjQn8EfAE8OMTwN/P6yPZsZL5n0X//9p0ahmZOWT1KeGnzO//ndkbHOu9uM3dCW12OiXRZJeo4yLpBbGRVIL4yKphXGR1MK4SGphXCS1MC6SWvwfldf/t/Dmb/QAAAAASUVORK5CYII=\n",
       "text/plain": [
        "<Figure size 432x288 with 1 Axes>"
       ]
@@ -80,7 +80,7 @@
     }
    ],
    "source": [
-    "plt.imshow(n['V'], vmin=-10, vmax=10, cmap='jet')\n",
+    "plt.imshow(n['U'], vmin=-10, vmax=10, cmap='jet')\n",
     "plt.show()"
    ]
   },