1_HPWpipe.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\"This script will concatenate uvf5 files in time based on certain parameters.\"\"\"\n",
    "# Python Standard Library Packages\n",
    "import os\n",
    "import glob\n",
    "import argparse\n",
    "\n",
    "# Community Developed Packages\n",
    "import numpy as np\n",
    "import astropy.units as u\n",
    "import astropy.coordinates as aco\n",
    "\n",
    "# HERA Collaboration Packages\n",
    "from pyuvdata import UVData"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "_StoreAction(option_strings=['-S', '--savepath'], dest='savepath', nargs=None, const=None, default=None, type=None, choices=None, help='Designate the path where the new hdf5 files will be saved.', metavar=None)"
      ]
     },
     "execution_count": 2,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "parser = argparse.ArgumentParser()\n",
    "parser.add_argument(\n",
    "    '-F',\n",
    "    '--files',\n",
    "    help='Designate the hdf5 files to be concatenated in time.',\n",
    "    nargs='*',\n",
    "    required=True)\n",
    "parser.add_argument(\n",
    "    '-P',\n",
    "    '--pols',\n",
    "    help='Designate which polarizations to analyze (e.g.: \"xx xy yx yy\").',\n",
    "    nargs='*',\n",
    "    required=True)\n",
    "parser.add_argument(\n",
    "    '-D',\n",
    "    '--day',\n",
    "    help='Designate which JD the hdf5 files come from.',\n",
    "    required=True)\n",
    "parser.add_argument(\n",
    "    '-E',\n",
    "    '--ext',\n",
    "    help='Designate which file extension (i.e., uvOCRS) the designated hdf5 files are.',\n",
    "    required=True)\n",
    "parser.add_argument(\n",
    "    '-L',\n",
    "    '--LSTrng',\n",
    "    help='Designate which LST in hours that will be analyzed (e.g.: \"6.0 7.0\").',\n",
    "    type=float,\n",
    "    nargs=2,\n",
    "    required=True)\n",
    "parser.add_argument(\n",
    "    '-X',\n",
    "    '--xants',\n",
    "    help='Designate which antenna numbers that should be excluded from analysis (e.g.: \"0 50 98\").',\n",
    "    type=int,\n",
    "    nargs='*',\n",
    "    default=list())\n",
    "parser.add_argument(\n",
    "    '-S',\n",
    "    '--savepath',\n",
    "    help='Designate the path where the new hdf5 files will be saved.',\n",
    "    required=True)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'Uncomment this code when running as .py:'"
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "\"\"\"Uncomment this code when running as .py:\"\"\"\n",
    "# args = parser.parse_args()\n",
    "# files = np.array(sorted(args.files))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "metadata": {},
   "outputs": [],
   "source": [
    "\"\"\"Uncomment this code when running as .ipynb:\"\"\"\n",
    "args = parser.parse_args(\n",
    "    \"-F /data4/paper/HERA19Golden/RawData/2457548/zen.2457548.?????.??.HH.uvcRP \\\n",
    "    -P xx xy yx yy\\\n",
    "    -D 2457548 \\\n",
    "    -E uvcRP \\\n",
    "    -L 10.5 23 \\\n",
    "    -X 22 43 80 81 \\\n",
    "    -S /data4/paper/fox/HERA19Golden/\".split())\n",
    "files = sorted(glob.glob(args.files[0]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Saving files to:\n",
      "/data4/paper/fox/HERA19Golden/uvcRP/2457548/LSThrs_10.5_to_23.0\n",
      "Polarizations: ['xx', 'xy', 'yx', 'yy']\n",
      "Extension: uvcRP\n",
      "LST Range: Hour 10.5 to Hour 23.0 \n",
      "Excluded Antennae: [22, 43, 80, 81]\n"
     ]
    }
   ],
   "source": [
    "\"\"\"Formatting command line arguments:\"\"\"\n",
    "pols = sorted(args.pols)\n",
    "day = args.day\n",
    "ext = args.ext\n",
    "LSTrng = args.LSTrng\n",
    "xants = sorted(args.xants)\n",
    "savepath = os.path.join(args.savepath, '{ext}/{day}/LSThrs_{LST0}_to_{LSTf}'.format(\n",
    "    ext=ext,\n",
    "    day=day,\n",
    "    LST0=LSTrng[0],\n",
    "    LSTf=LSTrng[-1]))\n",
    "os.system('mkdir -p {}'.format(savepath))\n",
    "print 'Saving files to:\\n{}'.format(savepath)\n",
    "print 'Polarizations: {}'.format(pols)\n",
    "print 'Extension: {}'.format(ext)\n",
    "print 'LST Range: Hour {} to Hour {} '.format(LSTrng[0], LSTrng[-1])\n",
    "print 'Excluded Antennae: {}'.format(xants)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "yy\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.16694.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.17390.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.18086.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.18782.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.19478.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.20173.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.20869.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.21565.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.22261.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.22957.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.23653.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.24349.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.25045.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.25741.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.26437.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.27133.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.27829.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.28525.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.29221.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.29917.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.30613.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.31309.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.32004.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.32700.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.33396.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.34092.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.34788.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.35484.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.36180.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.36876.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.37572.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.38268.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.38964.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.39660.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.40356.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.41052.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.41748.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.42444.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.43140.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.43836.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.44531.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.45227.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.45923.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.46619.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.47315.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.48011.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.48707.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.49403.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.50099.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.50795.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.51491.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.52187.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.52883.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.53579.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.54275.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.54971.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.55667.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.56363.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.57058.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.57754.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.58450.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.59146.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.59842.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.60538.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.61234.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.61930.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.62626.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.63322.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.64018.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.64714.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.65410.yy.HH.uvcRP\n",
      "/data4/paper/HERA19Golden/RawData/2457548/zen.2457548.66106.yy.HH.uvcRP\n"
     ]
    },
    {
     "ename": "ValueError",
     "evalue": "Required UVParameter _Nants_data has not been set.",
     "output_type": "error",
     "traceback": [
      "\u001b[0;31m---------------------------------------------------------------------------\u001b[0m",
      "\u001b[0;31mValueError\u001b[0m                                Traceback (most recent call last)",
      "\u001b[0;32m<ipython-input-6-fe394ed1c410>\u001b[0m in \u001b[0;36m<module>\u001b[0;34m()\u001b[0m\n\u001b[1;32m     42\u001b[0m             \u001b[0mant_str\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0;34m'cross'\u001b[0m\u001b[0;34m,\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     43\u001b[0m             time_range=[time[0], time[-1]])\n\u001b[0;32m---> 44\u001b[0;31m         \u001b[0muvd\u001b[0m \u001b[0;34m+=\u001b[0m \u001b[0muvdi\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m     45\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m     46\u001b[0m     \u001b[0;31m# Making new metadata\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/fortino/anaconda3/envs/hera/lib/python2.7/site-packages/pyuvdata-1.2.1-py2.7-linux-x86_64.egg/pyuvdata/uvdata.pyc\u001b[0m in \u001b[0;36m__iadd__\u001b[0;34m(self, other)\u001b[0m\n\u001b[1;32m    989\u001b[0m             \u001b[0mother\u001b[0m\u001b[0;34m:\u001b[0m \u001b[0mAnother\u001b[0m \u001b[0mUVData\u001b[0m \u001b[0mobject\u001b[0m \u001b[0mwhich\u001b[0m \u001b[0mwill\u001b[0m \u001b[0mbe\u001b[0m \u001b[0madded\u001b[0m \u001b[0mto\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    990\u001b[0m         \"\"\"\n\u001b[0;32m--> 991\u001b[0;31m         \u001b[0mself\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__add__\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mother\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0minplace\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mTrue\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    992\u001b[0m         \u001b[0;32mreturn\u001b[0m \u001b[0mself\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    993\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/fortino/anaconda3/envs/hera/lib/python2.7/site-packages/pyuvdata-1.2.1-py2.7-linux-x86_64.egg/pyuvdata/uvdata.pyc\u001b[0m in \u001b[0;36m__add__\u001b[0;34m(self, other, run_check, check_extra, run_check_acceptability, inplace)\u001b[0m\n\u001b[1;32m    789\u001b[0m         \u001b[0;32mif\u001b[0m \u001b[0;32mnot\u001b[0m \u001b[0misinstance\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mother\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mthis\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0m__class__\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    790\u001b[0m             \u001b[0;32mraise\u001b[0m \u001b[0mValueError\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0;34m'Only UVData objects can be added to a UVData object'\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0;32m--> 791\u001b[0;31m         \u001b[0mother\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mcheck\u001b[0m\u001b[0;34m(\u001b[0m\u001b[0mcheck_extra\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mcheck_extra\u001b[0m\u001b[0;34m,\u001b[0m \u001b[0mrun_check_acceptability\u001b[0m\u001b[0;34m=\u001b[0m\u001b[0mrun_check_acceptability\u001b[0m\u001b[0;34m)\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[0m\u001b[1;32m    792\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    793\u001b[0m         \u001b[0;31m# Check objects are compatible\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/fortino/anaconda3/envs/hera/lib/python2.7/site-packages/pyuvdata-1.2.1-py2.7-linux-x86_64.egg/pyuvdata/uvdata.pyc\u001b[0m in \u001b[0;36mcheck\u001b[0;34m(self, check_extra, run_check_acceptability)\u001b[0m\n\u001b[1;32m    357\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    358\u001b[0m         super(UVData, self).check(check_extra=check_extra,\n\u001b[0;32m--> 359\u001b[0;31m                                   run_check_acceptability=run_check_acceptability)\n\u001b[0m\u001b[1;32m    360\u001b[0m \u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    361\u001b[0m         \u001b[0;31m# Check internal consistency of numbers which don't explicitly correspond\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;32m/home/fortino/anaconda3/envs/hera/lib/python2.7/site-packages/pyuvdata-1.2.1-py2.7-linux-x86_64.egg/pyuvdata/uvbase.pyc\u001b[0m in \u001b[0;36mcheck\u001b[0;34m(self, check_extra, run_check_acceptability)\u001b[0m\n\u001b[1;32m    242\u001b[0m                 \u001b[0;32mif\u001b[0m \u001b[0mparam\u001b[0m\u001b[0;34m.\u001b[0m\u001b[0mrequired\u001b[0m \u001b[0;32mis\u001b[0m \u001b[0mTrue\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    243\u001b[0m                     raise ValueError('Required UVParameter ' + p\n\u001b[0;32m--> 244\u001b[0;31m                                      + ' has not been set.')\n\u001b[0m\u001b[1;32m    245\u001b[0m             \u001b[0;32melse\u001b[0m\u001b[0;34m:\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n\u001b[1;32m    246\u001b[0m                 \u001b[0;31m# Check parameter shape\u001b[0m\u001b[0;34m\u001b[0m\u001b[0;34m\u001b[0m\u001b[0m\n",
      "\u001b[0;31mValueError\u001b[0m: Required UVParameter _Nants_data has not been set."
     ]
    }
   ],
   "source": [
    "\"\"\"Looping through each polarization, finding correct files based on given LST range:\"\"\"\n",
    "pol_files = {pol: sorted([file for file in files if pol in file]) for pol in pols}\n",
    "for pol, dfiles in pol_files.items():\n",
    "    print pol\n",
    "    \n",
    "    # Finding the correct files based on the provided LST range\n",
    "    uvd = UVData()\n",
    "    files = []\n",
    "    times = []\n",
    "    for dfile in dfiles:\n",
    "        print dfile\n",
    "        # uvd.read_uvh5(dfile, read_data=False)\n",
    "        uvd.read_miriad(dfile)\n",
    "        LSTrads = np.unique(uvd.lst_array * u.rad)\n",
    "        LSThrs = aco.Angle(LSTrads).hour\n",
    "        LSTindices = np.where(np.logical_and(LSThrs >= LSTrng[0], LSThrs <= LSTrng[-1]))[0]\n",
    "\n",
    "        if LSTindices.size > 0:\n",
    "            JDtimes = np.take(np.unique(uvd.time_array), LSTindices)\n",
    "            files.append(dfile)\n",
    "            times.append(JDtimes.tolist())\n",
    "\n",
    "    # Loading in the correct data based on the provided LST range\n",
    "    uvd = UVData()\n",
    "#     uvd.read_uvh5(\n",
    "#         files[0],\n",
    "#         ant_str='cross',\n",
    "#         times=times[0])\n",
    "    uvd.read_miriad(\n",
    "        files[0],\n",
    "        ant_str='cross',\n",
    "        time_range=[times[0][0], times[0][-1]])\n",
    "\n",
    "    for file, time in zip(files[1:], times[1:]):\n",
    "        uvdi = UVData()\n",
    "#         uvdi.read_uvh5(\n",
    "#             file, \n",
    "#             ant_str='cross',\n",
    "#             times=time)\n",
    "        uvd.read_miriad(\n",
    "            file,\n",
    "            ant_str='cross',\n",
    "            time_range=[time[0], time[-1]])\n",
    "        uvd += uvdi\n",
    "\n",
    "    # Making new metadata\n",
    "    UVD0 = UVData()\n",
    "    UVDf = UVData()\n",
    "#     UVD0.read_uvh5(files[0], read_data=False)\n",
    "    UVD0.read_miriad_metadata(files[0])\n",
    "#     UVDf.read_uvh5(files[-1], read_data=False)\n",
    "    UVDf.read_miriad_metadata(files[-1])\n",
    "\n",
    "    JDt0 = np.round(np.unique(UVD0.time_array)[0], 5)\n",
    "    JDtf = np.round(np.unique(UVDf.time_array)[0], 5)\n",
    "    JD = str(int(JDt0))\n",
    "    JDt0 = '{:.5f}'.format(JDt0).split('.')[-1]\n",
    "    JDtf = '{:.5f}'.format(JDtf).split('.')[-1]\n",
    "    numfiles = len(files)\n",
    "\n",
    "    # Saving new metadata\n",
    "    uvd.extra_keywords['JDt0'] = JDt0\n",
    "    uvd.extra_keywords['JDtf'] = JDtf\n",
    "    uvd.extra_keywords['JD'] = JD\n",
    "    uvd.extra_keywords['numfiles'] = numfiles\n",
    "    uvd.extra_keywords['ext'] = ext\n",
    "    uvd.extra_keywords['LSTrng'] = LSTrng\n",
    "    uvd.extra_keywords['xants'] = xants\n",
    "    \n",
    "    # Naming new UVdata object\n",
    "    hdf5 = 'zen.{JD}.{JDt0}_{JDtf}.{pol}.HH.hdf5.{ext}'.format(\n",
    "        JD=JD,\n",
    "        JDt0=JDt0,\n",
    "        JDtf=JDtf,\n",
    "        pol=pol,\n",
    "        ext=ext)\n",
    "    hdf5 = os.path.join(savepath, hdf5)\n",
    "    print 'Writing:'\n",
    "    print hdf5\n",
    "    uvd.write_uvh5(hdf5, clobber=True)"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.15"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}