htcondor_demo.ipynb

{
 "cells": [
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "# HTCondor JobSubmission Demo"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Submiting a test job"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "from htcondor_utilities import run_testjob, check_htcondor_status, check_output, get_output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "A test case with the following parameters:\n",
      "Ehtdata:  torus.out0.05992.h5,torus.out0.05993.h5,torus.out0.05994.h5,torus.out0.05995.h5,torus.out0.05996.h5\n",
      "Rhigh-Ratio:  10:100:10\n",
      "Theta-Viewing-Angle:  30,50,60,70,80,90\n",
      "Rho0-Density-Normalization:  4.266338570441294e+17\n",
      "\n",
      "\n",
      "\n",
      "\n",
      "total jobs:  300\n",
      "ready to submit\n"
     ]
    }
   ],
   "source": [
    "run_testjob(submit=False)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Check the job status"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "JobID: 70, Total: 300, Submit Time: 2/27 09:31\n",
      "status report 61.33% done @ 02/29 10:21\n",
      "Done - 184:  ****************************************************************************************************************************************************************************************\n",
      "Run - 112:  ****************************************************************************************************************\n",
      "Idle - 0:  \n",
      "Hold - 4:  ****\n"
     ]
    }
   ],
   "source": [
    "status = check_htcondor_status()"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Check the output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Total output: 202\n",
      "['img_s05992_Rh100_i60.h5', 'img_s05992_Rh100_i70.h5', 'img_s05992_Rh100_i90.h5', 'img_s05992_Rh10_i60.h5', 'img_s05992_Rh30_i60.h5', 'img_s05992_Rh30_i70.h5', 'img_s05992_Rh30_i90.h5', 'img_s05992_Rh40_i30.h5', 'img_s05992_Rh40_i50.h5', 'img_s05992_Rh40_i60.h5', 'img_s05992_Rh40_i70.h5', 'img_s05992_Rh40_i80.h5', 'img_s05992_Rh40_i90.h5', 'img_s05992_Rh50_i50.h5', 'img_s05992_Rh50_i60.h5', 'img_s05992_Rh50_i70.h5', 'img_s05992_Rh50_i80.h5', 'img_s05992_Rh60_i60.h5', 'img_s05992_Rh60_i70.h5', 'img_s05992_Rh60_i90.h5', 'img_s05992_Rh70_i30.h5', 'img_s05992_Rh70_i50.h5', 'img_s05992_Rh70_i60.h5', 'img_s05992_Rh80_i50.h5', 'img_s05992_Rh80_i60.h5', 'img_s05992_Rh80_i70.h5', 'img_s05992_Rh80_i80.h5', 'img_s05992_Rh80_i90.h5', 'img_s05992_Rh90_i30.h5', 'img_s05992_Rh90_i50.h5', 'img_s05992_Rh90_i70.h5', 'img_s05992_Rh90_i90.h5', 'img_s05993_Rh100_i50.h5', 'img_s05993_Rh100_i60.h5', 'img_s05993_Rh100_i70.h5', 'img_s05993_Rh100_i80.h5', 'img_s05993_Rh10_i30.h5', 'img_s05993_Rh20_i60.h5', 'img_s05993_Rh20_i70.h5', 'img_s05993_Rh20_i80.h5', 'img_s05993_Rh30_i60.h5', 'img_s05993_Rh40_i80.h5', 'img_s05993_Rh50_i60.h5', 'img_s05993_Rh50_i70.h5', 'img_s05993_Rh50_i90.h5', 'img_s05993_Rh70_i70.h5', 'img_s05993_Rh70_i80.h5', 'img_s05993_Rh80_i90.h5', 'img_s05993_Rh90_i50.h5', 'img_s05993_Rh90_i90.h5', 'img_s05994_Rh100_i30.h5', 'img_s05994_Rh100_i50.h5', 'img_s05994_Rh100_i60.h5', 'img_s05994_Rh100_i80.h5', 'img_s05994_Rh100_i90.h5', 'img_s05994_Rh10_i30.h5', 'img_s05994_Rh10_i50.h5', 'img_s05994_Rh10_i60.h5', 'img_s05994_Rh10_i70.h5', 'img_s05994_Rh10_i80.h5', 'img_s05994_Rh10_i90.h5', 'img_s05994_Rh20_i30.h5', 'img_s05994_Rh20_i50.h5', 'img_s05994_Rh20_i60.h5', 'img_s05994_Rh20_i70.h5', 'img_s05994_Rh20_i80.h5', 'img_s05994_Rh20_i90.h5', 'img_s05994_Rh30_i30.h5', 'img_s05994_Rh30_i50.h5', 'img_s05994_Rh30_i60.h5', 'img_s05994_Rh30_i70.h5', 'img_s05994_Rh30_i80.h5', 'img_s05994_Rh30_i90.h5', 'img_s05994_Rh40_i30.h5', 'img_s05994_Rh40_i50.h5', 'img_s05994_Rh40_i60.h5', 'img_s05994_Rh40_i70.h5', 'img_s05994_Rh40_i90.h5', 'img_s05994_Rh50_i50.h5', 'img_s05994_Rh50_i60.h5', 'img_s05994_Rh50_i70.h5', 'img_s05994_Rh50_i80.h5', 'img_s05994_Rh50_i90.h5', 'img_s05994_Rh60_i30.h5', 'img_s05994_Rh60_i50.h5', 'img_s05994_Rh60_i60.h5', 'img_s05994_Rh60_i70.h5', 'img_s05994_Rh60_i80.h5', 'img_s05994_Rh60_i90.h5', 'img_s05994_Rh70_i30.h5', 'img_s05994_Rh70_i60.h5', 'img_s05994_Rh70_i70.h5', 'img_s05994_Rh70_i80.h5', 'img_s05994_Rh70_i90.h5', 'img_s05994_Rh80_i30.h5', 'img_s05994_Rh80_i60.h5', 'img_s05994_Rh80_i70.h5', 'img_s05994_Rh80_i80.h5', 'img_s05994_Rh80_i90.h5', 'img_s05994_Rh90_i30.h5', 'img_s05994_Rh90_i50.h5', 'img_s05994_Rh90_i70.h5', 'img_s05994_Rh90_i80.h5', 'img_s05994_Rh90_i90.h5', 'img_s05995_Rh100_i30.h5', 'img_s05995_Rh100_i50.h5', 'img_s05995_Rh100_i70.h5', 'img_s05995_Rh100_i80.h5', 'img_s05995_Rh100_i90.h5', 'img_s05995_Rh10_i30.h5', 'img_s05995_Rh10_i50.h5', 'img_s05995_Rh10_i60.h5', 'img_s05995_Rh10_i70.h5', 'img_s05995_Rh10_i80.h5', 'img_s05995_Rh10_i90.h5', 'img_s05995_Rh20_i30.h5', 'img_s05995_Rh20_i50.h5', 'img_s05995_Rh20_i60.h5', 'img_s05995_Rh20_i70.h5', 'img_s05995_Rh20_i80.h5', 'img_s05995_Rh30_i70.h5', 'img_s05995_Rh30_i80.h5', 'img_s05995_Rh30_i90.h5', 'img_s05995_Rh40_i30.h5', 'img_s05995_Rh40_i50.h5', 'img_s05995_Rh40_i70.h5', 'img_s05995_Rh40_i80.h5', 'img_s05995_Rh50_i50.h5', 'img_s05995_Rh50_i60.h5', 'img_s05995_Rh50_i80.h5', 'img_s05995_Rh60_i30.h5', 'img_s05995_Rh60_i60.h5', 'img_s05995_Rh60_i70.h5', 'img_s05995_Rh60_i80.h5', 'img_s05995_Rh60_i90.h5', 'img_s05995_Rh70_i30.h5', 'img_s05995_Rh70_i50.h5', 'img_s05995_Rh70_i70.h5', 'img_s05995_Rh70_i80.h5', 'img_s05995_Rh70_i90.h5', 'img_s05995_Rh80_i60.h5', 'img_s05995_Rh80_i90.h5', 'img_s05995_Rh90_i30.h5', 'img_s05995_Rh90_i50.h5', 'img_s05995_Rh90_i60.h5', 'img_s05995_Rh90_i70.h5', 'img_s05995_Rh90_i80.h5', 'img_s05996_Rh100_i30.h5', 'img_s05996_Rh100_i50.h5', 'img_s05996_Rh100_i60.h5', 'img_s05996_Rh100_i70.h5', 'img_s05996_Rh100_i80.h5', 'img_s05996_Rh100_i90.h5', 'img_s05996_Rh10_i30.h5', 'img_s05996_Rh10_i50.h5', 'img_s05996_Rh10_i60.h5', 'img_s05996_Rh10_i70.h5', 'img_s05996_Rh10_i80.h5', 'img_s05996_Rh10_i90.h5', 'img_s05996_Rh20_i30.h5', 'img_s05996_Rh20_i50.h5', 'img_s05996_Rh20_i70.h5', 'img_s05996_Rh20_i80.h5', 'img_s05996_Rh30_i30.h5', 'img_s05996_Rh30_i50.h5', 'img_s05996_Rh30_i60.h5', 'img_s05996_Rh30_i70.h5', 'img_s05996_Rh30_i80.h5', 'img_s05996_Rh30_i90.h5', 'img_s05996_Rh40_i30.h5', 'img_s05996_Rh40_i50.h5', 'img_s05996_Rh40_i60.h5', 'img_s05996_Rh40_i70.h5', 'img_s05996_Rh40_i80.h5', 'img_s05996_Rh40_i90.h5', 'img_s05996_Rh50_i30.h5', 'img_s05996_Rh50_i50.h5', 'img_s05996_Rh50_i60.h5', 'img_s05996_Rh50_i70.h5', 'img_s05996_Rh60_i30.h5', 'img_s05996_Rh60_i50.h5', 'img_s05996_Rh60_i70.h5', 'img_s05996_Rh60_i80.h5', 'img_s05996_Rh60_i90.h5', 'img_s05996_Rh70_i30.h5', 'img_s05996_Rh70_i50.h5', 'img_s05996_Rh70_i60.h5', 'img_s05996_Rh70_i70.h5', 'img_s05996_Rh70_i80.h5', 'img_s05996_Rh70_i90.h5', 'img_s05996_Rh80_i30.h5', 'img_s05996_Rh80_i50.h5', 'img_s05996_Rh80_i60.h5', 'img_s05996_Rh80_i70.h5', 'img_s05996_Rh80_i80.h5', 'img_s05996_Rh80_i90.h5', 'img_s05996_Rh90_i30.h5', 'img_s05996_Rh90_i50.h5', 'img_s05996_Rh90_i60.h5', 'img_s05996_Rh90_i70.h5', 'img_s05996_Rh90_i80.h5', 'img_s05996_Rh90_i90.h5']\n"
     ]
    }
   ],
   "source": [
    "h5list = check_output()\n",
    "print(f\"Total output: {len(h5list)}\")\n",
    "print(h5list)"
   ]
  },
  {
   "cell_type": "markdown",
   "metadata": {},
   "source": [
    "## Get the output"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "/home/cicuser/Projects/eht_gateway/img_s05992_Rh100_i60.h5\n",
      "File transfer successful!\n"
     ]
    }
   ],
   "source": [
    "# download the first output as an example\n",
    "# the file is downloaded to the current local dir\n",
    "get_output(h5list[0])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [],
   "source": [
    "# plot the file\n",
    "import h5py\n",
    "import matplotlib.pyplot as plt\n",
    "import numpy as np"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAcIAAAGiCAYAAAB07MYfAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjguMywgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/H5lhTAAAACXBIWXMAAA9hAAAPYQGoP6dpAABWVElEQVR4nO29e5AdZZ3//366+5wzZyYzJ5dJ5iIhBBfWSxB/BheS/WqCSIAvFxXKBaliSYkoC0EjUK6ILiEKURfQKlhkXakgu7qwVSteCksBISBfUGMMZQAvICEJkCEmJDOZycy5dD+/P55++nb6XGbmnDkz0+9XVVfmdPfpfrrPSb/P5/N8LkJKKUEIIYQkFKPVAyCEEEJaCYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiiaaoQPvnkkzj33HPR398PIQR++MMfhrZLKbFhwwb09/cjm81i9erVeP7550P75PN5XH311eju7kZHRwfOO+88vPrqq80cNiGEkATRVCEcGRnBiSeeiDvvvDN2+9e//nXcfvvtuPPOO7F161b09vbi9NNPx+HDh7191q9fjwcffBD3338/nnrqKQwPD+Occ86BbdvNHDohhJCEIKaq6LYQAg8++CA+/OEPA1DWYH9/P9avX49//ud/BqCsv56eHnzta1/Dpz71KQwODmLhwoX4z//8T1x44YUAgNdffx2LFy/GT3/6U5xxxhlTMXRCCCGzGKtVJ965cycGBgawZs0ab10mk8GqVavw9NNP41Of+hS2bduGYrEY2qe/vx/Lli3D008/XVEI8/k88vm899pxHLz55ptYsGABhBDNuyhCCCFNQUqJw4cPo7+/H4bRWGdmy4RwYGAAANDT0xNa39PTg127dnn7pNNpzJs3r2wf/f44Nm3ahJtuuqnBIyaEENJq9uzZg6OOOqqhx2yZEGqiFpqUsqbVVmuf66+/Htdcc433enBwEEcffTT+D/4vLKQmN2BCCCFTTglFPIWforOzs+HHbpkQ9vb2AlBWX19fn7d+3759npXY29uLQqGAgwcPhqzCffv2YeXKlRWPnclkkMlkytZbSMESFEJCCJlxuNEszZjealke4dKlS9Hb24tHHnnEW1coFPDEE094Ird8+XKkUqnQPnv37sVzzz1XVQgJIYSQemmqRTg8PIyXXnrJe71z5048++yzmD9/Po4++misX78et9xyC4477jgcd9xxuOWWW9De3o6LL74YAJDL5XDZZZfh2muvxYIFCzB//nxcd911OOGEE/DBD36wmUMnhBCSEJoqhL/97W9x6qmneq/1vN2ll16Ke++9F5/73OcwOjqKK6+8EgcPHsTJJ5+Mhx9+OOQD/sY3vgHLsvAP//APGB0dxWmnnYZ7770Xpmk2c+iEEEISwpTlEbaSoaEh5HI5rMaHOEdICCEzkJIsYgt+hMHBQXR1dTX02Kw1SgghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiioRASQghJNBRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkETTciE85phjIIQoW6666ioAwNq1a8u2nXLKKS0eNSGEkNmC1eoBbN26FbZte6+fe+45nH766fjoRz/qrTvzzDOxefNm73U6nZ7SMRJCCJm9tFwIFy5cGHr91a9+FW9961uxatUqb10mk0Fvb+9UD40QQkgCaLlrNEihUMB//dd/4eMf/ziEEN76LVu2YNGiRTj++ONx+eWXY9++fVWPk8/nMTQ0FFoIIYSQOKaVEP7whz/EoUOHsHbtWm/dWWedhe9973t47LHHcNttt2Hr1q34wAc+gHw+X/E4mzZtQi6X85bFixdPwegJIYTMRISUUrZ6EJozzjgD6XQaP/nJTyrus3fvXixZsgT3338/zj///Nh98vl8SCiHhoawePFirMaHYIlUw8dNCCGkuZRkEVvwIwwODqKrq6uhx275HKFm165dePTRR/GDH/yg6n59fX1YsmQJXnzxxYr7ZDIZZDKZRg+REELILGTauEY3b96MRYsW4eyzz66634EDB7Bnzx709fVN0cgIIYTMZqaFEDqOg82bN+PSSy+FZflG6vDwMK677jo888wzeOWVV7Blyxace+656O7uxkc+8pEWjpgQQshsYVq4Rh999FHs3r0bH//4x0PrTdPEjh07cN999+HQoUPo6+vDqaeeigceeACdnZ0tGi0hhJDZxLQQwjVr1iAuZiebzeLnP/95C0ZECCEkKUwL1yghhBDSKiiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiioRASQghJNBRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIomm5EG7YsAFCiNDS29vrbZdSYsOGDejv70c2m8Xq1avx/PPPt3DEhBBCZhMtF0IAeOc734m9e/d6y44dO7xtX//613H77bfjzjvvxNatW9Hb24vTTz8dhw8fbuGICSGEzBamhRBaloXe3l5vWbhwIQBlDX7zm9/EDTfcgPPPPx/Lli3Dd7/7XRw5cgTf//73WzxqQgghs4FpIYQvvvgi+vv7sXTpUlx00UV4+eWXAQA7d+7EwMAA1qxZ4+2byWSwatUqPP300xWPl8/nMTQ0FFoIIYSQOFouhCeffDLuu+8+/PznP8d//Md/YGBgACtXrsSBAwcwMDAAAOjp6Qm9p6enx9sWx6ZNm5DL5bxl8eLFTb0GQgghM5eWC+FZZ52FCy64ACeccAI++MEP4qGHHgIAfPe73/X2EUKE3iOlLFsX5Prrr8fg4KC37NmzpzmDJ4QQMuNpuRBG6ejowAknnIAXX3zRix6NWn/79u0rsxKDZDIZdHV1hRZCCCEkjmknhPl8Hn/4wx/Q19eHpUuXore3F4888oi3vVAo4IknnsDKlStbOEpCCCGzBavVA7juuutw7rnn4uijj8a+ffvwla98BUNDQ7j00kshhMD69etxyy234LjjjsNxxx2HW265Be3t7bj44otbPXRCCCGzgJYL4auvvoqPfexj2L9/PxYuXIhTTjkFv/rVr7BkyRIAwOc+9zmMjo7iyiuvxMGDB3HyySfj4YcfRmdnZ4tHTgghZDYgpJSy1YNoNkNDQ8jlcliND8ESqVYPhxBCyDgpySK24EcYHBxseNzHtJsjJIQQQqYSCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICUkaVTq3EJJEWl5ijRDSYOoRuvGI4ewvPkUSDoWQkJmCFi8pp9aqi56LwkhmGRRCQqYj1YSu1a5NCiOZZVAICWk1QvhiErT6gq+nM8ExUhTJDIRCSEgr0SISFbyZIIBxUBTJDIRCSGYWQetppjJTRW68RK1bQqYpTJ8gM4eZLCBC+EtLzm/E/zsl527xtRNSA1qEhDSLqXjwj0fQqomhdBo3pqpjoJVIph8UQjIzmI7WRHQ+rJ4xTjT1odkWnDCmTgwBCiKZVlAICRkvcUJWr7jVvV8LZi1oIZKEQiEk05/pFInYbMu0QQIojPGPUzqBexs3jmaKIwWRtBAKISH1MBEBHI8btIoA1i1q+hgTEEE4EsKMrIsIn3SMqtsbwmyICiYzDgohmd5Mh7nB8Y6h3mT4CuJXl/AFRE/EnceoIKxOvHjJuN0jwucJpSuAIWFspCjSOiRTDIWQkDgmI8ATEMCa4icMz9LzhE+LXfB8lQQQUCJoBsy+gNAIvT1ASBxj3KbC8PdviijSOiRTBIWQEP3AbbYlEhHAesUvJHyBv4UQYTdoQKyklL6wBa8nxiKUkesVQpStiz2Pvh7phK7FE8VGCCLFkEwBFEIyfWl2cefg8Sv93bBzGYE/46JOAyKpxS9o8QWFT0SExpHlwuc4YTELWnTSqZouIaPjCWK4Y4gLrHGPp6+vYVYixZA0GQohSSZTNfdYzQqMik1QAE0zLHz6fY5UohIUv4AIeuKnhaqS2Nl2jYH720Njtr2V5W+JiKsviLJcvMcLxZA0EQphUplOKQlTTQsqvlQUwKD4acvPNHzx0wIQFEDb9j+zoOWn94kQSouYAPHvt+uOZm2oIBLSBCiEJDm/tltgBdYlgNr6M01f/AD1mdhOuQAGrb+I+E1I9MYjSoFriJ6rljAKQ/jvoSCSaQSFkCiSIobNJk4Eqwmgqf5VRald12LUAtTiZ9u++AExeX7SE5uoKFUVyKibs5o4xW1z3x88RyVRDIlh8NwURNJCKITEZzqJYaOttymsCDMuATSC6QyuCNp2ZQGsIhj6vHEiFHJPjuNavHFVI7h9HKJYdk6KIWkRFEISZjqJYaNopQgGc/9MEyJl+QIYDYCJukErCWCV8meTnQ+MZTzl1vT6GBdqXeki1Y49G7+bZFpAIUwq1cp/RfPpkhxYU4tKIhgUQMOAsCwgZQGGaw0iMsfniqC0XQvQtQrVjk5NK61uAaxlddVb67RWge6YMQcFsaarlmJIphAKIalMpS4LM+lB1CxrMC4qtJoVmEpDWK4bVMpyEdTiZ9tKDDWGAOC+L2AZ1hS+iboZq8wBVqSSJRcn4BiHe5ZiSKYICiGZvVS0eCeb01ZdBENzgekURCql3KFuzl+sCJZKYSswWsmlmgDWcx31CEel+xU9fiVhjLuvFcSwLiiGZIqgEJLxM5MfQpNtc1RJBONcoZkMRDrl3y8tghrbgSyVgGLRE8dYtyliBLDiHN04PpeJWsu15iujgkgxJNMcCmGSmWi39GbTiDHFunWnSATTKYi2NsBy/3sFcv+EEMr1WSpBFgpKCPV7zUCNT1cEawpgvUIQV64u+N7o9nqb9MYEx8QeYzIwopQ0mRa0wfbZtGkT3vve96KzsxOLFi3Chz/8YfzpT38K7bN27VoIIULLKaec0qIRE49WCOhEzymM+h7K1farRwRNEyKTgchmgVTAEtQ1QKuJIAJWoG2Xi6B0Ii5HWV0ERSAxP1RH1YAwTYhUGiKTgdHWpsacTqvFSqnt0UIAweuvFEU6AbEShvCWGjtWWD8Nf8iRGUdLLcInnngCV111Fd773veiVCrhhhtuwJo1a/DCCy+go6PD2+/MM8/E5s2bvdfpdLoVwyUzhciDv773VJg3rFQqLRgtqkUwnYZoy4QtQUDNFzoOZMkGigXIQhFwHH8uUZdJq2UF1rL+ooJnuONy5yq9ces6o6YJ2IFSae55pW0ry1bPaeq5S8MKjLGClVbNQowhmPw/rlzH6HXTRUomQUuF8Gc/+1no9ebNm7Fo0SJs27YN73//+731mUwGvb29dR83n88jn897r4eGhiY/WFJOKx5Atc5ZxUKoWNUkvFPV93v7RC3BtozrDg1EeAbGGxVBr7OE46ZMxM0FjkcAtaXnJeobrhtWlWVzRkbdv8dfTk24pd+EZUGWShBmGkKLt22rMcdFjAaOUY1oJZyK6RXVgpwohmQStNQ1GmVwcBAAMH/+/ND6LVu2YNGiRTj++ONx+eWXY9++fVWPs2nTJuRyOW9ZvHhx08ZMWoB2+wXdf3p9aL8arY8mdO5xiiDgBsS4Igj4VqBOlZiICLrXLiwLRiYDI9sGY24OoqNdidPoKOzDh+GMjMA5cgSyVBy/69J1d8pSEbJQgDM6Bll03bpSKuG1LIiUBWGl4l3LEZdp0BUadIlKR4aWCUE3KZkgQpZ14GwNUkp86EMfwsGDB/HLX/7SW//AAw9gzpw5WLJkCXbu3IkvfelLKJVK2LZtGzKZTOyx4izCxYsXYzU+BEukmn4tM4qJPDzKogKb2CdwwseIF8GJPmRjUyS0O7Q9q5Llo8EntgNZKAD5PGSx5AcnBd2N0YLZdQmg4YqP5VlrciyvUjDcecf6L2zi8676ngjTCIi747d4inEzxxYhn0hVnEZEzZIZRUkWsQU/wuDgILq6uhp67GkTNbpu3Tr8/ve/x1NPPRVaf+GFF3p/L1u2DCeddBKWLFmChx56COeff37ssTKZTEWRJBHGEzka/bWvo/ka6ZZqogg2hEoi6J3QtVIdFRSDQlGJoN5WoWyaJ4I13L5eUEs6BRSV8Dmjo9WFE6gc4BJ6HVNJqBLSgbQBCCXqej4ShoBwCwBI29/Xe49jlH8m2gVrxMyPVqLinC5dpGT8TAshvPrqq/HjH/8YTz75JI466qiq+/b19WHJkiV48cUXp2h0pCrTOLR9MiIYnKcKRoj6KRKuCAaqxeh+gkoEbT8yVK+Ldo9wqSqCrigJ04TR3g7RnoUcG4McOQInnw+/xzDHec1m6FV5lGodwuhZcwbglML3KiiI1cTQ+9uEgF1fzmS1uUeKIRknLRVCKSWuvvpqPPjgg9iyZQuWLl1a8z0HDhzAnj170NfXNwUjTAj1WoXVEqMb8fBpoDUY2wGigmBXEw9vW0gEU250qOmPWTfWBZRLNJ8HtCVYpYWS70asZNEZMLJtEHM6AEfCGTkCZ+QI4NjeHGGlOqR1J+FXChoK3hYtjDWqz0jHAGBDSNdCdNRcImD6wTqeGIa/T0K3oyqVyoOaagXjTNMfY2Rm0FIhvOqqq/D9738fP/rRj9DZ2YmBgQEAQC6XQzabxfDwMDZs2IALLrgAfX19eOWVV/CFL3wB3d3d+MhHPtLKoZMg9bpIo8W8mzGOaq9j31KH+AojVDFGZLN+igQACAFpKiEUbjAJdL1QJ8bSC4pgFSsQwoA5pwOiox1ydBTO8IhrYap8RY9o2kVcCoj+EVMt5aHa9ZuudedIJcLBccZdm2NAIJCm4QqiZx0GxdBxA5AcR7mcgXIxrIS+rui10Sok46ClQvitb30LALB69erQ+s2bN2Pt2rUwTRM7duzAfffdh0OHDqGvrw+nnnoqHnjgAXR2drZgxLOYCVSZCQeh1HCRVovubDBlBbArPFDr6ajuHUOLYHtABF1x1CIIx1FWYKmkHvQ66T2YJ+hSSwSNTAZG7yJ1PMOAPTikBCkogHaMG1ENPOb6zMi/PhUFFHDH6HjzfSo4xwq7cysIovee4JhME5DCsw5DYmi6bmCdJxwUw2ou+MmUcCME08A1Wo1sNouf//znUzQaMl4xjCZDK7dYjYCPZhF8+Ee6QMSJYawIxj1stTVomkC2DcikQ+kRUkdMSglRUsWzVT9BRwmjFkGNdCqLoGEC0lFzgX2LgHwRpb1vQBgChiuAody96LVE3MCx+X3ViHMjR6x4z5XrWomx85uR75C0bUDI0GcjrJQvhjYgTDWXKAAvJQVjGL9lSKuQTIBpESxDZhD1/PpuxgOogkhVHE8Fa69SxGJFtDWYTvmWiuMAKQvSUonrMARQcq1BPQ8WdYlq92U1S1A6MDs7gbf0AH99E6UDb8LIZn0xDQigdx2uy9Jr4yQdX2S0y1C7MmOvzwxsr1BVR7u9NVIC0oaUfiJ/6LriflAFokyFIbxgGu0qlbbtiSGkhDCteMuwGpwnJBOE/gQSZpwCVmaZ6IdnI92f3jFj6l1GUyUiCe+i2jjq6LMnhFAimM0CpgGUVJCKJ4JutwjPGrR1pZhAjdHgHF4lERQGzFwX0LcIGNgP5/AwjPZ2PzevWPIscGGaSgCF4aVmeCkMUkIWC+757eoiCIS3x81lxlWI0Ski7jk80TUCwUOVaqHq+6BFO/Cd8YJpikV1z9waqJ7Y0/1JmgS/WaScWmJY5Zd308TQP0F9QTC1zl0rr04Lqmn6VWO0uzNlqYezZUIahlpfCiTIO/5xtFu0optSi2DXHOAtvZCvvwFZKCixC/QpDAkgEOhib6sIVSnHn0wfhxawkHVXwaoMCqJ7bV7kZ/B4ZecIBPboQJyAGIYKA6SUZShSE3ResdoMqQO6Rsmk8eYKQ4nRgUjSKcYrZq2xK1hF0ea3QVerKzyiLaPmBQEvqlGmLCWCplDeRJ08Xyr51iDgBcj4aQMRURDqHMacDmUJ7ntT1QeV0q1EExCXiEu3oou1kXgtoSJWY1kpOy2Gyl1aNndYyVWqUy3cex10k2IsD5FJA5YFoUXTkeGcREIaBC1CEs9ErMKAG63hFV1qEHWLqpUVxlA1GdtPnBeWpaxBt1waHEdZhloEhYCQEsJ2O0vYdsAd6ieQ1xTBnoUQhw6rRPlSKSSCulqL1xmiWPJdn1NNPW2fZKDSjAjkVlawDNU/SjSF2ykDcIVeW4WmG7GbTk3se0WrkNSAFiGpzDiiSKMdBJTLKxBlOJUEk9ujRCJKvV6ABtxWQ34JNWTcpHk932cYkCkL0jJC3SN0ykQwOjRkDZaNQXiVYrCoW80Djo4qd6cWwYBV6ozl48W0VVSrOKPvpztvWBZIU3asQN6hm2AvdIZHMEDGUOXlKna7IGQS0CIk1an28K2VnF1vQ9wJENfBIIRhhN2jejxAmQj674mUUMukw9agYahoUdO1CAFVSq0Y6ezgOGERjLmHIpsFFvehtLATolBURbO1COp2SggEkKgBj+seNZ1azYEdJVpl84ZA+Luh5we1wLldLQD3MwoWO2/LeFYjIY2CQkgaRlyIezAQolHU7mZe2z3qzSO6YhmKMDVNd37KDZDRgmaZyhoMWISiGJgb1I1sdfBI3A8FQ7n/jFwXnDYL1qFRJYJuJRphWZ5FqlseeUwXizBKtYLfbmBMxe9BZJ33HTJNv5i5dg070m97RTEkDYSuUVKbWjUm42pdRmpI1hXkUK0EWB2pDqGX2kUX3S9g9fnrwqkBIp0KB8i4+8iUBZnyrUEx6gD5QjgYR1uDXrWW8qLYxtwc7IVzYe4fgvPmIcgx1TJMtzOStg0UIlbmdBVBTY3viC60Ha1AFCoGoKNJQ3VGBQQAaQj3Pptu541i2EU6gcpIhGhoEZL6qPNBXGYV6uTpapZh1E02WQtSJ7lHOxxEA2liFqHTJXRaRFAI0ylI04BjudZgoajconHFrePaA0G5RO2/eQucNguyI6ty5uCKYCqlOlYUS6qRrne8aS6CmrKAoPJKN+Ho4khFHFP9UJC2DjySvoWtcZTgibZMZE6aIkgmDi1CUj9V6krWsgqjlUQ8GuE2rXYMw/BSEoCY1Ap3H49MxnfJBUQQlqmswZQJmAIiL5UQBh/U7t9x1qB25Rlzc3AcCWvnXtW13nZ8ERzLh4NrZooABgl+R4LfAe/HgRl2a2q3p2upS9NtxWQ7QKGgChboOUZHArDVeywLIpOBHB2rkJ/JcmukfmgRkvFRI/pP/Rko4hyN/NOW4ThdnePCccrO4c0B6ohSLX6G6rDup0tkVNWYYGUYANIyITMmpCUgDQFRVLluZW7RYHCLd92uCHZ2YvRve2DtPQgYAs7IqBJB03QryMxwEdTEVs8JBMRAWcBq8QOD1GehigboxHoUiuqe6ve7peYAqICmSp4GRpWScUAhJONnPGIIxIshUJ8YTtRVql2yUQsw8toTIh0gk7L8SNFAmoDMpOGkTDiWASmUEKoAmUr3IpBi4FpG+f/vWBTmWnAODcI5NOjXMTVNSD0nKIyZLYIaGfgxpAlahzoy17LCNVy1ezplqVSJQkHlEzq2X1XHkerzMYTq1RisZkTIBOC3h0yMeh/WsZGTomLU37gTpmOtAenNJcEQvvi5ixcdamihFUA6BbRlVNWYiDWohNCEkzLgWAKGLqnmtloSbgumkFs0MD5hmjAXzEOpw0Tu/+1S691ISmFZcPJ5NSdYq0D2TCNODKF+JHnWb8oCUhaMTMZvWwX4kbNSehG1/mHdIBmd7+nOLdaE84ikAhRCMnGqVAtRf1ZwkQI1xbDag62uTgTqQG5iuuEHw2hxBPy/TRMilYJMp9S4nIgQmgactKXcoqYAbKnSJmwnnItYJoLCy4s7fMoSOCkBZ+gwUCwqt2A6pSrJFAqVI2ZnC7pcnP5s3XlBGKZa0in1fXD8snzevS0WVdECfX/cyFzPKsxkmpqzSmY//OaQyVGvGAYQAUESVX7N1xLE6LkAuPlmgahNM9BdPtC1ISSIhgHZFnCLuiLoWSiGAZk24Jju/KAjVdqEPocRKDwdTA1wS4aZ/b1wLIGup1+BSFnquIbhlksrze4HeBXPgSwEImaFAZFt898TCWjygp3MwFyh/oHlulLrcpHSKiQxzOL/gWTKqEcMA1ZhOIJTPbxiRS+203o9wwnkogk1B+iVTTOD84WG715Lp1SyfExtTJmy1NygKaCS2hBxnYYDQdSDXBWeFuk0Rt/ajWK7u89YXv0Q0I12veoz06iEWiMJRmsGu07o4tp5t5iAjgRNp8P3NnQs9/ugP0/3OJCOqkTDJHsyQSiEpDHUEEPvdbRTfEAMY3FdXhMSQ8DPHzSEEsFAeyXvvKahUiOsyBi0cKcsOJaAk3ItyGBpMX3syJyi7qhg9CxEqcNE9y92oXT0IjWnmE6rcxaK6rpmqwgC5Sk3oVJ0Ek4+73XuEEKo9JVARK86hPSjcbVVH6w2o9MvtHeBViEZJ8wjJI0jrhizGwnpJVNLB3AMSMOfAxJCuEWvq8yTCbcoc73oxq/CgnTTFzyLQQr1E1AIL2pUuoE0UgiIaNFu1xqUwfzwiHBFX6sgGQOlnhyMkoRzaBBi/wF1zemUVzmm7vnOmU5AEMu+C6OjKm0FcK3ClOc2BaAq9cAGbBNCz6U6bv6mrgFrS18knQb0ZSSJghYhaTxxnc5RHjwTFA8REKWyX/RGwGqsmjwfqR6j63eahnrAWpYrcobbTd0VQdOAkBJSADAFpB6HG3EqhfDdogBkqJape6yAoAWT54cXZ9H2xqh/T1IpwDThjI5Vvl+zjUpWodeM11HuUV1NxrLKA2ZsWwXNuK7kcLUZvwKQ18rJS9UIlmsLVjCiVUh8KISkOUxADIEarlKjyjYg7HYNBsy4yfIwDL+zg+EKr9tfEACEI5VlaJmhIBsAcKxIubbouGIsWdk1ByO9Joxdb/gPaiGAYDWU2RwpGiTSqSJUJ1TXEXWkSh/RXesB3z3qSBVha9t+gXPtjtZtnnTEafTzme0/NMikoRCS5lFLDIGwGAbz/ILWna5GEkx/qGEd6vwzHaGpxQ3aOjS01Wf4nSRKDmDALaVm+SkXjgNpKo8qJFTkqJeLGBPgIwyItgyKfV1oe9NRDXd1vqGbMuHdh6Q+pAMiGLqXgVZMAABHlaBTBbuld+/UMQL3z/HnWYMNfuN7JtIqJGEohKS5VBPDstQHv7Zn6OHoik1ZtGmtwAhd8iySUK+sQMttqxToLeg4ELYrnClLzeWFokxVR3rdnV5HnQZzCUXKgjANGAvmoZS1MO/5IYz+n7d51wXT9FsrJcUa1ESExyuurXHcoBjbCYsh4OdZ2rb/PQmWwNMWoXbDRr8bMub7RogLhZA0n2gD16gYBq3CoBgGw+SBsr9FnBssiHaZ6e4SgJr/8yxE1y1qGIFGu+7YTFN1m3ALcBtF6QXLSB116l6btlJ0ori+lmKHAbF3P9r2u/OBpgnk8/HzV0khUKAccF2erqtT2nbYTRpA1yWVWiy1pe8Gy3iF1V0xFKJyDipLspEo/DaQqaNeMQziPuSCHSTKyqUFUySiuA9X6bX8EcoKdBdPBN1IUeFFIkJ1o2/LQBoGrBEbjqkCZaQByM5235UXLAHmuvqGl/Xi4PEmYNswXn5NbZrT4XWhT0y0aJAy74BfhEAGGiB7HT1ckfM+90BwjS58EPudAQLl84zQDyEAybz3pCoUQjK11COG0YTqaI3Q4PYYyzF0OtvxhEqaOoHe8GpcSrfbvOfuBCBsCWGrwBmZMSFTBsyiA2m4uYSGgN3ZBt1VQqeICMtSY+nshFGS6HrFAUolOEPDardi0bVQE/zfLuZzUgEwthLAYsl1j4bdpmquN5hEH6g+E+ciBeJ/GNE9SmJI8P9I0jKiYqgtpEpiqH/Rx1UO0dsClmMZpZIqixZowCsNw+sv6FgGpKFSJKSuhFKyIWwbMAzYHSlIQyB1xIHjZt7a7SnVKV272YyASOfmwEkJ5F4aAfp7AkEggeviA1mh3eaOrUrOFYrhllS6rqjjuztlpBZsJauwmntU7WDECjNJHhRC0hpi5g0rimFwPx3pGUULZUQQhRB+uL22BoVQ+YKWoRadO2gY/v8IKUMd6kvtJqwxB3ZaeDVHZTajjmVZMOZ0qPcZBoQj8doqQ7Vq2rPXt2x0GbDZXElmMjhu9RjD8L4DftNj1w3qVhnS7tFg4Xb9mYe+O3HRxfwRQiJQCElrqVcMg4SSsyPzhsEIUZ0UXyqp7gWAJ37KKlQCKC1DtVhKBY7hzkGJog1RkoADlNoMGDZgtxmAAJxcu1cz1AvucBzkj54PoyBgDI5AjuX9sZZKFMAaSB0V6kWIBqI9dccKwP8xIcIBSiEMN1+0VnAMrcLEQyEkraceMSwrYSbiH2BaALW14Ibdy1LJD77QYmgIzyXqpAScjOnXHHW7TIiSA+FIGLaEWXCjQdvVe+xsSlmFumB0RztEWwaHF2fgpNR1ibaMd11e7VIGayjiPj+d/+nIUIpNKN1GROYGI+k1weMLUUEEo8ciiYbfAjI9iNburEcMgzVBA5VbhBbBYGFt24YouVGJpgEnpUqrKbeoK2wZE3bG9ANoDOG6R9UclVGUSA07gADG5qmEe9meUZ3WM2n179wcAGDuHwE5POK3dIrLmyTxRDqXyOAPCG0VCsNrkxU7Lxwp3E1INfgtIdOPatGk0RqTwUR5BKwC10oLldwqKbdb0B2q0yqkIWBnBOw2E05bCk5az+dJCNuBkACkhDlmo+2AjVKbwOiiNEq5rOpl6KZi2PM6UegUKHUIyP6FkIWi/yAvFAOWKt1xlQiV4SvfWL4uInaVIoirHoufR6KhEJLpSS0xjKlRGnoA6khOXRTb7X4uSsqic1IG7IwBJ+XmGApVT9TOmkoMsxZkm6osI4rKPQoAMASsMRsdb5Q8MXTmdqjE+5SFfE87hCPhpICRY7tUPVPdsLdaWH8SmIhA6XsXjC4OzhVWO26lOeaycfExmHTYholMe1TbHsdv36QtwUjSdajaTPBBaQRaKxV1KgXc1kq64gwgDaDUJgDo/ED1l7BVwIwwpZdiYY2UkJUS+bkW8gvakAFgjORxZJEFawxIH3Zgjjp+l3shIEsleKMSBiAD5cWIj1SfT/w2B26bEPW6mrjWI4JuDqh3LAYzJRL+FCLTlwpd7lXbnioiEmmN5HWRMJV4iqI7V+gmyNsZPyXCSQmUsgbsNgN21oLdkYZMmSpoRgJCAtISgAFYIyVkDpVgZwQK8zIozc1ibL6AnQba/1rCwePTMHRUKeAljdfMbyNV51M9qzBQbq3MHepWnZF2heNQ8EiAGSOEd911F5YuXYq2tjYsX74cv/zlL1s9JNIKQh3OnZD7KxQ1qIUmIjjCkRB5WxXXFsoKdFJC5QdagJ0SKLUJFNsN2FkDdpsFJ2O5XemVMSIF4LgpGKnhEqxRB4VOE8XOFIodQH6+EtTDxzqqcLcanDtmtzapFuukEBWqWm7SOivwxKVN+CXaagQoBXIU2ZEi2cwIIXzggQewfv163HDDDdi+fTve97734ayzzsLu3btbPTTSKMb7AIrOGYbKrsU8dIMFu4slmGMlJWyGK2wpJV7SBOwMUGw3lBi2uWJoGTBsxyvFpucTAcA8UoJhS5Q6TOQXOMjPlTh8lIVFvwHkkVH/3FJClop+J4WJXHeSqCpi5dv8Eny+CNadqsKI3kQzI4Tw9ttvx2WXXYZPfOITePvb345vfvObWLx4Mb71rW/F7p/P5zE0NBRayDRmPGKgf7kHc8fixNC2/dy9QNK9NNQ8kHGkAOuIrVydBuCYgGMBTkr9W2pX3SOKHUoMZcrtWWhL2BnXjeqmXEjLgFFUFiMW5iGPHsWRPoHiHAEj1wUz1xUu+mzb7rxnQkRwsmJfhxtTBixutSIScDOeY/PHSeKY9kJYKBSwbds2rFmzJrR+zZo1ePrpp2Pfs2nTJuRyOW9ZvHjxVAyVNJpKrjG3w0NoXiiaWqGDVGQgYMXxrTBhS6QOF2CNOBAOAoLoJ9gX5wgUOg0U55goZU04pgGz4MAoSTgZfz+7zYSddsVQAP/3+OdRyDkotgs4C+YCAMzOTphdXd6YPGuF9S4nT5xrHH4aBgsYkFpMeyHcv38/bNtGT09PaH1PTw8GBgZi33P99ddjcHDQW/bs2TMVQyUToZIIRETQK5OlO5pHO1K4eHNDwZ52gYavouQ3dhWFEtKDBaSHHZgFCaME1XjX0HOHQLFDoDBHoNBpotRuwjEFjLwDYUMF2aQM5SbNCNhtAs5wCr3pIQgbGO11g3U6Orx5K6OzE8acORCWnjuc5S7SRlxXjWN4n7m3ImIB1nJ70ipMPDMmfSLugVcpcTaTySATiNYjM4w4EQxiBAIbdHPWALJUcsXSAWwAJiClUN0khIAQtuoyIVIwCjbSh4qQRsqfMJTa/QlIVwylUBVprJSANerAKOjEfEAKFWgDKWCOCLyenwtzTKB4dB7OnDSMQwKirQ3CNOCMHFFDb8tAOKr0m+5ROOvC95spJnGFtB0j8LIB93G2fR6kItNeCLu7u2GaZpn1t2/fvjIrkcwwJjA36P34Cf4bbM+k3WQAIA0IE0oM3XXCcQBHeOXWAJWRlh4SELYFwzZQKgkYJakiSVOAtAA76yfdSxNKDPMSQgKldhU8I2wJoyBgGTbsNonjj3oDdmYRZEfWy28UUqoAGkCF/Zum6rxuqzZEEMCM6E6hRaKSWFSz9OsNTNHHj74fCBxDpclIR6rPuuwYdba9ijsXSQzT3jWaTqexfPlyPPLII6H1jzzyCFauXNmiUZGmUc0aDPb8i+lW7yXZB+cEdQh90GXqdkMXJRuiUIQxWoR1uID04SJSwzbSIw7SwxLpwxKpYQlzVHWfcFJ+RGlhjgknowpzm2PKVSpNwBwVeOhPy2AUBP68pwdG3s1ZbEur6xFC1SXV5d/0dQoBI52CSFl+CbZKhcVbSXBMExHB4L+TEPtYi09XoJnwQekiTSrT3iIEgGuuuQaXXHIJTjrpJKxYsQLf/va3sXv3blxxxRWtHhqZDPrBU+Fh44lgQBzjXOLB1zKQquDt7wDaLJQAhPCtQeGuM6QEbAlRkiiVLBgZA2ZKwMmLUESprkDjpICidtFKwLDV+drfkMi/u4iR+WnMea4NojTkC7NpAm0ZiKIJWSiqcwe6YnguXtOECFg+ZcUDWmEtVugUUdd+Ze+boFhFfiRJR0KgvLBC1Vql3rHo9iQ+M0IIL7zwQhw4cAAbN27E3r17sWzZMvz0pz/FkiVLWj000gjqdUu5NSaVuPniIQNBM8I0K3YshyMB4UA6IlTBS+hAFilVObViCkZ7yq1FakCa8NyijqXmCwFAOK6r1FAJ+ZCAOQbMbR/FSKYNbfvN8LXpvy1Tjdd066HatlcBRbgVUYII0ww/3CNi31TG5b6usm/ZnF4dY4/Wjg0VyXaLrEcsw0nPDca6Yymas50ZIYQAcOWVV+LKK69s9TBIs/Dmm5yYh2agd1y04HIMXr3R0MqAO86RkHDCYlhyrUfbgWlLiKIDo82CnbXgpJQ70DClG2SjLEPooNQUAAjYGcAoAa+90o3FS/+Kwx19XnNYz3DRD1TTAERKuUhtGygWgVIJ0naUyLvl2OAoARCG8NMtPKLRkeMUl2rvqyZq491fd5WvtxVV7BhjIkFjLMTYfccL5wsTx4wRQpIg3IdctaTzilZh1J3qvRbl67TLUltjgCuSylUqHAfCduCkTciUbtvkDqCo3i/d5r5GSboFu4GOVyzswULkDFWKTaSUBShLbk9EPU+pH7iGoQQxlVKtokolFflq28qClQLQhccReOBHA0cqBZYE9wnfxMj+gXsX3b9SkEstwXBFsC7qFZ/IWDyRbSa0Cmc1FEIyfQi5EX0xBNz5Qm0FRK3CUI/CwLaoCOocRH0O3ZEi6rp0xVEIAUPYEI6ELAgIVxCl4UaJSnU+4UgIR0A4yiJMH5KYc/SbOLxrEex2C0ah5Kd6mAZguy2hHDegRzcMdvxgD2FZqk5pILhHXV+0J6OEyg+JFCgP3dcKFlLU9Rjcv5JVXq+1pKN8dQWgqDU4UYtLSmhLWJhm+bWNR/BrnYdWYWKgEJLpS8Al6gui/5CuNldYk0qdy8uiTgVgCxhu9wlZciANAZEx1b8AHBswikoM7QxgFIF9L3YjkwIOH5XG/DfHwlGuQh1TGhKiBO9/oc5zlKb0o11NUwmx208RpqkKgAPuNSOcTF4PAYtKRagaflCOMNQ6xIgXUJcFWHZL63WJViPOGgu6nNWJwvtX61NIMSQBKIRkehH38AlYKMHehBUDZxyoZHp3H01NkXQftsKdQwTg9jFUVpgoQuUgCgFRUgn1gAqYsWCgMAfIzxNID0m07Tcw2leC/Y4xzH3RgjniAAi4Yl23qye6Oh/SvR4Jx3/IG8rahA4EilqFaqBVLD8j8Gcg71IXprZtCEOEqt2MW7xqFUGYCDVckbJYqn6eRrsy6R6dtVAIyfSjDjEEbAipXJ7SQMU+rh5BN2mwMo0m6FJ1Uy6kcEXH1vs5EDAAqPEJd27RKFowijZSR9RArDGJtr8C4t0jkL/LARhRAi1tQLrzmFLNLQpbeueEbSv3qBYhO5wX54mgd02G+hFg6vGV5+eJ4PXqYgPB4woDIqXE3evzOAEBrDiXO1FrsJLgGGbYdat/AIz3HPUKGq3CREAhJNOTSmIIxFqHErZn8YWsQgSsHzNwPD03aMSsE8LrUgFHQgiVsiAdATUx6KZcQAmROVaCkzaRHrZhFlXt0fa/OnCeyMEalZCmK1jBZ7UrUFIIQBpASlmiIuSaDfwNlVoBGei159i+aOoC01K6QT+BOb1ISoYwDcDwE/plqQSnUAzf4zjGafWNWwTrEaeAS1fadnyD5kZbbcHvIq3CWQmFkMw8KliHEq5lKFWag4CphMtQGyrWp61YCSWcqK/yDQObHV0fVLlHIVU7p0KHQPZNB3Net3HwOBNdWROpA4FIUXcJuUaD7aIC7s+yBrPe/GBgne7NCITdptpl7EXU+nOoslSCLBRj5gDj507LRK9CRKqO9B2XCNYrLMHPyXA/37iCA4SMEwohmb5UDXgIRzZ6D+BSSa3Q82BSKgsIUO5IuGLpuQvhFxoss0ADoqWtsqArU6hWTEZRWYJCAnbGQOawo1ykow6sIyZee18Kb/2+cn2KoJUXtfrcB7rnutTb9fUCYdEDfLHTBAOGIuIlbQeyUPCS96MpCO4f8fe7EjHRqnWlMkzGqhKGNz84ZS2W6CKd1VAIyfSmlhgCEevQfai7zW8BqJQFt7g1TDfIxTT9VIpgGoWh/lXzj8ILmlEHDoyj5HaygAlxpAAAsOdmYI45sNsM5LtMdLxeQG5XCfkFJmTagpEvhlydZbiRocI0vNxBtz4cvHI2Zsx8aDBHMnhf3LQLadt+gj6grrEsOR9h4awlMONN12igO9Fzi4ZaL7l/N1OsKIazFgohmf7UegDFWIcA/AR0G4DtQBolVwxNJTZubiEMM5xjaJpu5RfhuVwhJWC5VmYw77CocgRFyW3w687zpQ87cNIGUkM2svssDL+1C5078l6ATWxum5SedRcvdjH3IJDqIUu2P2+oBTAQYCOECIthJSqJYETgKha+Do6tUYy30k2zCOaaklkDhZDMDOoRQ6CqIOqgGrj5eirB3YQQJZWa4FqDKNnKnep2gRCOmmeTWqj0gzAoXEUb5mgJTioFoyAhTcBJCaQHbcz7Yx6vrUqj848WYBfccWoLNHANwTnJYNpHdP7QTb9QBbv9oBnPLezEi23FGqzAxMUv+CMkOo/ZaHSptuCcYFlC/RQIFEVw1kEhJLOLCoLoVaZxm7dKww12cXMPYduBoBIB6ZhKcAwBYbsBJ5apBNEylWXnOIBwXa1Swjych91uQUqglFJl18a602j7ax6dr6Rhz8nA1O8DYq0LLXqi5PiiZ9vK2tPzhrYrfNriCxJIr6gqfEC8+MW4NcMuyDi3Z5MFMHgekYLfYDImn5KQCUAhJDOHSczRxFuH0uv+IFxB0ukMwpGAqbZLdw5O2DZgWSpwRgjlFnXFU8KBGMvDGkkjPy8Da0Q9rI08YIyWMO9Poxhe0o7c70chiiXPqvOiQ90AGamFL5g4H9ORoibjmOOL3qNa+8Ufr4mCFFcDVVuikhGjZPJQCEmiCFmHQMhCjAqiSsVwA1EMCelaIAIAUm7XCMe1GgvqeOaBYRhzUjDzDpyUoXIMsxZSA4PoGs4CpgF5eFilLuh5vJiuCZWiISdasaXu+qNVD1IhIGWKrLKyKNHJlm0jxIVCSBJFzSTwQMk25Xo0IIT0y52ZbkkyOwWZL6iHcaHgVzwZG0OblDj8rkXo2DWsgmhKDkS+CAwOw35rH8wDFuTwCGSxFDl57Qd7ywygVolgICK2rM0S3aKkQYwzaYiQmUtVEQxYGnFuSSnVnJwslVSbpEJBWXe2X9UFpZISx8MjyO4dg0y5UabFkhtYo0qqifZsE68yhmDeYtxSz3vjolanSgTdfozqnLQCSeOhEJLEUHMOzPET2KNi6K5UwqjbIgk3uEZHMVqugyWfR+rVAxg6tl29Nt32S7lOGK8MAKNjMLq6vC4PTaOW0FVKyai0jzACEaJTlcjuqGLgzb5XJNFQCEmikI6Mny+LEcSKYqjn9UolNw9ReMEbwhVDWShg3vYDOPTOnDpWOqXeOz8H2dkB5PMQKStspY63qkvsBVaw9AJ1VOsSwCg6UMWJiVRtBsEiAdL94dGsHEWSeCiEZOYwnod3DReaFsSy4IugKEbFUIhABwfbtwQB1cg3k1ZWoSOV0B0cQu7Pwzhy7DwVXVooQhwZgxg+AtHRDmFZEJlMuRhORBArCcNkK6Ho8UxFqkLNjiOENAcKIZn+TMSCAeqeT6pqJQbFEPD6HoYKYgMqzaJY9BPxAchSCcauN2AUHZTmd/g9Ex0HcnQUyGQg0ul44RuvGAa7I0StvomKYdAimwoRjAbkuAn03hg0tAZJg6EQkunNZC2acQRX1LQ6gnmMjp/v580ZunmGsKyQkLU9+wqO9Leh1JNTlWoMA1i4AKKzAzKfhzCN8kCeVgaFhOYCmyyCdfzIYXcJ0mwohGT2MxkxDFqFGl1uTSfEA6HIUdGWUX8Hukh0bXkRI2/JonTUArVtaBhy6LB/TNNsfWRk1C07FSLonStgDQoDIng/CGkyFEIyvZlObjAtekL4rY6ClV9KJfWvlBBtGZVqoYts2za6nngJxTkp2P0LvPeKbJsfEakf/o0ImtGiVs+xYgVwCsuXxeUo6rnBqDU4nb4PZNZAISTTn3offsGHf3SZLF4yt/uvFq9gvmHQRWoYEG1tSgzdjhYyX0Dbs68AhgGnfyFgmBC5Lhhdnf45gmI4mcW7dxWsy0r3RgvgVAXGRCzBoEuWATJkqqAQkplBXcnfk3MpCkNUdceVBc1Ec9scBygUPWEUpkqnkPmCqlFqqYoyxsuvAaaA/dY+tb9pwmhv91sjiZg5w/EQjX71L9Bb9LWG3LFTZQUGu3fo16FhVuhwT2uQNAkKIZlZ1BLESiJQgTJBCG+Mtya9nn5GKIoUUK48mS/4XeNNN1L0yChEJg3RloEzPAJj5+sQtoTdM9cTSqO93b8GYUCk6qyAGLzmcYifFy07VVZgHBFLsKIIEtJEKIRkZlKvhVhDGGNTJ4LvLzttJEjGNMvFMJ9XRbV1E2C3QLczcgRIpSHSaThDwzD+tAvGaBGyv1u1eAJUkIhleUnkoeoz0eupdF3VLD9EBbAFrYwqnK9qsBCtQdJEWHSbNJ8J5QDKcG5Zrf3qOmaMteRtkpFNMQ/lYC9DHcShG/UCbicKd5/Dh4FcF4RhQFoWUCiq+qS2rYQuZcE5cgTi5d0wFnYDmbRqCGzbkIUChGXBKRShq30rS8moKOjBaylL0HffM66i1XERnROl0ucYqB4TtFLLoAiSJkMhJNOTWsngk+hN6B+jsvtNOuXOEq+XoQ1IYbhtmvxxSLf8mo4olYNDEHNzyiq0LDhjebd8m+M2+lViJ/cOwGhvhzz2KBgDB9R7SyXPGpS27QuEFrZqwhe5xjIB1OMuu8AGpytUO54QEKbpt8UCwi5RUUX0CWkwdI2SmUmjH9pRYtyPobJsuomudonqOUFHJdjrjhXOoUEljCkLRlsGslSELJW8FkwiZUE6EvbwCOQf/gI5PwcxLweRbVNVagwBc04HjGyb2t8QEFbKd3maZnkEqNRNfu2wCxQodylXqz06mXnDYEBMbN1TI2T9VRRBWoNkCqBFSEg9BB/SgGvJOJA2lGVoGK6F4/YYdKQqu2Y7cA4PQ8zpgEinYDgOnHweEH5pNpGyIItKHJ0/vwyjsxPy6F51updfVS5VQImh7l6PwNxhQETqalxbz4+IyQgg4IpZTEWYYDFtfapKlWMogmSKoEVIyHiIWIghyxBQVpzeT4thoQh5ZFSlULS59UX1Pi4iZXmuUGd4BPjTThiDIzAWLoDIdQEpt3uFaUKk08rKdGSoK0PIBTrVIhi0KqVUXSpi94nJWyw7fwsCeEiioUVIyHiJBNmUWYYIGEOeGBaUeKVTMLJtkKYbPONuF6YJWH5vQyklnIF9SvQWzINIpSDHxiBHx/z3BcZQNr5K0aRRJjoPV23uttK+0fnNSuemCJIphkJIZj6NqBwzodNWEAPTVIW13dJrQghVn7RUUu9Jpd25PaEqz2hciw9SQrjvlaNjwOsDMDrnAPNyMCxLzTEeGfXnGSOpEe7Kei+idhANUNuSrPY+r40TRZBMTyiEZGYxhaJXtbpLdBzBfQ3VuV5K6adZuMhiSVmOqbSqM5ovAMWiN0+mRROGocTQAOCooBsMHoaR64Roa3NFs6jOUygo96x0ysZcu6PGOF2oUSGrJFyG6QtfJMqVIkimGy2bI3zllVdw2WWXYenSpchms3jrW9+KG2+8EQU3MEAj3ALHweXuu+9u0ahJS5mkCAaTy+tZYs8frTajO9QHCViCwrLKcvJksQTk88ryS6eAVMqb85Olkor2DAmjb1U5hwZh/3U/5JjbvqmjXVmgbRnV5DebVed05xurXlOs+1QgOpfnRabq91SqRKOv14h00qAIkmlOyyzCP/7xj3AcB//+7/+Ov/mbv8Fzzz2Hyy+/HCMjI7j11ltD+27evBlnnnmm9zqXy031cMlkaETO3ySYdDufOgVYui5NrwGvaxkK0/RyDL0O9zrnMJ2GSKcgAQgU1dyinifU59Viq3MbXaGUpRKEW+jbmDcXzpsH/aCdgBiGsO3A/KK7PZBrKAKiq+c5vajOatafexzvXth2fYEx1Y5LyBTRMiE888wzQ+J27LHH4k9/+hO+9a1vlQnh3Llz0dvbO9VDJNONqHUxFYzzfKGkeikhXVenEEL9HdjHfeEl3MM0VAUat4OFZxVKV8i9jheBnMWxvPpzcAjSPa+Y0wHn8HAooMYfYEyaRXC9DehqNrUqzwgrBVkqhg+jr6uaFcg8QTLNmFbpE4ODg5g/f37Z+nXr1qG7uxvvfe97cffdd8Nxqke65fN5DA0NhRbSYhr1wJtAlGPVeqL1nK9abU9Hxi6qc73ttWWSbmsmmKbfy9BRblJZKEDmlaCJTAbC7VShm9N6ifzFkp8uEXHJSrdqjZQScngEcCREyoLhJuYrC7Lou12D1+PY5a7OGLen6pBh+pV0PDEP3JdabaAogmQaMm2CZf7yl7/gjjvuwG233RZa/+UvfxmnnXYastksfvGLX+Daa6/F/v378cUvfrHisTZt2oSbbrqp2UMm46VRLtJIcnv9byt/8E7IbVqlZqmHey4JJTDaQtS9CYVte4EwcCSkoxLqRcqCyGTUOQpF5ep0A2FCgu740aahwtyOdM/peEn6kI67j6nyHFMpdy5P5Th61qlpuqLpAEVl6ek5SwBwRkcr34dA66jY7hH1BtkQ0gKElI39Rm7YsKGmCG3duhUnnXSS9/r111/HqlWrsGrVKnznO9+p+t7bbrsNGzduxODgYMV98vk88u4vbAAYGhrC4sWLsRofgiVSdV4JaRoNr2nZeMfGpOcV/QP5fxvCjwrVrlM9ZxcRd5FOQbRlgGJJpVi4ATje/GGtjhnjvScTziesUwCD56AIkglQkkVswY8wODiIrq6uhh674RbhunXrcNFFF1Xd55hjjvH+fv3113HqqadixYoV+Pa3v13z+KeccgqGhobwxhtvoKenJ3afTCaDTCYzrnGTKaTRwTN1dGQY/yErdKMY/4H8vx0jbCG6Zdkg/MAUb67ObeUk0imVZlEsAW63exWRKkNBLmXnmozo1PHZiEggTmyZNFqBZIbQcCHs7u5Gd3d3Xfu+9tprOPXUU7F8+XJs3rwZhlH7wbV9+3a0tbVh7ty5kxwpaSkVupM39hwTtI5iDzW+h3iscOrxOIZq5RTodi8Mw7UOg3Noji+IKUulWjjSFUQHUgoliKggRMDk7m9wPi9g+YUviVYgmfm0bI7w9ddfx+rVq3H00Ufj1ltvxV//+ldvm44Q/clPfoKBgQGsWLEC2WwWjz/+OG644QZ88pOfpMU3W4h2QmjKORoniPWfssp8pHS8Vk4whJdyIUwTMFHuLpWOCqgpFPxao5blB+JICQE1BxgVpUrRofXcCzXe8hSMih3kKYBkhtIyIXz44Yfx0ksv4aWXXsJRRx0V2qanLVOpFO666y5cc801cBwHxx57LDZu3IirrrqqFUMmzabZVmIkuGOqCeXq6fFoQXSvXecdaouvrLOEE27hJHQUqpTuPGL43mmL0b/2iLBVy/WLNvWtxwVdT7k2QqYZDQ+WmY4MDQ0hl8sxWGYmMlWJ+NOhXmkggV4ErltK6QlibJNduHN2uppMzD0L/Tev5uYNJsXXK37R8cz+RwppATMqWIaQhjJVVWmaEHBT32l90RCGOwZtJWq8nEHTC5IBzNB7vZxFANI/YM0gn/r6F1a5BxRAMgugEJLpz1QE1lQ8d51pBXHFpcd9qqCQ+MEvXmUZfWydYmF6b4R2eUbnBCsG0cQRM/ZQaoQ/0MjAKYBkZkMhJDOHFtcsrUogsKXxh9YCZce7UkMFsqNjiqk3Wvd5ZUVXLMWPzCYohGRmMZ3FsFmEilkjIHwBq7EBBQDqcpNSAMkshEJIZh6zRQxrFLWu/L5ykQoL5HjHUcWKpfCRBEAhJDOTmSiG4xGVieRXNsItS+EjCYRCSGYuM0EMGyEsMZ0gGhJARNEjBACFkMx0WhlRGsdUiEu1dkmEkHFDISRkMlCICJnxUAjJ7CBOkBptJVL0CJmVUAjJ7IXCRQipg9YUWCSEEEKmCRRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiiaakQHnPMMRBChJbPf/7zoX12796Nc889Fx0dHeju7sanP/1pFAqFFo2YEELIbMNq9QA2btyIyy+/3Hs9Z84c72/btnH22Wdj4cKFeOqpp3DgwAFceumlkFLijjvuaMVwCSGEzDJaLoSdnZ3o7e2N3fbwww/jhRdewJ49e9Df3w8AuO2227B27VrcfPPN6OrqmsqhEkIImYW0fI7wa1/7GhYsWIB3v/vduPnmm0Nuz2eeeQbLli3zRBAAzjjjDOTzeWzbtq3iMfP5PIaGhkILIYQQEkdLLcLPfOYzeM973oN58+bhN7/5Da6//nrs3LkT3/nOdwAAAwMD6OnpCb1n3rx5SKfTGBgYqHjcTZs24aabbmrq2AkhhMwOGm4RbtiwoSwAJrr89re/BQB89rOfxapVq/Cud70Ln/jEJ3D33XfjnnvuwYEDB7zjCSHKziGljF2vuf766zE4OOgte/bsafRlEkIImSU03CJct24dLrrooqr7HHPMMbHrTznlFADASy+9hAULFqC3txe//vWvQ/scPHgQxWKxzFIMkslkkMlkxjdwQgghiaThQtjd3Y3u7u4JvXf79u0AgL6+PgDAihUrcPPNN2Pv3r3euocffhiZTAbLly9vzIAJIYQkmpbNET7zzDP41a9+hVNPPRW5XA5bt27FZz/7WZx33nk4+uijAQBr1qzBO97xDlxyySX413/9V7z55pu47rrrcPnllzNilBBCSENomRBmMhk88MADuOmmm5DP57FkyRJcfvnl+NznPuftY5omHnroIVx55ZX4+7//e2SzWVx88cW49dZbWzVsQgghswwhpZStHkSzGRoaQi6Xw2p8CJZItXo4hBBCxklJFrEFP8Lg4GDDPYItzyMkhBBCWgmFkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiioRASQghJNBRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQtE8ItW7ZACBG7bN261dsvbvvdd9/dqmETQgiZZVitOvHKlSuxd+/e0LovfelLePTRR3HSSSeF1m/evBlnnnmm9zqXy03JGAkhhMx+WiaE6XQavb293utisYgf//jHWLduHYQQoX3nzp0b2pcQQghpFNNmjvDHP/4x9u/fj7Vr15ZtW7duHbq7u/He974Xd999NxzHqXqsfD6PoaGh0EIIIYTE0TKLMMo999yDM844A4sXLw6t//KXv4zTTjsN2WwWv/jFL3Dttddi//79+OIXv1jxWJs2bcJNN93U7CETQgiZBQgppWzkATds2FBThLZu3RqaB3z11VexZMkS/M///A8uuOCCqu+97bbbsHHjRgwODlbcJ5/PI5/Pe6+HhoawePFirMaHYIlUnVdCCCFkulCSRWzBjzA4OIiurq6GHrvhFuG6detw0UUXVd3nmGOOCb3evHkzFixYgPPOO6/m8U855RQMDQ3hjTfeQE9PT+w+mUwGmUym7jETQghJLg0Xwu7ubnR3d9e9v5QSmzdvxj/+4z8ilaptrW3fvh1tbW2YO3fuJEZJCCGEKFo+R/jYY49h586duOyyy8q2/eQnP8HAwABWrFiBbDaLxx9/HDfccAM++clP0uIjhBDSEFouhPfccw9WrlyJt7/97WXbUqkU7rrrLlxzzTVwHAfHHnssNm7ciKuuuqoFIyWEEDIbaXiwzHRkaGgIuVyOwTKEEDJDaWawzLTJIySEEEJaAYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiioRASQghJNBRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJhkJICCEk0VAICSGEJBoKISGEkERDISSEEJJoKISEEEISDYWQEEJIoqEQEkIISTQUQkIIIYmGQkgIISTRUAgJIYQkGgohIYSQREMhJIQQkmgohIQQQhINhZAQQkiioRASQghJNBRCQgghiaZpQnjzzTdj5cqVaG9vx9y5c2P32b17N84991x0dHSgu7sbn/70p1EoFEL77NixA6tWrUI2m8Vb3vIWbNy4EVLKZg2bEEJIwrCadeBCoYCPfvSjWLFiBe65556y7bZt4+yzz8bChQvx1FNP4cCBA7j00kshpcQdd9wBABgaGsLpp5+OU089FVu3bsWf//xnrF27Fh0dHbj22mubNXRCCCEJomlCeNNNNwEA7r333tjtDz/8MF544QXs2bMH/f39AIDbbrsNa9euxc0334yuri5873vfw9jYGO69915kMhksW7YMf/7zn3H77bfjmmuugRAi9tj5fB75fN57PTg4CAAooQjQmCSEkBlHCUUAaI5HUDaZzZs3y1wuV7b+S1/6knzXu94VWvfmm29KAPKxxx6TUkp5ySWXyPPOOy+0z+9+9zsJQL788ssVz3njjTdKKMnjwoULFy6zaPnLX/4yeWGK0DSLsBYDAwPo6ekJrZs3bx7S6TQGBga8fY455pjQPvo9AwMDWLp0aeyxr7/+elxzzTXe60OHDmHJkiXYvXs3crlcA6+CxDE0NITFixdjz5496OrqavVwEgHv+dTC+z31DA4O4uijj8b8+fMbfuxxCeGGDRs8l2cltm7dipNOOqmu48W5NqWUofXRfaRrFldyiwJAJpNBJpMpW5/L5filnUK6urp4v6cY3vOphfd76jGMxsd4jksI161bh4suuqjqPlELrhK9vb349a9/HVp38OBBFItFz+rr7e31rEPNvn37AKDMmiSEEEImwriEsLu7G93d3Q058YoVK3DzzTdj79696OvrA6ACaDKZDJYvX+7t84UvfAGFQgHpdNrbp7+/v27BJYQQQqrRtDzC3bt349lnn8Xu3bth2zaeffZZPPvssxgeHgYArFmzBu94xztwySWXYPv27fjFL36B6667Dpdffrnnarj44ouRyWSwdu1aPPfcc3jwwQdxyy23VI0YjSOTyeDGG2+MdZeSxsP7PfXwnk8tvN9TTzPvuZCyOdnpa9euxXe/+92y9Y8//jhWr14NQInllVdeicceewzZbBYXX3wxbr311tCF7tixA1dddRV+85vfYN68ebjiiivwL//yL+MSQkIIIaQSTRNCQgghZCbAWqOEEEISDYWQEEJIoqEQEkIISTQUQkIIIYlmVgkhWz+1nmOOOQZCiNDy+c9/PrRPPZ8BqZ+77roLS5cuRVtbG5YvX45f/vKXrR7SrGDDhg1l3+Xe3l5vu5QSGzZsQH9/P7LZLFavXo3nn3++hSOeeTz55JM499xz0d/fDyEEfvjDH4a213OP8/k8rr76anR3d6OjowPnnXceXn311XGNY1YJoW799E//9E+x23Xrp5GRETz11FO4//778b//+7+hlk669VN/fz+2bt2KO+64A7feeituv/32qbqMGc/GjRuxd+9eb/niF7/obavnMyD188ADD2D9+vW44YYbsH37drzvfe/DWWedhd27d7d6aLOCd77znaHv8o4dO7xtX//613H77bfjzjvvxNatW9Hb24vTTz8dhw8fbuGIZxYjIyM48cQTceedd8Zur+cer1+/Hg8++CDuv/9+PPXUUxgeHsY555wD27brH0jDy3hPAyp1vPjpT38qDcOQr732mrfuv//7v2Umk5GDg4NSSinvuusumcvl5NjYmLfPpk2bZH9/v3Qcp+ljn+ksWbJEfuMb36i4vZ7PgNTP3/3d38krrrgitO5tb3ub/PznP9+iEc0ebrzxRnniiSfGbnMcR/b29sqvfvWr3rqxsTGZy+Xk3XffPUUjnF0AkA8++KD3up57fOjQIZlKpeT999/v7fPaa69JwzDkz372s7rPPasswlo888wzWLZsmdf/EADOOOMM5PN5bNu2zdtn1apVoaT+M844A6+//jpeeeWVqR7yjORrX/saFixYgHe/+924+eabQ27Pej4DUh+FQgHbtm3DmjVrQuvXrFmDp59+ukWjml28+OKL6O/vx9KlS3HRRRfh5ZdfBgDs3LkTAwMDoXufyWSwatUq3vsGUc893rZtG4rFYmif/v5+LFu2bFyfQ8vaMLWCZrZ+IorPfOYzeM973oN58+bhN7/5Da6//nrs3LkT3/nOdwDU9xmQ+ti/fz9s2y67nz09PbyXDeDkk0/Gfffdh+OPPx5vvPEGvvKVr2DlypV4/vnnvfsbd+937drViuHOOuq5xwMDA0in05g3b17ZPuP5PzDtLcK4Cevo8tvf/rbu4zWr9dNsZjyfwWc/+1msWrUK73rXu/CJT3wCd999N+655x4cOHDAO149nwGpn7jvK+/l5DnrrLNwwQUX4IQTTsAHP/hBPPTQQwAQKh3Je998JnKPx/s5THuLkK2fWs9kPoNTTjkFAPDSSy9hwYIFdX0GpD66u7thmmbs95X3svF0dHTghBNOwIsvvogPf/jDAJRForvnALz3jURH6Fa7x729vSgUCjh48GDIKty3bx9WrlxZ97mmvUXY3d2Nt73tbVWXtra2uo61YsUKPPfcc9i7d6+3Lq7105NPPhma10p666fJfAbbt28HAO+LXM9nQOojnU5j+fLleOSRR0LrH3nkkXE9BEh95PN5/OEPf0BfXx+WLl2K3t7e0L0vFAp44okneO8bRD33ePny5UilUqF99u7di+eee258n8MkgnymHbt27ZLbt2+XN910k5wzZ47cvn273L59uzx8+LCUUspSqSSXLVsmTzvtNPm73/1OPvroo/Koo46S69at845x6NAh2dPTIz/2sY/JHTt2yB/84Aeyq6tL3nrrra26rBnD008/LW+//Xa5fft2+fLLL8sHHnhA9vf3y/POO8/bp57PgNTP/fffL1OplLznnnvkCy+8INevXy87OjrkK6+80uqhzXiuvfZauWXLFvnyyy/LX/3qV/Kcc86RnZ2d3r396le/KnO5nPzBD34gd+zYIT/2sY/Jvr4+OTQ01OKRzxwOHz7sPacBeM+PXbt2SSnru8dXXHGFPOqoo+Sjjz4qf/e738kPfOAD8sQTT5SlUqnuccwqIbz00kslgLLl8ccf9/bZtWuXPPvss2U2m5Xz58+X69atC6VKSCnl73//e/m+971PZjIZ2dvbKzds2MDUiTrYtm2bPPnkk2Uul5NtbW3yb//2b+WNN94oR0ZGQvvV8xmQ+vm3f/s3uWTJEplOp+V73vMe+cQTT7R6SLOCCy+8UPb19clUKiX7+/vl+eefL59//nlvu+M48sYbb5S9vb0yk8nI97///XLHjh0tHPHM4/HHH499Zl966aVSyvru8ejoqFy3bp2cP3++zGaz8pxzzpG7d+8e1zjYhokQQkiimfZzhIQQQkgzoRASQghJNBRCQgghiYZCSAghJNFQCAkhhCQaCiEhhJBEQyEkhBCSaCiEhBBCEg2FkBBCSKKhEBJCCEk0FEJCCCGJ5v8Hu0jSE5dxW5YAAAAASUVORK5CYII=",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "h5file = h5list[0]\n",
    "hfp = h5py.File(h5file,'r')\n",
    "dx_cgs = hfp['header']['camera']['dx'][()] * hfp['header']['units']['L_unit'][()]\n",
    "fovMuas = dx_cgs / hfp['header']['dsource'][()] * 2.06265e11\n",
    "cgsToJy = hfp['header']['scale'][()]\n",
    "unpol = np.copy(hfp['unpol']).transpose((1,0)) * cgsToJy\n",
    "hfp.close()\n",
    "\n",
    "# make plot and show\n",
    "ext = [ -fovMuas/2, fovMuas/2, -fovMuas/2, fovMuas/2 ]\n",
    "plt.imshow(unpol, origin='lower', extent=ext)\n",
    "plt.show()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "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.12.2"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}