rsignell-usgs/sandy_example.ipynb

## sandy_example.ipynb
{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# Use HSDS (data chunks in S3 objects) to access Model Output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:58:58.375407Z",
     "start_time": "2017-07-31T16:58:57.861042-04:00"
    }
   },
   "outputs": [],
   "source": [
    "%matplotlib inline\n",
    "import h5pyd\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "import matplotlib.image as mpimg\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "{'about': 'HSDS is a webservice for HSDS data',\n",
       " 'endpoint': 'http://149.165.157.109:5101',\n",
       " 'greeting': 'Welcome to HSDS!',\n",
       " 'hsds_version': '0.1',\n",
       " 'name': 'HSDS on JetStream',\n",
       " 'password': '***************',\n",
       " 'start_time': 1512254660,\n",
       " 'state': 'READY',\n",
       " 'username': 'rsignell'}"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "h5pyd.getServerInfo()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:58:58.755660Z",
     "start_time": "2017-07-31T16:58:58.393456-04:00"
    }
   },
   "outputs": [],
   "source": [
    "# Open the wind data \"file\"\n",
    "# server endpoint, username, password is found via a config file\n",
    "#f = h5pyd.File(\"/home/rsignell/Sandy_ocean_his_nc4.nc\", 'r')\n",
    "#f = h5pyd.File(\"/home/rsignell/sandy.nc\", 'r')\n",
    "f = h5pyd.File(\"/home/rsignell/sandy2.nc\", 'r')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:00.029849Z",
     "start_time": "2017-07-31T16:58:58.762683-04:00"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "('CPP_options', b'SANDY, ANA_BSFLUX, ANA_BTFLUX, ANA_FSOBC, ANA_M2OBC, ANA_SRFLUX, ANA_SSFLUX, ASSUMED_SHAPE, ATM_PRESS, ATM2OCN_FLUXE AVERAGES, COARE_TAYLOR_YELLAND, CURVGRID, DJ_GRADPS, DOUBLE_PRECISION, GLS_MIXING, KANTHA_CLAYSON, MASKING, MCT_LIB, MIX_GEO_TS, MIX_S_UV, MPI, NESTING, NONLINEAR, NONLIN_EOS, N2S2_HORAVG, !ONE_WAY, POWER_LAW, PROFILE, K_GSCHEME, RAMP_TIDES, RI_SPLINES, !RST_SINGLE, SALINITY, SOLAR_SOURCE, SOLVE3D, SPLINES_VDIFF, SPLINES_VVISC, SSH_TIDES, SWAN_COUPLING, TS_U3HADVECTION, TS_C4VADVECTION, TS_DIF2, UV_ADV, UV_COR, UV_U3HADVECTION, UV_C4VADVECTION, UV_LOGDRAG, UV_TIDES, UV_VIS2, VAR_RHO_2D, WAVES_OCEAN, WRF_COUPLING')\n",
      "('Conventions', b'CF-1.4, SGRID-0.3')\n",
      "('NLM_LBC', b'\\nEDGE:  WEST   SOUTH  EAST   NORTH  \\nzeta:  Clo    Cha    Cha    Cha    \\nubar:  Clo    Fla    Fla    Fla    \\nvbar:  Clo    Fla    Fla    Fla    \\nu:     Clo    RadNud RadNud RadNud \\nv:     Clo    RadNud RadNud RadNud \\ntemp:  Gra    RadNud RadNud Gra    \\nsalt:  Gra    RadNud RadNud Gra    \\ntke:   Clo    Gra    Gra    Gra')\n",
      "('_NCProperties', b'version=1|netcdflibversion=4.5.0|hdf5libversion=1.10.1')\n",
      "('_nc3_strict', 1)\n",
      "('ana_file', b'ROMS/Functionals/ana_btflux.h, ROMS/Functionals/ana_fsobc.h, ROMS/Functionals/ana_m2obc.h, ROMS/Functionals/ana_stflux.h')\n",
      "('avg_file', b'Sandy_ocean_avg.nc')\n",
      "('bry_file_01', b'Projects/Sandy/Sandy_bdy.nc')\n",
      "('code_dir', b'/peach/data2/tkalra/COAWST_regress_peach/COAWST_src2')\n",
      "('compiler_command', b'/share/apps/openmpi/bin/mpif90')\n",
      "('compiler_flags', b' -fastsse -Mipa=fast -tp k8-64 -Mfree')\n",
      "('compiler_system', b'pgi')\n",
      "('cpu', b'x86_64')\n",
      "('file', b'Sandy_ocean_his.nc')\n",
      "('format', b'netCDF-3 64bit offset file')\n",
      "('frc_file_01', b'Projects/Sandy/roms_namnarr_Sandy2012.nc')\n",
      "('frc_file_02', b'Projects/Sandy/tide_forc_Sandy.nc')\n",
      "('grd_file', b'Projects/Sandy/Sandy_roms_grid.nc')\n",
      "('header_dir', b'/peach/data2/tkalra/COAWST_regress_peach/COAWST_src2/Projects/Sandy')\n",
      "('header_file', b'sandy.h')\n",
      "('his_file', b'Sandy_ocean_his.nc')\n",
      "('history', b'ROMS/TOMS, Version 3.7, Thursday - February 23, 2017 -  7:31:24 AM')\n",
      "('ini_file', b'Projects/Sandy/Sandy_ini.nc')\n",
      "('os', b'Linux')\n",
      "('rst_file', b'Sandy_ocean_rst.nc')\n",
      "('svn_rev', b'1138M')\n",
      "('svn_url', b'https:://myroms.org/svn/src')\n",
      "('tiling', b'004x002')\n",
      "('title', b'Hurricane Sandy')\n",
      "('type', b'ROMS/TOMS history file')\n",
      "('var_info', b'ROMS/External/varinfo.dat')\n"
     ]
    }
   ],
   "source": [
    "[print(k) for k in f.attrs.items()];"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.167074Z",
     "start_time": "2017-07-31T16:59:00.033860-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "['AKv',\n",
       " 'Akk_bak',\n",
       " 'Akp_bak',\n",
       " 'Akt_bak',\n",
       " 'Akv_bak',\n",
       " 'Charnok_alpha',\n",
       " 'CrgBan_cw',\n",
       " 'Cs_r',\n",
       " 'Cs_w',\n",
       " 'Dwave',\n",
       " 'EminusP',\n",
       " 'FSobc_in',\n",
       " 'FSobc_out',\n",
       " 'Falpha',\n",
       " 'Fbeta',\n",
       " 'Fgamma',\n",
       " 'Hwave',\n",
       " 'Lm2CLM',\n",
       " 'Lm3CLM',\n",
       " 'LnudgeM2CLM',\n",
       " 'LnudgeM3CLM',\n",
       " 'LnudgeTCLM',\n",
       " 'LsshCLM',\n",
       " 'LtracerCLM',\n",
       " 'LtracerSponge',\n",
       " 'LtracerSrc',\n",
       " 'LuvSponge',\n",
       " 'LuvSrc',\n",
       " 'LwSrc',\n",
       " 'Lwave',\n",
       " 'Lwavep',\n",
       " 'M2nudg',\n",
       " 'M2obc_in',\n",
       " 'M2obc_out',\n",
       " 'M3nudg',\n",
       " 'M3obc_in',\n",
       " 'M3obc_out',\n",
       " 'N',\n",
       " 'Pair',\n",
       " 'Pwave_bot',\n",
       " 'Pwave_top',\n",
       " 'Tcline',\n",
       " 'Tnudg',\n",
       " 'Tnudg_SSS',\n",
       " 'Tobc_in',\n",
       " 'Tobc_out',\n",
       " 'Uwave_rms',\n",
       " 'Uwind',\n",
       " 'Uwind_eastward',\n",
       " 'Vstretching',\n",
       " 'Vtransform',\n",
       " 'Vwind',\n",
       " 'Vwind_northward',\n",
       " 'Znudg',\n",
       " 'Zob',\n",
       " 'Zos',\n",
       " 'Zos_hsig_alpha',\n",
       " 'angle',\n",
       " 'boundary',\n",
       " 'bustr',\n",
       " 'bvstr',\n",
       " 'dstart',\n",
       " 'dt',\n",
       " 'dtfast',\n",
       " 'el',\n",
       " 'eta_psi',\n",
       " 'eta_rho',\n",
       " 'eta_u',\n",
       " 'eta_v',\n",
       " 'evaporation',\n",
       " 'f',\n",
       " 'gamma2',\n",
       " 'gls',\n",
       " 'gls_Kmin',\n",
       " 'gls_Pmin',\n",
       " 'gls_c1',\n",
       " 'gls_c2',\n",
       " 'gls_c3m',\n",
       " 'gls_c3p',\n",
       " 'gls_cmu0',\n",
       " 'gls_m',\n",
       " 'gls_n',\n",
       " 'gls_p',\n",
       " 'gls_sigk',\n",
       " 'gls_sigp',\n",
       " 'grid',\n",
       " 'h',\n",
       " 'hc',\n",
       " 'lat_psi',\n",
       " 'lat_rho',\n",
       " 'lat_u',\n",
       " 'lat_v',\n",
       " 'latent',\n",
       " 'lon_psi',\n",
       " 'lon_rho',\n",
       " 'lon_u',\n",
       " 'lon_v',\n",
       " 'lwrad',\n",
       " 'mask_psi',\n",
       " 'mask_rho',\n",
       " 'mask_u',\n",
       " 'mask_v',\n",
       " 'nAVG',\n",
       " 'nHIS',\n",
       " 'nRST',\n",
       " 'ndefAVG',\n",
       " 'ndefHIS',\n",
       " 'ndtfast',\n",
       " 'nl_tnu2',\n",
       " 'nl_visc2',\n",
       " 'ntimes',\n",
       " 'ntsAVG',\n",
       " 'ocean_time',\n",
       " 'omega',\n",
       " 'pm',\n",
       " 'pn',\n",
       " 'rain',\n",
       " 'rdrg',\n",
       " 'rdrg2',\n",
       " 'rho0',\n",
       " 's_rho',\n",
       " 's_w',\n",
       " 'salt',\n",
       " 'sensible',\n",
       " 'shflux',\n",
       " 'spherical',\n",
       " 'ssflux',\n",
       " 'sustr',\n",
       " 'svstr',\n",
       " 'swrad',\n",
       " 'sz_alpha',\n",
       " 'temp',\n",
       " 'theta_b',\n",
       " 'theta_s',\n",
       " 'tke',\n",
       " 'tracer',\n",
       " 'u',\n",
       " 'ubar',\n",
       " 'v',\n",
       " 'vbar',\n",
       " 'w',\n",
       " 'xi_psi',\n",
       " 'xi_rho',\n",
       " 'xi_u',\n",
       " 'xi_v',\n",
       " 'xl',\n",
       " 'zeta']"
      ]
     },
     "execution_count": 5,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "list(f)  # list the datasets in the file"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.481363Z",
     "start_time": "2017-07-31T16:59:01.176093-04:00"
    }
   },
   "outputs": [],
   "source": [
    "dset = f['Hwave']"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.490385Z",
     "start_time": "2017-07-31T16:59:01.484368-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'d-66a40690-d860-11e7-a439-0242ac11000e'"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dset.id.id  # if this is an int, then you are using h5py!"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.505424Z",
     "start_time": "2017-07-31T16:59:01.493393-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(97, 64, 84)"
      ]
     },
     "execution_count": 8,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dset.shape"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 9,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.913809Z",
     "start_time": "2017-07-31T16:59:01.508433-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[<matplotlib.lines.Line2D at 0x7fe59403b5c0>]"
      ]
     },
     "execution_count": 9,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAXEAAAD8CAYAAACB3pQWAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJzt3Xd8VFXex/HPmfQE0kghpBJqEnoC\nBESUXixYUEFBfSy49rpre1x1V91d197XxQoWFFGKqCiKICVAkJICBEiDhJAQSO9znj8S92GRkkAm\nd+7M7/165RUymcx8rzf5eufMuecqrTVCCCHMyWJ0ACGEEGdOSlwIIUxMSlwIIUxMSlwIIUxMSlwI\nIUxMSlwIIUxMSlwIIUxMSlwIIUxMSlwIIUzM1RYPGhQUpGNiYmzx0EII4ZBSU1NLtNbBbf05m5R4\nTEwMmzdvtsVDCyGEQ1JK5Z7Jz8lwihBCmJiUuBBCmJiUuBBCmJiUuBBCmJiUuBBCmJiUuBBCmJiU\nuBBCmFir5okrpe4GbgYU8G+t9Us2TSWEnSmvbWDJ1gIOldeilMKiFL1COzEhPhQ3FzkWEsY5bYkr\npfrRXODDgHrgW6XU11rrLFuHE8Jo+4or+WBdDgtT91NV3/S774f5eTIrOZqrh0UR4ONuQELh7Fpz\nJB4HbNBaVwMopX4GLgWetWUwIYzUZNW8uWoPL/6QhYtSXDgwjOtHxjAgwh+tNU1Wzc+7i3l3bTb/\n/G4X//p5L3+7bAAXDAgzOrpwMq0p8TTgaaVUF6AGmArIOfXCYR0sq+WeBb+yYV8pFw/sxmMXxhPc\n2eM/31dK4eqiGBcXyri4UDILy3lo0Q5u/3gLP++O4PGLEvDxsMmKFkL8jtJan/5OSt0I3A5UAhlA\njdb63uPuMweYAxAVFZWYm3tGywAIYahvdhTyyJc7qGu08pdp/bh8SDhKqdP+XEOTlZd+2M0bq/bS\nvYsPH9+cTFc/zw5ILByFUipVa53U5p9rTYkf90TPAPu11m+c7D5JSUlaFsASZlJW08CTS9JZ9OsB\n+of78fKMQcQGd2rz46zbW8KcD1MJ8/Pks1tGyDi5aLUzLfFWva2ulApp+RwFXAZ80tYnEsJercws\nYvJLq1m8rYC7x/Vi0W0jz6jAAUb2COLf1yaRW1rN/7y/iaq6xnZOK8R/a+3cqC+UUhnAUuB2rfUR\nG2YSokPklFRxw/ubuPGDzfh4uPLFrSO5d0Lvs54yOKJHF16dOZjt+49yy7xU6hut7ZRYiN9r1bsv\nWutzbR1EiI5SXd/I6z/t4d+rs3F3tfDo1DiuGxmDu2v7zfeelNCVf1w+gD8u3M7fv9nJny+Kb7fH\nFuJY8ha6cBpaa5ZtL+SZ5ZkUltVy2eBwHprSlxBf27wBeUVSJOkF5by7NpuRPbowPj7UJs8jnJuc\naiacQmlVPde9t4k7P/mVQB93Fv5hBC9cNchmBf6bh6f2JaGbLw8s3EbB0RqbPpdwTlLiwuFtzT/K\nha+sYcO+wzx5cQJL7hhFUkxghzy3h6sLr109hIZGK3d/+iuNTTI+LtqXlLhwaB+n5HHlW+uxWBRf\n/GEk142MwcVy+nnf7al7kA9PX9qfTTlHeOXHPR363MLxyZi4cEhaa178IYtXVmZxXu9gXp4xCH9v\n4+ZsXzI4nDVZJbz2Yxbn9Q4iMbpjXgkIxydH4sLhWK2aJ5dm8MrKLK5IjOCd65IMLfDfPHFxPBEB\n3tz96VYqahuMjiMchJS4cCiNTVb+uHA776/L4cZR3fnH5QNwtZOlYjt7uvHiVYMoLKvlz4vTjY4j\nHIR9/HYL0Q5q6pu4ZV4qX2zZz30TevO/F8Rh6eDx79NJjA7grrG9+PLXAyzZVmB0HOEApMSFQzha\nXc81czfw465DPHVJP+4a16tVC1cZ4fYxPRgU6c9flqZTViPDKuLsSIkL0yuprGP6W+tJO1DOG1cP\nYVZytNGRTsnVxcJTl/SjtKqeF7/fbXQcYXJS4sLUtNY8+uUO8kqr+eCGYUzpb46LMvQL9+Oa4dF8\nuD6HjIJyo+MIE5MSF6a2ZFsB36UXcd+E3ozo0cXoOG3ywMQ++Hu78+fFaVitbVsSWojfSIkL0zpU\nUcvjS9IZFOnPzefGGh2nzfy83Xhocl825x5h0a8HjI4jTEpKXJhS8zBKGtX1TTx3xcAOPwuzvUxP\njGBwlD//+HanrD0uzoiUuDCl79KL+D6jiPsn9KZnyJldwMEeWCyK/70gnuKKOv69Zp/RcYQJSYkL\nU1qYmk83P09uMuEwyvESowOY2r8rb6/ex6GKWqPjCJOREhemU1nXyOqsEib162raYZTj/WlSXxqa\nrLz4fZbRUYTJSIkL0/lp5yHqG61MTuhqdJR2ExPkwzXDo1mwKY+sogqj4wgTkRIXpvNt+kGCOrl3\n2JrgHeWucb3wcXflb9/sNDqKMBEpcWEqtQ1N/LTzEBPiHWco5TeBPu7cMbYnP+48xJqsYqPjCJOQ\nEhemsiarhOr6Jib3c5yhlGNdf04M0V28+euyDLkKkGgVKXFhKt+mHcTX05URseY6O7O1PFxdeGRq\nHLuLKvl4Y57RcYQJSIkL02hosvJDZhHj40Jxd3XcX92J8aGM7NGFF77fzdHqeqPjCDvnuH8JwuFs\n2HeYspoGhx1K+Y1SiscujKe8poGXfpAph+LUpMSdSEOTlY3Zpew8WE51vflO8f4m7SBebi6M7h1s\ndBSbiwvzZcawKOZtyGXXQZlyKE5OLpTsBA4creHTjXl8uimf4oq6/9we1MmD5NhApvYPY0yfELzc\nXQxMeWr1jVaW7yhkQnwonm72m7M9PTCxD8t3FPLY4jQWzEm224tcCGNJiTswrTX/Wr2PZ7/diQbG\n9glhemIEjVZNXmk1e4sr+XlXMcu2F+Ll5sKkhFBmJUeTGB1gd4Xxy55ijlY3MG1QN6OjdJhAH3ce\nnNyXhxft4MtfD3DZkAijIwk7JCXuoJqsmseXpDF/Qx4X9A/j4al9iQjw/t39GluGWJbtKGTp1gK+\n2lpA366duX5kDNMTI+zmIsOLtxbg7+3Gub0cfyjlWFclRbJgUz7PLM9kXFwofl5uRkcSdsY+/kJF\nu/rtgsHzN+Rxy3mxvDpz8AkLHJovFTayZxDPXNqfDY+M4++X9cfFonho0Q4mvria5TsK0drYCxZU\n1zeyIr2Iqf3DHHpWyolYLOo/l3J7YcUuo+MIO+RcfxFOQGvNAwu38ePOIv4yLYGHp7T+iu8+Hq7M\nGBbFsjtH8e9rk3B1Udz20RYueX0t2/KP2jj5yX2fUURNQxPTBjrPUMqx+oX7MSs5mnkbckk7UGZ0\nHGFnpMQdzJJtBXy9vZD7J/bh2hExZ/QYSikmxIfyzd2j+ef0ARSU1XLJG2t59MsdlFV3/NXZF28t\nIMzPk6EOtlZKW9zfcim3vy7LMPyVkbAvUuIOpLCshse+SmNIlD+3jD77dbZdLIorkiJZef95XD8y\nhk825jHm+VW8sWoPFbUdU+alVfWs3l3MxQO7tfoVhSPy83Ljvgm9Scku5bv0g0bHEXZE3th0EFar\n5k8Lt9PQpHnhykHt+oakr6cbj1+UwBWJkfz92508++0u3ly1l9nJ0QyK9Ce4swddfDw4cLSG9IIy\nMgsrCA/wYnZyNMGdPc7quZfvKKTRqpk2KLydtsa8ZgyNZN76XJ5ensn5fUKcZqqlODUpcQfxUUou\na7JKeOqSfsQE+djkOeK7+fLhDcPYsb+MN1bt4c2f93KiV/ZBndw5XFXPWz/v5fIh4dwyuscZZapr\nbGL+hlx6hnQiLqxzO2yBubm6WHjswnhmvZPCe2tzuPX8HkZHEnZAStwBNDRZee2nPQzvHsg1w6Ns\n/nz9I/x4c1YiJZV1FB6tpbiylpKKekJ8PUjo5kdwZw/2FVcy95dsFqbuZ8nWAv41O4lRvYLa9DxP\nLs1g58EK3pqVaHfz1o0yqlcQ4+NCef2nPVyeGE5IZ0+jIwmDyZi4A1iZWURReR03jureoWUX1MmD\n/hF+jO0bypVDIzm/T8h/hk9igzvxzKX9+fmP5xMR4M0N729i+Y7CVj/2Z5vz+Tgljz+c18Ph10pp\nq0cviKOusYkXVuw2OoqwA60qcaXUvUqpdKVUmlLqE6WU/O/fjszfkEeYnydj+4YYHeV3wvy8+OyW\nEfSP8OP2j7cwb0PuaWdXpB0o43+/SmNkjy48MLF3ByU1j+5BPlw7IobPNueTWVhudBxhsNOWuFIq\nHLgLSNJa9wNcgBm2DiZaJ7ukil/2lDBzWJTdnF15PD9vN+bfOJwxfUJ47Ks05sxLpaj8xFd1X7un\nhBs/2ESQjzuvzhxst9tktDvH9qSzpxvPLM+UKYdOrrV/Ia6Al1LKFfAGCmwXSbTFxym5uFoUM4ZG\nGh3llLzcXXh7diKPTo1j9e5ixr/wM/PW57CvuJL6RisVtQ088uUOrpmbgo+7K+9cP5Qunc5uZosj\n8/d2565xvViTVcKq3XIpN2d22jc2tdYHlFLPAXlADbBCa73C5snEadU2NPF56n4mJoQS4mv/I1yu\nLhZuHh3LhPhQHvxiO48tTgfAosDTzYXahibmjI7lvgm9ZfpcK8xOjmbe+hye/jqTc3sGyasWJ3Xa\nEldKBQDTgO7AUeBzpdQsrfX84+43B5gDEBVl+xkSAr7eXsjR6gZmDY82OkqbxAT58MnNyWzdf5Ts\n4ipyS6sprqhlemIkidEBRsczDXdXCw9NieMP81P5ZFM+s5PN9Xsg2kdrphiOB7K11sUASqlFwEjg\nv0pca/028DZAUlKSDNJ1gI835hEb7MOIHua73qTFohgSFcCQKCntszEpIZTh3QN5YcUuLhoQhr+3\nu9GRRAdrzeuvPCBZKeWtmuevjQMybRtLnE5+aTWpuUe4IjFS5lA7MaUUj1+UQFlNA8/LlEOndNoS\n11qnAAuBLcCOlp9528a5xGl83TLn+sIBYQYnEUaL7+bLrORoPkrJJaNAphw6m1a9E6K1flxr3Vdr\n3U9rPVtrXXf6nxK2tGx7AQMj/YkMPPE64cK53DehN35ebjyxJF2mHDoZeTvbhHJKqkg7UM5FchQu\nWvh7u/PHSX3ZmFPKkm0yA9iZSImb0G9DKVP7S4mL/3fV0Ej6hfvy12WZHK6UF8vOQkrchJZuKyAx\nOoBu/l5GRxF2xMWi+Of0gZTXNPDwoh0yrOIkpMRNZs+hSnYerOACOQoXJxAX5ssDk3qzIqOIz1P3\nGx1HdAApcZP5enshSsEFMh4uTuKmUbEkxwby5JJ08kurjY4jbExK3GSWbS9gaEwgoSY4zV4Yw2JR\nPHfFQCxKcc+CrdQ1NhkdSdiQlLiJ7DxYTtahSpkbLk4rIsCbv13en9TcI9z32TasVhkfd1RyZR8T\nWby1ABeLkvFw0SoXDuhGwdEanlm+k+BOHjx+Ubyc3euApMRNwmrVLNlawLm9gmSJVtFqc0b34FB5\nHXN/ySa4swe3j+lpdCTRzqTETWJL3hEOHK3hgUlypRvRNo9MjaOkso5/freLhiYrd4/rJUfkDkRK\n3CQWby3A083ChHi53qRoG4tF8c8rBuLqYuGlH7IoOFrD05f2x03WH3cIUuIm0NBk5esdhYyPC6WT\nh+wy0XZuLhb+OX0A3fy9eGVlFkXldbxxzRB85PfJ9OR/xSbwy54SSqvqmTYo3OgowsSUUtw3oTd/\nv6w/a7KKuf69jVTWNRodS5wlKXETWLK1AD8vN87rHWx0FOEAZgyL4tWZQ9iSd5Rr30mhorbB6Eji\nLEiJ27ma+ia+Sz/I1P5dcXeV3SXaxwUDwnj96sFs31/GrHc2UlYjRW5W0gp27v11OVTXN3H5kAij\nowgHM7lfGG/OSiSjoIxb56fS0GQ1OpI4A1LiduxQRS2v/7SH8XEhJMUEGh1HOKAJ8aH87bIBrNt7\nmKeWZRgdR5wBeWvajj3/3W7qGpt49IJ4o6MIBzY9MYKdheXM/SWbPl19uXp4lNGRRBvIkbidSjtQ\nxmep+Vw/MobuQT5GxxEO7uGpcZzXO5g/L04jZd9ho+OINpASt0Naa/6yLIMAb3fuGNvL6DjCCbhY\nFK/MHExkoDf3LNgqb3SaiJS4nTlSVc+jX6WxMbuU+yc2X/xWiI7g5+XGi1cNoqi8lr/K+LhpyJi4\nHahvtHKkup5l2wt5+YfdVNU3cf3IGK5KijQ6mnAygyL9ue38nrz20x6m9OvKuLhQoyOJ05AS7wBH\nqupZt/cwa/eWsGHfYcqqG1BKoRTUNjRRUfv/Z82d2yuIxy6Mp3doZwMTC2d217he/JBZxEOLdrDi\nngACfNyNjiROQUrchgrLanjjp70s2JRPfZOVTh6uDOseSLcenlg1aA0erhYCfdwJ9HGnV0gnhnUP\nlBXmhKHcXS08f+VApr22lieWpvPyjMFGRxKnICVuAxW1DTy/Yjcfb8zDatVckRTJ9MQIBkT4ycpx\nwhQSuvlx+5ievLwyi6uGRjKyR5DRkcRJSIm3sz2HKpgzL5Xcw9VckRjB7WN6EhnobXQsIdrs1vN7\nsOjX/Ty+OJ3ld58rByB2SvZKO1q+o5Bpr62lvKaB+TcO5++XD5ACF6bl6ebC4xcmkHWokvfX5hgd\nR5yElHg70Frz6sosbvtoC727dmbpnaMY0aOL0bGEOGvj40MZ2zeEl37YTVF5rdFxxAlIiZ8lq7X5\nxJznv9/NZUPCWTBnBGF+XkbHEqLdPH5RPA1WzTPLM42OIk5ASvwsNDZZeWDhNt5bm8P/nBPDc9MH\nynKxwuFEd/HhD6NjWby1gPV75ZR8eyONc4bKaxu46cPNLNpygPsm9ObPF8ZjscjUQOGYbhvTk8hA\nLx5bnEZ9oyxZa0+kxM9A7uEqLntjHb9klfD0pf24S64eLhycp5sLT1yUwJ5Dlby7NtvoOOIYUuJt\ntG5vCZe8vpbiijo+vHEY1wyPNjqSEB1iXFwo4+NCefmHLAqO1hgdR7SQEm+lqrpGHl+cxtX/TiHQ\nx53Ft58jJ0AIp/P4RfFotCyQZUekxFth7Z4SJr20mg835HL9yBiW3jmKGFnjWzihyEBv7hjTk2/S\nDvJLVonRcQRS4qdUXtvAw4u2c83cFNxdLHx2ywieuDgBb3c50VU4r5vOjSUy0Iu/LsugUa7LabjT\nlrhSqo9SausxH+VKqXs6IpyRVu06xKQXV7NgUz63jI5l+d3nMlSucykEnm4uPDIljl1FFXy6Kd/o\nOE7vtIeUWutdwCAApZQLcAD40sa5DPVxSh6PfLmDXiGdeOPWkQyOCjA6khB2ZXK/rgzvHsgL3+/m\nooHd5OIlBmrrcMo4YK/WOtcWYezBRym5PPLlDsb0CWbpnaOkwIU4AaUUj10Yz5Hqel5dmWV0HKfW\n1hKfAXxiiyD2YN6GXB79Mo2xfUN4a3Yinm4uRkcSwm71C/fjqqRI3l+XQ3ZJldFxnFarS1wp5Q5c\nDHx+ku/PUUptVkptLi4ubq98HWbe+hwe+yqNcX1DeHPWEDxcpcCFOJ37J/bB3dXCcyt2GR3FabXl\nSHwKsEVrXXSib2qt39ZaJ2mtk4KDg9snXQd5f202jy1OZ3xcKG9IgQvRasGdPbhpVHe+3l5I2oEy\no+M4pbaU+EwccChl7pp9PLE0g0kJobxxjRS4EG110+hYArzdePY7ORo3QqtKXCnlDUwAFtk2Tsd6\nb202T32dyZR+XXnt6iGyAqEQZ8DX043bzu/J6t3FrNsrJwB1tFa1lta6WmvdRWvtMK+XPt+cz5NL\nM5ic0JVXZg6WS08JcRZmj4gmzM+TZ7/dhdba6DhOxSmb69u0Qh78Yjvn9gri5ZmDpMCFOEuebi7c\nM74XW/OP8l36Cd82EzbidO31S1YJd32ylUGR/vxrdqKMgQvRTi4fEkFssA/Prdglp+N3IKcq8T2H\nKrn1o1S6B/nw3vXDZA0UIdqRq4uFP03qy55DlXyxZb/RcZyG05T40ep6bvpgEx6uFt65Pgk/bzlN\nWIj2NikhlMFR/rzw/W5q6puMjuMUnKLEG5qs3PbRFgqO1vKv2YlEBHgbHUkIh6SU4qHJfSkqr+O9\ndXIFoI7gFCX+1LIM1u09zN8u609itKxEKIQtDY/t0nzm86q9HKmqNzqOw3P4El+RfpAP1udy46ju\nXJ4YYXQcIZzCnyb3pbKukTd/3mt0FIfn0CV+sKyWP32xnX7hvjw4ua/RcYRwGn26duaSQeHMW59L\nSWWd0XEcmsOWuNWque+zrdQ1WHl5xmA5G1OIDnb7mJ7UNjYxd42MjduSwzbb22v2sW7vYZ64OJ4e\nwZ2MjiOE0+kZ0okLB3Rj3vocGRu3IYcs8bQDZTy/YhdT+nXlyqRIo+MI4bTuHNuTqvom3l0rR+O2\n4nAlXlPfxF2f/koXHw/+dll/lFJGRxLCafUO7cyUfl15f20OZTUNRsdxSA5X4k8vz2BfcRXPXzkQ\nf293o+MI4fTuGNuTirpG3l+bY3QUh+RQJf5DRhHzN+Rx87ndOadnkNFxhBBAQjc/xseF8u7abKrq\nGo2O43AcpsSPVNXz4BfbiQvz5YFJfYyOI4Q4xm1jelBW08AnG/OMjuJwHKbEn/9+F0drGnjhyoGy\nMqEQdmZIVADJsYHMXZNNfaOscNieHKLE0wvK+Dglj9nJ0cSF+RodRwhxAree35OD5bV8tfWA0VEc\niulLXGvNE0vS8fd2597xvY2OI4Q4idG9gogP8+Wtn/ditcrVf9qL6Ut8ybYCNuUc4Y+T+sjyskLY\nMaUUt57fg33FVazIkKv/tBdTl3hVXSPPLM+kf7ifnNQjhAlM6deV6C7evLlqj1yLs52YusQXbdlP\nUXkdf74oHheLnNQjhL1zdbFw87mxbNtfxoZ9pUbHcQimLvFPN+WT0M2XoTGyRrgQZjE9MYIuPu68\nvVqWqW0Ppi3xtANlpBeUc9VQGUYRwkw83Vy4bmQMP+0qZtfBCqPjmJ5pS3zBpnw8XC1MGxhudBQh\nRBvNTo7Gy82Ft1fvMzqK6ZmyxGsbmvhq6wGm9g+TGSlCmFCAjztXDY1kybYDHCyrNTqOqZmyxL9J\nK6SitlFmpAhhYjeO6o5Vw3uyTO1ZMWWJf7oxn5gu3iTHyhuaQphVZKA3U/uH8VFKHuW1skztmTJd\niWeXVJGSXcqVQyNlrXAhTG7OubFU1jWycPN+o6OYlulK/MtfD2BRMH2IXLleCLPrH+FHUnQAH6zP\noUlOxT8jpivxFekHSYoJJMTX0+goQoh2cP05MeQermbVrkNGRzElU5V43uFqdh6sYGJ8qNFRhBDt\nZFJCV8L8PHl/XY7RUUzJVCW+IuMgABPjuxqcRAjRXtxcLMxKjmZNVglZRXLyT1uZrMSL6Nu1M1Fd\nvI2OIoRoRzOHReHuapGj8TNgmhI/XFnH5pxSJibIUbgQjibQx51LBnVj0ZYDlFXLdMO2ME2Jr9x5\nCKtGxsOFcFDXj+xOTUMTn23ONzqKqZimxFekFxHu70VCN7n8mhCOKL6bL0NjAvgoJVeu/NMGrSpx\npZS/UmqhUmqnUipTKTXC1sGOVV3fyJqsYibEh8oJPkI4sFnJ0eQcrmbNnhKjo5hGa4/EXwa+1Vr3\nBQYCmbaL9Hurd5dQ12hlYoIMpQjhyCb360pQJ3fmrc81OoppnLbElVK+wGjgHQCtdb3W+qitgx1r\nZWYRvp6uDJOLPwjh0DxcXbhqaCQ/7ixi/5Fqo+OYQmuOxGOBYuA9pdSvSqm5Simf4++klJqjlNqs\nlNpcXFzcriE35pQyPLYLri6mGcIXQpyhq4dHA/BxSp7BScyhNa3oCgwB3tRaDwaqgIeOv5PW+m2t\ndZLWOik4OLjdAh4qryX3cLUchQvhJML9vRgXF8qCTfnUNTYZHcfutabE9wP7tdYpLV8vpLnUO8TG\nnOaLqQ7tLiUuhLOYnRzN4ap6vtlx0Ogodu+0Ja61PgjkK6X6tNw0DsiwaapjbMouxcvNRaYWCuFE\nRvUMonuQD/M2yBucp9PaQeY7gY+UUtuBQcAztov03zbmHGFItD9uMh4uhNOwWBTXDI8iNfcI6QVl\nRsexa61qRq311pbx7gFa60u01kdsHQygrKaBnQfLGSrj4UI4nSsSI/F0szB/g7zBeSp2fXi7JfcI\nWiNvagrhhPy83bh4YDe++vWAXL7tFOy6xDfmlOJqUQyOCjA6ihDCALOTY6hpaOKLVLl828nYdYlv\nzimlX7gfXu4uRkcRQhigf4QfgyL9mbchF61lPZUTsdsSr21oYlt+GcNkaqEQTm12cjT7iqtYv/ew\n0VHskt2W+Pb9ZdQ3WeVNTSGc3AUDwgjwduOD9TlGR7FLdlvim1pO8kmKlvFwIZyZp5sLM4dF8X1G\nEfmlsp7K8ey2xDdml9I7tBMBPu5GRxFCGOzaETFYlJLLt52AXZZ4Y5OV1NwjMpQihACgq58nFwwI\nY8GmfCpkuuF/scsSzyysoLKukeGxXYyOIoSwEzec053KukY+3yzTDY9llyWekt38LvRwmZkihGgx\nMNKfpOgA3l+XQ5Ncvu0/7LTES4np4k2or6fRUYQQduSGUd3JK61mZWaR0VHsht2VuNWq2ZRTKvPD\nhRC/MzE+lHB/L+b+km10FLthdyW++1AFR6sbGNZdxsOFEP/N1cXC/5wTw8bsUrbkdcg6fHbP7kp8\nY3bz/HAZDxdCnMjMYVH4e7vx5qq9RkexC3ZX4in7Sunm50lEgJfRUYQQdsjHw5XrRsTwfUYRWUUV\nRscxnF2VuNaalOzm8XCllNFxhBB26rqRMXi5ufDmz3I0blclnl1SRUllncwPF0KcUqCPOzOGRbJk\nawH7jzj3qfh2VeIpLePhMjNFCHE6N50bC8DcNc49U8WuSnxjdilBnTyIDfIxOooQws6F+3txyeBw\nPt2UR1F5rdFxDGN3JT5cxsOFEK1059ieNFk1z6/YZXQUw9hNidc2NHFOzy5MTAg1OooQwiSiu/hw\n3YgYPk/dT2ZhudFxDGE3Je7p5sKz0wcybVC40VGEECZyx9ie+Hq68czyTKe8hJvdlLgQQpwJf293\n7hrXizVZJfy8u9joOB1OSlwIYXqzk6OJ7uLNM8szaWyyGh2nQ0mJCyFMz93VwsNT+rK7qNLprv4j\nJS6EcAiTEroyPi6E51bsIu/0bzqvAAAIIUlEQVSw85wAJCUuhHAISin+ekk/XC0WHlq03Wne5JQS\nF0I4jDA/Lx6e2pd1ew/z2eZ8o+N0CClxIYRDmTk0iuHdA3nq60ynOJNTSlwI4VAsFsXfLx9AfaOV\nvyzNMDqOzUmJCyEcTvcgH+4c25OvdxTy407Hvh6nlLgQwiHNGd2DniGdeOyrdKrrG42OYzNS4kII\nh+TuauGZS/tz4GgNL/+QZXQcm5ESF0I4rGHdA7kqKZK5v2STUeCYC2RJiQshHNrDU/vi7+XGvQu2\nUtvQZHScdiclLoRwaP7e7jx/5UB2FVXwpAPOVmlViSulcpRSO5RSW5VSm20dSggh2tP5fUL4w3k9\n+GRjHku3FRgdp1215Uh8jNZ6kNY6yWZphBDCRu6f2JshUf48vGgHuYerjI7TbmQ4RQjhFNxcLLwy\nczAuFsWt87c4zLTD1pa4BlYopVKVUnNOdAel1Byl1Gal1ObiYudbmF0IYf8iArx5acYgMg+W88fP\nHWORrNaW+Dla6yHAFOB2pdTo4++gtX5ba52ktU4KDg5u15BCCNFexvQJ4cHJffl6RyGv/7TH6Dhn\nrVUlrrUuaPl8CPgSGGbLUEIIYUu3jI5l2qBuPLdiN99nmPu0/NOWuFLKRynV+bd/AxOBNFsHE0II\nW1FK8Y/LB9A/3I+7P/2V1NwjRkc6Y605Eg8FflFKbQM2Al9rrb+1bSwhhLAtTzcX3rkuiZDOHlz/\n3kbSC8qMjnRGTlviWut9WuuBLR8JWuunOyKYEELYWoivJ/NvGk4nD1eufWcjew5VGh2pzWSKoRDC\nqUUEePPRTcNRCq6Zu8F0a6xIiQshnF5scCfm3zQchWL6W+tYkX7Q6EitJiUuhBBA366+LLnjHHqF\ndOKW+am8uWqvKeaRS4kLIUSLEF9PFtwyggv6h/GPb3cy650U8kurjY51SlLiQghxDE83F16dOZin\nLunHtvwyJr64mnd/yaaxyWp0tBOSEhdCiOMopZiVHM2Ke0czPDaQvyzLYPSzP/HGqj2UVtUbHe+/\nKFuM+SQlJenNm2XFWiGE+WmtWZl5iPfWZbN2z2HcXS1cOCCMWcnRDI70RynVLs+jlEo9k1ViXdvl\n2YUQwkEppRgfH8r4+FCyiir4YH0OX245wKItB4gL8+WGc2K4bEgELpb2KfM255MjcSGEaJvKukaW\nbC1g3oZcMgvL6R3aiQcn92Vs35AzPjI/0yNxGRMXQog26uThytXDo1h+1yhev3oI9Y1WbvxgMzPe\n3tDh1/GU4RQhhDhDSikuGBDGxIRQPt2YR3pBOZ5uLh2aQUpcCCHOkpuLhdkjYgx5bhlOEUIIE5MS\nF0IIE5MSF0IIE5MSF0IIE5MSF0IIE5MSF0IIE5MSF0IIE5MSF0IIE7PJ2ilKqWIg9wx/PAgoacc4\nZuLM2w7Ovf2y7c7rt+2P1loHt/WHbVLiZ0MptflMFoFxBM687eDc2y/b7pzbDme//TKcIoQQJiYl\nLoQQJmaPJf620QEM5MzbDs69/bLtzuustt/uxsSFEEK0nj0eiQshhGgluylxpdRkpdQupdQepdRD\nRuexNaVUpFLqJ6VUplIqXSl1d8vtgUqp75VSWS2fA4zOaitKKRel1K9KqWUtX3dXSqW0bPsCpZS7\n0RltQSnlr5RaqJTa2bL/RzjZfr+35Xc+TSn1iVLK05H3vVLqXaXUIaVU2jG3nXB/q2avtPTgdqXU\nkNM9vl2UuFLKBXgdmALEAzOVUvHGprK5RuB+rXUckAzc3rLNDwErtda9gJUtXzuqu4HMY77+B/Bi\ny7YfAW40JJXtvQx8q7XuCwyk+b+BU+x3pVQ4cBeQpLXuB7gAM3Dsff8+MPm42062v6cAvVo+5gBv\nnu7B7aLEgWHAHq31Pq11PfApMM3gTDaltS7UWm9p+XcFzX/I4TRv9wctd/sAuMSYhLallIoALgDm\ntnytgLHAwpa7OOS2K6V8gdHAOwBa63qt9VGcZL+3cAW8lFKugDdQiAPve631aqD0uJtPtr+nAR/q\nZhsAf6VU2Kke315KPBzIP+br/S23OQWlVAwwGEgBQrXWhdBc9ECIccls6iXgT4C15esuwFGtdWPL\n1476OxALFAPvtQwlzVVK+eAk+11rfQB4DsijubzLgFScY98f62T7u81daC8lrk5wm1NMm1FKdQK+\nAO7RWpcbnacjKKUuBA5prVOPvfkEd3XE3wFXYAjwptZ6MFCFgw6dnEjL2O80oDvQDfCheQjheI64\n71ujzX8H9lLi+4HIY76OAAoMytJhlFJuNBf4R1rrRS03F/328qnl8yGj8tnQOcDFSqkcmofOxtJ8\nZO7f8hIbHPd3YD+wX2ud0vL1QppL3Rn2O8B4IFtrXay1bgAWASNxjn1/rJPt7zZ3ob2U+CagV8s7\n1O40v9GxxOBMNtUyBvwOkKm1fuGYby0Brmv593XA4o7OZmta64e11hFa6xia9/WPWutrgJ+A6S13\nc9RtPwjkK6X6tNw0DsjACfZ7izwgWSnl3fI38Nv2O/y+P87J9vcS4NqWWSrJQNlvwy4npbW2iw9g\nKrAb2As8anSeDtjeUTS/TNoObG35mErz2PBKIKvlc6DRWW383+F8YFnLv2OBjcAe4HPAw+h8Ntrm\nQcDmln3/FRDgTPsdeBLYCaQB8wAPR973wCc0j/830HykfePJ9jfNwymvt/TgDppn8Zzy8eWMTSGE\nMDF7GU4RQghxBqTEhRDCxKTEhRDCxKTEhRDCxKTEhRDCxKTEhRDCxKTEhRDCxKTEhRDCxP4PMkaF\ngfD82AsAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fe594124b38>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.plot(dset[:,30,40])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:01.923834Z",
     "start_time": "2017-07-31T16:59:01.916817-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "[64, 64, 84]"
      ]
     },
     "execution_count": 10,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dset.chunks"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:02.405581Z",
     "start_time": "2017-07-31T16:59:01.927844-04:00"
    }
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "CPU times: user 8 ms, sys: 0 ns, total: 8 ms\n",
      "Wall time: 15.7 ms\n"
     ]
    }
   ],
   "source": [
    "%time data =  dset[55,:,:]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:02.423628Z",
     "start_time": "2017-07-31T16:59:02.416611-04:00"
    }
   },
   "outputs": [],
   "source": [
    "data = np.ma.masked_equal(data, dset.fillvalue)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "metadata": {
    "ExecuteTime": {
     "end_time": "2017-07-31T20:59:02.688332Z",
     "start_time": "2017-07-31T16:59:02.427639-04:00"
    }
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7fe588209eb8>"
      ]
     },
     "execution_count": 13,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUIAAAD8CAYAAAACGq0tAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAEFlJREFUeJzt3X+s3XV9x/HniwIiqIPWgpWSFRPC\nMGaAu2E4FjNBHDoj/IGLzizN0qT/uA2niYJLlpjsD00Wf/yxYBpRm8XxQ9RBiLM2FbK4LMWLVC0W\nLCLD2kK7TuaPZSrlvT/Ot3p6aXtP7z3fc+7t5/lIbs75fu/33O8753z76ufz+X6/n5OqQpJadsq0\nC5CkaTMIJTXPIJTUPINQUvMMQknNMwglNc8glNQ8g1BS8wxCSc07dZI7e/nKFbXugtPG9wdPe83o\n2/5q5/j2O58TqetYhusd99/TCfnet8+cdglaoP/j5/yyfpH5tptoEK674DQe3HLB2P7eKa+YHXnb\n55++aGz7nc+J1HUsw/WO++/pxPzxKy+bdglaoO21baTt7BpLat5EW4STNNkW4O6Rtx13Xbb0pMWz\nRSipeQahpOYZhJKad9KOEU5S3+N0jgNOlmeJ22OLUFLzRgrCJGcnuTvJo0l2JXldkpVJtibZ3T2e\n03exktSHUbvGnwC+UlU3JjkdOBP4ILCtqj6c5GbgZuADPdV5VOPoMs699OV4f/NELpM5Fru50tIz\nb4swycuA1wO3AVTVL6vqWeB6YHO32Wbghr6KlKQ+jdI1fhVwAPhMkoeTfCrJWcB5VbUPoHs892gv\nTrIxyWyS2QMHD42tcEkal1GC8FTgtcCtVXU58HMG3eCRVNWmqpqpqpnVq1YssExJ6s8oY4R7gD1V\ntb1bvptBED6TZE1V7UuyBtjfV5HDjjdO18eY4WI5JigtffO2CKvqaeCHSS7uVl0DfBe4F1jfrVsP\n3NNLhZLUs1HPGv8V8LnujPETwF8wCNG7kmwAngLe3k+JktSvkYKwqnYAM0f51TXjLWd+C+1qjrvL\nO8zu78lly94dRyx7p8nJzztLJDXPIJTUvCU56UKf3VhJmssWoaTmGYSSmmcQSmrekhkjXCrjgl4K\nIy+XaY8tQknNMwglNW/JdI2Hu6QL7SbbrdXx2OXVsdgilNQ8g1BS8wxCSc1bMmOEw/r+AiUtP8Of\n+7WnOOObxssWoaTmGYSSmrcku8bH4yUybRr+3Lc+P/rwiN1ojcIWoaTmGYSSmrfsusY6uYx6FcCR\nXVy7uxovW4SSmmcQSmqeQSipeY4RaqqGL4txdhhNiy1CSc0bqUWY5Engp8Ah4LmqmkmyErgTWAc8\nCfxpVf24nzIlqT8n0iJ8Q1VdVlUz3fLNwLaqugjY1i1L0rKzmK7x9cDm7vlm4IbFlyNJkzdqEBbw\n1SQPJdnYrTuvqvYBdI/nHu2FSTYmmU0ye+DgocVXLEljNupZ46uqam+Sc4GtSR4ddQdVtQnYBDBz\n6Rm1gBolqVcjtQiram/3uB/4EnAF8EySNQDd4/6+ipSkPs0bhEnOSvLSw8+BNwE7gXuB9d1m64F7\n+ipSkvo0Stf4POBLSQ5v/89V9ZUk3wDuSrIBeArvhJe0TM0bhFX1BHDpUdYfBK7poyhJmiTvLJHU\nPINQUvMMQknNc/YZHdfwjDBb9u7o9e9L02KLUFLzDEJJzbNrrBcY/kKlLXtH+x7p43Vxtz7/+V8/\n93uGtRTZIpTUPINQUvPsGusFhruvW/aO9pq5Z5Q9G6zlxBahpOYZhJKaZxBKap5jhHqBY13uciJ3\nlgxv6yUzWupsEUpqnkEoqXl2jTWy4TtOAJ5/erS7TqSlzhahpOYZhJKaZxBKap5jhHqBI8f+fnOr\n3NzLYEa9/U5a6mwRSmqeQSipeXaNdVxOqqoW2CKU1LyRgzDJiiQPJ7mvW74wyfYku5PcmeT0/sqU\npP6cSNf4JmAX8LJu+SPAx6rqjiSfBDYAt465Pk3Zsc4gw5GTrw53ocFutJaXkVqESdYCfwJ8qlsO\ncDVwd7fJZuCGPgqUpL6N2jX+OPB+4PlueRXwbFU91y3vAc4/2guTbEwym2T2wMFDiypWkvowbxAm\neSuwv6oeGl59lE3raK+vqk1VNVNVM6tXrVhgmZLUn1HGCK8C3pbkLcAZDMYIPw6cneTUrlW4FvA+\ng8Z4aY1OFvO2CKvqlqpaW1XrgHcAX6uqdwH3Azd2m60H7umtSknq0WKuI/wA8N4kjzMYM7xtPCVJ\n0mSd0J0lVfUA8ED3/AngivGXpKXkyMlYj+z+2h3WycI7SyQ1zyCU1DyDUFLzDEJJzTMIJTXPIJTU\nPCdm1XENzz5zvO8oGZ6JRlpubBFKap5BKKl5BqGk5hmEkppnEEpqnkEoqXlePqORHTkTjbPP6ORh\ni1BS8wxCSc2za6yR2RXWycoWoaTmGYSSmmcQSmqeY4Qa2Za9O475O2ef0XJmi1BS8wxCSc0zCCU1\nzyCU1Lx5gzDJGUkeTPKtJI8k+VC3/sIk25PsTnJnktP7L1eSxm+Us8a/AK6uqp8lOQ34epJ/Bd4L\nfKyq7kjySWADcGuPtWpCPAOs1szbIqyBn3WLp3U/BVwN3N2t3wzc0EuFktSzkcYIk6xIsgPYD2wF\nvg88W1XPdZvsAc4/xms3JplNMnvg4KFx1CxJYzVSEFbVoaq6DFgLXAFccrTNjvHaTVU1U1Uzq1et\nWHilktSTEzprXFXPAg8AVwJnJzk8xrgWOM633krS0jXKWePVSc7unr8YeCOwC7gfuLHbbD1wT19F\nSlKfRjlrvAbYnGQFg+C8q6ruS/Jd4I4kfw88DNzWY52S1Jt5g7Cqvg1cfpT1TzAYL9QyNPf7R47k\nBKxqi3eWSGqeQSipeQahpOY5MetJZvTb4xwHlA6zRSipeQahpObZNV7mnClGWjxbhJKaZxBKap5d\n42XI7rA0XrYIJTXPIJTUPINQUvMcI1wiHPeTpscWoaTmGYSSmmfXuGd2eaWlzxahpOYZhJKaZxBK\nap5jhGPgOKC0vNkilNQ8g1BS8+wanwC7wNLJad4WYZILktyfZFeSR5Lc1K1fmWRrkt3d4zn9lytJ\n4zdK1/g54H1VdQlwJfDuJK8Gbga2VdVFwLZuWZKWnXmDsKr2VdU3u+c/BXYB5wPXA5u7zTYDN/RV\npCT16YROliRZB1wObAfOq6p9MAhL4NxxFydJkzByECZ5CfAF4D1V9ZMTeN3GJLNJZg8cPLSQGiWp\nVyMFYZLTGITg56rqi93qZ5Ks6X6/Bth/tNdW1aaqmqmqmdWrVoyjZkkaq3kvn0kS4DZgV1V9dOhX\n9wLrgQ93j/f0UuEUebmM1IZRriO8Cvhz4DtJdnTrPsggAO9KsgF4Cnh7PyVKUr/mDcKq+jqQY/z6\nmvGWI0mT18ydJXZzJR2L9xpLap5BKKl5BqGk5p20Y4SOCUoalS1CSc0zCCU1b1l3je3+ShoHW4SS\nmmcQSmqeQSipeQahpOYZhJKaZxBKap5BKKl5BqGk5hmEkpq37O4s8W4SSeNmi1BS8wxCSc0zCCU1\nzyCU1DyDUFLzDEJJzTMIJTVv3iBM8ukk+5PsHFq3MsnWJLu7x3P6LVOS+jNKi/CzwHVz1t0MbKuq\ni4Bt3bIkLUvzBmFV/Rvw33NWXw9s7p5vBm4Yc12SNDELHSM8r6r2AXSP546vJEmarN5PliTZmGQ2\nyeyBg4f63p0knbCFBuEzSdYAdI/7j7VhVW2qqpmqmlm9asUCdydJ/VloEN4LrO+erwfuGU85kjR5\no1w+czvwH8DFSfYk2QB8GLg2yW7g2m5ZkpaleecjrKp3HuNX14y5FkmaimUxMespr9g9tPT2qdUh\n6eTkLXaSmmcQSmqeQSipectijPD5py/69fMte3+z3i9ykjQOtgglNc8glNQ8g1BS8wxCSc0zCCU1\nzyCU1DyDUFLzDEJJzTMIJTXPIJTUPINQUvMMQknNWxaTLhzLlr07jlh2EgZJC2GLUFLzDEJJzTMI\nJTXPIJTUPINQUvMMQknNMwglNW9RQZjkuiSPJXk8yc3jKkqSJmnBQZhkBfCPwJuBVwPvTPLqcRUm\nSZOymBbhFcDjVfVEVf0SuAO4fjxlSdLkLCYIzwd+OLS8p1snScvKYoIwR1lXL9go2ZhkNsnsgYOH\nFrE7SerHYoJwD3DB0PJaYO/cjapqU1XNVNXM6lUrFrE7SepHql7QiBvthcmpwPeAa4AfAd8A/qyq\nHjnOaw4A/wm8HPivBe14fJZCDWAdc1nHkazjSCdax29X1er5NlrwNFxV9VySvwS2ACuATx8vBLvX\nrAZIMltVMwvd9zgshRqswzqsY2nUsaj5CKvqy8CXx1SLJE2Fd5ZIat60gnDTlPY7bCnUANYxl3Uc\nyTqO1EsdCz5ZIkknC7vGkpo30SCc1iQNST6dZH+SnUPrVibZmmR393jOBOq4IMn9SXYleSTJTdOo\nJckZSR5M8q2ujg916y9Msr2r484kp/dZx1A9K5I8nOS+adWR5Mkk30myI8lst24ax8jZSe5O8mh3\nnLxuCsfHxd37cPjnJ0neM6X342+6Y3Rnktu7Y3fsx8fEgnDKkzR8FrhuzrqbgW1VdRGwrVvu23PA\n+6rqEuBK4N3dezDpWn4BXF1VlwKXAdcluRL4CPCxro4fAxt6ruOwm4BdQ8vTquMNVXXZ0OUZ0zhG\nPgF8pap+B7iUwfsy0Tqq6rHufbgM+D3gf4EvTbqOJOcDfw3MVNVrGFym9w76OD6qaiI/wOuALUPL\ntwC3THD/64CdQ8uPAWu652uAxyZVy1AN9wDXTrMW4Ezgm8DvM7hQ9dSjfV497n8tg39UVwP3Mbh1\ncxp1PAm8fM66iX4uwMuAH9CN3U+rjjn7fhPw71N6Pw7PZ7CSwaV+9wF/3MfxMcmu8VKbpOG8qtoH\n0D2eO8mdJ1kHXA5sn0YtXXd0B7Af2Ap8H3i2qp7rNpnU5/Nx4P3A893yqinVUcBXkzyUZGO3btKf\ny6uAA8BnuqGCTyU5awp1DHsHcHv3fKJ1VNWPgH8AngL2Af8DPEQPx8ckg3CkSRpakOQlwBeA91TV\nT6ZRQ1UdqkHXZy2DKdUuOdpmfdaQ5K3A/qp6aHj1pOvoXFVVr2UwdPPuJK+fwD7nOhV4LXBrVV0O\n/JzJdMePqht7exvw+Snt/xwGU/tdCLwSOIvB5zPXoo+PSQbhSJM0TNAzSdYAdI/7J7HTJKcxCMHP\nVdUXp1kLQFU9CzzAYMzy7O4ecpjM53MV8LYkTzKYz/JqBi3ESddBVe3tHvczGA+7gsl/LnuAPVW1\nvVu+m0EwTuv4eDPwzap6pluedB1vBH5QVQeq6lfAF4E/oIfjY5JB+A3gou6Mz+kMmtz3TnD/c90L\nrO+er2cwXterJAFuA3ZV1UenVUuS1UnO7p6/mMEBtwu4H7hxUnVU1S1Vtbaq1jE4Hr5WVe+adB1J\nzkry0sPPGYyL7WTCn0tVPQ38MMnF3aprgO9Ouo4h7+Q33WKmUMdTwJVJzuz+7Rx+P8Z/fExq0LUb\n2HwLgxlrvg/87QT3ezuDMYZfMfhfdwODsahtwO7uceUE6vhDBs34bwM7up+3TLoW4HeBh7s6dgJ/\n161/FfAg8DiD7tCLJvgZ/RFw3zTq6Pb3re7nkcPH5pSOkcuA2e6z+RfgnCnVcSZwEPitoXXTqOND\nwKPdcfpPwIv6OD68s0RS87yzRFLzDEJJzTMIJTXPIJTUPINQUvMMQknNMwglNc8glNS8/wcOyIYm\n+htQ9QAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fe594025828>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(data, origin=\"lower\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "0.0"
      ]
     },
     "execution_count": 14,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "dset.fillvalue"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = np.ma.masked_equal(data, 1.e37)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.image.AxesImage at 0x7fe5881876d8>"
      ]
     },
     "execution_count": 16,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAUIAAAD8CAYAAAACGq0tAAAABHNCSVQICAgIfAhkiAAAAAlwSFlz\nAAALEgAACxIB0t1+/AAAADl0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uIDIuMS4wLCBo\ndHRwOi8vbWF0cGxvdGxpYi5vcmcvpW3flQAAIABJREFUeJztnXusZXd137/rPO9jXh7PtRnbwNgJ\nIc6jwWRESJyixA4wpChQiURQFFmNI1cVSUmTKhgitUoVtZAmQKRUpBaQWFXKow6pEQ0mlgOqElVO\nhkASgwETcLFjx3OxPeOZO/dxzt6rf5wzZ3/X2mf/Zt/XuTPe6yONZu/z2+e319nnzG/2d+31EFVF\nEARBk2nttQFBEAR7TSyEQRA0nlgIgyBoPLEQBkHQeGIhDIKg8cRCGARB44mFMAiCxhMLYRAEjScW\nwiAIGk9nlic7cuSIHjt2bJanDIJt87Wv/eNk2+RhiTvQ7xNqxhIH1mb7GWFSc4q8vRP27g3rK89g\nuLZy0Q8w04Xw2LFjOHny5CxPGQTb5pZb//NkmxcF7VhBpbSrbvEwY63tLyxSNzU2cZjkiWPJxPWD\neygcK2yqy5f+9/trHRfSOAiCxjPTO8IguBw5d01vss13ern712PvCBNjQrc2qbscdzdn7gJp09/Z\nGcnrx/Lp26WxxB2nZMV2K/Njxft4vtbQzic57afsT9AaFG9srw7NWOfLj47+fm6t3lz1ThkEQfD8\nJRbCIAgaTyyEQRA0nvARBoHj1h//T2Z/+JK5yTb7BbVjHXxmzN1imH2p2AasX7DkI5SpY6mnvyU/\noPHvaeWY5DL1dQBgF6e2nJFdep+xw35Q9iWW/IzkPxRy/UmW8Fvm2wsnijvCIAgaT62FUEQOicg9\nIvIVEXlYRH5YRA6LyP0i8sj47yt229ggCILdoK40/h0A96nqm0SkB2ABwLsAPKCq7xaROwHcCeAd\nu2RnEMyMlaN9s5/1OIi6eD3v2velpDFLYK0ZGFySvBVjpXCTxBi/TzMvV3m7Wrra+VPpNNPP6+fM\nS/K9mLPVITsG7jiS2yYkCQAmUrmeZL7oHaGIHADwKgAfAgBV3VDV0wDeAODu8WF3A3hjrTMGQRBc\nYtSRxjcAWAbw+yLyBRH5oIgsArhaVZ8EgPHfV017s4jcISInReTk8vLyjhkeBEGwU9RZCDsAXg7g\nA6p6E4AVjGRwLVT1LlU9rqrHl5aWtmhmEATB7lHHR/g4gMdV9cHx/j0YLYRPichRVX1SRI4COLVb\nRgbBLGlveMdascl+QJ9iZ3yELsUOpuhCPb+VDzmp9P1tNY3OTS9ko3BxCe8jZN+fT8Wr8At6OzhV\nMRXio5w5507V4uIVvpDFZL+eQ/aid4Sq+o8AHhORl45fuhXAlwF8EsBt49duA3BvrTMGQRBcYtR9\navyLAP5w/MT4GwD+JUaL6MdF5HYA3wLw07tjYhAEwe5SayFU1S8COD5l6NadNScI9p6/uOffmf2b\n/vV7J9tcSzCzUTYutMbqOJbKtcNnvBZkqWkyRNyELCdTxWPdmDlbKjwnlTFSYaPPCmmZLBZnIslt\nk51SCgXiCjaRWRIEQbAtYiEMgqDxRNGFILgIh79cFPc88x1FAQYvC1cpkrZUtJUyJEpZJ4SRfz7z\no0KS+qIIJovFF3/gmrCpsbpKM5G5wpkgraHT4YPpRWbHs9B8dN12sXVK3BEGQdB4YiEMgqDxxEIY\nBEHjCR9hEDhOfO+77AtX7Ztszj9dOOR6Z63TarhQxMiUQms4fMYUafUhMuzgcyEng+lOslKoTsrX\nlxqkMTOU8hd6HyFnhbSmb49OxZVj3FhLpx7nz9XuFZP61qo7Xn0mCILg+U4shEEQNJ6QxkFjeM0P\n/cfJ9sYVVrsO54t7gpV/esSMmbCVbPo2AHTPFQcOF6ze25gjudfl9AtnJGdZOM2Y9cmQHg2kYl28\n7kxVe63K4kj0HU4VTBAfMsOkekDT+1gm+yKwQzqu3bf3dK5m7kWJO8IgCBpPLIRBEDSeWAiDIGg8\n4SMMLjtOvNT2CNPFIu0tW+yZsY2Dhbfo2R/ZP9ke7DeHmUZM2k4UGzX+wurKK8afB1eMldLt4Iq0\nKvnBuHhpiUQFGGUjN3Orw++r6KEMAEJhPD6kx4zxdXSFU7ngaivR6Io/W+6+F74+yWtVg7gjDIKg\n8cRCGARB4wlpHFx23PfV99Q+9sZfe99ke/1IoV2zed9od/N2eGkMzqrou/lJDkunGOv0h+awNsm/\nbGiNMrK5ro01j/Pz5xzCsmHjWxKRNbYPyoCzR1yWDB1Y6mdS1QPaf5hE0dYd71kSBEHwfCcWwiAI\nGk9I4+Cy4Ef/+X+ZbK8ftFJt9Ughf9aWrEYaXlvoLp0vtGurb9NCpFXxxBSA1KxSyse1O3nlGG/P\n9QbmuPmulcqMkaQsk529qbGMnsIOc3sfNBwW13XYogIS7sl2PijG1Fd3pd3cvOyeGnPBVffEt7pA\nhT1VUvW22xc/xtgXBEHQcGIhDIKg8cRCGARB4wkfYXDJ8CM/81tmf/Vw8f/0xkuLn+rGAeuzGuyf\n7gcEgNZ84XPr9oqxTsce124Xc7R9qEeLxshf1m5ZP2CX9rttZwf7DxMVYOba1T7CHNP9ey1nL/sI\nB7nzpw6KFJr1YfU/f+MXFVcehvxu/pNohZNQfYgMhR6VmllV+PXKBVx521fZ2VymSdwRBkHQeGrd\nEYrIowDOYhQyOlTV4yJyGMDHABwD8CiAn1HVZ3fHzCAIgt1jM9L4x1X127R/J4AHVPXdInLneP8d\n098aBBdnfb8VKKtXkcTbT0VP97vQl4Vivztnw1H6lLnR6xTbPkxlvlu8b7GzYcYWaL/TKs7Vb1k7\nurTfcVVbu7SfkRBrO3GZJURaXtHYd+gqm65nxT/r1cyWKD3fKYpSnB3Y4rTnWTYPai4NvuWKaVfc\noterC1TkuZ2EM3aEzM9dgQfuHZ13vDS+cO7dzyx5A4C7x9t3A3jjNuYKgiDYM+ouhArgT0Xk8yJy\nx/i1q1X1SQAY/33VtDeKyB0iclJETi4vL2/f4iAIgh2mrjS+WVWfEJGrANwvIl+pewJVvQvAXQBw\n/PjxurniQRAEM6PWQqiqT4z/PiUifwzgFQCeEpGjqvqkiBwFcGoX7QwawGC/a86zjyqx7CucSrLP\n+vfmFgofnk9Z29cvxg7NrRbbvfP2OPIDLrbX7Rjtz7WK+bu+exPRb1k75qTYH5BPb02tD2+dKsQO\nnO8vo3gRDqXJXPzJerv4Z72e23/iK+3CL9hzoTp9GjvfLnyJKy7Fbo32B25sUBFa41o0G19oKnzG\nFGZ1kUXs/sy7zgfZbpXmSnFRaSwiiyKy/8I2gNcAeAjAJwHcNj7sNgD31jtlEATBpUWdO8KrAfyx\njAIUOwD+h6reJyJ/BeDjInI7gG8B+OndMzMIgmD3uOhCqKrfAPADU15/GsCtu2FU0EwyG82BnIqZ\nKlWL6blipgcW1ibb+/tW1i7NnZtsH+oV0tjL336rmHOhZcNnWA73E9K4J8UcLIX9HIyXtSyVMxdy\nMtDin2vdMBt+DwCczYr+LufcBT89XCi2N+aL7e68Oe659WKO1Y6df40ydDZa9FnarrgrFaDN2wlh\nmujzzP1SvDTGxK4ozBoEQVCLWAiDIGg8sRAGQdB4ovpMkOQ1vbdMtv904yM7Pv/Nbyoqzgy/3/qR\n8rnCj8QVpef6vqpzsX+Q/IBAtV/wio4NnzlI+973NycUnkO+Pu8H7NH7umL9mDwnp9V5H96KFmEr\na7nt0ZyhOF8b9arZeB/koFNc45Xc+gjPZIWP8JnuYrHdWzTHne4VPsPnBnN2bI3GyF+43rZhQoNW\n8bl9aI1SpWw1/abtcRwZ5DIJgQufs2YVmrgjDIKg8cRCGARB4wlpHCSpK4dfe/DnKsc+c+bDk+1X\nvf43zdiZGymT4oCVeLpI1WIoe2TfnA192d8r9l8wd9aMHeoWkpfDYva318xxLHO9rK2Sw6kQGT/G\nYTG8fTa3oSleyjIsvVulkqjTaTnZzDI6NZaCK+v0XXZKh+boUKHas65Q7flWIcvXfQUbqj7TWieZ\nnCrg6ou2dqJ5UxAEwaaIhTAIgsYT0jgoceLwz0+273vmg7Xew/IXAE58951Tjzt31GUiHKHCCgdc\nMdN5kppUTGGuY+XYUr/IHjnQsU+NjVxtVctfloW+Z0nVcT5bhOVw21Us5SfF57WQhWu5fdzJsrlU\n1IGm5IIMG+7Jc1UBV6Dc34ThObOEpuQ5Ok5Oz3WKa7BA/VJSNnnW6dhsQL1ZbMKP6Yecu97IRV/j\neGocBEFQi1gIgyBoPLEQBkHQeMJHGJRgv+DrrvmFyfann/jd+nN85d2Tbe5XvPadNj1gSH7BzqLL\nGKFGTNxs6cq5FXPc/m4RCuMLkXJVGfa5eV8Zh6Psb7nslLbNQrmAb7x0Ni+yLHzGCMP+t4WWDQXi\nKX3R1tOU+XGOqsicz2wGii/oypj+ys5+HmOf3nmX4bJB13joYlr4fTyf7/Pc61CRWVdMlzNN1ofk\nL9yw14M/tjMR2h/bGD7CIAiCesRCGARB4wlpHKTp+mz2zZP1qFDoopWkOl9IpoUF1yuEMkhesPhc\nsT33nDmOJR5nMwA2xMVnUjAsfw+0bNZJVZGEs2oLDnCRhFK/EbrnSElXphwWU8zBYTfeHZAl+gnz\nHAMna4cU7sJzbib0heVwj9wSw5Y9V5+KsWrX20hyeKGwaThw9q7R2LydQ/uj66MhjYMgCOoRC2EQ\nBI0nFsIgCBpP+AiDNHnhHzux9K/M0H3L/63eFJwK5cIcOnOFH8lXlTlq/IK2qgzDPrc+hpVjHAbi\nGzRxqEqe2/sDDq3J6d4hVSlmf3u1csw0YfKFU8leH1qzSPvcROrZoS2c6tP2qlh11Uyf06ISDofF\nrA237yf2vltftYYxPsI5ut4bzu+6UIwNF1xP7MWxzT71roK4IwyCoPHEQhgEQeMJaRwk+fRjvzPZ\nrqooczHWD5HUmbcZBvPUf2Rf18pVlsPz7WLMh4vYkBArn6pkopfGZ7NCFvqewdyvmOXqopOuXH1m\nUVxvZKmWgsZeks2lwqlUfoY/59O5lcbLwwOTbd+XZJ2ux1mx4T/8uVeGhRvBZ48MsnblWJZPv7cS\nn8mTqILDY/w+6djrkfcpE8atZFl/ZEfdyJ+4IwyCoPHUXghFpC0iXxCRT433rxeRB0XkERH5mIj0\nLjZHEATBpchmpPHbATwM4MJ993sAvE9VPyoivwfgdgAf2GH7gksJ13fxxHf96mT7vq/ZXiQ/9Nbf\nnmxvfCdlEcy7p4dUTOHqBZsxstgppKctqmrlNUs6X0iA33eQnuT6niJMz8nY/ZRpwgUZvDQ+RPsL\nvqgqwRZuOGnJsnCxZEeh87r0dHwAmwmz0v32ZHvZSeNHB0cm208NDjm7ivnPtAtXQWtov3eWw+tD\nu4TwfpYX83lp3G5Rq1Y3tsHtPEnbSssel7eL/VKrz+7Yxp2UxiJyHYB/BuCD430BcAuAe8aH3A3g\njfVOGQRBcGlRVxq/H8CvovjP7EoAp1X1wn9LjwO4dtobReQOETkpIieXl5e3ZWwQBMFucNGFUERe\nD+CUqn6eX55y6NTHQKp6l6oeV9XjS0tLWzQzCIJg96jjI7wZwE+JyE8CmMPIR/h+AIdEpDO+K7wO\nwBO7Z2ZwSeAqedz31fdMtn/i5t8wY2s3FSEdQ6o4Iz3rO1ugopyHe7YAalUf4m7HzsGhJN5/eLhT\nNHY61C4KuvrmSlwstQs7x5X0viUK4znYso6pLqh6i/OnnqYMnafJb3fe+fCYdtsWoD1I3sW+FP90\nF9z3clW7CJG5Rq0v9MrWP0y2H2vbbB32Hw7y6go53NfYh+AMKXxmmBU2csgNAOR5tfOO/YJZRvdq\nPhamxeEz7vucGx2rrR3KLFHVd6rqdap6DMCbAfyZqr4VwGcBvGl82G0A7q11xiAIgkuM7cQRvgPA\nL4vI1zHyGX5oZ0wKgiCYLZvKLFHVzwH43Hj7GwBesfMmBZcqctZKNQ6fGXzHlWaMWmog71KfjK4N\nbznQLySvD31hmcsZHT7cog/uJ2zn4BAXLnDgg1tMMQV3f7CiRYhslhVS67GhlXucxfFMts+Mcb8R\nLvDgJTpLeR/iczW5B9o1C452YG08TAVSN9rnzBi7Bzgjx1/vxXbxWb7tepFwwQQOmVlz/UY2SHoP\n3XVkOax0vTVz9210CXz4zPCCNI7MkiAIgnrEQhgEQeOJhTAIgsYT1WeC+vSsnye/ovCDDfZZJ41W\nRD14f1OvVfiYfOUYTvni9825yi69RDob++14e91VpUk1VGLf2fmsmOPZ4YI57sygSEs7s2HDStj+\ng9SH+eq+TSvkgqu+rzEolY59f0Pn8cy08JOuqytUSz7JPHEfZP2z9nrnHSq060KKOXymJfOogv2H\n6htMkS+Qt3XD2Us9j31k88RHWPNWL+4IgyBoPLEQBkHQeEIaB7X59DffWzn2yrf8ttlnBcz1RX2P\nXK5ksprZSm5nhoW0Ynm22LVVXzi0pu2KmT49LOT7twf7i7kzK9vODgspu565iipU9fPcoJDGp9ft\nHGdWizk2XEgIV9m57uCZybaXxpz94sNnukJVdmh7qC6EheSvl8YreTG24voml6X4+LwtO0ef3Ag9\nN7bQKb6nYSI7pUs9THxlGibn8Bmnc4X6HHvvyKTPcUjjIAiCesRCGARB4wlpHOwI/ukcK2AheZMN\n7IEsQ7lPBgB06YnyWUpV4ewLjy+YcKC9NvW4ZzZsn4+n1grZvDKwdqxRsdHz68XYyqo9brhB/5yc\n2tPF4hqwnDzaO2OOe1Hnmcn2kmsJOifTCzRkrsDDAPzk2boK1hL9XZguFYUt9Vuht/n+MatU8GHY\nSUnjROFaLtwwKObwhRWs/8W5XMYP9OOpcRAEQU1iIQyCoPHEQhgEQeMJH2GwZU4cfdtku/Wq682Y\ncT9xlRBXkJOb/Zx3PsIO+ZEOd4sJnxnayi4tqjhzVccWG12QIrTmxf2iqZHPCnlspWhk9Ox5Gxaz\nPqACo+yzcsVGW9RMqNe3oS8vOvzsZPvlBx+bbH//3GPmuBso66Tr7lP6UvjfBhQy40NkMu5/DAtn\ngvjKN6YvM4UrDVyYTZfON2jbMfY7cr9pn1F0flh8lnU3R4tCa1ptCrNxfY21W8yf9e38rfnILAmC\nINgUsRAGQdB4QhoHtXntwZ8z+63DhZzsPWfDIVqcWcGqZWj/710dFBJpo2el5lpWjJ2jYgfcqxiw\nfT9OZ84OyjThkJCXzD9ljutTSMvX+7bJ2NNrRagNZ8b4EJADvULWXjNvw2L4fN/bL/qG3NC1mSUH\nW8Vn8WEx69R/5GzO0rU6M2PghkxvFpeOMQeWxjZ7h1nJCxeGl8YsvX0WEbORyDrhNiMcFdPyfUl6\nxXfre6AMx72YozBrEARBTWIhDIKg8cRCGARB4wkfYVCbz5z5cOXYie99l9nvvrjokbvKLjcbAZGs\nPDLXLnxW7G/yqWEcTnMms2ExVT2Py/2Pi6ovP3jQ9lfOD1CYBsVjeF8lV8HZ71L7ltqFL/Bwu5jf\n34mcz4s5By68ZSPnoqrVZPS2NV/sln2crlJPm3yEA+qb7I3kCj9r6tIiyX/Ifte+q1LDBXm9r7VN\n4TNtCp/JWtbevMuVaayN+SCqzwRBEGyKWAiDIGg8IY2DnaFV/X+qcthDz8qb/f1CTh7q2WorBzrF\nGEvZZwdW/voKKAxnN7As7LetVJunTIr5tpe8xdhCu7CpXDi1uqLK2ZyyVejUZ12YCld68dK1Cp8h\nwmSumQdnePj3tWmMQ438HLzfc5VpOFyJz9Vx8pel8pz7LnrUK7ldkWUCAC3lUB0zhPzC7yzhejFz\n1ToqCILgecxFF0IRmRORvxSRvxGRL4nIr49fv15EHhSRR0TkYyLSu9hcQRAElyJ1pPE6gFtU9ZyI\ndAH8uYh8GsAvA3ifqn5URH4PwO0APrCLtgYz4sSh24sdauEp+2wx0+zIgcn26g0Hzdj6FSRb5gtJ\nc/iILYrwon1FMYKDXSuNmdMkh88MbKtMzkDx2QxV2Q3cWwOw8sxLY95fJLnunxr3Wda6p6QsQ7ll\np8/gYLm93xVm5aKzfO4eqrNpeu75cj8h3wcV90WZq1zA+y03P7sH2Eb/pH9I13Su465jp7h2HZLJ\nw7a1g7NJtG2/Z71QAGOnMkt0xIWSwN3xHwVwC4B7xq/fDeCN9U4ZBEFwaVHLRygibRH5IoBTAO4H\n8PcATqtO6vE8DuDaivfeISInReTk8vLyTtgcBEGwo9RaCFU1U9WXAbgOwCsA3DjtsIr33qWqx1X1\n+NLS0rRDgiAI9pRNhc+o6mkR+RyAVwI4JCKd8V3hdQCe2AX7gj1ADhSNjLIXXDHZPn+tDVs5v1T4\nfdaOWGfM+aOF72jxmsIveP2hZ8xx3MjorPP9cVgM9z/2/XJXqcjnyoZ9Zsc+Qg7F8L2LTaZD21Vl\nYX9WuwiD8RkR7EvsSHW4CPvRFtrWV7mPMlJ8n2HOXFkEbYudo0vxOe2a4SMAMDCNnagYLVxVILIr\nd/dS5rNR2JH3Mw4p3Mr3RubMlXarsL/VcsVXaV/dWN6u/7mBek+Nl0Tk0Hh7HsBPAHgYwGcBvGl8\n2G0A7t3UmYMgCC4R6twRHgVwt4i0MVo4P66qnxKRLwP4qIj8BoAvAPjQLtoZBEGwa1x0IVTVvwVw\n05TXv4GRvzB4nrFxrPDlnntRIVdXr3Q9iQ8X22vX2RCIK64uigxcuVgUGdhwWSBPrhZhN8PchUdw\noQWSw5k7jt/XbvlCAtPn23D9Rni/PXQFUTuFzecqMicAK699SMgiSWAOwfFzsLRccX2MObQmnU1C\nRSJc/AiHyPjQohUqoMDba7l1N/geJlU2st704TPnW8WcXr5rRciTL9DB0jh33/vk3FGYNQiCoB6x\nEAZB0HhiIQyCoPFE9ZnnGa+79hcn23qgKFia77P+pmy+CIHIFuzP4Pwx8gsuFf9Xrh1xzXOOFaEe\n+xdtIdIOVQrhBk28DdjgU++zYl8g+428r4j9gn0X+uKrzFTB4TneBzkg/+E6+ce8fy+jc7Vc5Zhe\ni31z1fcfttlUdercBvncfJpb2/QutmPsq/O+viq/oPfvMSUbwUVbuQ+znYNDo1Yz+5sYVvRDLvkS\n6TfmmzdJu17lnsl5NnV0EATB85BYCIMgaDwhjS9zTtz4TrOfvfiqyfZwkaTJPlf9Y56yCBasrNg4\nUOyvHS7kyPBFVv7qkI5bc5KXpOwq7BjD0sdLTYalts/86Laqx3w4zQW8DGeJ10kURGXZ5o+bp5CZ\nnpfoFFrD1WfaNYuvAsBaTteRbmF8dgf3J/YZHXzshpOrVWExXuazPyNz5+Y5zpO8PpvZrCHOIlob\nVv8++Pvrduw1FWrO4kNuhheKAUdh1iAIgnrEQhgEQeMJaXwZ8rprfmGyPfjOo2ZscIDk8CIVKphz\n2RjcQmPeyQqqrcB5//qMffIMSmwfuKd22bCQXUIZAD7Qn/tQtN2Tvg5JIc6p9091U8kDg2z6E88q\nyTwNlsD8FHrOFXDl4q6LrvDrPGWWdLgPh8sQ4SfK67mVjKbNJfWv9FkmG5RPkzu3hJfDZn7wE1q6\nPk5ec3bKmrPxbF5I3nMkh1eG9reT6jPD17tb83sa5r5P7A4XXQiCIHi+EwthEASNJxbCIAgaT/gI\nLxFee/DnzL7MFf4V6bvqHze8oNg+YH00VX5B9gmO9gvP2nDBj1FVD3IpebdLzj2Khy5MI6NMEO5r\n7AtocgOevGapEAeHwvgKNlIRnuOP47CbZBgP+aw6zn/FvXt90VYOz+Fsj7xueRRHxpV54MNgqFJP\n4l6nK9UFUc25fAUbmt8Xjz2fFb5A9h8OfBgP2e/Dc/g7M9veLvq9ZJmdf7O/pbgjDIKg8cRCGARB\n4wlpvMu8duFnzT73BpZuIR30u15kjsv6xVeT93xWSLE/XHBhMSSHOZg/67vwFhpzERAwCQaJvhDI\nEvJDSPJSqAc69bMU6pLn3NfC9dnlkJyE5OXMhPIchazl/hq9VnX2iO9ZUpVB48NnmFIx07yQnV2h\n4gnuXKlsFS6K0HVfH2ehsBz2ITK8f94VbeWiFKbnjOszwy6BVBEK/l4y37Oa3Bu5c3Xohd9mRZFX\nT9wRBkHQeGIhDIKg8cRCGARB4wkf4Q7w6s6bzX6b/YAvvs6M5fsL51w216Ft50Ppkb+mb/+/ynrk\nN3G+v2GFX9AV/zB+wbzrCl5WuFXEhyRwhpn7L9X6E+l9YuewVUOsbysXCv8ZJvx75Af0P+gqN6ZP\nsROpSC9DtV/Qh8iYMB7n3+tqcSz7/nwVHPad9V14S9ai75OO81VjOBzFF2Y1hVSdK5ELn7KNK7lN\nj2O/4Dn3wzpDcVqcVudT6rhZ1tD5CNlXOSDf38bQhQnRHL4wq14I56qZaRd3hEEQNJ5YCIMgaDwh\njTfBq9s/M9lu9Yvb/s61tgKM7i9SNbIDVjpk8xRSQHI4c9VhcoptcC0d7JgrCFMlh8tzkL3+V2DU\nKmWIqNUZwiEyvqgJzWlkckl3p0JaeIfCLXx/CrIjc9Lbyvnq//czkoVZyx63RhePQz28pGPZ7DNG\nOPSjS5VSfFZFXws57O2YYzlc85+ul/ksjX2/ESu3qU+LrzCTqCrD/Ud4ey2z9vI1HbjQGu5rw/J3\n6CoJ8e9D/c9osMPSWEReKCKfFZGHReRLIvL28euHReR+EXlk/PcV9U4ZBEFwaVFHGg8B/Iqq3gjg\nlQDeJiLfA+BOAA+o6ksAPDDeD4IguOy46P21qj4J4Mnx9lkReRjAtQDeAODHxofdDeBzAN6xK1Ze\nIrQPHSp2riy2s/22akG+UNzacxYIYJ8Os4zNXZh/1uUxawfvlzJGSKnwg7qy/KVk9kTtS07GUCf3\njBwuVWTgAyu2p+1XYJ6ulqRxi7bt+zJzXLEnTqLz08nctRwdtEgmcg8NJ/d6FRkoo/3iWFOcwV38\nYZvbbVYXU/BPg6vw2Tp138fnWs1s9ghL3tLTYJK5G3R9vDRep31fPJefKK8PKLLCZ4/wd+h70Bwa\n9dcxBT8SbOphiYgcA3ATgAdp2sqKAAASvUlEQVQBXD1eJC8slldVvzMIguDSpfZCKCL7APwRgF9S\n1ec28b47ROSkiJxcXl7eio1BEAS7Sq2FUES6GC2Cf6iqnxi//JSIHB2PHwVwatp7VfUuVT2uqseX\nlpZ2wuYgCIId5aI+QhnFJnwIwMOq+l4a+iSA2wC8e/z3vbti4R5y4vDP2xeuLQqiZvsLZxyHxAA2\nLCb3PjzOGDHZHS7DgP2Aboz9faXKMRX/tblICeNTKbWtNRVnaMCHyPD7fGZJxfyl6Bl289Qspqnu\nZFnCEM4YYZt8uAX7n9T5Ozn0pUXHDcT5fyn0JW+7jB/ynbGP0Bd3ZZ9b3/kZN6g5FFet8SE4HOKT\nCuPx2MZR5OtzfkAuapuqKsO+Pt+7mH2y3kc4pCKrmakw46vP8Pfp/caba95UJxjpZgA/C+DvROSL\n49fehdEC+HERuR3AtwD89KbOHARBcIlQ56nxn6O6Y+KtO2tOEATB7GlMZonJCpkvEsNl3mZ+SK8I\nFcivtwUTsn3FGBdC8AUTkmExXDAhkd2Rd6qzNhKtaW0US1UIi98vqQiZOpionzlF8taUJsbeegUZ\ncje1sBwupRhMn4+LrQK2GKv6GBw+F38uN4cpAuAuiOnV2y7s7amVv0MKBdrwGRf0g2F53U5IeU8q\nM4Yl74Ypqlrdb8TD15jtH7g5+FpluXcjkJvCZI9U/z78t77ZTjCRaxwEQeOJhTAIgsYTC2EQBI3n\neesj9MVSOy9+4WRbFwsfYb5g04eyPoe+eN8fPco3vj4f+pIYq0h7K/nfUk6OhH/PhMnw2Gb8e5uL\nPNgeCR8h+xk14Qc0zaHazvghH1c5BVpU6cX3NW5V9Nn1PquceyM7nxj78dg3N3QVZjicpmRHRchM\nqg9zqfBrIrSGfXMctuJ9iSn4fXVDZPx1ZNhd60NijI/Q/3Y2eYsXd4RBEDSeWAiDIGg8l7U0fnXL\nxnC3FooqMO0bbJ/gfF8hh7NFCoPxoS+96fIXsJKXw1tKmR80Ze5DXXbivx5SCF4V6fTIl/JxlTvV\nlOZIzM8y1xRm9cVPqA9HSa9WfJbSYTzm5s9ZWyXkWIts9H2Nq+RZ5uJ4OAyk7XpAc48Ulqs+AyWj\nD5DsT1zRJ9lTlsbV2Ri21/D0bY+fo0pep0JffHhOVQZQSkJv150Td4RBEDSeWAiDIGg8sRAGQdB4\nLmsfIfsEAaB1dVHmKztoxzIKkxkuUIhM16UPJcNiaJt8O6UUOJ6ybq7PJlLgjOss9b7N5hlNPXe1\nD89M7xoImwouNFZyZ+nUzfILqfCf6iFD3mI7XCgG+fv4OMD6BVumMk8qnc+O+bTA6veR70yqP3TK\nL7hVqnyBqcouKf8h49Pokr5KOjZVYSZnH2TNykVVxB1hEASNJxbCIAgaz2UtjaVve6rmlDGSzdtC\nkCyHuYcwh8EAri9wIiwmKX9TlV0qQkJSc5QLp1bPLxXnLoWcpHr4SMV2ykj3S2oNWdKwnHTXlD5L\nySY61shQXx0mIVfN/DzmjuMxL5s5nEaMp8CHz5BUa/k4HqrsQhLP9y7eSITFbLbYKFCuTLMVSe3l\nb2oOHkmFz/BYltlrkJErRRNz5Hxctr17urgjDIKg8cRCGARB47nspPFr+m+dbLeuf6EZy/dVZ4zk\nnDHCRRFc9siQpHGpAOr2HkxtCpZnvp9wlfwFXLZHSv4mskL4babfiH86zk/3vPbmY/mJnptDWtWy\n2XzslCtC+DurfsRuxvwTXx5rux4g5uRUwNUVIzBPlN0YS3abnWLNqCrw4EkVf+Axn+GyE9I4BRdT\nGCbkb27GvOSlqA4jf50dvD9093Q2aOSixB1hEASNJxbCIAgaTyyEQRA0nsvOR9g6sG+ynS/axkvZ\nXPFx8r57JF+RMTJ0fYd9n2Am5ZurflPNw1IZIj4khHxupek5ZIZDR7zrLFH0pVUxR+kz01jL+5Go\nQGzeYX+WM5ezPdwUPGWin5ILkfEhT3Ru0/7YHZfzcc7nlpNfkKvFOF8ih9N4/x5NYcJxWv67NTtb\nc0qzz9D7BDtUPDZVEcZkfiTOlfvGThVZIb4nsSb6FWtVWEzCRyh+bJPEHWEQBI0nFsIgCBrPZSeN\nTd/hOWt+VYgM4Aqpkkz2/YRZa5bCT6o0Qip7xB9aNZbQHyU7TCVSPza9wEG57zBtu/k5118rtktj\nieyXFifOd7TyOH/djM2mOmr1caX2GmQXyycfZqMJt4fJbkhI6BbP6bNCzHFUwLXti8DW87kki5QS\nuZfo9L6yfJ/eT9jL31Sx1Lyi+G2qOIi/jsqhMDzm5S9fOttietNuhYveEYrIh0XklIg8RK8dFpH7\nReSR8d9XbOqsQRAElxB1pPEfADjhXrsTwAOq+hIAD4z3gyAILksuuhCq6v8B8Ix7+Q0A7h5v3w3g\njTtsVxAEwczYqo/walV9EgBU9UkRuWoHbUrTK+JbtOMe3ZMf0DfFMf2Eq6rIAOnKMRXum2RjJH9s\nYqyKugVsRmPTfVil8JPKCjNuUvaxpXyEPv2uoqFSyR/Ec5bmp8+SsjcRJmT8VMYvWP1NlPxZHGpj\nfITex5bIW+QpOKzEpZ5JRaWbpI0pf5ivPkPX1PsnmVQKnPGZervYH8zfdcnvmvD98T5/Nu/L5vCZ\nS70wq4jcISInReTk8vLybp8uCIJg02x1IXxKRI4CwPjvU1UHqupdqnpcVY8vLS1VHRYEQbBnbFUa\nfxLAbQDePf773h2z6GJQIUt18ldNH5HEWF0lkQppqTm2E50lSnKPpaA71oq/6go2Zn7/QpWCTLkK\nElknJhzHf5ZE6IuYwqyJ41iu+opB9kg+c2LMUVXgtpQtkahuY9RedSFZkxmTUntb9L8YRTr0566Q\nyj4Dhbe9bObwmYQ0Nob40+YVkjcljV34zODMqGizDutJ5jrhMx8B8H8BvFREHheR2zFaAF8tIo8A\nePV4PwiC4LLkoneEqvqWiqFbd9iWIAiCPeGyyyzRdrUuNPIplaVgBuqfu1IOp2RyzUNL8jfxlDSZ\nMbIV6roAEh+m1FeFn8wnJHSq+IMtkkDvKT1dTkkwLrjKr1dnKXiZZeR26il3oqiDuZBcSLb0O+Xj\n7NiWQg6SxV3dC+R2koQdRvKW3AMVUrZ0vRMFEyrel8r0iqILQRAE2yQWwiAIGk8shEEQNJ7Lzkd4\n39d+c7L92pf/BzM2nF+cbItteWxdNKnwlkRyQN3jkuE5qYZKfNz0U118sGLM9+rVVJhGTV+oVO64\nY1P+PZl+nD+2KlMFcH7GhP/QZLj46jOpoq3kM0yG8fC5fBMpHkv5AZH4XuqScuwm3KlVIT6bsoMn\nrcoQgfPpJcJijO/WH5fyH26SuCMMgqDxxEIYBEHjueyksWFo74dbFCkveUKv1pW/nprHbqm3SeJU\nJdWZUi1VEsEf6K9PFVvNbuACoFWZGW7+UghORaiKl9eV73H7Rhr7jAhTiMNnXPBxidAXlsM+rITm\nTxVTSJIqPGGOq3cCTRSG2LqRfAKaLhEi48ekZviMJEKeNkvcEQZB0HhiIQyCoPHEQhgEQeO5rH2E\nMhjaffYROp+B9RkmHt1vpfrMTpSYKZ2gen72jXg/VaVrJ1HMNMlWPxtXhMF0fyGQDn2p8uX68Jl0\nrFHFmLeD5/fhM2xXVcoeYL8MN2Yq5NStMOOpeWzdtMvabsCd8BNvIj3OHJuqdZvxdqTYBUEQbItY\nCIMgaDyXtzReWzf77fXiXllyq02kQi+Uwme2GlqTmrMGJYlbN6ul7pypuqN1C9VuBlZ/yrLQn7ta\nNpv/pisyVfy5UnaY95WKgSbm59AX7snhs3U4BMdNocYHwPb5jB8e85OgFrW/s5oaOvVvJPm+VGHW\n1L+zqqpGJRdRdWHWzRJ3hEEQNJ5YCIMgaDyXtTT+9Dffa/ZPfN+vTbZbB7v+8Ampogv1ZUXN92wx\ngSPZkqKujTWfyKYKvyb7ktQ7dbXEBWxvj0QBAtOnJZWBUirMWsMmWJVYymAwT5Sn21c6t8vckYpi\nDeoe3drjErI5Qe2IBl8IYQvfeyqzKWlH6mlwlaRORE+ENA6CINgmsRAGQdB4YiEMgqDxXNY+whKU\nadJyPVurfA0lH0fdbJKa76lNItPBF1VN+fe0qhdw2x+XMMWEvkzfLtnhYbvYf7OJyjFV2UBJf2fC\nDuMSS1zv1BxJn1XC98efO6cqNaVCtSZ0xGW4VMy3qfCnRAFTM03quC34/kqZH1vI2kplp7TCRxgE\nQbA9YiEMgqDxPL+kMesKH75At9WtjDIdEpnnJUlqzpWwo27+91bkmJ/C2V/VU6MkJ2v2PUmGbKSq\nx1bNl+rV69UT9//lrI3EHCUXgC+MMBlwU6TkXoWNpV41HM5R+izFdrtmgYDS56wIrdmUNK7ye2zC\nEPO2VKhR6lQ1w2JSr+9k4ZNt3RGKyAkR+aqIfF1E7tyeKUEQBHvDlhdCEWkD+K8AXgfgewC8RUS+\nZ6cMC4IgmBXbuSN8BYCvq+o3VHUDwEcBvGFnzAqCIJgd2/ERXgvgMdp/HMAPbc+c7SFrG5Ptzqp9\nnt5ZL5xFnWe32QR1h1g9krj8dVPzMncgfWx2j6Wq2wzmd7ZRT2qs5N5jO/Y5f+cWfK2tDTdWNUfC\n3o2DNefYZYaLLsVuR+zYvDMtn/dVjmdy2vIUnZ2IU5vOdu4Ip12OkqUicoeInBSRk8vLy9s4XRAE\nwe6wnYXwcQAvpP3rADzhD1LVu1T1uKoeX1pa2sbpgiAIdgfxxSVrv1GkA+BrAG4F8A8A/grAv1DV\nLyXeswzg/wE4AuDbWzrxznEp2ACEHZ6wwxJ2WDZrx4tV9aJ3YFv2EarqUER+AcBnMHJHfTi1CI7f\nswQAInJSVY9v9dw7waVgQ9gRdoQdl4Yd2wqoVtU/AfAnO2RLEATBnhApdkEQNJ69Wgjv2qPzMpeC\nDUDY4Qk7LGGHZVfs2PLDkiAIgucLIY2DIGg8M10I96pIg4h8WEROichD9NphEblfRB4Z/33FDOx4\noYh8VkQeFpEvicjb98IWEZkTkb8Ukb8Z2/Hr49evF5EHx3Z8TER6u2kH2dMWkS+IyKf2yg4ReVRE\n/k5EvigiJ8ev7cVv5JCI3CMiXxn/Tn54D34fLx1fhwt/nhORX9qj6/Fvx7/Rh0TkI+Pf7o7/Pma2\nEO5xkYY/AHDCvXYngAdU9SUAHhjv7zZDAL+iqjcCeCWAt42vwaxtWQdwi6r+AICXATghIq8E8B4A\n7xvb8SyA23fZjgu8HcDDtL9Xdvy4qr6MwjP24jfyOwDuU9XvBvADGF2Xmdqhql8dX4eXAfhBAOcB\n/PGs7RCRawH8GwDHVfX7MArTezN24/ehqjP5A+CHAXyG9t8J4J0zPP8xAA/R/lcBHB1vHwXw1VnZ\nQjbcC+DVe2kLgAUAf41Rnvi3AXSmfV+7eP7rMPpHdQuAT2GUurkXdjwK4Ih7babfC4ADAL6Jse9+\nr+xw534NgL/Yo+txoZ7BYYxC/T4F4LW78fuYpTSeVqTh2hme33O1qj4JAOO/r5rlyUXkGICbADy4\nF7aM5egXAZwCcD+AvwdwWlUvNH6Z1ffzfgC/iqLE55V7ZIcC+FMR+byI3DF+bdbfyw0AlgH8/thV\n8EERWdwDO5g3A/jIeHumdqjqPwD4LQDfAvAkgDMAPo9d+H3MciGsVaShCYjIPgB/BOCXVPW5vbBB\nVTMdSZ/rMCqpduO0w3bTBhF5PYBTqvp5fnnWdoy5WVVfjpHr5m0i8qoZnNPTAfByAB9Q1ZsArGA2\ncnwqY9/bTwH4n3t0/iswKu13PYBrACxi9P14tv37mOVCWKtIwwx5SkSOAsD471OzOKmIdDFaBP9Q\nVT+xl7YAgKqeBvA5jHyWh8Y55MBsvp+bAfyUiDyKUT3LWzC6Q5y1HVDVJ8Z/n8LIH/YKzP57eRzA\n46r64Hj/HowWxr36fbwOwF+r6lPj/Vnb8RMAvqmqy6o6APAJAD+CXfh9zHIh/CsALxk/8elhdMv9\nyRme3/NJALeNt2/DyF+3q4iIAPgQgIdV9b17ZYuILInIofH2PEY/uIcBfBbAm2Zlh6q+U1WvU9Vj\nGP0e/kxV3zprO0RkUUT2X9jGyC/2EGb8vajqPwJ4TEReOn7pVgBfnrUdxFtQyGLsgR3fAvBKEVkY\n/9u5cD12/vcxK6fr2LH5kxhVrPl7AL82w/N+BCMfwwCj/3Vvx8gX9QCAR8Z/H56BHT+K0W383wL4\n4vjPT87aFgD/BMAXxnY8BODfj1+/AcBfAvg6RnKoP8Pv6McAfGov7Bif72/Gf7504be5R7+RlwE4\nOf5u/heAK/bIjgUATwM4SK/thR2/DuAr49/pfwfQ343fR2SWBEHQeCKzJAiCxhMLYRAEjScWwiAI\nGk8shEEQNJ5YCIMgaDyxEAZB0HhiIQyCoPHEQhgEQeP5/zIV7kA1OaEIAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<matplotlib.figure.Figure at 0x7fe5881c4c18>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.imshow(data, origin=\"lower\")"
   ]
  }
 ],
 "metadata": {
  "_draft": {
   "nbviewer_url": "https://gist.github.com/6810bb6e29a9118f9569abfc31e67637"
  },
  "gist": {
   "data": {
    "description": "hsds_examples/nrel/notebooks/Sandy_example.ipynb",
    "public": true
   },
   "id": "6810bb6e29a9118f9569abfc31e67637"
  },
  "kernelspec": {
   "display_name": "Python [conda env:hsds]",
   "language": "python",
   "name": "conda-env-hsds-py"
  },
  "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.6.3"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 1
}