joezuntz/h5py-example.ipynb

## h5py-example.ipynb
{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "0cca298f-dc13-4dd2-8bae-93e9eb739ca3",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Populating the interactive namespace from numpy and matplotlib\n"
     ]
    }
   ],
   "source": [
    "import h5py\n",
    "%pylab inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "b1809fc2-56cb-498e-bb61-efd30c54ac61",
   "metadata": {},
   "outputs": [],
   "source": [
    "filename = \"/global/projecta/projectdirs/lsst/groups/WL/users/zuntz/data/cluster_mag/cluster_mag_background.hdf5\"\n",
    "f = h5py.File(filename, \"r\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "d0a2bc67-2f14-4d72-b667-1d9d16e822bc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<KeysViewHDF5 ['provenance', 'sample']>"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "# list of data sets and groups in this file:\n",
    "f.keys()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "a21f3e44-a34f-4504-afbd-6a081f00923e",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<KeysViewHDF5 ['dec', 'ra']>\n"
     ]
    }
   ],
   "source": [
    "# picks out a single group\n",
    "g = f[\"sample\"]\n",
    "print(g.keys())"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "7ea6a7e0-3f96-44fb-a7fa-16e31301152b",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "<HDF5 dataset \"ra\": shape (257299825,), type \"<f4\">\n"
     ]
    }
   ],
   "source": [
    "print(g[\"ra\"])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "433c4018-fcdc-47dd-befb-7b676704b66b",
   "metadata": {},
   "outputs": [],
   "source": [
    "ra_chunk = g[\"ra\"][:5000]\n",
    "dec_chunk = g[\"dec\"][:5000]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "d08f0641-9d94-4ff4-9654-8b6827fb24fc",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "(0, 4)"
      ]
     },
     "execution_count": 7,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAD4CAYAAADy46FuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2de6xnV3XfP6vYM40SGAMzZC7jMePIozquh7qei2XUpkXEocFENmr5A9QUE4W6DeJRIde8WhKX0Doiik3aRmCSYGhMoDFN4/ghpyZBUaXEZmzGM4OwawcG8GBkTwhjEJXJNKt/3HPsM3v265zf89zf9yP9dH+/c/Zj3TNz9/qtvR7b3B0hhBAixd9atABCCCGWGykKIYQQWaQohBBCZJGiEEIIkUWKQgghRJYzFi3AELZv3+579uxZtBhCCDEq7r///uPuvqNvv1Eqij179nDgwIFFiyGEEKPCzL42pJ+2noQQQmSRohBCCJFFikIIIUQWKQohhBBZpCiEEEJkkaIQQgiRRYpCCCFEllEqisPHTixaBCGEWBlGqSgA9rz7jlN+CiGEmA02xoOLtq7t9acff2TRYgghxKgws/vdfb1vv1GW8IDTLYmj179mQZIIIcTmZpRbT/t2bePo9a955gXaihJCiFmxqbae9rz7DlkWQgiRYOjW02gtihSyKIQQYrqMUlG07Hn3HUlfhRSGEEJMh1EqijaPovVRxJRCdwtKSkMIIYYzSh/F+vq6H7/suui90KKQz0IIITZYKR8FcFrUU8yCSFkbQggh6hmtokihyCchhJguo916Sp2ZnbImpDyEEKvOym09QTzJLuXgDttqS0oIIeoYtaJoFUJoLdRsP7V9pTCEECLPqBUF9Mub6FoTUhBCCFHHRIrCzD5gZofM7KCZ/ZGZvbi5/nwz+/3m3n1mdmGi/7lmdq+ZPWJmnzGzLTXzHj524pSFPrQgQksjvNf+VHKeEEKUmdSi+JC7v9TdLwJuB97fXH8vcNDdXwq8Efhwov+vADe4+17gr4Cf7zN5axnkFvxQacQyueXoFkKINBMpCnd/qvPxh4E2hOoC4HNNm4eAPWb2o92+ZmbAK4Fbm0ufAF5bM29b66nWcd1FikEIIfoxsY/CzD5oZt8A/jnPWhQPAv+0uX8J8BLg7KDrC4HvuPvJ5vNjwK7MPFeb2QEzO/Dkk0+eoiDCZLtwyym8n1Ik8l0IIcTpFBWFmd1jZkcirysB3P197r4buAV4a9PteuD5ZnYQeBvwReBkOHRkumRSh7vf5O7r7r6+Y8eOU+7FFEYsEqq9XqoPJWUhhBDPMrWEOzN7CXCHu18YXDfgq8BLu1tVzfUngZ3uftLMXg78krv/k9Jcba2nXOG/ISGzIcryFkJsJoYm3E2kKMxsr7s/0rx/G/CP3f11ZnYW8H13/4GZ/UvgJ9z9jZH+vwd81t0/bWYfAQ65+2+U5s2dmV2q85Rb+KUYhBCbmUVlZl/fbEMdAl4FvKO5/uPAl8zsIeDVneuY2Z1tGC3wLuCdZvYoGz6L3+ozeRgiGyqJmI8i1z+lXLQVJYRYZUZZ66m1KEILIGdN5MqPD83uFkKIMbGQradFsXVtr69ddWNysY8VA8wplRyxflIgQogxsnJFAbu1mvr4IvpWlO1aJ9qCEkKsIqNUFG3CHZyeK9H1NeS2k/pYBd3xFEIrhFg1RqkouoRbTjElkdsu6rPg63wLIcQqMmpFUVOOI6Y8Ugt+TmnE5pFVIYRYBUatKCC/WKfuxc6iCMNjh/o9hBBiszHKqKf2KNTQMkiFvnbv11STzVkcOb+HIqKEEMvMSoXHhmdm53ImSpnY3fal6937Ug5CiLGxcuGxLd1opDAqKbQaUltEqbO3S/MOuSeEEGNjU1kUNQl14fVJLYLQtyHrQgixrKyURXH42AkgrSBSWdhdahb2nAXSnSfMsRBCiM3EKBVFN+EOymdmd9u190OndZjl3f0Ze587+yKcT1tRQogxM0pF0VoUsSzr7rU+VkOsX6xvKloq5+eQpSGEGDOjVBT7dm07zRro/myJLdAxBVI6wyIcL+cLid2TRSGEGDNnLFqASQnrN6VKerT3U31jn1PXWnKly2PyCCHEGBmlojh87ARrxBfxvsX+WmJKpUZJlOZURJQQYuyMcusJ6stp5JzPuTLlbbuYgzscP7YNFpM15hQXQohlZ5QWReujCMtyQHr7pyWnFML3qQzsXCnzknXRlUnVaIUQo8DdR/fav3+/53jJu25/5mf3fexnrF+qfaxPrG9t+5wsQggxbYADPmDNHe3WU0t3+ygXgVSKaArJhd628+ZkCuXL3Zc1IYRYZkZbwuP4Zdf1KvudCostFf6rXchzEVN9SowIIcSsWLkSHrULbLe0RpiRnXJqhwt8jSUSc1jX+Ctic8jRLYRYJkbrzE4ROohzdZ9iNZpSDuZaCyCWy9HtFyqhcC5tRQkhlo2JLAoz+4CZHTKzg2b2R2b24ub6883s95t795nZhYn+N5vZV5v+B83soknk6RKLiGqvQ3qrqNumdC0kl/EdztkqqFRJECGEWBYm8lGY2fPc/anm/duBC9z9X5vZh4Dvuft1ZnY+8F/d/Scj/W8Gbnf3W/vM2z3hLpdJ3Sf8tKb8R6n/EKugpq2sDCHENFiIj6JVEg0/DLRa5wLgc02bh4A9Zvajk8wVEi7MqSS23LZTrF24FVUjR9uvlSFlqYTXahWAlIQQYqEMiantvoAPAt8AjgA7mmv/Efi15v0lwElgf6TvzcDDwCHgBmBrZp6rgQPAgXPOOcfdy7kIYV5ELLdhFjkMsXyMrgxD51a+hRBiEphVHoWZ3WNmRyKvKxtF8z533w3cAry16XY98HwzOwi8DfhioyxC3gOcD7wMeAHwroxCu8nd1919fceOHcCpUUy1W1AhNd/+Y6QsmJilktoCk6UghBgDxagnd7+scqxPAXcAv+gbW1I/B2BmBny1eYVjP968fdrMPg5cUzkXUHcGdu5MiFr/QKp/n0qzbftp+BtU+kMIMU8mjXra2/l4BfBQc/0sM9vSXH8z8Kd+qj+j7b/W/DTgtWxsX/Wia1WEC2gpGztnTXQVRCmxrzYKKjZnH2IRU4qWEkLMmkkT7q5vtqEOAa8C3tFc/3HgS2b2EPDqznXM7M42jBa4xcwOA4eB7cAv95k8V5CvtkLrkGinMP+i1ipJlfKIKahY29g9WRRCiFkzyhIeW9f2+tpVNwJlH0PN9tCswk9L45aS8UrlRYQQog8rVcJj365t2WzqcGsmVZqj2ydGn22dkgWTshRSc8d+vz4htUIIMTWGhEot+rVl53nung8XrQ1FrSk/Pim5kuO50ua1Y6T6CCFEFwaGxy580R/y2rLzvGReREguv2JI3z7t++ZM5NrXnqshhSGESLFSiiJ3cFHtQjlrC2IaY/dVAjoISQiRY6iiGKWPIkcqX6LkIxgSZpqKYqo94Kg0Vi6qK0bof1HorBBiKgzRLot+9TkKddYM2V4qbRv12cLqbsHJkhBC5GAVLYrUN+bUmRCz/oYdS/yLtUllaKcS6GLj5hICZUkIIabJqBVF7hS61DZMaRHtVqKtCXeNyZTbeip9bseokbUrb04BCSHERAwxQxb9Km095RiyLVUTWlsb/lqqcjtEpnDraZJILyHE5oWBW0+jzMzuHlzUUsqArj1UaJoHCZXkG3rYUTh2aY6wjawOIVaToZnZC7cOhrxqnNm5b/g1365rzrfIzRebozbpryR3zb0aC0cIsVqgPIrh2culhbr2cKQ+c/ahpGxySqkUQSXFIcTqMFRRjNqZHZKqj9T+LNVVqslTiDFplFGpYmwu6qmVK1YPKrwW/g59j30VQqwoQ7TLol85iyL2eUiZixpHdThWTY5ETq7UGNPavuqOL+e2EKsHq7j1FFvMJ63PFI4ba1Prn+gz35BtoEnaS1kIsXoMVRTFo1CXne6WSi6KKLX1lDqbIrVllSsBktrCSs1dOnei1CYmV+x3zm1ZCSFEiVGHx3YZEvJZGwoLpyqi2gOGwjHCcXIyDAnBjckZklOMfX4nIcT4WNnw2FlE7pS2k8L7tVFJ3b59/SQp+Wp9ErExtfUkxGrBKvoohtAne7kUdprLtO6rDFKfa/sOCd2NtZ/GOEKI5WSoohjl1tPWtb3+9OOPPPM55ado7+W2doZmKQ/Zoin5N4ac550Kky3JFso/5BkKIcbF0K2nUSqKsIRHyambos8CmFuAh5biiC3WNYoit8jH2tScwTG0lIgQYjyspKLokloMu/f7Lsix8Wu+tQ+tF5VScjGHc8pZ3SUma66Mefd9TgkJIcbJSjqzQ4ds7R57TTJdzoEcu993/JoxJk0MzMk7iW9ECDFOWLSPwsyuAT4E7HD342ZmwIeBy4HvA29y9wci/fYDNwM/BNwJvMMLQq2vr/vxy65L3s99K85920/lTMS+zXfbl6yNWn/AMmz75CyIZZBPCDGchVoUwG7gbuBrwPbm2uXAXYABlwL3JvreB7y8aXcX8OrSfLGop1yJjD6hrLFxY5bLUCaJKipZUENk7BO1NYtQZCHE/GCRFoWZ3Qp8APgDYN03LIqPAp93999t2jwMvMLdH+/0WwP+xN3Pbz6/oWnzr3LzdX0UQ85xGBrlNOlZDimrBPLf4HOJczk/Rel8ilKRQ0U+CbG5WJhFAVwBfLh5f5RnLYrbgX/Yafc5NpRIt+86cE/n808AtyfmuRo4ABzYsvO8uX2zndY8Q/wV3W/2uTpNpdpUJbn6IqtCiHHCLBPugHuAI5HXlcC9wDY/XVHcEVEU+4NxXxZRFH9Ykqe09ZS6V9qKKvUPx+g7Tkmumvlr24VKpo/TvmZuZXYLMT5mqiiSnWEf8ESjII4CJ4GvAzuBjwJv6LR9GFgL+q8BD3U+vwH4aGneGkVREwHU95v3LPwUuQW9r4KYJAKsZo5JFaMQYrEsRFGcNtipFsVrONWZfV+izxea+60z+/LSPJMchRprW7o2z0Uxt6WUu969l9qOmoYFIAUhxHgZqihmecLdncBXgEeBjwFvaW+Y2cFOu18AfrNp9xeNspiIVFnw8JS4tm14PecELiW55e6388ROnusjQ4pumfPQaR46v2PPInWt+zP2vGLthBCbh02RmZ1aYEvZ2jFqM5m77ft8TsnWna+Uw5DK3yideVEqJT5JTocyuIVYfoZGPY3yzOzDx05Ev+l2P4d0F8jab+btYtm+uqQWxtTn0HpILaip36PbvhTW2v0cWhfhta5c4Rw11kG3b03IrhBihAzZr1r0a//+/dn99ppInbBdzBGca1uKsir5EUqy5saOyRz7HVLtc3MPjYLq3pcfQ4jlhEWX8Jgnsa2n1JZJTYJZjNK2Tk0pkD5bObXbSLlv7bEtr9p5u2Pk5ql9DtqKEmL5WMmigCGTfhtOtc19Oy/lJ9TkMoT9+n7T7zN2qk9q/JqxS1ZTrI0QYv6wDOGx83qFimKS7ZtU/z4KIHW97wJbK3tJSdQu1DkZw9+tZp5JnqsQYvastKJIUfIV1PRLtc8tgkMtlxrZSv1KSmqof6Tkt4mNnfKhCCEWw0opir61nmILaM02Uo2i6bM4D7F8aiyI1NhDHc99LKGh1pyUhhDzZ6iiGGV47L5d26pOkEsl2PUhlbDWvs+FzsbGiMmUkj8WmhubL5QvVf01J2MpUTAcu2aOGOHzFEIsP6OPeipF1+QURSwiqibRLIxQivWviVyqJfU71CT3lZIRaxPwcmPF5pMiEGL5WKmEu5ZYGGfXksglzEE8SaxPYl5uUeyrJGJzdn+H7pjh9Zjsqd8zNW94LczizimJMJEwFVqrRDwhRsqQ/apFv7bsPO+UfbdJk7xK/Ut+gpKTNzZeH9ly70uRSTVjhG37RDhN+tyFEPMDObNPpSbsM3YttwBOK9IpJ9Mk4a65+7WKpI9ctXOH16QghFgMK6coWlLfcvssrkMWrqGLXZ8Fu2TlTEOWof1TUWR9FLQQYr6slKLIHVxU+tY/zQWrxsKYVAnlFELfLbBSiG+NPDGZan9XKQohFstKKYrU1lOtr2HSBSu2wNYukn39E7PKO0hZBH3H6I6Vui+EWA6GKorRRz11f3bLZ+eIldjuE5EzyTkNtX3D6KfYvZTMtfkM3QilIXKV+ipEVojNwSgVxb5d2555H4ZxQn7h6i5ufc54CMNuUyGpuZyH2LjdsbvvS+dcpOaLyTQtcnkauTBlIcS4GWXC3da1vb521Y2nXc+VBG/vh+SSymJ9+mQfl+YukSr3HUu4y5VC75vk17dEed+S6kKIxTA04W6UimJ9fd2PX3Zd8n534cottiklEetfqyCGZGTnFuY+i/20lNLQMZWVLcRys1LnUYTO7FyI7LTI5TdMEqUUjjM0xyPsM0/ktBZiHLBKUU+l8NhJwmFzIaS59qXFvObzLJLS5pljIYUhxHKzkooiF5Y5acJd7tv+LKyWWVsLs1zEpSCEGAdDFcVofRQHDhyoinQKyVVcTRXay43VpdbxXCtbbr4+ck5CjUzySQgxDhZaPdbMrjEzN7PtzWczs183s0fN7JCZXZzo93kze9jMDjavF/WZN+aIrnE6x5za4Xipiqs1FWXDCqopp3lM0ZSqucbkzMk1JE8klGmS+0KITcAQM6T7AnYDdwNfA7Y31y4H7gIMuBS4N9H388B63zm7PopZlJ6oye6uuVfjiO4jf+l+n601IcTqwcCtpzOmoGtuAK4F/qBz7Urgk41gf25mZ5nZmrs/PoX5OHzsRDZfAtIHFJW2akrbT7nQ1/D8hppv+7lkwdK2VSwpr1ZmIYSoZaKtJzO7Ajjm7g8Gt3YB3+h8fqy5FuPjzbbTvzczq5m3exRqbKEtLaqxBbzvgT/dfuGCXsqyTmV1xw4K6usDqJFZCCH6UHRmm9k9wM7IrfcB7wVe5e4nzOwoG9tIx83sDuA/ufv/bsb4HHCtu98fjL3L3Y+Z2XOBzwK/4+6fTMhxNXA1wDnnnLPf3vAbWblrMpbDz6WEvBKpTOXYvH2znmO/W6rftDOk5bAWYnMwM2e2u1/m7heGL+ArwLnAg42SOBt4wMx2smFB7O4MczbwzcjYx5qf3wU+BVySkeMmd1939/Vv/WDLM9fbb9B9tlpCR3PsXmycWP2iXCRSrv2QKK2ccz1WcylmofRFSkIIMXjryd0Pu/uL3H2Pu+9hQzlc7O7fAm4D3thEP10KnAj9E2Z2RidK6kzgZ4AjNXO3RQFzUUItoeUQuxdaAClFktrSim0T5RRYd9yU/H23kHK+jpTsNUqktnyJEGLzMqvqsXeyYXE8CnwMeEt7w8wONm+3Aneb2SHgIHCsaVvk8LET2X3/0oJcu/8fsxBCUoqnNGbN/LF2uf4l30zK+kg5/nPjCiFWiCGhUot+dY9CDemTkZ2636feUm6+0mE+pTFK4b21v0+JmLwKqRVi88EqlfBoiwLmFu2akhg1CqH2fqzMR62CGqLc+rbrM4aUhBCbk5VSFLmigO371Df3VJ9Y25qCgCWrIdVnEuVQomYcKQMhVo+VVBQ1lVuHVj+tLQA4SWHAaY7ZZ/sr1W+aRQ6FEMvHUEUxyqNQ+9A3cgjyDt1uyGspO7x2zpiTuW+0USrKqSbnIhYIIIQQLaNWFOFiWCq813cRD/vGwl1zY6byKFK5FEMiompIyVijTIQQYuHbSENekxQFzG315PwH4XhDI6v6FO4bGvVUM7YQYvVglbae2qKAkC/uF1Kq8dTHMqnZGopZD7migil5Y59lBQgh5sUoFUWbmV2iVL21bRNbgGMZ16mFOxy/z55/bc2nSchlZQshRIlRn3BXS8xpW2obXuv2Kzmxp1GUTzWWhBDTZqEn3C0DJauh+zPXP2YlhGVBakt0pL7JlxzsuTlSv6csBCHErNiUFkXMcqi1Jrpthpb3jl1PlQ2v3ZoKx6/ZspJFIoToslIWxeFjJ4D82dMhNYtmrvx4bd9c/5oS5rE+sUU/53+RkhBCTJUhoVKLfoWZ2eH79nMs1LVvDai+BQVr7tWUGymNo5BXIURfWKUSHm1RQPfTi/G112LkciRS/Uq5GH1KcKT6TVICRAghahmqKEa59dSemZ0KQ+0TDtrHX1A7fup8h1w0VU1ORo0cQggxbUbrzD5+2XXRhTfm8O1erwl/DUn5QMJxwzGGONBz1Pge5J8QQqQY6swepaLYurbX1666MakUIP7tvU9+Qy6KqR0/tBZKTvWcYkkpMC36QohpMVRRLNzfMOQVq/XU58Ci1P2Yk7nmbImacWpkGOrvEEKIGhjooxitRfH044+ccq12C2fSLO3YGDXzlshZFLOYTwixeqysRZFjWmGuNf2GhMl27w8JkRVCiD6wSlFPMVLJdzFfQOx8iFSkUmrcltozHbrzxnwZpaitmgKHQggxC0avKGqc111yUUq5OVKZ1yVF0reSbOx9TH4hhJgXZyxagFkRftPvLvalMyFK0VF9Iqb6KCadNyGEWEZG6cwOiwL2cWS37WrCX7vXSsqllj59FB4rhJgmCy0KaGbXmJmb2fbm8/lm9mdm9rSZXZPpd66Z3Wtmj5jZZ8xsS+2cqTLg4f1Uu1Q2dNumxhfQZxGvKTTYp1qtEELMi4kVhZntBn4K+Hrn8reBtwO/Wuj+K8AN7r4X+Cvg52vmPHzsRLFaa8o5nSMVfhpTHLnPrTzdV06R5RzdQgixaKZhUdwAXAs8s4fl7k+4+xeAv051MjMDXgnc2lz6BPDa2klzC2p30a3ZXoopmvBaLgs8FcUU+iliFk2sjRBCLBMTObPN7ArgmLs/uLHu9+KFwHfc/WTz+TFgV2auq4GrAc455xwgX18pJOdgbq91lUPMCT1kayjnGM8piT5RWUIIMVNKiRbAPcCRyOtK4F5gW9PuKLA96PtLwDWJcXcAj3Y+7wYO1yR/xM6jCCmVBw+v9ymfkZtraGJfdy4l1wkhZgHzPo8C2Ac80SiIo8BJNvwUOzttcorCgOPAGc3nlwN318y9Zed5z/zipQW+9pCgUGHUKpo+lA5aEkKIWTJ3RXHaQD0tiub+7wGvb95/BHhLzVypEh6TLLY1hfxKbWpO0pNCEEIsiqGKYmp5FGZ2FFh39+NmthM4ADwP+Bvge8AF7v6Umd0JvNndv2lmPwZ8GngB8EXgZ9396dJcYVHAWAnx3BkS4bUUqbyLXBRTOK7KhwshloWVP48iJHZoUXi/pc+5EalxUs7vWoUihBCzRtVjK0j5LWJtasZp34fnd/eRRQgh5gWrVD328LETyXu5/IowhyFmOXTHyeVJlOaKJeB138uaEEKMhU1XFDC11ZPKYUjlScQys/sWCMzlTgghxFgYpY8iLAoYMq1v7LUlzGOJe11qCgwKIcSsGeqj2DQWRc25DznHc4zc9Zj1UVIYuTGFEGJZGaVFkTszG/IRSkPKgufGHDquEELMm4WWGZ83+3ZtA06tvBpaErXf5mPtciXMa06h695TNVghxNgZpaJo6S7iYTRRLpopHCMs9Z0rYZ6LeCpVrRVCiFEyJKZ20a8tO8+ryoVIldnIFd/rk0ehXAghxJhgYB7FKJ3Z+3Zt43jnc9dxHAtHLYXIwrDDjWLXZUkIITYbo3Rmx87MjtVj6pLLkxiyyE96ToUQQsyblav1FEY91TLLRV2+CSHEMrNSUU8wvWiiSSKXcselCiHEZmGUimLfrm3VIavhYl8bItu2zVkJbdSVQmCFEJuZUW49lUp4QPyMivB9+xnKCXWlNkIIseys3NZTl9gWUc0hQ91EvTCXottOloMQYpXZFIoiJOc3CKOicolzuUgpIYRYFUarKGKZ1OHPFLF2sWu15TqEEGIzM0ofRTc8tnQmdZhjURseO8SXIYQQy8xK+SjaooApYpZC7H7uPOvSuEIIsSqMUlF0qYlWivkZ+i76UhJCiFVl9Ipi2sgXIYQQpzJKH0VNHkUXldYQQogF+yjM7BozczPb3nw+38z+zMyeNrNrMv1uNrOvmtnB5nVR7Zx9v/nLUhBCiGFMrCjMbDfwU8DXO5e/Dbwd+NWKIf6tu1/UvA4OkSH0RYRhrUPKiQshhNhgGhbFDcC1wDN7WO7+hLt/AfjrKYwfJXbORKzUeO2WkxSIEELEmUhRmNkVwDF3f3CCYT5oZofM7AYz25qZ62ozO2BmB5588kmg35kSuVyKXD8hhFh1iorCzO4xsyOR15XA+4D3TzD/e4DzgZcBLwDelWro7je5+7q7r+/YseO0+7nEuVSSnawIIYQoU1QU7n6Zu18YvoCvAOcCD5rZUeBs4AEz21k7ubs/3hzl+jTwceCSmn6Hj52oWuRTW1KgSCghhKhl8NaTux929xe5+x533wM8Blzs7t+qHcPM1pqfBrwWODJElpwzO3WmhKrBCiFEHTNJuDOznWb2GPBO4N+Z2WNm9rzm3p1m9uKm6S1mdhg4DGwHfrlm/O7BRaVT6CBdqkMWhRBClNl0CXdSAkIIEWeligJ2iTmqpSSEEGJ6jF5RKIJJCCFmy2gVRe25EkIIISZjtIqihCwMIYSYDqNVFDELohsBJQtDCCGmw2gVBZweGjv0UCIhhBBpRq0oWqQYhBBidmwKRSGEEGJ2jFpRyJIQQojZM1pFoagmIYSYD6NVFLImhBBiPoxWUQghhJgPo1UU2noSQoj5MFpFoa0nIYSYD6NVFEIIIeaDFIUQQogsUhRCCCGySFEIIYTIIkUhhBAiixSFEEKILFIUQgghskhRCCGEyCJFIYQQIou5+6Jl6I2ZfRd4eNFyVLAdOL5oISoYg5xjkBEk57SRnNPl77j7c/t2OmMWksyBh919fdFClDCzA5JzOoxBRpCc00ZyThczOzCkn7aehBBCZJGiEEIIkWWsiuKmRQtQieScHmOQESTntJGc02WQnKN0ZgshhJgfY7UohBBCzAkpCiGEEFmWWlGY2U+b2cNm9qiZvTtyf6uZfaa5f6+Z7VlCGd9kZk+a2cHm9eZ5y9jI8dtm9oSZHUncNzP79eb3OGRmF89bxkaOkpyvMLMTnef5/gXIuNvM/sTMvmxmXzKzd0TaLPx5Vsq5DM/zb5vZfWb2YCPndZE2y/C3XiPnsvy9P8fMvmhmt5qKmucAAAOHSURBVEfu9X+W7r6UL+A5wF8APwZsAR4ELgjavAX4SPP+9cBnllDGNwH/ZQme5z8CLgaOJO5fDtwFGHApcO+SyvkK4PYFP8s14OLm/XOB/xP5d1/486yUcxmepwE/0rw/E7gXuDRos9C/9R5yLsvf+zuBT8X+bYc8y2W2KC4BHnX3r7j7D4BPA1cGba4EPtG8vxX4STOzJZNxKXD3PwW+nWlyJfBJ3+DPgbPMbG0+0j1LhZwLx90fd/cHmvffBb4M7AqaLfx5Vsq5cJpn9L3m45nNK4yyWfTfeq2cC8fMzgZeA/xmoknvZ7nMimIX8I3O58c4/T/5M23c/SRwAnjhXKQL5m+IyQjwz5rth1vNbPd8ROtN7e+yDLy8Mf/vMrO/u0hBGrP977Px7bLLUj3PjJywBM+z2So5CDwB/C93Tz7PBf2tA1VywuL/3m8ErgX+JnG/97NcZkUR03Ch9q5pM0tq5v9DYI+7vxS4h2c1+bKx6GdZywPAS9z97wH/GfifixLEzH4E+Czwb9z9qfB2pMtCnmdBzqV4nu7+/9z9IuBs4BIzuzBoshTPs0LOhf69m9nPAE+4+/25ZpFr2We5zIriMaCrjc8GvplqY2ZnANuY77ZFUUZ3/0t3f7r5+DFg/5xk60vN81447v5Ua/67+53AmWa2fd5ymNmZbCy+t7j7/4g0WYrnWZJzWZ5nR57vAJ8Hfjq4tei/9VNIybkEf+//ALjCzI6ysRX+SjP7naBN72e5zIriC8BeMzvXzLaw4XS5LWhzG3BV8/51wB9746FZFhmDfekr2NgnXkZuA97YROtcCpxw98cXLVSIme1s91PN7BI2/g//5ZxlMOC3gC+7+68lmi38edbIuSTPc4eZndW8/yHgMuChoNmi/9ar5Fz037u7v8fdz3b3PWysR3/s7j8bNOv9LJe2eqy7nzSztwJ3sxFd9Nvu/iUz+w/AAXe/jY0/gv9mZo+yoRFfv4Qyvt3MrgBONjK+aZ4ytpjZ77IR4bLdzB4DfpENZxzu/hHgTjYidR4Fvg/83JLK+TrgF8zsJPB/gdfPe8Fg41vbvwAON/vVAO8FzunIuQzPs0bOZXiea8AnzOw5bCiq/+7uty/T33oPOZfi7z1k0mepEh5CCCGyLPPWkxBCiCVAikIIIUQWKQohhBBZpCiEEEJkkaIQQgiRRYpCCCFEFikKIYQQWf4/cwtzuttE+SQAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plot(ra_chunk, dec_chunk, ',')\n",
    "xlim(0, 4)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "id": "963a8b30-d767-42b7-afd9-e6909c0591af",
   "metadata": {},
   "outputs": [],
   "source": [
    "# would get the full column:\n",
    "# g[\"ra\"][:]"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.7.4"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}
	{
	"cells": [
	{
	"cell_type": "code",
	"execution_count": 1,
	"id": "0cca298f-dc13-4dd2-8bae-93e9eb739ca3",
	"metadata": {},
	"outputs": [
	{
	"name": "stdout",
	"output_type": "stream",
	"text": [
	"Populating the interactive namespace from numpy and matplotlib\n"
	]
	}
	],
	"source": [
	"import h5py\n",
	"%pylab inline"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 2,
	"id": "b1809fc2-56cb-498e-bb61-efd30c54ac61",
	"metadata": {},
	"outputs": [],
	"source": [
	"filename = \"/global/projecta/projectdirs/lsst/groups/WL/users/zuntz/data/cluster_mag/cluster_mag_background.hdf5\"\n",
	"f = h5py.File(filename, \"r\")"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 3,
	"id": "d0a2bc67-2f14-4d72-b667-1d9d16e822bc",
	"metadata": {},
	"outputs": [
	{
	"data": {
	"text/plain": [
	"<KeysViewHDF5 ['provenance', 'sample']>"
	]
	},
	"execution_count": 3,
	"metadata": {},
	"output_type": "execute_result"
	}
	],
	"source": [
	"# list of data sets and groups in this file:\n",
	"f.keys()"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 4,
	"id": "a21f3e44-a34f-4504-afbd-6a081f00923e",
	"metadata": {},
	"outputs": [
	{
	"name": "stdout",
	"output_type": "stream",
	"text": [
	"<KeysViewHDF5 ['dec', 'ra']>\n"
	]
	}
	],
	"source": [
	"# picks out a single group\n",
	"g = f[\"sample\"]\n",
	"print(g.keys())"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 5,
	"id": "7ea6a7e0-3f96-44fb-a7fa-16e31301152b",
	"metadata": {},
	"outputs": [
	{
	"name": "stdout",
	"output_type": "stream",
	"text": [
	"<HDF5 dataset \"ra\": shape (257299825,), type \"<f4\">\n"
	]
	}
	],
	"source": [
	"print(g[\"ra\"])"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 6,
	"id": "433c4018-fcdc-47dd-befb-7b676704b66b",
	"metadata": {},
	"outputs": [],
	"source": [
	"ra_chunk = g[\"ra\"][:5000]\n",
	"dec_chunk = g[\"dec\"][:5000]"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 7,
	"id": "d08f0641-9d94-4ff4-9654-8b6827fb24fc",
	"metadata": {},
	"outputs": [
	{
	"data": {
	"text/plain": [
	"(0, 4)"
	]
	},
	"execution_count": 7,
	"metadata": {},
	"output_type": "execute_result"
	},
	{
	"data": {
	"image/png": "iVBORw0KGgoAAAANSUhEUgAAAYoAAAD4CAYAAADy46FuAAAABHNCSVQICAgIfAhkiAAAAAlwSFlzAAALEgAACxIB0t1+/AAAADh0RVh0U29mdHdhcmUAbWF0cGxvdGxpYiB2ZXJzaW9uMy4xLjEsIGh0dHA6Ly9tYXRwbG90bGliLm9yZy8QZhcZAAAgAElEQVR4nO2de6xnV3XfP6vYM40SGAMzZC7jMePIozquh7qei2XUpkXEocFENmr5A9QUE4W6DeJRIde8WhKX0Doiik3aRmCSYGhMoDFN4/ghpyZBUaXEZmzGM4OwawcG8GBkTwhjEJXJNKt/3HPsM3v265zf89zf9yP9dH+/c/Zj3TNz9/qtvR7b3B0hhBAixd9atABCCCGWGykKIYQQWaQohBBCZJGiEEIIkUWKQgghRJYzFi3AELZv3+579uxZtBhCCDEq7r///uPuvqNvv1Eqij179nDgwIFFiyGEEKPCzL42pJ+2noQQQmSRohBCCJFFikIIIUQWKQohhBBZpCiEEEJkkaIQQgiRRYpCCCFEllEqisPHTixaBCGEWBlGqSgA9rz7jlN+CiGEmA02xoOLtq7t9acff2TRYgghxKgws/vdfb1vv1GW8IDTLYmj179mQZIIIcTmZpRbT/t2bePo9a955gXaihJCiFmxqbae9rz7DlkWQgiRYOjW02gtihSyKIQQYrqMUlG07Hn3HUlfhRSGEEJMh1EqijaPovVRxJRCdwtKSkMIIYYzSh/F+vq6H7/suui90KKQz0IIITZYKR8FcFrUU8yCSFkbQggh6hmtokihyCchhJguo916Sp2ZnbImpDyEEKvOym09QTzJLuXgDttqS0oIIeoYtaJoFUJoLdRsP7V9pTCEECLPqBUF9Mub6FoTUhBCCFHHRIrCzD5gZofM7KCZ/ZGZvbi5/nwz+/3m3n1mdmGi/7lmdq+ZPWJmnzGzLTXzHj524pSFPrQgQksjvNf+VHKeEEKUmdSi+JC7v9TdLwJuB97fXH8vcNDdXwq8Efhwov+vADe4+17gr4Cf7zN5axnkFvxQacQyueXoFkKINBMpCnd/qvPxh4E2hOoC4HNNm4eAPWb2o92+ZmbAK4Fbm0ufAF5bM29b66nWcd1FikEIIfoxsY/CzD5oZt8A/jnPWhQPAv+0uX8J8BLg7KDrC4HvuPvJ5vNjwK7MPFeb2QEzO/Dkk0+eoiDCZLtwyym8n1Ik8l0IIcTpFBWFmd1jZkcirysB3P197r4buAV4a9PteuD5ZnYQeBvwReBkOHRkumRSh7vf5O7r7r6+Y8eOU+7FFEYsEqq9XqoPJWUhhBDPMrWEOzN7CXCHu18YXDfgq8BLu1tVzfUngZ3uftLMXg78krv/k9Jcba2nXOG/ISGzIcryFkJsJoYm3E2kKMxsr7s/0rx/G/CP3f11ZnYW8H13/4GZ/UvgJ9z9jZH+vwd81t0/bWYfAQ65+2+U5s2dmV2q85Rb+KUYhBCbmUVlZl/fbEMdAl4FvKO5/uPAl8zsIeDVneuY2Z1tGC3wLuCdZvYoGz6L3+ozeRgiGyqJmI8i1z+lXLQVJYRYZUZZ66m1KEILIGdN5MqPD83uFkKIMbGQradFsXVtr69ddWNysY8VA8wplRyxflIgQogxsnJFAbu1mvr4IvpWlO1aJ9qCEkKsIqNUFG3CHZyeK9H1NeS2k/pYBd3xFEIrhFg1RqkouoRbTjElkdsu6rPg63wLIcQqMmpFUVOOI6Y8Ugt+TmnE5pFVIYRYBUatKCC/WKfuxc6iCMNjh/o9hBBiszHKqKf2KNTQMkiFvnbv11STzVkcOb+HIqKEEMvMSoXHhmdm53ImSpnY3fal6937Ug5CiLGxcuGxLd1opDAqKbQaUltEqbO3S/MOuSeEEGNjU1kUNQl14fVJLYLQtyHrQgixrKyURXH42AkgrSBSWdhdahb2nAXSnSfMsRBCiM3EKBVFN+EOymdmd9u190OndZjl3f0Ze587+yKcT1tRQogxM0pF0VoUsSzr7rU+VkOsX6xvKloq5+eQpSGEGDOjVBT7dm07zRro/myJLdAxBVI6wyIcL+cLid2TRSGEGDNnLFqASQnrN6VKerT3U31jn1PXWnKly2PyCCHEGBmlojh87ARrxBfxvsX+WmJKpUZJlOZURJQQYuyMcusJ6stp5JzPuTLlbbuYgzscP7YNFpM15hQXQohlZ5QWReujCMtyQHr7pyWnFML3qQzsXCnzknXRlUnVaIUQo8DdR/fav3+/53jJu25/5mf3fexnrF+qfaxPrG9t+5wsQggxbYADPmDNHe3WU0t3+ygXgVSKaArJhd628+ZkCuXL3Zc1IYRYZkZbwuP4Zdf1KvudCostFf6rXchzEVN9SowIIcSsWLkSHrULbLe0RpiRnXJqhwt8jSUSc1jX+Ctic8jRLYRYJkbrzE4ROohzdZ9iNZpSDuZaCyCWy9HtFyqhcC5tRQkhlo2JLAoz+4CZHTKzg2b2R2b24ub6883s95t795nZhYn+N5vZV5v+B83soknk6RKLiGqvQ3qrqNumdC0kl/EdztkqqFRJECGEWBYm8lGY2fPc/anm/duBC9z9X5vZh4Dvuft1ZnY+8F/d/Scj/W8Gbnf3W/vM2z3hLpdJ3Sf8tKb8R6n/EKugpq2sDCHENFiIj6JVEg0/DLRa5wLgc02bh4A9Zvajk8wVEi7MqSS23LZTrF24FVUjR9uvlSFlqYTXahWAlIQQYqEMiantvoAPAt8AjgA7mmv/Efi15v0lwElgf6TvzcDDwCHgBmBrZp6rgQPAgXPOOcfdy7kIYV5ELLdhFjkMsXyMrgxD51a+hRBiEphVHoWZ3WNmRyKvKxtF8z533w3cAry16XY98HwzOwi8DfhioyxC3gOcD7wMeAHwroxCu8nd1919fceOHcCpUUy1W1AhNd/+Y6QsmJilktoCk6UghBgDxagnd7+scqxPAXcAv+gbW1I/B2BmBny1eYVjP968fdrMPg5cUzkXUHcGdu5MiFr/QKp/n0qzbftp+BtU+kMIMU8mjXra2/l4BfBQc/0sM9vSXH8z8Kd+qj+j7b/W/DTgtWxsX/Wia1WEC2gpGztnTXQVRCmxrzYKKjZnH2IRU4qWEkLMmkkT7q5vtqEOAa8C3tFc/3HgS2b2EPDqznXM7M42jBa4xcwOA4eB7cAv95k8V5CvtkLrkGinMP+i1ipJlfKIKahY29g9WRRCiFkzyhIeW9f2+tpVNwJlH0PN9tCswk9L45aS8UrlRYQQog8rVcJj365t2WzqcGsmVZqj2ydGn22dkgWTshRSc8d+vz4htUIIMTWGhEot+rVl53nung8XrQ1FrSk/Pim5kuO50ua1Y6T6CCFEFwaGxy580R/y2rLzvGReREguv2JI3z7t++ZM5NrXnqshhSGESLFSiiJ3cFHtQjlrC2IaY/dVAjoISQiRY6iiGKWPIkcqX6LkIxgSZpqKYqo94Kg0Vi6qK0bof1HorBBiKgzRLot+9TkKddYM2V4qbRv12cLqbsHJkhBC5GAVLYrUN+bUmRCz/oYdS/yLtUllaKcS6GLj5hICZUkIIabJqBVF7hS61DZMaRHtVqKtCXeNyZTbeip9bseokbUrb04BCSHERAwxQxb9Km095RiyLVUTWlsb/lqqcjtEpnDraZJILyHE5oWBW0+jzMzuHlzUUsqArj1UaJoHCZXkG3rYUTh2aY6wjawOIVaToZnZC7cOhrxqnNm5b/g1365rzrfIzRebozbpryR3zb0aC0cIsVqgPIrh2culhbr2cKQ+c/ahpGxySqkUQSXFIcTqMFRRjNqZHZKqj9T+LNVVqslTiDFplFGpYmwu6qmVK1YPKrwW/g59j30VQqwoQ7TLol85iyL2eUiZixpHdThWTY5ETq7UGNPavuqOL+e2EKsHq7j1FFvMJ63PFI4ba1Prn+gz35BtoEnaS1kIsXoMVRTFo1CXne6WSi6KKLX1lDqbIrVllSsBktrCSs1dOnei1CYmV+x3zm1ZCSFEiVGHx3YZEvJZGwoLpyqi2gOGwjHCcXIyDAnBjckZklOMfX4nIcT4WNnw2FlE7pS2k8L7tVFJ3b59/SQp+Wp9ErExtfUkxGrBKvoohtAne7kUdprLtO6rDFKfa/sOCd2NtZ/GOEKI5WSoohjl1tPWtb3+9OOPPPM55ado7+W2doZmKQ/Zoin5N4ac550Kky3JFso/5BkKIcbF0K2nUSqKsIRHyambos8CmFuAh5biiC3WNYoit8jH2tScwTG0lIgQYjyspKLokloMu/f7Lsix8Wu+tQ+tF5VScjGHc8pZ3SUma66Mefd9TgkJIcbJSjqzQ4ds7R57TTJdzoEcu993/JoxJk0MzMk7iW9ECDFOWLSPwsyuAT4E7HD342ZmwIeBy4HvA29y9wci/fYDNwM/BNwJvMMLQq2vr/vxy65L3s99K85920/lTMS+zXfbl6yNWn/AMmz75CyIZZBPCDGchVoUwG7gbuBrwPbm2uXAXYABlwL3JvreB7y8aXcX8OrSfLGop1yJjD6hrLFxY5bLUCaJKipZUENk7BO1NYtQZCHE/GCRFoWZ3Qp8APgDYN03LIqPAp93999t2jwMvMLdH+/0WwP+xN3Pbz6/oWnzr3LzdX0UQ85xGBrlNOlZDimrBPLf4HOJczk/Rel8ilKRQ0U+CbG5WJhFAVwBfLh5f5RnLYrbgX/Yafc5NpRIt+86cE/n808AtyfmuRo4ABzYsvO8uX2zndY8Q/wV3W/2uTpNpdpUJbn6IqtCiHHCLBPugHuAI5HXlcC9wDY/XVHcEVEU+4NxXxZRFH9Ykqe09ZS6V9qKKvUPx+g7Tkmumvlr24VKpo/TvmZuZXYLMT5mqiiSnWEf8ESjII4CJ4GvAzuBjwJv6LR9GFgL+q8BD3U+vwH4aGneGkVREwHU95v3LPwUuQW9r4KYJAKsZo5JFaMQYrEsRFGcNtipFsVrONWZfV+izxea+60z+/LSPJMchRprW7o2z0Uxt6WUu969l9qOmoYFIAUhxHgZqihmecLdncBXgEeBjwFvaW+Y2cFOu18AfrNp9xeNspiIVFnw8JS4tm14PecELiW55e6388ROnusjQ4pumfPQaR46v2PPInWt+zP2vGLthBCbh02RmZ1aYEvZ2jFqM5m77ft8TsnWna+Uw5DK3yideVEqJT5JTocyuIVYfoZGPY3yzOzDx05Ev+l2P4d0F8jab+btYtm+uqQWxtTn0HpILaip36PbvhTW2v0cWhfhta5c4Rw11kG3b03IrhBihAzZr1r0a//+/dn99ppInbBdzBGca1uKsir5EUqy5saOyRz7HVLtc3MPjYLq3pcfQ4jlhEWX8Jgnsa2n1JZJTYJZjNK2Tk0pkD5bObXbSLlv7bEtr9p5u2Pk5ql9DtqKEmL5WMmigCGTfhtOtc19Oy/lJ9TkMoT9+n7T7zN2qk9q/JqxS1ZTrI0QYv6wDOGx83qFimKS7ZtU/z4KIHW97wJbK3tJSdQu1DkZw9+tZp5JnqsQYvastKJIUfIV1PRLtc8tgkMtlxrZSv1KSmqof6Tkt4mNnfKhCCEWw0opir61nmILaM02Uo2i6bM4D7F8aiyI1NhDHc99LKGh1pyUhhDzZ6iiGGV47L5d26pOkEsl2PUhlbDWvs+FzsbGiMmUkj8WmhubL5QvVf01J2MpUTAcu2aOGOHzFEIsP6OPeipF1+QURSwiqibRLIxQivWviVyqJfU71CT3lZIRaxPwcmPF5pMiEGL5WKmEu5ZYGGfXksglzEE8SaxPYl5uUeyrJGJzdn+H7pjh9Zjsqd8zNW94LczizimJMJEwFVqrRDwhRsqQ/apFv7bsPO+UfbdJk7xK/Ut+gpKTNzZeH9ly70uRSTVjhG37RDhN+tyFEPMDObNPpSbsM3YttwBOK9IpJ9Mk4a65+7WKpI9ctXOH16QghFgMK6coWlLfcvssrkMWrqGLXZ8Fu2TlTEOWof1TUWR9FLQQYr6slKLIHVxU+tY/zQWrxsKYVAnlFELfLbBSiG+NPDGZan9XKQohFstKKYrU1lOtr2HSBSu2wNYukn39E7PKO0hZBH3H6I6Vui+EWA6GKorRRz11f3bLZ+eIldjuE5EzyTkNtX3D6KfYvZTMtfkM3QilIXKV+ipEVojNwSgVxb5d2555H4ZxQn7h6i5ufc54CMNuUyGpuZyH2LjdsbvvS+dcpOaLyTQtcnkauTBlIcS4GWXC3da1vb521Y2nXc+VBG/vh+SSymJ9+mQfl+YukSr3HUu4y5VC75vk17dEed+S6kKIxTA04W6UimJ9fd2PX3Zd8n534cottiklEetfqyCGZGTnFuY+i/20lNLQMZWVLcRys1LnUYTO7FyI7LTI5TdMEqUUjjM0xyPsM0/ktBZiHLBKUU+l8NhJwmFzIaS59qXFvObzLJLS5pljIYUhxHKzkooiF5Y5acJd7tv+LKyWWVsLs1zEpSCEGAdDFcVofRQHDhyoinQKyVVcTRXay43VpdbxXCtbbr4+ck5CjUzySQgxDhZaPdbMrjEzN7PtzWczs183s0fN7JCZXZzo93kze9jMDjavF/WZN+aIrnE6x5za4Xipiqs1FWXDCqopp3lM0ZSqucbkzMk1JE8klGmS+0KITcAQM6T7AnYDdwNfA7Y31y4H7gIMuBS4N9H388B63zm7PopZlJ6oye6uuVfjiO4jf+l+n601IcTqwcCtpzOmoGtuAK4F/qBz7Urgk41gf25mZ5nZmrs/PoX5OHzsRDZfAtIHFJW2akrbT7nQ1/D8hppv+7lkwdK2VSwpr1ZmIYSoZaKtJzO7Ajjm7g8Gt3YB3+h8fqy5FuPjzbbTvzczq5m3exRqbKEtLaqxBbzvgT/dfuGCXsqyTmV1xw4K6usDqJFZCCH6UHRmm9k9wM7IrfcB7wVe5e4nzOwoG9tIx83sDuA/ufv/bsb4HHCtu98fjL3L3Y+Z2XOBzwK/4+6fTMhxNXA1wDnnnLPf3vAbWblrMpbDz6WEvBKpTOXYvH2znmO/W6rftDOk5bAWYnMwM2e2u1/m7heGL+ArwLnAg42SOBt4wMx2smFB7O4MczbwzcjYx5qf3wU+BVySkeMmd1939/Vv/WDLM9fbb9B9tlpCR3PsXmycWP2iXCRSrv2QKK2ccz1WcylmofRFSkIIMXjryd0Pu/uL3H2Pu+9hQzlc7O7fAm4D3thEP10KnAj9E2Z2RidK6kzgZ4AjNXO3RQFzUUItoeUQuxdaAClFktrSim0T5RRYd9yU/H23kHK+jpTsNUqktnyJEGLzMqvqsXeyYXE8CnwMeEt7w8wONm+3Aneb2SHgIHCsaVvk8LET2X3/0oJcu/8fsxBCUoqnNGbN/LF2uf4l30zK+kg5/nPjCiFWiCGhUot+dY9CDemTkZ2636feUm6+0mE+pTFK4b21v0+JmLwKqRVi88EqlfBoiwLmFu2akhg1CqH2fqzMR62CGqLc+rbrM4aUhBCbk5VSFLmigO371Df3VJ9Y25qCgCWrIdVnEuVQomYcKQMhVo+VVBQ1lVuHVj+tLQA4SWHAaY7ZZ/sr1W+aRQ6FEMvHUEUxyqNQ+9A3cgjyDt1uyGspO7x2zpiTuW+0USrKqSbnIhYIIIQQLaNWFOFiWCq813cRD/vGwl1zY6byKFK5FEMiompIyVijTIQQYuHbSENekxQFzG315PwH4XhDI6v6FO4bGvVUM7YQYvVglbae2qKAkC/uF1Kq8dTHMqnZGopZD7migil5Y59lBQgh5sUoFUWbmV2iVL21bRNbgGMZ16mFOxy/z55/bc2nSchlZQshRIlRn3BXS8xpW2obXuv2Kzmxp1GUTzWWhBDTZqEn3C0DJauh+zPXP2YlhGVBakt0pL7JlxzsuTlSv6csBCHErNiUFkXMcqi1Jrpthpb3jl1PlQ2v3ZoKx6/ZspJFIoToslIWxeFjJ4D82dMhNYtmrvx4bd9c/5oS5rE+sUU/53+RkhBCTJUhoVKLfoWZ2eH79nMs1LVvDai+BQVr7tWUGymNo5BXIURfWKUSHm1RQPfTi/G112LkciRS/Uq5GH1KcKT6TVICRAghahmqKEa59dSemZ0KQ+0TDtrHX1A7fup8h1w0VU1ORo0cQggxbUbrzD5+2XXRhTfm8O1erwl/DUn5QMJxwzGGONBz1Pge5J8QQqQY6swepaLYurbX1666MakUIP7tvU9+Qy6KqR0/tBZKTvWcYkkpMC36QohpMVRRLNzfMOQVq/XU58Ci1P2Yk7nmbImacWpkGOrvEEKIGhjooxitRfH044+ccq12C2fSLO3YGDXzlshZFLOYTwixeqysRZFjWmGuNf2GhMl27w8JkRVCiD6wSlFPMVLJdzFfQOx8iFSkUmrcltozHbrzxnwZpaitmgKHQggxC0avKGqc111yUUq5OVKZ1yVF0reSbOx9TH4hhJgXZyxagFkRftPvLvalMyFK0VF9Iqb6KCadNyGEWEZG6cwOiwL2cWS37WrCX7vXSsqllj59FB4rhJgmCy0KaGbXmJmb2fbm8/lm9mdm9rSZXZPpd66Z3Wtmj5jZZ8xsS+2cqTLg4f1Uu1Q2dNumxhfQZxGvKTTYp1qtEELMi4kVhZntBn4K+Hrn8reBtwO/Wuj+K8AN7r4X+Cvg52vmPHzsRLFaa8o5nSMVfhpTHLnPrTzdV06R5RzdQgixaKZhUdwAXAs8s4fl7k+4+xeAv051MjMDXgnc2lz6BPDa2klzC2p30a3ZXoopmvBaLgs8FcUU+iliFk2sjRBCLBMTObPN7ArgmLs/uLHu9+KFwHfc/WTz+TFgV2auq4GrAc455xwgX18pJOdgbq91lUPMCT1kayjnGM8piT5RWUIIMVNKiRbAPcCRyOtK4F5gW9PuKLA96PtLwDWJcXcAj3Y+7wYO1yR/xM6jCCmVBw+v9ymfkZtraGJfdy4l1wkhZgHzPo8C2Ac80SiIo8BJNvwUOzttcorCgOPAGc3nlwN318y9Zed5z/zipQW+9pCgUGHUKpo+lA5aEkKIWTJ3RXHaQD0tiub+7wGvb95/BHhLzVypEh6TLLY1hfxKbWpO0pNCEEIsiqGKYmp5FGZ2FFh39+NmthM4ADwP+Bvge8AF7v6Umd0JvNndv2lmPwZ8GngB8EXgZ9396dJcYVHAWAnx3BkS4bUUqbyLXBRTOK7KhwshloWVP48iJHZoUXi/pc+5EalxUs7vWoUihBCzRtVjK0j5LWJtasZp34fnd/eRRQgh5gWrVD328LETyXu5/IowhyFmOXTHyeVJlOaKJeB138uaEEKMhU1XFDC11ZPKYUjlScQys/sWCMzlTgghxFgYpY8iLAoYMq1v7LUlzGOJe11qCgwKIcSsGeqj2DQWRc25DznHc4zc9Zj1UVIYuTGFEGJZGaVFkTszG/IRSkPKgufGHDquEELMm4WWGZ83+3ZtA06tvBpaErXf5mPtciXMa06h695TNVghxNgZpaJo6S7iYTRRLpopHCMs9Z0rYZ6LeCpVrRVCiFEyJKZ20a8tO8+ryoVIldnIFd/rk0ehXAghxJhgYB7FKJ3Z+3Zt43jnc9dxHAtHLYXIwrDDjWLXZUkIITYbo3Rmx87MjtVj6pLLkxiyyE96ToUQQsyblav1FEY91TLLRV2+CSHEMrNSUU8wvWiiSSKXcselCiHEZmGUimLfrm3VIavhYl8bItu2zVkJbdSVQmCFEJuZUW49lUp4QPyMivB9+xnKCXWlNkIIseys3NZTl9gWUc0hQ91EvTCXottOloMQYpXZFIoiJOc3CKOicolzuUgpIYRYFUarKGKZ1OHPFLF2sWu15TqEEGIzM0ofRTc8tnQmdZhjURseO8SXIYQQy8xK+SjaooApYpZC7H7uPOvSuEIIsSqMUlF0qYlWivkZ+i76UhJCiFVl9Ipi2sgXIYQQpzJKH0VNHkUXldYQQogF+yjM7BozczPb3nw+38z+zMyeNrNrMv1uNrOvmtnB5nVR7Zx9v/nLUhBCiGFMrCjMbDfwU8DXO5e/Dbwd+NWKIf6tu1/UvA4OkSH0RYRhrUPKiQshhNhgGhbFDcC1wDN7WO7+hLt/AfjrKYwfJXbORKzUeO2WkxSIEELEmUhRmNkVwDF3f3CCYT5oZofM7AYz25qZ62ozO2BmB5588kmg35kSuVyKXD8hhFh1iorCzO4xsyOR15XA+4D3TzD/e4DzgZcBLwDelWro7je5+7q7r+/YseO0+7nEuVSSnawIIYQoU1QU7n6Zu18YvoCvAOcCD5rZUeBs4AEz21k7ubs/3hzl+jTwceCSmn6Hj52oWuRTW1KgSCghhKhl8NaTux929xe5+x533wM8Blzs7t+qHcPM1pqfBrwWODJElpwzO3WmhKrBCiFEHTNJuDOznWb2GPBO4N+Z2WNm9rzm3p1m9uKm6S1mdhg4DGwHfrlm/O7BRaVT6CBdqkMWhRBClNl0CXdSAkIIEWeligJ2iTmqpSSEEGJ6jF5RKIJJCCFmy2gVRe25EkIIISZjtIqihCwMIYSYDqNVFDELohsBJQtDCCGmw2gVBZweGjv0UCIhhBBpRq0oWqQYhBBidmwKRSGEEGJ2jFpRyJIQQojZM1pFoagmIYSYD6NVFLImhBBiPoxWUQghhJgPo1UU2noSQoj5MFpFoa0nIYSYD6NVFEIIIeaDFIUQQogsUhRCCCGySFEIIYTIIkUhhBAiixSFEEKILFIUQgghskhRCCGEyCJFIYQQIou5+6Jl6I2ZfRd4eNFyVLAdOL5oISoYg5xjkBEk57SRnNPl77j7c/t2OmMWksyBh919fdFClDCzA5JzOoxBRpCc00ZyThczOzCkn7aehBBCZJGiEEIIkWWsiuKmRQtQieScHmOQESTntJGc02WQnKN0ZgshhJgfY7UohBBCzAkpCiGEEFmWWlGY2U+b2cNm9qiZvTtyf6uZfaa5f6+Z7VlCGd9kZk+a2cHm9eZ5y9jI8dtm9oSZHUncNzP79eb3OGRmF89bxkaOkpyvMLMTnef5/gXIuNvM/sTMvmxmXzKzd0TaLPx5Vsq5DM/zb5vZfWb2YCPndZE2y/C3XiPnsvy9P8fMvmhmt5qKmucAAAOHSURBVEfu9X+W7r6UL+A5wF8APwZsAR4ELgjavAX4SPP+9cBnllDGNwH/ZQme5z8CLgaOJO5fDtwFGHApcO+SyvkK4PYFP8s14OLm/XOB/xP5d1/486yUcxmepwE/0rw/E7gXuDRos9C/9R5yLsvf+zuBT8X+bYc8y2W2KC4BHnX3r7j7D4BPA1cGba4EPtG8vxX4STOzJZNxKXD3PwW+nWlyJfBJ3+DPgbPMbG0+0j1LhZwLx90fd/cHmvffBb4M7AqaLfx5Vsq5cJpn9L3m45nNK4yyWfTfeq2cC8fMzgZeA/xmoknvZ7nMimIX8I3O58c4/T/5M23c/SRwAnjhXKQL5m+IyQjwz5rth1vNbPd8ROtN7e+yDLy8Mf/vMrO/u0hBGrP977Px7bLLUj3PjJywBM+z2So5CDwB/C93Tz7PBf2tA1VywuL/3m8ErgX+JnG/97NcZkUR03Ch9q5pM0tq5v9DYI+7vxS4h2c1+bKx6GdZywPAS9z97wH/GfifixLEzH4E+Czwb9z9qfB2pMtCnmdBzqV4nu7+/9z9IuBs4BIzuzBoshTPs0LOhf69m9nPAE+4+/25ZpFr2We5zIriMaCrjc8GvplqY2ZnANuY77ZFUUZ3/0t3f7r5+DFg/5xk60vN81447v5Ua/67+53AmWa2fd5ymNmZbCy+t7j7/4g0WYrnWZJzWZ5nR57vAJ8Hfjq4tei/9VNIybkEf+//ALjCzI6ysRX+SjP7naBN72e5zIriC8BeMzvXzLaw4XS5LWhzG3BV8/51wB9746FZFhmDfekr2NgnXkZuA97YROtcCpxw98cXLVSIme1s91PN7BI2/g//5ZxlMOC3gC+7+68lmi38edbIuSTPc4eZndW8/yHgMuChoNmi/9ar5Fz037u7v8fdz3b3PWysR3/s7j8bNOv9LJe2eqy7nzSztwJ3sxFd9Nvu/iUz+w/AAXe/jY0/gv9mZo+yoRFfv4Qyvt3MrgBONjK+aZ4ytpjZ77IR4bLdzB4DfpENZxzu/hHgTjYidR4Fvg/83JLK+TrgF8zsJPB/gdfPe8Fg41vbvwAON/vVAO8FzunIuQzPs0bOZXiea8AnzOw5bCiq/+7uty/T33oPOZfi7z1k0mepEh5CCCGyLPPWkxBCiCVAikIIIUQWKQohhBBZpCiEEEJkkaIQQgiRRYpCCCFEFikKIYQQWf4/cwtzuttE+SQAAAAASUVORK5CYII=\n",
	"text/plain": [
	"<Figure size 432x288 with 1 Axes>"
	]
	},
	"metadata": {
	"needs_background": "light"
	},
	"output_type": "display_data"
	}
	],
	"source": [
	"plot(ra_chunk, dec_chunk, ',')\n",
	"xlim(0, 4)"
	]
	},
	{
	"cell_type": "code",
	"execution_count": 8,
	"id": "963a8b30-d767-42b7-afd9-e6909c0591af",
	"metadata": {},
	"outputs": [],
	"source": [
	"# would get the full column:\n",
	"# g[\"ra\"][:]"
	]
	}
	],
	"metadata": {
	"kernelspec": {
	"display_name": "Python 3",
	"language": "python",
	"name": "python3"
	},
	"language_info": {
	"codemirror_mode": {
	"name": "ipython",
	"version": 3
	},
	"file_extension": ".py",
	"mimetype": "text/x-python",
	"name": "python",
	"nbconvert_exporter": "python",
	"pygments_lexer": "ipython3",
	"version": "3.7.4"
	}
	},
	"nbformat": 4,
	"nbformat_minor": 5
	}