demo_storm_coordinates.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "e7df11fd-aee9-4a57-ba89-e3ade0eaf613",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# demo storm coordinates\n",
    "import pandas as pd\n",
    "import numpy as np\n",
    "import geojson\n",
    "import matplotlib.pyplot as plt\n",
    "%run -i storm_coords.py"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "7e23bcdc-04ca-4360-b68d-58da77b26d01",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "111.19492664455889 180.0\n"
     ]
    }
   ],
   "source": [
    "### test \n",
    "lat1 = 47.\n",
    "lon1 = -62.1\n",
    "lat2 = 46.\n",
    "lon2 = -62.1\n",
    "\n",
    "d, az = dist_bearing( lon1, lat1, lon2, lat2)\n",
    "print(d, az)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b3631470-f5c4-4fda-b819-24b54af90bb3",
   "metadata": {},
   "source": [
    "#### Explore the geojson file a bit"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "4cfed371-f9f1-43de-b70b-14f474b70ca9",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "{\"geometry\": {\"coordinates\": [-84.0, 16.0], \"type\": \"Point\"}, \"properties\": {\"max_wind_speed_mph\": 17.26, \"minimum_sea_level_pressure_mb\": 0.0, \"radius_to_max_wind_nmi\": 0.0, \"storm_class\": \"DB\", \"time_utc\": \"2023-08-24T18:00:00\"}, \"type\": \"Feature\"}\n",
      "{\"coordinates\": [-84.0, 16.0], \"type\": \"Point\"}\n",
      "-84.0\n",
      "{'time_utc': '2023-08-24T18:00:00', 'max_wind_speed_mph': 17.26, 'minimum_sea_level_pressure_mb': 0.0, 'radius_to_max_wind_nmi': 0.0, 'storm_class': 'DB'}\n",
      "2023-08-24T18:00:00\n",
      "65\n"
     ]
    }
   ],
   "source": [
    "with open('/vortexfs1/home/csherwood/proj/NOPP_Idalia/idalia_track.json') as f:\n",
    "    gj = geojson.load(f)\n",
    "features = gj['features']\n",
    "print(features[0])\n",
    "print(features[0]['geometry'])\n",
    "print(features[0]['geometry']['coordinates'][0])\n",
    "print(features[0]['properties'])\n",
    "print(features[0]['properties']['time_utc'])\n",
    "\n",
    "print(len(features))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "cdb1e38d-6789-4672-917c-9ae361446c0a",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "mph2ms = 0.44704\n",
    "nmi2km = 1.852\n",
    "time_utc = []\n",
    "lat = []\n",
    "lon = []\n",
    "max_wind_speed_ms = []\n",
    "radius_to_max_wind_km = []\n",
    "for feature in features:\n",
    "    lon.append( feature['geometry']['coordinates'][0])\n",
    "    lat.append( feature['geometry']['coordinates'][1])\n",
    "    time_utc.append( feature['properties']['time_utc'])\n",
    "    max_wind_speed_ms.append( feature['properties']['max_wind_speed_mph']*mph2ms )\n",
    "    radius_to_max_wind_km.append( feature['properties']['radius_to_max_wind_nmi']*nmi2km )\n",
    "\n",
    "# convert time stamps to datetime64\n",
    "dt = pd.to_datetime(time_utc)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "a7516b36-20e0-4a6c-9e16-43f70fa24a13",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# Put in a dataframe and drop duplicates.\n",
    "# dfo = pd.DataFrame( {'time_utc': time_utc}, index = dt)\n",
    "dfo = pd.DataFrame({'time_utc': time_utc,'lat': lat, 'lon': lon,\n",
    "                    'max_wind_speed_ms': max_wind_speed_ms, \n",
    "                    'radius_to_max_wind_km': radius_to_max_wind_km}, \n",
    "                   index = dt).drop_duplicates()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 20,
   "id": "5a68a836-c4c3-4b5c-ac05-9eecac2ea4e7",
   "metadata": {
    "tags": []
   },
   "outputs": [],
   "source": [
    "# compute distance and bearing between unique points\n",
    "%run -i 'storm_coords.py'\n",
    "lon = np.array( dfo['lon'].values)\n",
    "lat = np.array( dfo['lat'].values)\n",
    "\n",
    "dist = np.ones_like( lon )\n",
    "brg = np.ones_like( dist )\n",
    "for i in range(len( lon )-1):\n",
    "    dist[i],brg[i] = dist_bearing(lon[i], lat[i], lon[i+1], lat[i+1])\n",
    "\n",
    "# make all the lists the same length\n",
    "dist[i+1]=dist[i]\n",
    "brg[i+1]=brg[i]\n",
    "\n",
    "# add to the dataframe\n",
    "dfo['distance_to_next_km']=dist\n",
    "dfo['bearing_to_next_geo']=brg\n",
    "dfo.to_csv('idalia_track_dist_brg.csv')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 19,
   "id": "ced863b0-c491-4897-bfd0-af3af96ac112",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "[1692900000 1692921600 1692943200 1692964800 1692986400 1693008000\n",
      " 1693029600 1693051200 1693072800 1693094400 1693116000 1693137600\n",
      " 1693159200 1693180800 1693202400 1693224000 1693245600 1693267200\n",
      " 1693288800 1693310400 1693332000 1693353600 1693375200 1693396800\n",
      " 1693418400 1693440000 1693461600 1693483200 1693504800 1693526400\n",
      " 1693548000 1693569600 1693591200 1693612800 1693634400 1693656000\n",
      " 1693677600]\n",
      "37\n",
      "1693328400\n",
      "[19] [1693310400]\n",
      "[1693310400 1693332000]\n",
      "[0. 1.]\n",
      "0.8333333333333333\n"
     ]
    }
   ],
   "source": [
    "# example interpolation\n",
    "import calendar\n",
    "def toTimestamp(d):\n",
    "    return calendar.timegm(d.timetuple())\n",
    "\n",
    "t = np.array( [toTimestamp( d ) for d in dfo.index ] )\n",
    "print(t)\n",
    "print(len(t))\n",
    "dti = calendar.timegm( pd.to_datetime( '2023-08-29 1700').timetuple() )\n",
    "print(dti)\n",
    "tindex=np.argwhere(t<=dti)[-1]\n",
    "print(tindex, t[tindex])\n",
    "tarray = np.squeeze( np.array( (t[tindex],t[tindex+1])) )\n",
    "one = np.squeeze( np.array( (0., 1.) ) )\n",
    "print(tarray)\n",
    "print(one)\n",
    "result = np.interp(dti, tarray, one)\n",
    "print(result)"
   ]
  },
  {
   "cell_type": "markdown",
   "id": "b49dd265-c8a2-489a-8d81-53c58e933ec8",
   "metadata": {},
   "source": [
    "* load track geojson\n",
    "* convert to lists of points with time, lat, lon, radius, and max winds\n",
    "* calculate distance and bearing between points\n",
    "* load list of buoy locations with time, lat, long\n",
    "* determine storm location, bearing, and translation speed for buoy time\n",
    "* determine storm coordinates (with and without radial scaling) for buoy location"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "fd47fa34-7c57-4833-b66c-0ddc683a6fe1",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/plain": [
       "<matplotlib.collections.PathCollection at 0x2aac0509d490>"
      ]
     },
     "execution_count": 11,
     "metadata": {},
     "output_type": "execute_result"
    },
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAiwAAAGdCAYAAAAxCSikAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjcuMSwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy/bCgiHAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA/NklEQVR4nO3deXwU9f3H8ffsQjYJSRYCJiEHCAgiIofchwgolwd4ghdISz0Br7YqVn+Ftpqq1XoWbaugrRxVQagIEuRWUK4goAICEoSEO9kkkGt3fn9QozG7m2xIdifJ6/l47OPRne/Mdz8zrN13vjPzHcM0TVMAAAAWZgt1AQAAABUhsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMtrEOoCqovH49GhQ4cUHR0twzBCXQ4AAKgE0zSVm5urxMRE2Wy+x1HqTGA5dOiQUlJSQl0GAACoggMHDig5Odlne50JLNHR0ZLO7HBMTEyIqwEAAJXhcrmUkpJS+jvuS50JLD+cBoqJiSGwAABQy1R0OQcX3QIAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMsjsAAAAMurMzPdAqi61V/v0ytp63TUladzohvpniG9NahDm1CXBQClCCxAPffY3I+1cNNXpe+PuvI1eeZCXdHlfD1zyxUhrAwAfkRgAULkeN4pPfXBcn26a788HlPnNz9HU0YNVIfk+KDVsH53Rpmw8lMfpe/UVV0v0IALWgWtnrORV1CoP368Qst37VWx26M2zWI1ZegAdW/h++mvAGoPrmEBQuDgiRwNeeqf+vjL3corKNKpomJt2X9Io1+apSVbdwWtjpc//sxv+6tp64JUydk5nndKl7z4D33w5ddyFRTqdHGxtmce1q1vvas5m74MdXkAqgGBBQiB+95aqKISt9e2x+YukcfjCUodR3Pz/bYfq6DdKu5//0OdKir22jZt8XKVBOl4Aqg5BBYgyEpKPNqZecxne1GJW2nbvw1KLQnOKL/t8TH+261i04FDPts8pql3NmwNYjUAagLXsABBlldUVOE6B45l13whku4f3k+3v/auz/b7hvcLSh1ny2Oaftu/z84JUiW+HTiZowfnLdL2Q4dlSnI0sOuGLh31fyMGh7o0oFZghAUIspjwMNkMw+86PVoH50LRbq2TdVOfzl7bbujZUb3btghKHWcrzG73296zZVKQKvHusCtPV0yfqW3/CyuSVFji1jsbt2r8v94LaW1AbUFgAYLMZrPp0gta+2xvFh2pzucmBq2ex68drH/dO1rdWyUpOdapbucm6a27b9TUG4YErYazNfKi9j7bohxhGtK+bRCrKe/xRWkqcnu/jmbddwe0I/NwkCsCah/DNCsYS60lXC6XnE6ncnJyFBMTE+pyAL9KSjy69q9va9/Rk2WWR4Q11LwHblNKs8ahKayW8ng8uvHNOdr+sx/+MLtdc395kzokxIWosjM6pb6kQh8XWUvSiA7t9ML1VwaxIsA6Kvv7zTUsQAg0aGDTf387Xsu279bsT7eqsKREgy5so18M6CabjYHPQNlsNr3/q1u0Zs93enPdJp0qLlb/1i10V7+eCmsQ+v+bc3v8/114utj7HU4AfhT6/5KBeuzyjm11ecfQnq6oSy5pc64uaXNuqMsoJ6WJU/uOn/TZPrid71OEAM7gTzkAqGEPX3aJz7aosDDd2KVjEKsBaicCCwDUsMHnt9HDl11S7u6w2MgIfXDnbZwGBCqBi24BIEiKSko0a+NWHcnLV79WLdWvTctQlwSEHBfdAoDFhDVooPG9u4W6DKBWYhwSAABYHiMsAFDPHMnJ0/MfrdHXB48oJiJcEwZ118AObUJdFuAXgQUA6pHF6d/o4VmL9dOLFyfNXKjurZI0857RIasLqAiBBXVW9qkC/XnZKn22d78kQ/1bt9TDlw9Q48jwUJcGhMSpoiI9MnuJvN1psXHfQb269DNNHNo36HUBlcE1LKiTdh89rv5/fV3zt36lw7n5Opybp/e37lD/F17X3mMnQl0eEBJ//+QLv0+2nv3Z1iBWAwSGwII66Y5Z81TsKf+wuWK3R7+aNT8EFQGhtyvzmN/23NOFQaoECByBBXXOgZM5ynTl+Ww/mONSZo4riBUB1pAc6/TbHh7WMEiVAIELKLBMnz5dnTp1UkxMjGJiYtSnTx8tXrxYklRcXKxHHnlEF110kRo1aqTExESNGzdOhw4d8tvnzJkzZRhGuVdBQUHV9wr12rdH/f8VKUnfcloI9dC9Q/v4bb+qa/sgVQIELqDAkpycrD//+c/auHGjNm7cqMGDB2vUqFHasWOHTp06pc2bN+uJJ57Q5s2bNW/ePO3atUsjR46ssN+YmBhlZmaWeYWHc2EkquaChLgK12l3TrMgVAJYS+PIcD0wvJ/XtpRYp6aMHBTkioDKC+guoauvvrrM+yeffFLTp0/X+vXrNWHCBKWlpZVpf/nll9WzZ09lZGSoRYsWPvs1DEMJCQmBlAL4lBATrZZNGmv/yWyv7a2bNlF8TFRwiwIs4leDe6pb6yT95cM12n/spCIdYbq+Z0fdMagHzzSCpVX5tma32613331X+fn56tPH+zBjTk6ODMNQ48aN/faVl5enli1byu12q0uXLvrjH/+orl27+t2msLBQhYU/XiDmcnFNAn701tjrdcX0t3WquLjM8siwhpp52w0hqgqwhq7nJumdSTeFugxlnnQpt6BQrc9pqgYNCEvwL+CHH27btk19+vRRQUGBoqKiNGvWLF1xxRXl1isoKFD//v3Vvn17/fvf//bZ3/r16/Xtt9/qoosuksvl0osvvqiPPvpIW7duVdu2bX1uN3XqVE2bNq3cch5+iB8UlZTob2s+19JvvpUhaegFbXVP/54Ka8D0Q0AopW3branvLVPO6TPXKtoNQ1defIH+dOMQRnnqoco+/DDgwFJUVKSMjAxlZ2fr/fff1z//+U+tWrVKHTp0KF2nuLhYN954ozIyMrRy5cqAAoTH49HFF1+sAQMG6KWXXvK5nrcRlpSUFAILAFjY6q/36d4ZH3htu6R9K03/5TVBrQehV2NPaw4LC9N5550nSerevbs2bNigF198Ua+//rqkM2Fl9OjR2rdvn5YvXx5weLDZbOrRo4d2797tdz2HwyGHwxFo+QCAEPrT/E98tq35Zp8OZ+cqvnF0ECtCbXHWY2+maZaOdPwQVnbv3q1ly5apadOmVeovPT1dzZs3P9vSAAAWcyg712/73HVfBqmSH+3LOqFVW/foq/2HFeBJBwRRQCMsjz32mEaMGKGUlBTl5uZqzpw5WrlypZYsWaKSkhLdcMMN2rx5sz788EO53W5lZWVJkmJjYxUWFiZJGjdunJKSkpSamipJmjZtmnr37q22bdvK5XLppZdeUnp6ul599dVq3lUAgNWVmOVnqK4pR7Lz9NibH2nz7oOly85LbKonf3mF2iYx9YHVBBRYDh8+rLFjxyozM1NOp1OdOnXSkiVLNGTIEH333XdauHChJKlLly5ltluxYoUGDhwoScrIyChzUVV2drbuvPNOZWVlyel0qmvXrlq9erV69ux5dnsGALCcplGROp53ymf7tT0uDEodxW637vrre/r+WHaZ5XuzTuiO5/+j+VPHq0l0ZFBqQeUEfNGtVVX2oh0AQOgs2LhDv/vPUq9tF6UkaPbkm4NSR9qmXXrkn4u8ttkMQ/dc3VcTRvCHczDU2EW3AABU1ajuFyr7VIFeWLxWxe4fT//0aJ2sf/zq+qDVsfnbg2pgs6nEy0NSPaapzd9+rwmq3sBy3JWvpRt3KTv/tNolnaMBnVurod1erZ9RE77Yc0AzVm3S6cIi9T//XI0f0D0k8+YQWAAAQXX7gG4a27+r1u7cr5zTBerfrqWaRAX39IujoV2mvJ9gsBmGwhtW78/je2u+1NNzlstjmrIbZ4JSQmy0Xp18nVolxFbrZ1UXj8ej2/42V19mZJUu27jvoKZ/8rnm3neLzosP/Maas8EMPQCAoLPZbBpwQStdffEFQQ8rkjTk4nZye7wHFo9paki3dtX2Wel7DuqpWZ/I7TFlmiod1TmanafJr8xXiTt4FxoHYtq8ZWXCyg8Ki0s07m9zg14PgQUAUO9ceG6Crux1QbnlNsNQlzaJuuxi3zOtB2rW8i2y24xyy90eU4eOu/Tpjn3V9lnV6b+bv/HZ5jpdqJVf7QliNZwSAgDUU1PHDVX7lDjNXrFFh467FBsdoesv6aTxw3pU67UlOw8c9TmaY7cZ2nPouC7t1KZSfS3Zukuvf7Jex1z5io2K1B2De+qqi8sHL2++PXxcf16wUrsyjyqsYQNd0eV8TRrax+fjSopK3H7727Dnew3sULm6qwOBBQBQL9ltNt162cW69bKL5fGYsnkZBakOzZyN9P2xbHm7J9ftMRVbydunf/vOIi3euqv0/clTBXp0zhJ9/OVuvTx+pN9tP9z8tR6ds6TMsjdXbtT8L7br40cnKDI8rNw2hiGvNf8gsUlw78jllBAAoN6rqbAiSaP6Xujzhz+sgV2XdT2vwj4+27W/TFj5qRVf7VHaNt+Psykp8eh3//nYa9vJUwX6zayPvLZ1buF7xnmbYWhM785+Kq5+BBYAAGrQFT0v0MDOZ06d2IwzwchuM2QzDE0dN1TRkeEV9vHK0s/8tr/+yec+2+au3+rzlJQkfbprv9flz992lc+7pR4deWnQb23mlBBqVE7Rad380Vx9dfKIpDMJuV9iS828/AbZa8H8AwBwthrYbXr2zqu0dOMuLVy3QydyT6lDi3iNGdRF7VPiKtWHv9mBJemEn/Z9R0/63dbtZS4aSYpzRmnZ7+7Qk/M/0epv9qnE7VGruFj99qoB6nVei4qLrmYEFtSY00VF6j77VRX/5D8Gj6Q1h/ar97uvacNNE0NXHAAEkd1m04ie7TWiZ/sqbZ/UOEYHT7h8tif4ecJ1x5R4aZ3vvv1dYNw4MlzP3nplpWqsaZwSQo25b/WHZcLKTx09na+ZOzYFuSIAqJ0evOISv+0PjOjvs23kxRfI0cB3KBnWqfpu4a5JBBbUmNUH/c8t8Pr2L4JUCQDUbhe1SNAvB3b32nZL3y7q2SbF57Y2m03TJ1xbev3MT7U6p4n+dOOwaquzJnFKCDWmosfEn3YXB6kSAKj9HrriEl3Ztb3+umitDmW7FO+M0v3D+5855VOBnm1StPb3d+uvi9dq876DCg9rqFv6dtao7sF5OnZ1ILCgxsRFRCnzVK7P9guanBPEagCg9ju/+Tl67VfXVmnbmMhw/f76y6u5ouDhlBBqzB96+/8P48VLrw5SJQCA2o7AghozpGVbTejg/Zzr032HKy4yKsgVAQBqK8M0/U28W3u4XC45nU7l5OQoJia40wXDv7yiIv1u3cfan5utLs2a64keg5iDBQAgqfK/31zDghoXFRbG6R8AwFnhlBAAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALA8AgsAALC8BqEuANZS5M7UkRO/kdudqbCGFyqu8TOy2yNCXRYAoJ4LaIRl+vTp6tSpk2JiYhQTE6M+ffpo8eLFpe2maWrq1KlKTExURESEBg4cqB07dlTY7/vvv68OHTrI4XCoQ4cOmj9/fuB7grOWeeJ+7c/sptOFK1RU8o3yTr+vvZltdDJ3RqhLAwDUcwEFluTkZP35z3/Wxo0btXHjRg0ePFijRo0qDSXPPPOMnn/+eb3yyivasGGDEhISNGTIEOXm5vrsc926dRozZozGjh2rrVu3auzYsRo9erQ+//zzs9szBCQn/z3lnXrXa9uxnN+pyH0kyBUBAPAjwzRN82w6iI2N1bPPPqtf/vKXSkxM1AMPPKBHHnlEklRYWKj4+Hg9/fTTuuuuu7xuP2bMGLlcrjIjNcOHD1eTJk00e/bsStfhcrnkdDqVk5OjmJiYs9mlemnvoU5ye475bA8P66uUuPeCWBEAoD6o7O93lS+6dbvdmjNnjvLz89WnTx/t27dPWVlZGjp0aOk6DodDl156qT777DOf/axbt67MNpI0bNgwv9tIZ8KQy+Uq80LVuT0n/bYXl+wMUiUAAJQXcGDZtm2boqKi5HA4dPfdd2v+/Pnq0KGDsrKyJEnx8fFl1o+Pjy9t8yYrKyvgbSQpNTVVTqez9JWSkhLorqAM/9dfG0ajINUBAEB5AQeW888/X+np6Vq/fr3uuece3X777frqq69K2w3DKLO+aZrllv1cVbaZMmWKcnJySl8HDhwIcE/wU43CB/ttj435TZAqAQCgvIBvaw4LC9N5550nSerevbs2bNigF198sfS6laysLDVv3rx0/SNHjpQbQfmphISEcqMpFW0jnTnd5HA4Ai0fPsQ3ma69me0lFZRra2BvIWejG4JfFAAA/3PWE8eZpqnCwkK1atVKCQkJSktLK20rKirSqlWr1LdvX5/b9+nTp8w2krR06VK/26D62e1hat18h8LDuuvHr0UDRUVco1bN14eyNAAAAhtheeyxxzRixAilpKQoNzdXc+bM0cqVK7VkyRIZhqEHHnhATz31lNq2bau2bdvqqaeeUmRkpG655ZbSPsaNG6ekpCSlpqZKku6//34NGDBATz/9tEaNGqUFCxZo2bJlWrt2bfXuKSpkt0coJW5hqMsAAKCcgALL4cOHNXbsWGVmZsrpdKpTp05asmSJhgwZIkl6+OGHdfr0ad177706efKkevXqpaVLlyo6Orq0j4yMDNlsPw7s9O3bV3PmzNHjjz+uJ554Qm3atNHcuXPVq1evatpFAABQ2531PCxWwTwsAADUPjU+DwsAAECwEFgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVgAAIDlEVjqGM/pT+TJf0ce94lQlwIAQLVpEOoCUD08uS9K+dMlec4syJ0mj72FFPuRbPawkNYGAMDZYoSlDvDkvS3lv6rSsPIDd4Z0fGAoSgIAoFoRWOqCvOd8t3mOyVO4IXi1AABQAwgsdcJp/81504NTBgAANYTAAgAALI/AUidE+m+OnhycMgAAqCEElrogeorvNlu8bGFdg1cLAAA1gMBSB9gajZGifq1y/5z286SmK0NQEQAA1Yt5WOoIW9RdUtRdZ+4IcmdKjsGy2aNCXRYAANWCwFLH2Bw9Ql0CAADVjlNCAADA8ggsAADA8ggsAADA8ggsAADA8ggsAADA8gIKLKmpqerRo4eio6MVFxena665Rjt37iyzjmEYXl/PPvusz35nzpzpdZuCgoKq7RUAAKhTAgosq1at0sSJE7V+/XqlpaWppKREQ4cOVX5+fuk6mZmZZV5vvvmmDMPQ9ddf77fvmJiYctuGh4dXba8AAECdEtA8LEuWLCnzfsaMGYqLi9OmTZs0YMAASVJCQkKZdRYsWKBBgwapdevWfvs2DKPctgAAANJZXsOSk5MjSYqNjfXafvjwYS1atEgTJkyosK+8vDy1bNlSycnJuuqqq7Rlyxa/6xcWFsrlcpV5AQCAuqnKgcU0TT300EPq37+/Onbs6HWdt956S9HR0bruuuv89tW+fXvNnDlTCxcu1OzZsxUeHq5+/fpp9+7dPrdJTU2V0+ksfaWkpFR1VwAAgMUZpmmaVdlw4sSJWrRokdauXavk5GSv67Rv315DhgzRyy+/HFDfHo9HF198sQYMGKCXXnrJ6zqFhYUqLCwsfe9yuZSSkqKcnBzFxMQE9HkAACA0XC6XnE5nhb/fVXqW0OTJk7Vw4UKtXr3aZ1hZs2aNdu7cqblz5wbcv81mU48ePfyOsDgcDjkcjoD7BgAAtU9Ap4RM09SkSZM0b948LV++XK1atfK57htvvKFu3bqpc+fOARdlmqbS09PVvHnzgLcFAAB1T0AjLBMnTtSsWbO0YMECRUdHKysrS5LkdDoVERFRup7L5dK7776r5557zms/48aNU1JSklJTUyVJ06ZNU+/evdW2bVu5XC699NJLSk9P16uvvlrV/QIAAHVIQIFl+vTpkqSBAweWWT5jxgyNHz++9P2cOXNkmqZuvvlmr/1kZGTIZvtxcCc7O1t33nmnsrKy5HQ61bVrV61evVo9e/YMpDwAAFBHVfmiW6up7EU7AADAOmr0oltYw86c5fr02BsqcucrqkGchjZ/WHER54W6LAAAqh2BpZaate9uHS38tvR9YVGeZu+/W50aj9KghMkhrAwAgOrH05proc+OzCgTVn7qy+wFOl6wP8gVAQBQswgstVD6yXl+29OyfD8ZGwCA2ojAUgsVmwV+23OLDwepEgAAgoPAUgvZjYZ+2yPsziBVAgBAcBBYaqH2MZf5bR+U8EBwCgEAIEgILLXQoLgHFGmP9drWslFPJUV6f3o2AAC1FYGlFrLb7bqj7X/UufF1CjMiZciuCLtTg+Me1DUpT4W6PAAAqh0z3QIAgJCp7O83IywAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyCCwAAMDyGoS6AFTNphM7tCPnW13UuJ26Nrkg1OUAAFCjCCy1zKYTO/TkV6/LLY8k6f2DaWpgNNDUC+/VRY3bhbg6AABqBqeEapGjp0/oD19NLw0rPygxS/T49peUU5QbosoAAKhZBJZa5Nmdb/ptf37XW0GqBACA4CKw1CL7Th30274z97vgFAIAQJARWGoRQ8ZZtQMAUFsFFFhSU1PVo0cPRUdHKy4uTtdcc4127txZZp3x48fLMIwyr969e1fY9/vvv68OHTrI4XCoQ4cOmj9/fmB7Ug90jDnPb3v3Jh2CVAkAAMEVUGBZtWqVJk6cqPXr1ystLU0lJSUaOnSo8vPzy6w3fPhwZWZmlr4++ugjv/2uW7dOY8aM0dixY7V161aNHTtWo0eP1ueffx74HtVhD7X7hc9RFJtsurftrUGuCACA4DBM0zSruvHRo0cVFxenVatWacCAAZLOjLBkZ2frgw8+qHQ/Y8aMkcvl0uLFi0uXDR8+XE2aNNHs2bMr1YfL5ZLT6VROTo5iYmIC2o/a5PtTh/Xol88rt+THkOhsEKXnOj+scyJiQ1gZAACBq+zv91nNw5KTkyNJio0t+0O5cuVKxcXFqXHjxrr00kv15JNPKi4uzmc/69at04MPPlhm2bBhw/TCCy+cTXl1UnJkvP7d+2m5ivK0/9QhtYpMUVRYRKjLAgCgRlU5sJimqYceekj9+/dXx44dS5ePGDFCN954o1q2bKl9+/bpiSee0ODBg7Vp0yY5HA6vfWVlZSk+Pr7Msvj4eGVlZfn8/MLCQhUWFpa+d7lcVd2VWikmLEoXhTFRHACgfqhyYJk0aZK+/PJLrV27tszyMWPGlP7vjh07qnv37mrZsqUWLVqk6667zmd/hlH22gzTNMst+6nU1FRNmzatitUDAIDapEq3NU+ePFkLFy7UihUrlJyc7Hfd5s2bq2XLltq9e7fPdRISEsqNphw5cqTcqMtPTZkyRTk5OaWvAwcOBLYTAACg1ggosJimqUmTJmnevHlavny5WrVqVeE2x48f14EDB9S8eXOf6/Tp00dpaWllli1dulR9+/b1uY3D4VBMTEyZFwAAqJsCOiU0ceJEzZo1SwsWLFB0dHTpqIjT6VRERITy8vI0depUXX/99WrevLm+++47PfbYY2rWrJmuvfba0n7GjRunpKQkpaamSpLuv/9+DRgwQE8//bRGjRqlBQsWaNmyZeVONwEAgPopoBGW6dOnKycnRwMHDlTz5s1LX3PnzpUk2e12bdu2TaNGjVK7du10++23q127dlq3bp2io6NL+8nIyFBmZmbp+759+2rOnDmaMWOGOnXqpJkzZ2ru3Lnq1atXNe0mAACozc5qHhYrqS/zsAAAUJdU9vebZwkBAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLI7AAAADLaxDqAhC4p7cu0+y9G1XgLlFkg4Yae15PPdhxUKjLAgCgxhBYapnLF7+ijPyTpe/zSoo0/Zu1+uTQTn049O4QVgYAQM3hlFAt8trXa8qElZ/a5Tqq2Xs2BbkiAACCg8BSi7y563O/7a98vTpIlQAAEFwEllrktLvYb3t+cWGQKgEAILgILLVIY0eE3/Zm4VFBqgQAgOAisNQij3Ua6rd9atcrglQJAADBRWCpRUakdNBlzdt5bbs6paP6J7QOckUAAAQHtzXXMtP7jdGmYxn6v00f6VhhvuLCo/Rkt6vUqWlSqEsDAKDGEFhqoW7NWmjRMOZcAQDUH5wSAgAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlkdgAQAAlsfTmi3O7Xbr9rT39Gnmfpn/W9Y8MlrvDB+t1s6mIa0NAIBgYYTF4i55/+9a+5OwIkmZp3J12bw3dDAvJ2R1AQAQTAEFltTUVPXo0UPR0dGKi4vTNddco507d5a2FxcX65FHHtFFF12kRo0aKTExUePGjdOhQ4f89jtz5kwZhlHuVVBQULW9qiMWf7dTh/JzvbaZku78ZH5wCwIAIEQCCiyrVq3SxIkTtX79eqWlpamkpERDhw5Vfn6+JOnUqVPavHmznnjiCW3evFnz5s3Trl27NHLkyAr7jomJUWZmZplXeHh41faqjnh+y1q/7V+fOBKkSgAACK2ArmFZsmRJmfczZsxQXFycNm3apAEDBsjpdCotLa3MOi+//LJ69uypjIwMtWjRwmffhmEoISEhkHLqvIKSEr/tniDVAQBAqJ3VNSw5OWeuoYiNjfW7jmEYaty4sd++8vLy1LJlSyUnJ+uqq67Sli1b/K5fWFgol8tV5lXX9E88129747D6PQIFAKg/qhxYTNPUQw89pP79+6tjx45e1ykoKNCjjz6qW265RTExMT77at++vWbOnKmFCxdq9uzZCg8PV79+/bR7926f26SmpsrpdJa+UlJSqrorlvWn3pfL8NM+pfvAYJUCAEBIGaZpmhWvVt7EiRO1aNEirV27VsnJyeXai4uLdeONNyojI0MrV670G1h+zuPx6OKLL9aAAQP00ksveV2nsLBQhYWFpe9dLpdSUlKUk5MT0GdZ3edZGbplyVy5f/bPdOv5nfVk32EhqgoAgOrhcrnkdDor/P2u0jwskydP1sKFC7V69WqfYWX06NHat2+fli9fHnCAsNls6tGjh98RFofDIYfDEXDttU2vhBbaM/63mrljk1Yc3KPkRk493mOQIsLCQl0aAABBE1BgMU1TkydP1vz587Vy5Uq1atWq3Do/hJXdu3drxYoVato08MnNTNNUenq6LrroooC3ravGX9hN4y/sFuoyAAAIiYACy8SJEzVr1iwtWLBA0dHRysrKkiQ5nU5FRESopKREN9xwgzZv3qwPP/xQbre7dJ3Y2FiF/W9UYNy4cUpKSlJqaqokadq0aerdu7fatm0rl8ull156Senp6Xr11Verc18BAEAtFVBgmT59uiRp4MCBZZbPmDFD48eP1/fff6+FCxdKkrp06VJmnRUrVpRul5GRIZvtx+t9s7OzdeeddyorK0tOp1Ndu3bV6tWr1bNnzwB3BwAA1EVVvujWaip70Q4AALCOyv5+8ywhAABgeQQWAABgeVW6rRnBsXr/Pj396Rp9l31SjgYNNKJNO/3ukksVyS3NAIB6hhEWi3p1w3qNXzBPXx87qtMlJcouKNDsHV+q/4x/KK+oKNTlAQAQVAQWC8otLNDz6z712pZdWKDfpi3x2gYAQF1FYLGg6Ru/kL9bt1Z8tzdotQAAYAUEFgs6mJvrt73Y7Q5SJQAAWAOBxYIuTmjutz2Ki24BAPUMgcWCxnbqojC73W87AAD1CYHFgmw2m94edYPshlGurV9KC/2m7yUhqAoAgNBhHhaL6pmcrK13TdKLX6zTxkMHFRPm0L09eqlHUnKoSwMAIOgILDXg8MGTythzRI2bRum8DokyvIyUVEZkWJim9L+0mqsDAKD2IbBUo9zsU3rusff0+cqv9cN9ySmtz9Fvnh6tdh0ZGQEAoKp4WnM1OLD/uHZsO6D/vL5CWfuOyXR7SttsNkOOiDC9tvABxSU2rnSfq/fv04e7dymqYUNN6NpdSTyBGgBQB1X295sRlrNw+lSR/vyHD/TZ6l0/LowOl+1UkWxFJZIkj8dUYUGR/vvOOk347YgK+zx26pSunv22Dufnly6buXWLRpzXTq9ecXW17wMAALUBdwmdhb889V+tX7u77ELDkKeRQ56GP96W7HGb2rL+20r1ec3cf5cJKz9Y/O0uPf3p6rOqFwCA2orAUkVZmdlavfxreTxezqiZpjyOhj++NyRHeMPy6/3MlsxDOuRnltu3t26pSqkAANR6BJYq2vV1pu9Gw5Aa/HhoDUkDRnSqsM+0vf5HYU6XlMjj8fhdBwCAuohrWKooslEF0+P/b+DFZjPUok2chl3XvcI+m0REVriOzUbGBADUP/z6VVHni89VjDPCZ7tRVKKIRg6NvK2vnv33XQqPrPj5P7d27CR/M7a0i21ahUoBAKj9CCxV1LChXQ8+eqUMmyGb/ceYYbMZim/u1Kwlv9b7G36vu6ZcpagY38HmpyLDwnRP955e22yGoReGXVkttQMAUNswD8tZ+mrb95r778+0LT1DjvCG6tq1hdyHT+r7nZlqEh+joTf1Vf+ruwZ0Kuc/O7bpuXWf6uipfNkMQx3PidMzQ4arXdNmNbgnAAAEX2V/vwks1Shtzjr99YF/yWYz5HZ7ZLMZ8nhMDbyuh3776niuPwEA4Gcq+/vNL2g1yT6Wq5d+845M05T7fzPd/nDL88p5G/Tph9ySDABAVRFYqsnqBZvkLvF+y7HNZujjWZ8F3KfH4+E2ZgAAxG3N1Sb7qEs2u03uEne5No/H1Mkjrkr39bfVn+uf6zYov6hYktSiiVOpI4eqewseoAgAqJ8YYanAqfxCnTieV3qax5eW7RO9hhVJstttanVhUqU+7/eLPtGLqz4rDSuSlHEyR7e99a42Znxf+cIBAKhDGGHx4dD3J/T6y8u0/tPd8nhMNYltpNG39tH1N/WSYZSfLaXPiM6KjXcq+1iuPD8LN26PRyMnDKrwM/MKCjV385de20xJjy1M09JJv6jS/gAAUJsxwuLFsaMu3XfnzNKwIkknT+Tr9ZeX6R+vfuJ1mzBHQ/1pzmQ1bhYtSbI3sMswDNkb2PXQC+PUrkvLCj931qat8nfL1v6T2YHuCgAAdQIjLF68N+tz5bpOe32w4XtzPtd1Y3qq2Tnlb71q1SFJMzb8UesWb9X+bw7J2Sxal17TvTTEVCS3sKjCdTweD7dHAwDqHX75vFiz6hvvT2GWZHpMfbFuj89twxwN1bZLS7Xr0lLndz1XzqZRlf7cKy883297lCOMsAIAqJcYYfHi59eglGv3catxzvE8/WXyTG38ZEfpssRW5+i3r/5C7bu1qvBz28efozbNYrXn2Amv7Xf161FhHwAA1EX8ue5Fr77nyW73/RjCi7uXDx+maer/bn1Fm1d+XWZ51v5jmnLDCzqccbxSn/3uhJvVpllsmWWGpFu6ddad/bw/ZwgAgLqOERYvbriljz75eLtMs6TMqSHDkC4f3kmJybHltklfs1O7tuwvt9zjMVVUWKKFb67UHVOvr/CzG4WF6aN7bteuw0f10Ve7FB3u0E3dOqlRWMVPewYAoK4isHiRnBKrv7w6Vi8+85F278ySJIWFNdDV13XThHsGe91m+/rdsjeweZ3t1uP2aOvanQHV0C7+HLWLPyfw4gEAqIMILD6cf0Gi/jbjVzr4/Qnl5RYopUVTRTZy+FzfER4m08eFujKk8EhGSAAAqKqArmFJTU1Vjx49FB0drbi4OF1zzTXaubPsyIFpmpo6daoSExMVERGhgQMHaseOHT56/NH777+vDh06yOFwqEOHDpo/f35ge1JDkpJjdf4FiX7DiiT1u6qrzzuLJOnSa7pLkvbsO6Lpb6zUMy8u0QeLtij/VGG11gsAQF0UUGBZtWqVJk6cqPXr1ystLU0lJSUaOnSo8vPzS9d55pln9Pzzz+uVV17Rhg0blJCQoCFDhig3N9dnv+vWrdOYMWM0duxYbd26VWPHjtXo0aP1+eefV33PgiypdZyuveuyM29+cr2uzWaodYdkXT66t16fsUq/nDhT736wUUuWbddfX03TTb94Xbv3HA5N0QAA1BKGaZr+Jlf16+jRo4qLi9OqVas0YMAAmaapxMREPfDAA3rkkUckSYWFhYqPj9fTTz+tu+66y2s/Y8aMkcvl0uLFi0uXDR8+XE2aNNHs2bMrVYvL5ZLT6VROTo5iYspP6hYMpmlq8b/Wat70ZTq494iiGkdq+K39dNMDI7Rp+wE98acPym1jsxlqGhulOTPuUgM7N20BAOqXyv5+n9U1LDk5OZKk2Ngzd83s27dPWVlZGjp0aOk6DodDl156qT777DOfgWXdunV68MEHyywbNmyYXnjhhbMpL+gMw9AV4y7RFeMuKTcj7byFm2WzGeVOG3k8po4ey9UXG/eqb6/zgl0yAAC1QpUDi2maeuihh9S/f3917NhRkpSVdeaOmvj4+DLrxsfHa//+8rf8/iArK8vrNj/0501hYaEKC3+8/sPlcgW8D9Xh+wMn9MG7X2jr5v2KiAzT4CEdNfzqLgoPb1h2vYMnfF/j0sDQ3DVb9conn8uUqXPPbaovc44q42S2Gths6tu6hf5wxeVqGhUZhD0CAMB6qhxYJk2apC+//FJr164t1/bzpxmbpun1Ccdns01qaqqmTZsWQMXVb+uW/Zry4Cy53R553GfCyNc7DurjRVv13Ktjy1yoGxcXo2Mn8vXzM3Aeu3Qq3q61ezNkmlJxpLS55OiZRkMqcru1bOcerfn2O308cbyaO0NzugsAgFCq0kUTkydP1sKFC7VixQolJyeXLk9ISJCkciMjR44cKTeC8lMJCQkBbzNlyhTl5OSUvg4cOFCVXakyt9ujP0/7QCXFP4YVSZIp7f32sOb867My64+6omu5sCJJhbF2mXZDpil5JBU5/9fws6xW6HbrvvcWVe9OAABQSwQUWEzT1KRJkzRv3jwtX75crVqVnaK+VatWSkhIUFpaWumyoqIirVq1Sn379vXZb58+fcpsI0lLly71u43D4VBMTEyZVzBt//KAjh3J9RpCPB5TSz5ML7NsyKAOGn75mVNndtuZNGLYDZVEGKXhxB2pM//bx8DStkO+T5EBAFCXBXRKaOLEiZo1a5YWLFig6Ojo0lERp9OpiIgIGYahBx54QE899ZTatm2rtm3b6qmnnlJkZKRuueWW0n7GjRunpKQkpaamSpLuv/9+DRgwQE8//bRGjRqlBQsWaNmyZV5PN1lF9sl8v+2unNNl3ttshh59cISGDOqgtBVfyeU6rfgkp97evK10HY/d/2eaUrmLeQEAqA8CCizTp0+XJA0cOLDM8hkzZmj8+PGSpIcfflinT5/Wvffeq5MnT6pXr15aunSpoqOjS9fPyMgo86Pbt29fzZkzR48//rieeOIJtWnTRnPnzlWvXr2quFs1r1XrOJ9thiG1bFV+Wn3DMNS967nq3vVcSVJhcYne2/GNThUWS5LshVJJdLnNSjW02wkrAIB66azmYbGSUMzD8vB9/9bWLfvLXsPyQ9sTIzVkRKcK+3j5g7WauXSDfvhXyI/XmRN1Xk4LXd/5Qj01cmj5BgAAaqnK/n7z5/pZeGzatTqv7ZkLje12mwzDkGFIN4zpqaRmjZT53dEK+7j76j4a2u18SZLNMNTouM6c+/mZLkkJ+tNVl1dn+QAA1BqMsJwl0zS1dfN+bf/ygOw2Q3vX79K6hZtKn9p8QY/WeuD529SiXXO//ew+eExrt++VaUq9L2ihLzIPasXufYoMa6gJfbqpe4tkv9sDAFAbVfb3m8BSTUzT1GM3vqQvP91ZZoI4m92mRjERem3VE4qNd/rpAQCA+odTQkH29ca9Sl/zTfmp990e5eec0oczVoWoMgAAaj8CSzXZsupr2Xw8vNDjMbXhkx1BrggAgLqDwFJNbHb/k6g0aFjBJCsAAMAnAks16T2skzxuj9c2wzDU76quQa4IAIC6g8BSTVp1SNJlo3vp589rtNltik+J1fBb+4WmMAAA6gACSzV68K9jNfbhqxXTpJGkM6eBBl3fQ88t+q2inJEhrg4AgNqL25prgNvtUe6JPEVEhcsRERbSWgAAsLLK/n4H9CwhVI7dblPjc0IbmgAAqEs4JQQAACyPwAIAACyPwAIAACyPwAIAACyPwAIAACyPwAIAACyPwBICuSfzlbX/mIqLSkJdCgAAtQLzsATR93sO67XH/6NNK76STCmqcaSuvesyjbl/uOw+nvQMAAAILEFz9NBJPXTls8p3nZb+N7dwXvYp/euZ/+p4VrYmP3NLaAsEAMDC+LM+SD54/RPlu06Xf6KzKX309hpl7T8WmsIAAKgFCCxB8tnireXDyv8Ykjau2BHcggAAqEUILEFS0TMmTU+deAYlAAA1gsASJL2GXCSbjwtrTVPqNqhDkCsCAKD2ILAEyXX3XK7wSEe50GIY0pCb+iixVVyIKgMAwPoILEESn9JUz/331+rQs3XpsvBGDo2ePEz3P3drCCsDAMD6DLOiiytqCZfLJafTqZycHMXExIS6HL+OZWYrP+eU4ls0U3hkWKjLAQAgZCr7+808LCHQrHljNWveuNzyooJirfnvZqWv+UZ2u029h3dWj8s7MqkcAKDeI7BYxIkjOXrk2r/q+28Py2a3yTCkj2d9pov6tNUfZk1iJAYAUK/xp7tFvPjQv3Vo31FJksftkbvkzJwtOz7/Vv965r+hLA0AgJAjsFjAscxsfbFsu9eJ5TweU4v/tUYlxe4QVAYAgDUQWCzgyPfHS58v5M3pvELl5ZwKXkEAAFgMgcUCzkls4rc9PDJMUc7IIFUDAID1EFgs4JykWHUffKHXmXBtNkPDb+uvBg3tIagMAABrILBYxP3P36b4lKaSzoSUH8LL+d1aadyjV4eyNAAAQo7bmi2iWfPG+tuKx7Xqg41KX/217A3s6jOis3oP6yR7A0ZXAAD1GzPdAgCAkKns73fAp4RWr16tq6++WomJiTIMQx988EGZdsMwvL6effZZn33OnDnT6zYFBQWBlgcAAOqggANLfn6+OnfurFdeecVre2ZmZpnXm2++KcMwdP311/vtNyYmpty24eHhgZZXbxzak6X0FduV9d2RUJcCAECNC/galhEjRmjEiBE+2xMSEsq8X7BggQYNGqTWrVv72OIMwzDKbYvyDu8/qmd/8Yq2rvyqdNnFQzrptzMmqllibAgrAwCg5tToXUKHDx/WokWLNGHChArXzcvLU8uWLZWcnKyrrrpKW7Zs8bt+YWGhXC5XmVdddzrvtB669P+0bc03ZZZvXbFdvxk0VUUFRSGqDACAmlWjgeWtt95SdHS0rrvuOr/rtW/fXjNnztTChQs1e/ZshYeHq1+/ftq9e7fPbVJTU+V0OktfKSkp1V2+5XzyzlodOXCs3BT+7hKPDu7O1Kp314WoMgAAalaNBpY333xTt956a4XXovTu3Vu33XabOnfurEsuuUT/+c9/1K5dO7388ss+t5kyZYpycnJKXwcOHKju8i0nfcU2GYbhtc3ewKb0FduDXBEAAMFRY/OwrFmzRjt37tTcuXMD3tZms6lHjx5+R1gcDoccDsfZlFjrNAhrIMMwZHp98JChhg2ZVgcAUDfV2AjLG2+8oW7duqlz584Bb2uaptLT09W8efMaqKz26jeqp9cnOkuSu8Stftf2DHJFAAAER8CBJS8vT+np6UpPT5ck7du3T+np6crIyChdx+Vy6d1339WvfvUrr32MGzdOU6ZMKX0/bdo0ffzxx9q7d6/S09M1YcIEpaen6+677w60vDqt76geuuiSC2Szlf1nM2yGug/trG5DAw+HAADUBgGfQ9i4caMGDRpU+v6hhx6SJN1+++2aOXOmJGnOnDkyTVM333yz1z4yMjLK/OhmZ2frzjvvVFZWlpxOp7p27arVq1erZ09GDH7K3sCupxb/Tv+a9q4W/T1N+TmnFN2kka6+Z5huffz6ckEGAIC6gqn5aymPx6PTeQWKiAonqAAAaq3K/n5zlWYtZbPZ1CgmMtRlAAAQFPxpDgAALI/AAgAALI/AAgAALI9rWOqo/V9/r1lPvq/PFm6U6THVY3gX3fq763Ve11ahLg0AgIBxl1Ad9G36Pj3Q/wmVFBXLXXJmojmb3Sab3aZnl/2fOva/IMQVAgBwRmV/vzklVAe9/uu3VVz4Y1iRJI/bI3eJW6/c92YIKwMAoGoILHWM63iu0lds9zqFv+kxtSf9O2XuPRyCygAAqDoCSx1TcKqwwnVO5xUEoRIAAKoPgaWOaZrYRE0Tm/hsb+SMVHI7HioJAKhdCCx1jN1u1y2PXe+z/cbfjFRYeFgQKwIA4OwRWOqgq+8ZqvF/vEkNwxuWLmvQ0K4xD4/SzVOuDWFlAABUDbc112F52flKX7FdpsdUp0s7yNmM4wIAsBYefghFNW6k/tf2CnUZAACcNU4JAQAAyyOwAAAAyyOwAAAAy+Malnom55hLS95coZ0bdisyOlKDbu6niy/vJMMwQl0aAAA+EVjqkV2b9ujhIX/QKddpyZRsdkMfz1yhQTf31yNvT5Ldbg91iQAAeMUpoXrC7XZr2vV/0encApkeU6Zplj4cccXstVr8z+UhrhAAAN8ILPVE+vLtOpJxzOtDEQ1D+u9rH4egKgAAKofAUk8c3n/MZ5tpSof3Hw1iNQAABIbAUk8kton32WbYDCW29t0OAECoEVjqiU6XdlBim3jZ7OX/yU2PqZETR4SgKgAAKofAUk/YbDZNnf+woptESYZks9tkb3DmrqAr7rhMw8YPDG2BAAD4wW3N9Uirji301rcv65N/r9E3X+xWo5gz87Bc0Lsd87AAACyNpzUDAICQqezvN6eEAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5RFYAACA5dWZZwn98IQBl8sV4koAAEBl/fC7XdGTgupMYMnNzZUkpaSkhLgSAAAQqNzcXDmdTp/tdebhhx6PR4cOHVJ0dHSdffKwy+VSSkqKDhw4wAMeg4xjH1oc/9Dh2IdWfTj+pmkqNzdXiYmJstl8X6lSZ0ZYbDabkpOTQ11GUMTExNTZL67VcexDi+MfOhz70Krrx9/fyMoPuOgWAABYHoEFAABYHoGlFnE4HPr9738vh8MR6lLqHY59aHH8Q4djH1oc/x/VmYtuAQBA3cUICwAAsDwCCwAAsDwCCwAAsDwCCwAAsDwCSy2wa9cujRo1Ss2aNVNMTIz69eunFStWlFnHMIxyr9deey1EFdctlTn+GRkZuvrqq9WoUSM1a9ZM9913n4qKikJUcd2xcuVKr99twzC0YcOG0vX4/le/yh57vvs1a9GiRerVq5ciIiLUrFkzXXfddWXa69N3v87MdFuXXXnllWrXrp2WL1+uiIgIvfDCC7rqqqu0Z88eJSQklK43Y8YMDR8+vPR9ZWYORMUqOv5ut1tXXnmlzjnnHK1du1bHjx/X7bffLtM09fLLL4e6/Fqtb9++yszMLLPsiSee0LJly9S9e/cyy/n+V6/KHHu++zXr/fff1x133KGnnnpKgwcPlmma2rZtW7n16s1334SlHT161JRkrl69unSZy+UyJZnLli0rXSbJnD9/fggqrNsqc/w/+ugj02azmQcPHixdZ/bs2abD4TBzcnKCXnNdVlRUZMbFxZl/+MMfyizn+1/zvB17vvs1p7i42ExKSjL/+c9/+l2vPn33OSVkcU2bNtUFF1ygt99+W/n5+SopKdHrr7+u+Ph4devWrcy6kyZNUrNmzdSjRw+99tpr8ng8Iaq67qjM8V+3bp06duyoxMTE0u2GDRumwsJCbdq0KVSl10kLFy7UsWPHNH78+HJtfP9rlrdjz3e/5mzevFkHDx6UzWZT165d1bx5c40YMUI7duwot259+e5zSsjiDMNQWlqaRo0apejoaNlsNsXHx2vJkiVq3Lhx6Xp//OMfddlllykiIkKffPKJfv3rX+vYsWN6/PHHQ1d8HVCZ45+VlaX4+Pgy2zVp0kRhYWHKysoKQdV11xtvvKFhw4YpJSWlzHK+/zXP27Hnu19z9u7dK0maOnWqnn/+eZ177rl67rnndOmll2rXrl2KjY2VVM+++6Ee4qmvfv/735uS/L42bNhgejwec+TIkeaIESPMtWvXmps2bTLvueceMykpyTx06JDP/v/yl7+YMTExQdyj2qU6j/8dd9xhDh06tNxnNGzY0Jw9e3awd61WqOzx/6kDBw6YNpvNfO+99yrsn++/b9V57PnuB66yx/+dd94xJZmvv/566bYFBQVms2bNzNdee81n/3X5u88IS4hMmjRJN910k991zj33XC1fvlwffvihTp48Wfpo8b/97W9KS0vTW2+9pUcffdTrtr1795bL5dLhw4fL/QWE6j3+CQkJ+vzzz8tse/LkSRUXF3Psfajs8f+pGTNmqGnTpho5cmSF/fP99606jz3f/cBV9vjn5uZKkjp06FC63OFwqHXr1srIyPC5bV3+7hNYQqRZs2Zq1qxZheudOnVKkmSzlb3cyGaz+T1PuWXLFoWHh5c5bYQfVefx79Onj5588kllZmaqefPmkqSlS5fK4XCUu84IZ1T2+P/ANE3NmDFD48aNU8OGDStcn++/b9V57PnuB66yx79bt25yOBzauXOn+vfvL0kqLi7Wd999p5YtW/rcrk5/90M9xAP/jh49ajZt2tS87rrrzPT0dHPnzp3mb37zG7Nhw4Zmenq6aZqmuXDhQvPvf/+7uW3bNvPbb781//GPf5gxMTHmfffdF+Lqa7/KHP+SkhKzY8eO5mWXXWZu3rzZXLZsmZmcnGxOmjQpxNXXHcuWLTMlmV999VW5Nr7/Ncvfsee7X7Puv/9+Mykpyfz444/Nb775xpwwYYIZFxdnnjhxwjTN+vfdJ7DUAhs2bDCHDh1qxsbGmtHR0Wbv3r3Njz76qLR98eLFZpcuXcyoqCgzMjLS7Nixo/nCCy+YxcXFIay67qjo+Jumae7fv9+88sorzYiICDM2NtacNGmSWVBQEKKK656bb77Z7Nu3r9c2vv81y9+xN02++zWpqKjI/PWvf23GxcWZ0dHR5uWXX25u3769tL2+ffcN0zTNUI/yAAAA+MM8LAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPIILAAAwPL+HwCHHJwu40chAAAAAElFTkSuQmCC\n",
      "text/plain": [
       "<Figure size 640x480 with 1 Axes>"
      ]
     },
     "metadata": {},
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.scatter(lon,lat,25,max_wind_speed_ms)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "d673dc2d-2263-4eec-90b5-9a9873fbb052",
   "metadata": {
    "tags": []
   },
   "outputs": [
    {
     "data": {
      "text/html": [
       "<div>\n",
       "<style scoped>\n",
       "    .dataframe tbody tr th:only-of-type {\n",
       "        vertical-align: middle;\n",
       "    }\n",
       "\n",
       "    .dataframe tbody tr th {\n",
       "        vertical-align: top;\n",
       "    }\n",
       "\n",
       "    .dataframe thead th {\n",
       "        text-align: right;\n",
       "    }\n",
       "</style>\n",
       "<table border=\"1\" class=\"dataframe\">\n",
       "  <thead>\n",
       "    <tr style=\"text-align: right;\">\n",
       "      <th></th>\n",
       "      <th>LONGITUDE</th>\n",
       "      <th>LATITUDE</th>\n",
       "      <th>STORMNAME</th>\n",
       "      <th>DTG</th>\n",
       "      <th>YEAR</th>\n",
       "      <th>MONTH</th>\n",
       "      <th>DAY</th>\n",
       "      <th>HHMM</th>\n",
       "      <th>MSLP</th>\n",
       "      <th>BASIN</th>\n",
       "      <th>STORMNUM</th>\n",
       "      <th>STORMTYPE</th>\n",
       "      <th>INTENSITY</th>\n",
       "      <th>SS</th>\n",
       "      <th>LAT</th>\n",
       "      <th>LON</th>\n",
       "    </tr>\n",
       "  </thead>\n",
       "  <tbody>\n",
       "    <tr>\n",
       "      <th>0</th>\n",
       "      <td>-86.1</td>\n",
       "      <td>21.3</td>\n",
       "      <td>TEN</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>26.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>1006.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TD</td>\n",
       "      <td>25.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>21.3</td>\n",
       "      <td>-86.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>1</th>\n",
       "      <td>-86.3</td>\n",
       "      <td>21.2</td>\n",
       "      <td>TEN</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>27.0</td>\n",
       "      <td>0</td>\n",
       "      <td>1005.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TD</td>\n",
       "      <td>30.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>21.2</td>\n",
       "      <td>-86.3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>2</th>\n",
       "      <td>-86.8</td>\n",
       "      <td>20.5</td>\n",
       "      <td>TEN</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>27.0</td>\n",
       "      <td>600</td>\n",
       "      <td>1001.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TD</td>\n",
       "      <td>30.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>20.5</td>\n",
       "      <td>-86.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>3</th>\n",
       "      <td>-86.3</td>\n",
       "      <td>19.8</td>\n",
       "      <td>TEN</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>27.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>1001.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TD</td>\n",
       "      <td>30.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>19.8</td>\n",
       "      <td>-86.3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>4</th>\n",
       "      <td>-85.7</td>\n",
       "      <td>20.0</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>27.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>996.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>35.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>20.0</td>\n",
       "      <td>-85.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>5</th>\n",
       "      <td>-85.4</td>\n",
       "      <td>19.8</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>28.0</td>\n",
       "      <td>0</td>\n",
       "      <td>992.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>45.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>19.8</td>\n",
       "      <td>-85.4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>6</th>\n",
       "      <td>-85.2</td>\n",
       "      <td>19.8</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>28.0</td>\n",
       "      <td>600</td>\n",
       "      <td>988.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>55.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>19.8</td>\n",
       "      <td>-85.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>7</th>\n",
       "      <td>-85.2</td>\n",
       "      <td>20.6</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>28.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>989.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>55.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>20.6</td>\n",
       "      <td>-85.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>8</th>\n",
       "      <td>-85.2</td>\n",
       "      <td>21.2</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>28.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>987.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>60.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>21.2</td>\n",
       "      <td>-85.2</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>9</th>\n",
       "      <td>-85.1</td>\n",
       "      <td>21.8</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>29.0</td>\n",
       "      <td>0</td>\n",
       "      <td>983.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>60.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>21.8</td>\n",
       "      <td>-85.1</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>10</th>\n",
       "      <td>-85.0</td>\n",
       "      <td>22.6</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>29.0</td>\n",
       "      <td>600</td>\n",
       "      <td>981.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>65.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>22.6</td>\n",
       "      <td>-85.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>11</th>\n",
       "      <td>-84.8</td>\n",
       "      <td>23.8</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>29.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>977.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>70.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>23.8</td>\n",
       "      <td>-84.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>12</th>\n",
       "      <td>-84.8</td>\n",
       "      <td>25.3</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>29.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>974.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>80.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>25.3</td>\n",
       "      <td>-84.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>13</th>\n",
       "      <td>-84.7</td>\n",
       "      <td>26.9</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>30.0</td>\n",
       "      <td>0</td>\n",
       "      <td>965.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>90.0</td>\n",
       "      <td>2.0</td>\n",
       "      <td>26.9</td>\n",
       "      <td>-84.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>14</th>\n",
       "      <td>-84.5</td>\n",
       "      <td>28.3</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>30.0</td>\n",
       "      <td>600</td>\n",
       "      <td>944.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>105.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>28.3</td>\n",
       "      <td>-84.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>15</th>\n",
       "      <td>-83.6</td>\n",
       "      <td>30.0</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>30.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>955.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>100.0</td>\n",
       "      <td>3.0</td>\n",
       "      <td>30.0</td>\n",
       "      <td>-83.6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>16</th>\n",
       "      <td>-82.5</td>\n",
       "      <td>31.5</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>30.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>981.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>HU</td>\n",
       "      <td>65.0</td>\n",
       "      <td>1.0</td>\n",
       "      <td>31.5</td>\n",
       "      <td>-82.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>17</th>\n",
       "      <td>-80.9</td>\n",
       "      <td>32.7</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>31.0</td>\n",
       "      <td>0</td>\n",
       "      <td>985.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>55.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.7</td>\n",
       "      <td>-80.9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>18</th>\n",
       "      <td>-78.8</td>\n",
       "      <td>33.6</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>31.0</td>\n",
       "      <td>600</td>\n",
       "      <td>991.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>33.6</td>\n",
       "      <td>-78.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>19</th>\n",
       "      <td>-77.0</td>\n",
       "      <td>33.7</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>31.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>990.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>TS</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>33.7</td>\n",
       "      <td>-77.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>20</th>\n",
       "      <td>-74.4</td>\n",
       "      <td>33.6</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023083e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>8</td>\n",
       "      <td>31.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>994.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>55.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>33.6</td>\n",
       "      <td>-74.4</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>21</th>\n",
       "      <td>-72.3</td>\n",
       "      <td>33.4</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>1.0</td>\n",
       "      <td>0</td>\n",
       "      <td>996.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>55.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>33.4</td>\n",
       "      <td>-72.3</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>22</th>\n",
       "      <td>-70.5</td>\n",
       "      <td>32.8</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>1.0</td>\n",
       "      <td>600</td>\n",
       "      <td>998.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.8</td>\n",
       "      <td>-70.5</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>23</th>\n",
       "      <td>-68.8</td>\n",
       "      <td>32.3</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>998.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.3</td>\n",
       "      <td>-68.8</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>24</th>\n",
       "      <td>-67.0</td>\n",
       "      <td>32.1</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>1.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>1000.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>45.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>32.1</td>\n",
       "      <td>-67.0</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>25</th>\n",
       "      <td>-65.6</td>\n",
       "      <td>31.9</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>2.0</td>\n",
       "      <td>0</td>\n",
       "      <td>1000.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>45.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.9</td>\n",
       "      <td>-65.6</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>26</th>\n",
       "      <td>-64.7</td>\n",
       "      <td>31.7</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>2.0</td>\n",
       "      <td>600</td>\n",
       "      <td>999.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.7</td>\n",
       "      <td>-64.7</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>27</th>\n",
       "      <td>-63.9</td>\n",
       "      <td>31.7</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>2.0</td>\n",
       "      <td>1200</td>\n",
       "      <td>999.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.7</td>\n",
       "      <td>-63.9</td>\n",
       "    </tr>\n",
       "    <tr>\n",
       "      <th>28</th>\n",
       "      <td>-63.2</td>\n",
       "      <td>31.9</td>\n",
       "      <td>IDALIA</td>\n",
       "      <td>2.023090e+09</td>\n",
       "      <td>2023.0</td>\n",
       "      <td>9</td>\n",
       "      <td>2.0</td>\n",
       "      <td>1800</td>\n",
       "      <td>997.0</td>\n",
       "      <td>al</td>\n",
       "      <td>10.0</td>\n",
       "      <td>EX</td>\n",
       "      <td>50.0</td>\n",
       "      <td>0.0</td>\n",
       "      <td>31.9</td>\n",
       "      <td>-63.2</td>\n",
       "    </tr>\n",
       "  </tbody>\n",
       "</table>\n",
       "</div>"
      ],
      "text/plain": [
       "    LONGITUDE  LATITUDE STORMNAME           DTG    YEAR  MONTH   DAY  HHMM  \\\n",
       "0       -86.1      21.3       TEN  2.023083e+09  2023.0      8  26.0  1800   \n",
       "1       -86.3      21.2       TEN  2.023083e+09  2023.0      8  27.0     0   \n",
       "2       -86.8      20.5       TEN  2.023083e+09  2023.0      8  27.0   600   \n",
       "3       -86.3      19.8       TEN  2.023083e+09  2023.0      8  27.0  1200   \n",
       "4       -85.7      20.0    IDALIA  2.023083e+09  2023.0      8  27.0  1800   \n",
       "5       -85.4      19.8    IDALIA  2.023083e+09  2023.0      8  28.0     0   \n",
       "6       -85.2      19.8    IDALIA  2.023083e+09  2023.0      8  28.0   600   \n",
       "7       -85.2      20.6    IDALIA  2.023083e+09  2023.0      8  28.0  1200   \n",
       "8       -85.2      21.2    IDALIA  2.023083e+09  2023.0      8  28.0  1800   \n",
       "9       -85.1      21.8    IDALIA  2.023083e+09  2023.0      8  29.0     0   \n",
       "10      -85.0      22.6    IDALIA  2.023083e+09  2023.0      8  29.0   600   \n",
       "11      -84.8      23.8    IDALIA  2.023083e+09  2023.0      8  29.0  1200   \n",
       "12      -84.8      25.3    IDALIA  2.023083e+09  2023.0      8  29.0  1800   \n",
       "13      -84.7      26.9    IDALIA  2.023083e+09  2023.0      8  30.0     0   \n",
       "14      -84.5      28.3    IDALIA  2.023083e+09  2023.0      8  30.0   600   \n",
       "15      -83.6      30.0    IDALIA  2.023083e+09  2023.0      8  30.0  1200   \n",
       "16      -82.5      31.5    IDALIA  2.023083e+09  2023.0      8  30.0  1800   \n",
       "17      -80.9      32.7    IDALIA  2.023083e+09  2023.0      8  31.0     0   \n",
       "18      -78.8      33.6    IDALIA  2.023083e+09  2023.0      8  31.0   600   \n",
       "19      -77.0      33.7    IDALIA  2.023083e+09  2023.0      8  31.0  1200   \n",
       "20      -74.4      33.6    IDALIA  2.023083e+09  2023.0      8  31.0  1800   \n",
       "21      -72.3      33.4    IDALIA  2.023090e+09  2023.0      9   1.0     0   \n",
       "22      -70.5      32.8    IDALIA  2.023090e+09  2023.0      9   1.0   600   \n",
       "23      -68.8      32.3    IDALIA  2.023090e+09  2023.0      9   1.0  1200   \n",
       "24      -67.0      32.1    IDALIA  2.023090e+09  2023.0      9   1.0  1800   \n",
       "25      -65.6      31.9    IDALIA  2.023090e+09  2023.0      9   2.0     0   \n",
       "26      -64.7      31.7    IDALIA  2.023090e+09  2023.0      9   2.0   600   \n",
       "27      -63.9      31.7    IDALIA  2.023090e+09  2023.0      9   2.0  1200   \n",
       "28      -63.2      31.9    IDALIA  2.023090e+09  2023.0      9   2.0  1800   \n",
       "\n",
       "      MSLP BASIN  STORMNUM STORMTYPE  INTENSITY   SS   LAT   LON  \n",
       "0   1006.0    al      10.0        TD       25.0  0.0  21.3 -86.1  \n",
       "1   1005.0    al      10.0        TD       30.0  0.0  21.2 -86.3  \n",
       "2   1001.0    al      10.0        TD       30.0  0.0  20.5 -86.8  \n",
       "3   1001.0    al      10.0        TD       30.0  0.0  19.8 -86.3  \n",
       "4    996.0    al      10.0        TS       35.0  0.0  20.0 -85.7  \n",
       "5    992.0    al      10.0        TS       45.0  0.0  19.8 -85.4  \n",
       "6    988.0    al      10.0        TS       55.0  0.0  19.8 -85.2  \n",
       "7    989.0    al      10.0        TS       55.0  0.0  20.6 -85.2  \n",
       "8    987.0    al      10.0        TS       60.0  0.0  21.2 -85.2  \n",
       "9    983.0    al      10.0        TS       60.0  0.0  21.8 -85.1  \n",
       "10   981.0    al      10.0        HU       65.0  1.0  22.6 -85.0  \n",
       "11   977.0    al      10.0        HU       70.0  1.0  23.8 -84.8  \n",
       "12   974.0    al      10.0        HU       80.0  1.0  25.3 -84.8  \n",
       "13   965.0    al      10.0        HU       90.0  2.0  26.9 -84.7  \n",
       "14   944.0    al      10.0        HU      105.0  3.0  28.3 -84.5  \n",
       "15   955.0    al      10.0        HU      100.0  3.0  30.0 -83.6  \n",
       "16   981.0    al      10.0        HU       65.0  1.0  31.5 -82.5  \n",
       "17   985.0    al      10.0        TS       55.0  0.0  32.7 -80.9  \n",
       "18   991.0    al      10.0        TS       50.0  0.0  33.6 -78.8  \n",
       "19   990.0    al      10.0        TS       50.0  0.0  33.7 -77.0  \n",
       "20   994.0    al      10.0        EX       55.0  0.0  33.6 -74.4  \n",
       "21   996.0    al      10.0        EX       55.0  0.0  33.4 -72.3  \n",
       "22   998.0    al      10.0        EX       50.0  0.0  32.8 -70.5  \n",
       "23   998.0    al      10.0        EX       50.0  0.0  32.3 -68.8  \n",
       "24  1000.0    al      10.0        EX       45.0  0.0  32.1 -67.0  \n",
       "25  1000.0    al      10.0        EX       45.0  0.0  31.9 -65.6  \n",
       "26   999.0    al      10.0        EX       50.0  0.0  31.7 -64.7  \n",
       "27   999.0    al      10.0        EX       50.0  0.0  31.7 -63.9  \n",
       "28   997.0    al      10.0        EX       50.0  0.0  31.9 -63.2  "
      ]
     },
     "execution_count": 12,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "df = pd.read_csv('/vortexfs1/home/csherwood/proj/NOPP_Idalia/idalia_track.csv')\n",
    "df"
   ]
  }
 ],
 "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.9.16"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}