calipatria_client.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "id": "568f9a78",
   "metadata": {},
   "outputs": [],
   "source": [
    "from obspy import UTCDateTime\n",
    "from obspy import Stream\n",
    "from obspy import read\n",
    "from obspy.clients.fdsn import Client\n",
    "import os\n",
    "import json"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "id": "13391847",
   "metadata": {},
   "outputs": [],
   "source": [
    "t = UTCDateTime(2021, 6, 5, 10, 55, 58) # largest event time\n",
    "\n",
    "start_time = UTCDateTime(\"2021-06-05 00:00:00\")\n",
    "end_time = UTCDateTime(\"2021-06-06 00:00:00\")\n",
    "lon = -115.635\n",
    "lat = 33.140\n",
    "depth = 5.8\n",
    "\n",
    "# chan_priority_list=[\"HH[ZNE12]\", \"BH[ZNE12]\", \"EH[ZNE12]\", \"HN[ZNE12]\"]\n",
    "chan_priority_list=[\"HHZ\", \"HHN\", \"HHE\", \"HH1\", \"HH2\",\n",
    "                    \"BHZ\", \"BHN\", \"BHE\", \"BH1\", \"BH2\",\n",
    "                    \"EHZ\", \"EHN\", \"EHE\", \"EH1\", \"EH2\",\n",
    "                    \"HNZ\", \"HNN\", \"HNE\", \"HN1\", \"HN2\"]\n",
    "\n",
    "client1 = Client(\"IRIS\")\n",
    "client2 = Client(\"SCEDC\")\n",
    "\n",
    "IRIS_inventory = client1.get_stations(\n",
    "    starttime=start_time, endtime=end_time,\n",
    "    longitude=lon, latitude=lat, maxradius=0.5)\n",
    "\n",
    "SCEDC_inventory = client2.get_stations(\n",
    "    starttime=start_time, endtime=end_time,\n",
    "    longitude=lon, latitude=lat, maxradius=0.5)\n",
    "\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "id": "0213f269",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "Inventory created at 2022-09-22T00:22:25.680600Z\n",
      "\tCreated by: IRIS WEB SERVICE: fdsnws-station | version: 1.1.50\n",
      "\t\t    http://service.iris.edu/fdsnws/station/1/query?starttime=2021-06-05...\n",
      "\tSending institution: IRIS-DMC (IRIS-DMC)\n",
      "\tContains:\n",
      "\t\tNetworks (9):\n",
      "\t\t\tAZ, BC, CI, EN, NP, PB, SB, SY, ZY\n",
      "\t\tStations (96):\n",
      "\t\t\tAZ.TONN (Tonner Ranch, Anza, CA, USA)\n",
      "\t\t\tAZ.USGCB (US Gypsum Company, site B (borehole), CA, USA)\n",
      "\t\t\tBC.CEA0 (Planta de bombeo La Rosita)\n",
      "\t\t\tBC.ICBC (Instituto Cultural de Baja California)\n",
      "\t\t\tCI.BOM (Bombay Beach)\n",
      "\t\t\tCI.CLI2 (Calipatria 2)\n",
      "\t\t\tCI.COK2 (Cook Ranch 2)\n",
      "\t\t\tCI.CRR (Carrizo)\n",
      "\t\t\tCI.DRE (Desert Research Extended Center)\n",
      "\t\t\tCI.ERR (Elmore's Ranch)\n",
      "\t\t\tCI.FRK (Frink)\n",
      "\t\t\tCI.IMP (Imperial)\n",
      "\t\t\tCI.NSS2 (North Shore Salton Sea 2)\n",
      "\t\t\tCI.Q0027 (1001 West Kramer Road)\n",
      "\t\t\tCI.Q0082 (Netquake Q0082)\n",
      "\t\t\tCI.RXH (Rock Hill)\n",
      "\t\t\tCI.SAL (Salton City)\n",
      "\t\t\tCI.SGL (Signal Mtn)\n",
      "\t\t\tCI.SLB (Salton Beach)\n",
      "\t\t\tCI.SLV (Salvation Pass)\n",
      "\t\t\tCI.SNR (Schaffner Ranch)\n",
      "\t\t\tCI.SWS (Sam W. Stewart)\n",
      "\t\t\tCI.WI2 (Wister, Niland, CA, USA)\n",
      "\t\t\tCI.WMD (Westmorland)\n",
      "\t\t\tCI.WWF (Wister Waterfowl)\n",
      "\t\t\tEN.ELM (Elmore, , , )\n",
      "\t\t\tEN.ENG (English, , , )\n",
      "\t\t\tEN.HAT (Hatfield, , , )\n",
      "\t\t\tEN.LIN (Lindsey, , , )\n",
      "\t\t\tEN.OBS (Obsidian Butte, , , )\n",
      "\t\t\tEN.RED (Red Hill Marina, , , )\n",
      "\t\t\tEN.SIM (Simpson, , , )\n",
      "\t\t\tEN.YOU (Young, Calipatria, CA, USA)\n",
      "\t\t\tNP.117 (CA:El Centro;ECA9 Irr Dist SS)\n",
      "\t\t\tNP.286 (CA:Imp Vly;Superstition Mtn Pk)\n",
      "\t\t\tNP.412 (CA:El Centro;ECA10 Reg Hosp)\n",
      "\t\t\tNP.464 (CA:El Centro;Meadows Union Sch)\n",
      "\t\t\tNP.5028 (CA:Imperial;ECA7;Imp Vly Coll)\n",
      "\t\t\tNP.5051 (CA:Imp Vly;Parachute Test Site)\n",
      "\t\t\tNP.5053 (CA:Calexico;FS)\n",
      "\t\t\tNP.5055 (CA:Holtville;PO)\n",
      "\t\t\tNP.5056 (CA:Brawley;ECA1  Borchard Rnch)\n",
      "\t\t\tNP.5057 (CA:Holtville;ECA3 Pine Union S)\n",
      "\t\t\tNP.5058 (CA:El Centro;ECA11 McCabe Sch)\n",
      "\t\t\tNP.5060 (CA:Brawley;Arpt Hangar)\n",
      "\t\t\tNP.5061 (CA:Calipatria;FS)\n",
      "\t\t\tNP.5158 (CA:Imperial;ECA6 Huston Rd)\n",
      "\t\t\tNP.5165 (CA:El Centro;Differential Arr)\n",
      "\t\t\tNP.5210 (Imperial Wildlife)\n",
      "\t\t\tNP.5274 (CA:Imp Vly;Superstition Mtn Bs)\n",
      "\t\t\tNP.5375 (CA:Ocotillo Wells;Veh Rec Area)\n",
      "\t\t\tNP.5413 (CA:Imperial;Imp Co Landfill)\n",
      "\t\t\tNP.5490 (CA:Calexico;TB Hay Farm)\n",
      "\t\t\tNP.931 (CA:El Centro; ECA 12)\n",
      "\t\t\tNP.952 (CA:El Centro;ECA5 James Rd)\n",
      "\t\t\tNP.955 (CA:El Centro;ECA4 Anderson Rd)\n",
      "\t\t\tNP.SAAR (Shape Acceleration Array)\n",
      "\t\t\tPB.DHL1 (Durmid Hill1, CA, USA)\n",
      "\t\t\tPB.GLRS (GLRS_SCGN_CS1999, CA)\n",
      "\t\t\tPB.P494 (WestsideESCS2005, CA)\n",
      "\t\t\tPB.P505 (ImperialSpCS2006 , CA)\n",
      "\t\t\tPB.P506 (Ramer_LakeCS2005, CA)\n",
      "\t\t\tPB.SLMS (SLMS_SCGN_CS1999, CA)\n",
      "\t\t\tPB.USGC (USGC_SCGN_CS1998, CA)\n",
      "\t\t\tSB.WLA (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA01 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA02 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA03 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA04 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA05 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA06 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA09 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA10 (Wildlife Liquefaction Array)\n",
      "\t\t\tSY.BOM (BOM synthetic)\n",
      "\t\t\tSY.BTC (BTC synthetic)\n",
      "\t\t\tSY.CRR (CRR synthetic)\n",
      "\t\t\tSY.DRC (DRC synthetic)\n",
      "\t\t\tSY.DRE (DRE synthetic)\n",
      "\t\t\tSY.ERR (ERR synthetic)\n",
      "\t\t\tSY.FRK (FRK synthetic)\n",
      "\t\t\tSY.IMP (IMP synthetic)\n",
      "\t\t\tSY.NSS (NSS synthetic)\n",
      "\t\t\tSY.NSS2 (NSS2 synthetic)\n",
      "\t\t\tSY.PLR (PLR synthetic)\n",
      "\t\t\tSY.RXH (RXH synthetic)\n",
      "\t\t\tSY.SAL (SAL synthetic)\n",
      "\t\t\tSY.SGL (SGL synthetic)\n",
      "\t\t\tSY.SLB (SLB synthetic)\n",
      "\t\t\tSY.SMTC (SMTC synthetic)\n",
      "\t\t\tSY.SNR (SNR synthetic)\n",
      "\t\t\tSY.SSW (SSW synthetic)\n",
      "\t\t\tSY.SWS (SWS synthetic)\n",
      "\t\t\tSY.WES (WES synthetic)\n",
      "\t\t\tSY.WMD (WMD synthetic)\n",
      "\t\t\tSY.WWF (WWF synthetic)\n",
      "\t\t\tZY.COON (Bruce Coon Residence)\n",
      "\t\tChannels (0):\n",
      "\n",
      "Inventory created at 2022-09-22T00:22:29.000000Z\n",
      "\tCreated by: SCEDC WEB SERVICE: fdsnws-station | version: 1.1\n",
      "\t\t    http://service.scedc.caltech.edu/fdsnws/station/1/query?start=2021-...\n",
      "\tSending institution: SCEDC (SCEDC)\n",
      "\tContains:\n",
      "\t\tNetworks (8):\n",
      "\t\t\tAZ, BC, CE, CI, EN, NP, PB, SB\n",
      "\t\tStations (77):\n",
      "\t\t\tAZ.TONN (Tonner Ranch, Anza, CA, USA)\n",
      "\t\t\tBC.CEA0 (Planta de bombeo La Rosita)\n",
      "\t\t\tCE.01711 (El Centro - Imperial & Ross)\n",
      "\t\t\tCE.01794 (El Centro - Meloland Geotechnical Array)\n",
      "\t\t\tCE.11023 (Niland - Fire Station)\n",
      "\t\t\tCE.11217 (North Shore - Hwy 111 & Parkside Dr)\n",
      "\t\t\tCE.11369 (Westmorland - Fire Station)\n",
      "\t\t\tCE.11591 (North Shore - Durmid)\n",
      "\t\t\tCE.11627 (Bombay Beach - Bertram)\n",
      "\t\t\tCE.11628 (Salton City)\n",
      "\t\t\tCE.11684 (Frink)\n",
      "\t\t\tCE.12626 (Desert Shores)\n",
      "\t\t\tCI.BOM (Bombay Beach)\n",
      "\t\t\tCI.CLI2 (Calipatria 2)\n",
      "\t\t\tCI.COK2 (Cook Ranch 2)\n",
      "\t\t\tCI.CRR (Carrizo)\n",
      "\t\t\tCI.DRE (Desert Research Extended Center)\n",
      "\t\t\tCI.ERR (Elmore's Ranch)\n",
      "\t\t\tCI.FRK (Frink)\n",
      "\t\t\tCI.IMP (Imperial)\n",
      "\t\t\tCI.IRG5 (Time Code)\n",
      "\t\t\tCI.NSS2 (North Shore Salton Sea 2)\n",
      "\t\t\tCI.Q0027 (1001 West Kramer Road)\n",
      "\t\t\tCI.Q0082 (Netquake Q0082)\n",
      "\t\t\tCI.RXH (Rock Hill)\n",
      "\t\t\tCI.SAL (Salton City)\n",
      "\t\t\tCI.SGL (Signal Mtn)\n",
      "\t\t\tCI.SLB (Salton Beach)\n",
      "\t\t\tCI.SLV (Salvation Pass)\n",
      "\t\t\tCI.SNR (Schaffner Ranch)\n",
      "\t\t\tCI.SWS (Sam W. Stewart)\n",
      "\t\t\tCI.WMD (Westmorland)\n",
      "\t\t\tCI.WWF (Wister Waterfowl)\n",
      "\t\t\tEN.ELM (Elmore, , , )\n",
      "\t\t\tEN.ENG (English, , , )\n",
      "\t\t\tEN.HAT (Hatfield, , , )\n",
      "\t\t\tEN.LIN (Lindsey, , , )\n",
      "\t\t\tEN.OBS (Obsidian Butte, , , )\n",
      "\t\t\tEN.RED (Red Hill Marina, , , )\n",
      "\t\t\tEN.SIM (Simpson, , , )\n",
      "\t\t\tEN.YOU (Young, Calipatria, CA, USA)\n",
      "\t\t\tNP.117 (CA:El Centro;ECA9 Irr Dist SS)\n",
      "\t\t\tNP.286 (Imperial Valley, Superstition Mtn Peak)\n",
      "\t\t\tNP.412 (CA:El Centro;ECA10 Reg Hosp)\n",
      "\t\t\tNP.464 (CA:El Centro;Meadows Union Sch)\n",
      "\t\t\tNP.5028 (CA:Imperial;ECA7;Imp Vly Coll)\n",
      "\t\t\tNP.5051 (CA:Imp Vly;Parachute Test Site)\n",
      "\t\t\tNP.5053 (CA:Calexico;FS)\n",
      "\t\t\tNP.5055 (CA:Holtville;PO)\n",
      "\t\t\tNP.5056 (CA:Brawley;ECA1  Borchard Rnch)\n",
      "\t\t\tNP.5057 (CA:Holtville;ECA3 Pine Union S)\n",
      "\t\t\tNP.5058 (CA:El Centro;ECA11 McCabe Sch)\n",
      "\t\t\tNP.5060 (CA:Brawley;Arpt Hangar)\n",
      "\t\t\tNP.5061 (CA:Calipatria;FS)\n",
      "\t\t\tNP.5062 (Salton Sea Wildlife (SSW))\n",
      "\t\t\tNP.5158 (CA:Imperial;ECA6 Huston Rd)\n",
      "\t\t\tNP.5165 (CA:El Centro;Differential Arr)\n",
      "\t\t\tNP.5271 (5271, Bombay Beach, CA (BBB))\n",
      "\t\t\tNP.5274 (CA:Imp Vly;Superstition Mtn Bs)\n",
      "\t\t\tNP.5375 (CA:Ocotillo Wells;Veh Rec Area)\n",
      "\t\t\tNP.5413 (CA:Imperial;Imp Co Landfill)\n",
      "\t\t\tNP.5490 (CA:Calexico;TB Hay Farm)\n",
      "\t\t\tNP.931 (CA:El Centro; ECA 12)\n",
      "\t\t\tNP.952 (CA:El Centro;ECA5 James Rd)\n",
      "\t\t\tNP.955 (CA:El Centro;ECA4 Anderson Rd)\n",
      "\t\t\tPB.DHL2 (durmidhillcs2005, CA, USA)\n",
      "\t\t\tPB.SCS1 (saltncity1cs2006, Salton Sea, CA, USA)\n",
      "\t\t\tPB.SCS2 (saltncity2cs2006, Salton Sea, CA, USA)\n",
      "\t\t\tSB.WLA (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA01 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA02 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA03 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA04 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA05 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA06 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA09 (Wildlife Liquefaction Array)\n",
      "\t\t\tSB.WLA10 (Wildlife Liquefaction Array)\n",
      "\t\tChannels (0):\n",
      "\n"
     ]
    }
   ],
   "source": [
    "print(IRIS_inventory)\n",
    "print(SCEDC_inventory)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e90fe10e",
   "metadata": {},
   "outputs": [],
   "source": [
    "for ev in IRIS_inventory:\n",
    "    print(ev.code)\n",
    "    for st in ev:\n",
    "        print(st.channels)\n",
    "        print(st.code)\n",
    "        print(st.latitude)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 4,
   "id": "d0896b38",
   "metadata": {},
   "outputs": [],
   "source": [
    "IRIS_contents = IRIS_inventory.get_contents()\n",
    "SCEDC_contents = SCEDC_inventory.get_contents()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "9c122056",
   "metadata": {},
   "outputs": [],
   "source": [
    "print(IRIS_contents[\"networks\"])\n",
    "print(len(IRIS_contents[\"stations\"]))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 5,
   "id": "f631e773",
   "metadata": {},
   "outputs": [],
   "source": [
    "IRIS_waveforms = []\n",
    "SCEDC_waveforms = []"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "17e260c9",
   "metadata": {},
   "outputs": [],
   "source": [
    "l = IRIS_contents['stations'][0].split()\n",
    "print(l)\n",
    "net = l[0].split('.')[0]\n",
    "sta = l[0].split('.')[1]"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "c2de83b8",
   "metadata": {},
   "outputs": [],
   "source": [
    "st = client1.get_waveforms(net, sta, \"*\", \"*\", start_time, end_time)\n",
    "st.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "3a945738",
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "start_time_file = (f\"{start_time.year}{str(start_time.month).zfill(2)}{str(start_time.day).zfill(2)}T{str(start_time.hour).zfill(2)}{str(start_time.minute).zfill(2)}{str(start_time.second).zfill(2)}Z\")\n",
    "end_time_file = (f\"{end_time.year}{str(end_time.month).zfill(2)}{str(end_time.day).zfill(2)}T{str(end_time.hour).zfill(2)}{str(end_time.minute).zfill(2)}{str(end_time.second).zfill(2)}Z\")\n",
    "for s in st:\n",
    "#     print(s.stats.station, s.stats.channel)\n",
    "    if s.stats.channel in chan_priority_list:\n",
    "        s.write(f\"waveforms/{s.stats.network}.{s.stats.station}.{s.stats.channel}__{start_time_file}__{end_time_file}.mseed\", format=\"MSEED\")\n",
    "#     else:\n",
    "#         st.remove(s)   \n",
    "#         print(st)\n",
    "# st[0].stats.channel"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "e7b2d9a2",
   "metadata": {},
   "outputs": [],
   "source": [
    "l = IRIS_contents['stations'][1].split()\n",
    "print(l)\n",
    "net = l[0].split('.')[0]\n",
    "sta = l[0].split('.')[1]\n",
    "\n",
    "st = client1.get_waveforms(net, sta, \"*\", \"*\", start_time, end_time)\n",
    "st.plot()"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "999c8fd9",
   "metadata": {},
   "outputs": [],
   "source": [
    "for i in range(len(IRIS_waveforms)):\n",
    "    print(IRIS_waveforms[i])"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "id": "b7f775ba",
   "metadata": {},
   "outputs": [],
   "source": [
    "start_time_file = (f\"{start_time.year}{str(start_time.month).zfill(2)}{str(start_time.day).zfill(2)}T{str(start_time.hour).zfill(2)}{str(start_time.minute).zfill(2)}{str(start_time.second).zfill(2)}Z\")\n",
    "end_time_file = (f\"{end_time.year}{str(end_time.month).zfill(2)}{str(end_time.day).zfill(2)}T{str(end_time.hour).zfill(2)}{str(end_time.minute).zfill(2)}{str(end_time.second).zfill(2)}Z\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "id": "55f1beb9",
   "metadata": {
    "scrolled": true
   },
   "outputs": [],
   "source": [
    "for i in IRIS_contents[\"stations\"]:\n",
    "    print(i)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "id": "ed54c117",
   "metadata": {},
   "outputs": [],
   "source": [
    "out_station_json_file='station_list.json'\n",
    "\n",
    "station_list = {}\n",
    "\n",
    "for ev in IRIS_inventory:\n",
    "    for st in ev:\n",
    "        station_list[st.code] = {\"latitude\": st.latitude, \"longitude\": st.longitude}\n",
    "        \n",
    "with open(out_station_json_file, 'w') as fp:\n",
    "    json.dump(station_list, fp)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 10,
   "id": "677c7ee2",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No data available for BC.CEA0\n",
      "No data available for BC.ICBC\n",
      "No data available for CI.BOM\n",
      "No data available for CI.CLI2\n",
      "No data available for CI.COK2\n",
      "No data available for CI.CRR\n",
      "No data available for CI.DRE\n",
      "No data available for CI.ERR\n",
      "No data available for CI.FRK\n",
      "No data available for CI.IMP\n",
      "No data available for CI.NSS2\n",
      "No data available for CI.Q0027\n",
      "No data available for CI.Q0082\n",
      "No data available for CI.RXH\n",
      "No data available for CI.SAL\n",
      "No data available for CI.SGL\n",
      "No data available for CI.SLB\n",
      "No data available for CI.SLV\n",
      "No data available for CI.SNR\n",
      "No data available for CI.WI2\n",
      "No data available for CI.WMD\n",
      "No data available for CI.WWF\n",
      "No data available for EN.ELM\n",
      "No data available for EN.ENG\n",
      "No data available for EN.HAT\n",
      "No data available for EN.LIN\n",
      "No data available for EN.OBS\n",
      "No data available for EN.RED\n",
      "No data available for EN.SIM\n",
      "No data available for EN.YOU\n",
      "No data available for NP.117\n",
      "No data available for NP.412\n",
      "No data available for NP.464\n",
      "No data available for NP.5028\n",
      "No data available for NP.5051\n",
      "No data available for NP.5053\n",
      "No data available for NP.5055\n",
      "No data available for NP.5056\n",
      "No data available for NP.5057\n",
      "No data available for NP.5058\n",
      "No data available for NP.5060\n",
      "No data available for NP.5061\n",
      "No data available for NP.5158\n",
      "No data available for NP.5165\n",
      "No data available for NP.5210\n",
      "No data available for NP.5375\n",
      "No data available for NP.5413\n",
      "No data available for NP.5490\n",
      "No data available for NP.931\n",
      "No data available for NP.952\n",
      "No data available for NP.955\n",
      "No data available for NP.SAAR\n",
      "No data available for PB.DHL1\n",
      "No data available for PB.SLMS\n",
      "No data available for SB.WLA09\n",
      "No data available for SY.DRC\n",
      "No data available for SY.PLR\n",
      "No data available for SY.SMTC\n",
      "No data available for ZY.COON\n"
     ]
    }
   ],
   "source": [
    "os.mkdir(\"waveforms_IRIS\")\n",
    "\n",
    "st_list = []\n",
    "for i in IRIS_contents[\"stations\"]:\n",
    "    curr = i.split()\n",
    "    \n",
    "    network = curr[0].split('.')[0]\n",
    "    station = curr[0].split('.')[1]\n",
    "    try:\n",
    "        wave_form = client1.get_waveforms(network, station, \"*\", \"*\", start_time, end_time)\n",
    "\n",
    "    except:\n",
    "        print(f\"No data available for {network}.{station}\")\n",
    "        \n",
    "    for s in wave_form:\n",
    "\n",
    "        if s.stats.channel in chan_priority_list:\n",
    "            s.write(f\"waveforms_IRIS/{s.stats.network}.{s.stats.station}.{s.stats.channel}__{start_time_file}__{end_time_file}.mseed\", format=\"MSEED\")\n",
    "            \n",
    "            if s.stats.station not in st_list:\n",
    "                st_list.append(s.stats.station)\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 11,
   "id": "2382b48b",
   "metadata": {
    "scrolled": true
   },
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "54\n"
     ]
    }
   ],
   "source": [
    "waveform_IRIS_files = os.listdir(\"waveforms_IRIS\")\n",
    "print(len(waveform_IRIS_files))"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 12,
   "id": "8f42f2eb",
   "metadata": {},
   "outputs": [],
   "source": [
    "net_list = []\n",
    "sta_list = []"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 13,
   "id": "086ef2c7",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAvsAAADtCAYAAAAhkbrhAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAABAtUlEQVR4nO3deXxU9b3/8dcnIQRISAIIBGRH2RRxwV161brVrS51l1u0t1btr1rb2mLrRlvbq9dabWtbbVWquLS2ttartliXq6yKCiIom4DsAbKTBLJ8fn+ck+lJMpMMmBAyeT8fj/Ng5vv5nu/5njNnwmfOfM93zN0REREREZHUk9beHRARERERkbahZF9EREREJEUp2RcRERERSVFK9kVEREREUpSSfRERERGRFKVkX0REREQkRSnZFxERERFJUUr2RURERERSlJJ9EREREZEUpWRfRERERCRFKdkXaUdmlmlmj5nZejMrMbM3zGx8ozpTzWyrmRWa2T1mZmH5aDP7XzPbFsZnmFmvyHrXm9lCM6sxs6lJ9OVIM1tkZhVm9n9mNrRR/CozW2FmO8zsIzMbtSdtmdl0M9tpZuXhsiQs/36kbKeZVUee/zas8wUzWxn24flG+9vXzF4Mt7nMzD7fwv7GPa7JHItG7XQPj32ZmX1qZpc1ik8JX9/S8LXu2kxbn2n/zKyLmS02s48blXvYZv3x/H5zx2ZPmFk/M/uTmW0Jj+kLZjYkEk8zs/vNrDisc1MkdqyZvW5mRWa2ycx+FT1OZjbNzJaaWZ2ZXZpEXxIexzA+1czWha/ZQjPL25O2LHi/VkWO68th+W8jZdVmtit67M3sbDOba8F7fr2Z3dFomy+aWUEYn29mxzbTvzsbvVfKG8VvNLO1YVvTzSyzmbbi7k8Y62FmD4X92mJm32mmnTQz+4mZbQjPhXvNGry/Trbg/VVuwftrWCS2X3geFVrwfroiEhsWnsvlkeUKRKR57q5Fi5Z2WoAs4DZgEJAOfAtYFYmfCawFRgADgKXA1WHsKOA/gdywnT8Bj0TWPQ84G3gOmNpCPzKBdcDVQDfgbuDNSPwc4EPgIMCAkUDvPWxrehL9mQpMb1TWDygCvhDu7wzg8Uj8T8DvgB7hvhcCvRK039xxbbb/cdq6B3gJyAGOA4qBUWFsfNiPieHr9DrwwwTtfOb9A74JzAI+blTuQH4bn8vDgRuA/YCuwM+BVyPx64H3wv0cDWwCTg5jZ4T7lAX0Ad4EbouseyVwKjAPuLSFfrR0HL8B/B8wNDyXDwa67WFbbyTRn98CdzYquwz4fHiuDQE+AiZH4gcBXSLvvc2AJWj/TuC3CWKnAuuBA8Lz5jngp830NeH+AHcBM8PzfDCwAjg9Qd2vAouA/kBvYA7wtTC2X3junknwN+9mYG5k3SeBP4TH5iBgCzAujA0DqtryPNaiJRWXdu+AFi1a/r2ESVId0Cd8/jTwvUj8auC1BOueCiyOUz6dlpPr04GPIs+zgApgaPj8beDUJPehpbaS6U+8ZP9rwMuR5yOByjApyAZ2AgMi8TeB/0zQfsLj2lL/47S1CTg68vxx4Pbw8U+B30RiJwOfJGjnM+1fmFgtBc6iHZL9OPtzIFAWeT4XuCTy/IfAownW/SrwQpzyN2g5uW7uOKaHr9eBSe5DwrZ2oz9Nkv04de4Cfhmn3MLX04GsBOveSeJk/2fA3ZHnxwIbmulHwv0B3gW+EHl+C/Bkgrp/Aa6LPL8MmB0+PgeYH4mlhcf0wPD5dmBsJP4QcFf4eBhK9rVo2e1Fw3hE9i3HAlvcfXv4fBywOBJfRHC1K57jgCV7uN0G23H3HcAqYJyZpQOHAePDIQerzez26NfyybYVqXOzmW03szlm9rk97OMqoIbg6vyBQIm7b4rUjx0rMzvBzIoTtUXD49ps/8MhIP8bPu4F5CfbVhgbbmbdW3P/QncDPwF2xGkb4L1wWMV0M+uToE5ranw+tte5HD2Og4DuwEXhUJRlZnbtHrZV75cWDAd7xcwO2cM+N9nf8ByrAv4X+EV4HmJml5vZB43WvzQc9vK+mV3QKGaNHg+sH7ZkZr82s183qt/c/jRuK9Hr11LdeH87EsUbr9vVzDaGQ5PuN7MezfRBRNCYfZF9hpnlElzF+kGkOBsojTwvDcsar3sowRCKH+3h5htvJ7qt/kAXgmEHBwP/AVwCTNmDtgAeIBhWMAB4EHjBzAZ/xj42u013n+Xuec20Fe1fS239t7ufHalb6+4VSbZVGilvbI/3LxzTPYpgCEQ8nyMYtnIowXCORxPUaxXh6/nfwK2R4mTP5S8QDJv5+R5uvrljtT/BcKqRBFeJLwDuMLOT9qAtgO8SDF8aArwCvGRm8V7bhMzsawTvhT9Ey8NzrCfwJYLhT/XlT7l7NAn/EzCGYMjR94DHzGxiGJsJXG5mo8J+fTss7xG2db27Xx9pq7n9mQl808zywjH2U+rbiWMm8DUzG2Bm+wFfj9SdC4wys7PMLINgGE9mJD4T+J4F98KMJ3iN6mPbCIbEDSY4pw8l+PZMRJqhZF9kH2Bm3YC/AS+6ezQRKycYI1svJyyLrjsceAH4irsndTXUzJZEbnAbEmc70W1Vhs/vdvdid/+UIEk/cw/awt3fd/cid9/l7k8S/Od/ahLdbq7dZreZRFvRurvTVjmQ3ujqYnNt5UTKW+pTtK2EMTNLA34B3OjuHqdd3P0td692960EHwrPtGZuFK6X6KbPFtbpDfwD+Im7/ysSSuZcPhJ4DDjP3bckub3GfUzmXJ7m7pXh++UJ/n0u705buPvb7l4etnVPWH5UMv0Ot3c2cDtwtrtXNo6H75G/AN82s7Hx2nD3pe6+2d1r3H0mwRC1c8PYP4H7gJeB5QRj52sIxsHHa6u5/fkxsIbg/oKXgWeBDQl27fdhnbeBdwjuadkQbmMbcFHY3maCb1uWRNq6gSC5X0vwofTpyLrl7v6eu9e6+1qC4X7nJ+iDiISU7Iu0MzPrAjwDbAQaz3CxlOAmz3oTiHzdb2b5BFfgfuTuf0t2m+5+kLtnh8unjbdjZlkEVz+XuntR2LcG3d6TthJ0py7JbjdudwTBNw6fENwsmBsej3oNjlVzbTWqm3T/w2OzOdm2wtjqeIndZ9i/HOBwgm9INhPchHmAmW1OMMSh/ngnGoYV3b/61zWpq9XhVeCXgOfd/ZfN7R9Nz+WxwN8JbpSel8z2EvSxueO4HNjVuNt72FY8yZ7LhMPXHgHOdfeVLVTvQnDFPRkN+uDu97n7SHcfSDB0aqG71+5uW+6+w92vcfcB7j6W4D6CBfFWCpPxqe4+2N2HE3y4WBCJv+Luh7l7H4JvfwYQTACAu29194vcvZ+7Hwn0SrSdxvsqIgm01c0AWrRoSW4huJL5TyAjTuwsgqtpwwnGhn/Iv2eNySX4z/v2BO12IZhN5gmC/1C7AekJ6tbPQDMlfPxTGs6gcxfB2OGewECCccxT9rCtCwlueu1CMByoFBjWqI3mZuM5neDK3xM0nBnlWYJhUN0Jrmw2NxtPc8e12f7Haet/IsfmGILZeEaHsfEENxweHr5er9LybDy7tX8EyWp+ZLmA4MNBPv8e7zyB4ObUXsBTwEttcB53JRiC8UiC+NcJkra+BEOONgKfD2ODgU+BLydYNyM8f98EJoeP0/bwOD4JPBy+tqPDfpy0u20BeQTfSGWG+34TwQe/nEZtxJuN51CgADglzjaHEsyi1S1s+/8RvEf6J+jjueG5lUZwA3gJcEwY6w6MDc+DUQR/L85L0E6z+0NwBT4/PI9OJbjaPjBBW30IhkkZwbCb1cARjfa/C8FMPU8B90diI8PzNAO4lOAKf48wdhTB/StG8HfoVeDXrX0ua9GSaku7d0CLls68hP+xO8HwgvLIMilS5xaCsapFBNM8Wlj+5XDd6HrlkfXuDOPRZUozfTkS+CDsy5tEZp8J//P/XZhIbACm1fdjD9qaFbZTDMwnTPgard8k2Q/LzyS4WbaC4Cpwr0isL8FV5QqCK7inRGKTosemueOaRP+/T8MZWroTJJDlBB8SLm+0nSnhMSsjmIkoMxJbAlzxWfev0fZOJDIbD0ECuJzgxt3NBAlrvzY4l/8jPMd2NDonh4TxNOD+8HUvAL4VWfcOgqu00fWWROLTaXoun9hMX5o7jnkEs8WUEXzgu66F/YrbVvh6LAj7Wgi8BhweZ/14yf5jQG2j/X058jdhLkGCXwS8BfxHZN0rGh2bZ8LtlxGcs9EZj3qH59iOcF+vjdO33yazP8BJBOd3BcE9BJMatRX7u0XwAWNlWPdj4IJGdf8c7t92giFo0ffE5QTfBOwgmCL14EaxNWFsHfBLILu1z2UtWlJtqU8aREREREQkxWjMvoiIiIhIilKyLyIiIiKSopTsi4iIiIikKCX7IiIiIiIpSsm+iIiIiEiKUrIvIiIiIpKilOyLiIiIiKQoJfsiIiIiIilKyb6IiIiISIpSsi8iIiIikqKU7IuIiIiIpCgl+yIiIiIiKUrJvoiIiIhIilKyLyIiIiKSopTsi4iIiIikKCX7IiIiIiIpSsm+iIiIiEiKUrIvIiIiIpKilOyLiIiIiKQoJfsiIiIiIilKyb6IiIiISIpSsi8iIiIikqKU7IuIiIiIpCgl+yIiIiIiKUrJvoiIiIhIilKyLyIiIiKSopTsi4iIiIikKCX7IiIiIiIpSsm+iIiIiEiKUrIvIiIiIpKilOyLiIiIiKQoJfsiIiIiIimqS3t3oCMwMwMGAmXt3RcRERERkVBPYKO7e6IKSvaTMxBY396dEBERERFpZBCwIVFQyX5yygDWrVtHTk5Oe/dFRERERDq50tJSBg8eDC2MPFGyvxtycnKU7IuIiIhIh6Fkvw2UlZWxYMECSktL27srkqTu3btzxBFH0KdPn/buioiIiEirUbLfyh577DEeeughampq2rsrspvMjPPPP5+pU6eSlqaJqkRERKTjU7Lfit544w0efPBBJk+ezKWXXkq/fv0IJvKRfZm7U1xczEsvvcQDDzzA4MGDmTx5cnt3S0REROQz0+XLVvTyyy8zbtw4brzxRvr3769Ev4MwM3r16sUVV1zBKaecwssvv9zeXRIRkb2soqKCo446ih/84Aft3RWRVqVkvxWtXr2a8ePHt3c35DOYMGECa9asae9uiIjIXjZjxgzeeecdfvKTn7R3V0RalZL9VlRbW0tGRkZ7d0M+g4yMDN1vISLSCdXW1rZ3F0TahJL9veDNN9/kmGOOITc3lz59+nDSSSexevVq7rzzTjIyMsjOzqZnz54cccQRzJo1K24bf/jDHzj++OMblPXv35+vfe1rsecFBQWkpaVRWFjI9OnTOeOMM5rt18knn8zAgQOpq6trUL5161bOOussevTowejRo3n11VebrFtTU8P48eMZM2ZMs9uYPn06gwYNIicnh6uuuopdu3bFYsOGDaNHjx5kZ2eTnZ3Ntddem7CdYcOGMW/evAZl1157LXfeeScQ3C9hZtxyyy0N6nTr1i12pX7KlClkZmbGtpednc2UKVOa7b+IiHQO6enp7d0FkTahZL+NlZSUcN5553HLLbdQVFTE2rVrueGGG2J/VL7yla9QXl5OaWkp1113HRdffHHcdk444QQWLFhAVVUVACtWrCA3N5c5c+bE6syePZuxY8fSu3fvFvu1YcMG3nrrLSoqKnjttdcaxL7+9a8zcOBAtm3bxt13381FF11EUVFRgzq/+tWvyM3NbXYbixcv5lvf+hZ/+9vfWLduHWvWrOHHP/5xgzqvvfYa5eXllJeX89vf/rbFfjcnJyeH3/zmN2zfvj1hnWnTpsW2V15ezvTp0z/TNkVEJDUo2ZdUpWS/jS1fvpwePXrwxS9+kbS0NLKzszn//PMZMmRIg3pmxuWXX86mTZsoKSlp0s7IkSPp06cP77zzDgCzZs3i7LPPpqamhsLCwljZpEmTkurXU089xfHHH88ll1zCk08+GSsvLy/n+eef54c//CE9evTgvPPO4+CDD+aFF16I1dmyZQsPP/xwk6vo8bZxySWXMHHiRHJzc7ntttuYMWNGUv3bEwMGDOCss87iZz/7WZttQ0REUpOmXJZUpTO7jY0aNYqKigquueYa/vGPfyT8oa3a2lqefPJJRo8enfCK+fHHHx8b5jN79myOP/54jjnmGGbPnh0rO+GEE5Lq14wZM7jkkku45JJLeO6555p8YzBgwIBY3QkTJrBkyZLY8+9973t8//vfJysrq9ltLF26tMENyxMmTGD16tVUVlbGys477zz69+/P+eefz9q1a5Pqe3Nuu+02fv3rX8c+AImIiCRDV/YlVSnZb2O5ubm8+eabVFZWMmXKFPr27cuVV15JWVkZEPwIV15eHtnZ2XzjG99odhaAE044oUmyf9xxxzFr1iwqKyt57733krqy/+GHH7JkyRIuvPBCTjzxRLp37x67cl9eXk5OTk6D+jk5OZSXlwMwd+5cli9fzhVXXNHidhq3Vf+4vq2nnnqKNWvWsGLFCoYMGcJ5552Huyds79RTTyUvLy+2PPbYY03qjBkzhjPPPDPh1f1p06Y1aOPee+9tcT9ERCT1KdmXVKVkfy84+OCDeeKJJ9i8eTNz5sxhzpw53HXXXQBcddVVFBcXU1FRweuvv87VV1/N+++/H7edE044gTlz5rB161aqq6vJz8+PJftvv/02/fv3Z+jQoS32Z8aMGZx00kn069ePtLQ0vvSlL8WG8mRnZzf59qG0tJTs7Gzq6uq44YYbeOCBB5L6DYHGbdU/zs7OBuC4446jW7du5OTkcN9997FixQpWr16dsL1XXnmF4uLi2HLVVVfFrXf77bcnvLp/xx13NGjjO9/5Tov7ISIiqU/JvqQq/YLuXnbEEUdwwQUX8OGHHzJx4sRYuZlx7LHHcuCBB/Lqq69y2GGHNVn30EMPpaamhocffpjjjjsOgHHjxrFq1SpeffXVpK7quztPP/0027ZtIz8/H4DKykqqqqooLCzkwAMPpKSkhM2bN8fiixYt4r/+678oLS3lvffe45xzzgFg165dlJaWkp+fzyeffEKPHj0abGvcuHEsXrw49nzRokUMHz6c7t27N+mXmbXaj5CNGTOGM844g/vuu69V2hMRkdSnZF9Sla7st7GPP/6Yn//852zcuBEIbth94YUXOOqoo5rUfeedd1i6dCljx46N21Z6ejrHHHMMv/jFL2LTcJoZhx12GA899FCT8fp1dXVUVVXFll27dvHmm2+ybds2PvjgAxYuXMjChQtZtmwZw4YN489//jPZ2dmce+653HHHHVRWVvL3v/+dDz/8kHPOOYfc3Fw2bNgQW+/3v/89w4cPZ+HChXET+Msvv5w//elPvPfee5SUlHDXXXdx5ZVXAvDpp58yd+5cqqur2bFjBzfffDNDhw5l2LBhn+Vwx9x+++08+OCDmjNfRESSomRfUpWS/TbWs2dP5syZwxFHHEFWVhannHIKZ511FlOnTgXgkUceITs7m6ysLC666CJuv/12zjrrLACefPJJDjrooAbtTZo0iYKCgtiVfQiGwxQUFDRJ9l955RW6d+8eWw4//HCefPJJLr74YkaOHEl+fn5sueaaa2Iz5fz6179m3bp19OnTh+985zv86U9/olevXphZg3V69+5Neno6+fn5savy2dnZvPXWWwCMHz+en/3sZ5xzzjkMGjSIwYMHx36GvKysjGuuuYa8vDyGDRvGypUref7552OzIfzkJz/hC1/4wh4f97Fjx3L66ac3+ZGUO+64o8E8+6eddtoeb0NERFKHZuORVGXN3RApATPLAUpKSkqa3LwadeGFFzJp0iS++c1v7rW+Set67rnn+O///m/efvvt9u6KiIjsRS+88ALnnnsuEHwz3lpDS0XaSmlpaf0MjrnuHn+6R3RlX0RERKTBMJ7Gvywv0pEp2W9l+qakY9PrJyLSOUWT/cZDQEU6MiX7rSgrK4uioqL27oZ8BoWFhS3+WJiIiKQeXdmXVNUmyb6ZTTOzpWZWZ2aXNoodY2bzzKzczNab2cWR2JFmtsjMKszs/8xsaCTW3cxmmFmZmX1qZpc1andK2F6pmT1mZl0jsZFmNjts9z0zm9AW+33UUUcxa9Ystm/f3hbNSxurrKxk5syZcWdKEhGR1Ba9QVdX9iWVtNU8+yuAG4EfRQvNbADwF+Aa4J9ALpAXxjKB54A7gKeAacATwOfC1acBvYH9gYOBl8zsXXdfbmbjgfuA08Jt/w24Fbg9XPdp4EXg88B/AX81s1Hu3qrzMl544YW8+OKLXHHFFZx22mn07dtXd/d3AO5OcXExr732Gtu2bWPatGnt3SUREdnLNIxHUlWbzsZjZm8Av3X3Z8Ln9wDV7v6DOHVPB+5397Hh8yxgKzDW3dea2SbgPHefH8YfB1a6+w/N7KdAnrtfF8ZOBn7v7iPMbDTwNtDX3XeF8bXAZHd/M8n9SGo2HoD169fzxBNPMHfuXEpKSpJpXvYBPXr04Mgjj+Syyy5L+DsHIiKSumbNmhX7ccrCwkJ69erVzj0SaV6ys/Hs7V/QPRKYbWZLCK7S/wu4wd2LgHFA7OdW3X2Hma0CxplZKZAfjQOLgPrxFuMIvimIxoabWfcwtqw+0Q99ABwExE32w28ZMiNFPZPdwUGDBnHLLbckW11ERET2AbqyL6lqb48x2R+4AjgfOIDgw8b9YSwbaPyppDQszwZq3b0iTizeuqWR8ubaTeQWoCSyrG+mroiIiHRw0WG3ukFXUsneTvYrgcfcfbm77yAY039mGCsHGo+RyQnLy4F0M+sRJxZv3ZxIeXPtJvJTgvsJ6pdBze+WiIiIdGTRH9HSlX1JJXs72f+w0fPoz9MtBcbHAsGY/ZHA0nCYz+ZoHJgALIm3bhhb7e6VYWy0mWVE4odE1m3C3Xe6e2n9ApQls3MiIiLS8SnZl1TSVlNvZphZt7D9DDPrZmZpwHTgKjMbEY6nv4VglhyAN4DscArNTILZdBa4+9owPgO4zcx6mtkxwLnAH8PYU8DFZna4meUCPwjr4+7LgGXAVDPLNLPrgVpgTlvsu4iIiHRsSvYllbTVlf3fEQzZmQQ8Hj7+nLu/AvwcmE0wDj4NuAmCq+nABcC3gGLgeGBypM3bCcbPbwKeBa4PE3ncfTHwbeCFsN11wF2RdS8Hzgjb/SpwQWtPuykiIiKpQcm+pJI2nXozVezO1JsiIiLS8bz99tscffTRAKxYsYIDDjignXsk0rxkp97ULz6JiIiIROjKvqQSJfsiIiIiEUr2JZUo2RcRERGJULIvqUTJvoiIiEiEkn1JJUr2RURERCL0C7qSSpTsi4iIiEToyr6kEiX7IiIiIhFK9iWVKNkXERERiVCyL6lEyb6IiIhIhJJ9SSVK9kVEREQidIOupBIl+yIiIiIRurIvqUTJvoiIiEiEkn1JJUr2RURERCKU7EsqUbIvIiIiEqFkX1KJkn0RERGRiGeeeaa9uyDSapTsi4iIiEQ8+eST7d0FkVajZF9EREREJEW1S7JvZoeb2RwzKzWzT8zsqkhsqpltNbNCM7vHzCwSO9LMFplZhZn9n5kNjcS6m9kMMyszs0/N7LJG25xiZuvDbT5mZl33zt6KiIiIiLSP9rqy/zjwIpAHfAn4hZmNMrMzgeuAo4GDgLOBqwDMLBN4DngA6A3MA56ItDktLN8fuBT4jZmNCtcdD9wHnAcMBoYBt7bh/omIiIiItLv2SvaHAk+5e527vwd8BIwGJgO/dvdP3H0TcC9wZbjOiUC5uz/q7lXAD4GJkav7k4Fp7l7q7nOAvxMk/QCXA3909wXuXgL8KNJuE2aWaWY59QvQsxX3XURERERkr2ivZP9XwGQz62JmRxFcbZ8PjAMWR+otIrjCT+OYu+8AVgHjzKwXkJ/sumFsuJl1T9C/W4CSyLJ+d3dQRERERKS9tVey/w/gK0AVMAf4gbsXANlAaaReaVhGnFg0ng3UuntFkuuWRsrj+SmQG1kGJbVXIiIiIiL7kC57e4Nm1gd4AfgywVCbscA/zGwJUA7kRKrnhGXEiUXj5UC6mfWIJPzNrZsTKW/C3XcCOyN9Tnb3RERERET2Ge1xZX8EUOLuf3X3Wnf/EHgD+BywFBgfqTsBWBI+bhAzsyxgJLDU3YuAzcmuG8ZWu3tla+2UiIiIpA53b+8uiLSK9kj2lwM9zewcC4wBTiYYUz8DuM7MhptZPvCtsAyCDwTZ4RSamQSz6Sxw97VhfAZwm5n1NLNjgHOBP4axp4CLwyk/c4EfRNoVERERaWDXrl3t3QWRVrHXk/1wNpxLCGbEKQVmAr9y93+4+4vAQ8A7BDP0vAQ8Fq63E7iA4ANAMXA8wQw89W4nuJl2E/AscL27LwvXXQx8m2D40HpgHXBXW+6niIiIdFxVVVXt3QWRVmH6mqpl4fSbJSUlJeTkNL5tQERERDq6t99+m6OPPjr2fNOmTeTn57e4XnV1NRkZGW3ZNZG4SktLyc3NBch198aT2MS012w8IiIiIvusZK7sv/nmm+Tk5HD//fe3fYdE9pCSfREREelUKisrufnmm3nrrbcS1kkm2f/a175GVVUVN910U2t2T6RVKdkXERGRTuWXv/wl9957L5/73OcS1qmsbHnCvry8vFbslUjbULIvIiIincq6detij8vL4/7kTlJX9vfbb79W65NIW1GyLyIiIp1KeFMjAD179uTPf/5zkzpVVVVUV1c3e4Vfyb50BEr2RUREpFNpPLPeZZdd1qTOlVdeybhx4+jVqxfFxcVx24km+3V1da3St8mTJ3PhhRfqR72k1SjZFxERkU4lemUfoKampkmdjRs3snLlSnbu3MlLL70Ut50+ffrEHif6QLA7ampqmDFjBs899xxLliz5zO2JgJJ9ERER6UQefvhhHn/88Sblr732WsJ1Pvnkk7jlXbt2jT3etm1b3Druzvnnn89VV10FQEFBAU8//XTcDxhRzz//PAAPPfQQEydOZOPGjc3Wj/rGN77B5ZdfTm1tbdLrdGRr165l6tSpbN68OWGdnTt3cvrpp3PjjTe2eOzby65du3jooYdYuXJl6zbs7lpaWIAcwEtKSlxERGRfVVRU5BdddJG//PLL7d2VfdKGDRsc2O3lpJNO8nnz5nm/fv38sssu84svvtjPOOMMP+GEE2J1Zs+eHXebW7ZsidVZs2aNn3HGGQ74t7/97SZ1q6urG2w3+vy0005Lah+Liopi67z44oteU1PT4jq33Xab33XXXXFj69ev91NPPdX/8Y9/JLX9trBz506fMmWK/+53v4sbP+ussxzwAw88MGEbL7/8cuy4XHDBBW3V1QaKior84IMP9ltvvTWp+tOnT4/1MRklJSX19XO8uTy2uaAWJfsi+6KPP/7YCwoKvK6urr270il9+umnXlpamnT9NWvWeHl5eYOyuro6nzlzpr/00ku+YcOGWHlVVZU/+OCDPmvWLH/22We9qqpqt/tXW1vr//rXv3zLli27vW4iNTU1Xl1d3Wrt7Y66ujqfO3eul5WVNYlt3rzZH374YS8uLvZdu3b5z372swaJYry2Kisrvaamxmtra1vc9q5du3zTpk1J93Xjxo0+f/78pOpWVVX52rVrk247GatXr/aqqiq/4447fOLEib5mzRpft26dn3baaX7vvfd6fn7+HiX7ySz333+/z5s3z5944gk/8cQT/fHHH/fq6mrfuHFjrM5FF13UYJ1Fixb56tWr/cUXX/RDDjnE//nPfzaI//jHP27wfMeOHf7Xv/7VP/74Y6+rq/O6ujovLS1t8D459NBDY/WPPPJI7969u0+bNs0XLFjgO3bsiNWrq6vz7du3e2FhYYP+1CsoKPC1a9f69773vVi8/pypra314uLi2L9z5szx7373u37RRRf53Llzvby8PNa/6upqr66u9pkzZ/q2bduSeh3Ly8t9+vTpvnnzZv/Sl77U4Bhs377dlyxZ4lOnTvXXX3/d3d3NLBYvKyuL+/dp3rx5DdqJ935qrLq62v/5z3/6qlWrYu+dsrKypN47W7du9fvvvz+2vTlz5vjdd9/tp5xyihcUFMTq1dTUxN439R9aAC8vL/fnn3/eH3zwwQbbq62t9Z07d7p78sm+uesGkJaYWQ5QctRRR5GRkYGZUVhYyKZNmxg9ejTvvvsuBx98MJmZmWzdupXy8nKGDx9ev27s3+XLl7PffvvRu3dv0tLSSEtLY+3ataxdu5b+/fuzZcsWjj76aNLS0jAzlixZwtChQ8nKysLM2LJlC1VVVQwZMqRBu2bGRx99RP/+/cnJyYmVrV27lszMTPr379+g7qpVq8jJyaFPnz6xMiDWv7y8vCbtv//++4wZM4auXbtiZqxZs4asrCz69u3boN7KlSvJy8ujd+/eDcoBFi5cyOjRo8nMzGzQPsD8+fM54ogjYvXrl+bKJ06c2KAMgp87nzhxYpNtz58/nyOPPLJBubvzzjvvNKlfV1fHe++9x+GHH97kOGzdupWSkhKGDRsWK6usrGTVqlUcdNBBTV7zdevWkZ6eTn5+fpO2PvzwQ4YMGUJWVlaT/i5YsIAJEyaQnp7eoLx+36P7Uv/vrl27+PDDD5kwYUKD9ioqKli5ciUHH3xwg/KSkhI2btzI6NGjG7RTWFhIUVERI0aMaLCNgoICKisrGTJkSINygGXLltGvX7/YONhon99//33GjRsXe+9E161/verfD2lpacyaNYsJEyZQVVXFypUrmTBhAh988EGs/wsWLGjw/jz66KNZsGABxx13HN26dWPevHn069ePVatWceyxxzY5frW1tbz//vscfvjhDfpaXV3N4sWLOfTQQxuUV1ZWsmLFiibHr/71NTMGDBjQILZ69Wp69OjR5P0BsHjxYkaMGEH37t2bxN555x0OO+yw2N+A+vjbb7/d4H0QPRcavw/MjLlz5zJx4kTS09NJS0tjwYIFjB49mm7dulFSUkJhYSHbtm2jtraWcePG0bt3bwoKCigsLGTkyJG8++67HH744WRmZrJ+/XrKysooKSnhkEMO4d133yWe0aNHs3z5co455hiWLl3KgQce2OS1OvbYY5k/f/4e38jYs2dPevXqxaeffsr48eNZvHgxEAylyMvLo6CggEMPPZSFCxc2WXfEiBGxoRjdunXjsMMOY+7cuQ3qpKWlxfo2fPhw1q9fzzHHHNPgh5dGjBhBRUUFXbt2ZciQIXTr1o1//etfAGRmZtK7d2+GDh3KvHnzOOigg1iyZAlHHHEEH374IWPGjKFbt27Mnz8/1l5eXl5srPegQYMoKyvj6KOPZubMmc0ei6FDh7J27doWj1l+fn6zQxuihg0bxpo1axLGR44cSXFxMQceeCDp6enMnj2bo446ii5durB8+fImw1iGDBnCp59+So8ePaioqIiVH3300SxZsqTBlJfdu3dvMutN/frxnkePm0hr6NmzJ2VlZbu93sCBA3driFUbyXX30kRBJftJqE/227sfIiIiIiKNNJvsd9mbPenoZsyYQffu3amrq8Pd2bFjBxkZGXTv3p1t27aRn59PXV0d5eXlpKenx67c+b+HA1FUVERubm7sCrK7s3nzZqqrqxk+fDiFhYXst99+sXUKCwtjd/u7O+Xl5XTt2pXMzMwGbbs7RUVF7Lfffg3aLi0tJTs7m4yMjAb9KC4uJicnh/T09AbtFBUV0atXL9LS0mJ16//dvn07ffv2jdUtKSmhZ8+edOnSpUnbubm5pKenJ2wjWpao/frywsLCBus03t/GX1cVFxc3Ka+rq6OkpCRWXt9WXV0dpaWl9O7du0H92tpaysrKmpTXv+5du3alR48esbKdO3eyc+fO2JXtaP9LS0vJysqia9euDcqjr2/91droh+/o8Ui0742PY01NDRUVFfTq1atB/V27dsX6F61fVVVFbW0tPXv2bFBeUVGBmZGVldWgfMeOHXTp0iXhuZ3o3Nm2bRv9+vWL+7pv27aNvn37UldXF1u2b9/OgAEDKC0tJT09nZ49e7Jlyxb69u3LkiVLmDFjBsXFxUyePJlJkyZRU1NDQUEBOTk57Ny5k4KCAsaMGcO2bdviHqdEr2/j41dff9euXVRVVSV8fXv06BG7Ua++vKSkhKysLDIyMmKv6e68Dxpvp6X3R+PztLCwkH79+jU5pnV1dRQVFVFRUcHGjRsZPHgw/fr1o7a2NnauZmZmUlBQwH777UdNTQ1FRUX07t2bbdu2UVBQwFNPPcX27du55pprOOSQQ1i7di25ubn07duXsrIy8vLyWL9+PV27dmX+/Pmxtk499VQAXnnlFWbPns2qVasAGDNmDGeddRaHHnookydPJuroo4/m7LPPJiMjg759+7Jr1y5ycnLYunUrffr04ZNPPqF///5UVVUxcOBACgsLWb58OS+88AJdunThuOOOY//996d379706NGDVatW0b9/f/bff3/MjMWLF/P+++9TWlrKpEmTGDx4MGvXrmXYsGF07dqV6upqdu7cyaOPPkr37t254oorGDBgANu2bYu9PwoLC9myZQsjR46krq6OvLw86urq2Lx5M1lZWXTp0oXc3NzYObxp0yZee+01qqqqOP744xk8eDCrV6+mX79+dOnShaqqKkpKSnj33XcZO3Ysn3zyCUuWLGHChAmcc845rF69muzsbLp3787VV18NQP/+/dlvv/0azODy1a9+lXHjxpGbm8uiRYvIzs5mwYIFbN68mZycHHbt2sXIkSPJycmhb9++5ObmMnDgQN5++202bdrEli1bOOCAA9iyZQsnn3wyI0aMoKysDDOL/V9Y/39WTU0Na9as4emnnwbgiiuuoHfv3qxevTrW7pw5c+jTpw/Dhg0DiN2wevzxx3PaaadRV1dHZWUlO3fuJDMzM9ZudXU1gwYNoqioiAMOOIBly5aRl5dHly5dWLRoEaWlpQwYMIDVq1ezYcMGtm7dyscffwzABRdcwEsvvZTUD2S1hcbfTjSWlpZGXl4ehYWFu9324YcfHnsNevXqxahRo1ixYkXcb7bqHXXUUYwaNYrt27ezdOlSKioq2Lp1ayx+1llnUVFRwaBBg3jttdfYsGFDg/UnT57Mxo0bef/992N97tKlS5MbXkeOHMnIkSPZvHkz/fv3Z8CAAaxdu5bevXtz0kknMXfuXPr06UPPnj0pKSmhsrKSiooKysvLyc7OZuPGjaxfv55NmzY1ee1OPvnkBjdTH3DAAQ1uaB03bhzLly/n7LPP5qCDDmLs2LGsWLGCadOmNWjnV7/6FQsXLqSqqor8/HxGjx7N9u3bmTdvHu7O0KFD2bhxI8uWLSMjI4O8vDxef/113J2LL744NoIjPz+f/v37895775GZmcnYsWNZs2YNL730EnPmzGnQr6VLlwJw0kknMXHiRPLz83nggQfiniMZGRlUV1cDcNttt5GRkUFlZSXz589n06ZNTJo0iREjRjB16tSEr3dMc2N8tGjMvohIW4l3z0VNTY0/99xzvnHjxnboUcdTW1vrK1eujB3Lt956y7Ozs/1HP/pRO/esZfXjvdtCYWGhv/vuu7Hn69at89raWr/pppv2aDz+Bx980Gz8nHPOcfdgDHV5eblv2rQpNs561qxZsXqRMdZ+ww03NDgWCxYs8FtvvbVBu48++mjs8VNPPeXu3uJY86VLl8bWufHGGx1o9gbRZ555JlY/3nuyrq7OZ8+e3ex2a2tr/eabb/Yf/vCHzfZtT0ybNs27devmgA8dOjRWvmvXrtjNx/X9nzRpUsJ27rnnnli9N954o9X7GU9tba0/88wz/vHHH8eeR++baCz62rsH9wzUj8+PR2P2W1H9MJ6SkpImP8QhIiKyL6murm7wjZI09Nhjj8W+EUnWJ598EruPCeDll1+mrKyMiy++GIATTjihwb0djd17771kZWVx3XXXxb7JvfDCC5v8cu+zzz4ba3PIkCGsWbOGtLRglvQ777yTO+64I6n+PvrooxQXF/Otb32LhQsXxu4Xiae6uppTTz2VvLw8/va3vyXVfnuYOXMmw4cP58ADD2wSqz+m48aNS/j7BDU1NUydOpUjjjgi7o+o7QuysrJi97ckk5+XlpbWf+Oc6xqz/9ko2RcREUkd0QkGkrFq1SpGjhwJwIABA2I3ZCaTZDb2i1/8gmnTpvGvf/2Lww47rEFs6dKlsckeNmzYwMCBA2Pb+MMf/sB//ud/7la/O4sXXniBr3/96zz++OOceOKJ7d2dPXbBBRfw17/+FVCyv0fMrC8wHTgJWAdc7+6vJrmukn0REZEUEb2Cmoz33nuPgoICrr32Wh555BFOPvlk4N/JfvQDQDLcPe4HjpqaGjIyMujWrRulpaVkZGSwdOlS3njjDb761a/qG5sUt337dq666iq+8pWv8MUvfrHF+kr2GzGzPxHMqHMjcBrwKDDS3YuSWFfJvoiISIoYPHgw69evjz2/+uqrefTRR+PWHTVqFIsXL27wa7n16hP2zMzMVrsJuLi4mNra2tjkHCKJJJvsp+29LrUfM8sGvgjc7u4V7v434EPgnAT1M80sp34Beu693oqIiEhbiibS7s4jjzwSt96ll17KRx99FDfRj9q5c2er9S0vL0+JvrSqTpHsAwcCJe6+KVK2CDgoQf1bCL4FqF/WJ6gnIiIiHUz9Dz8258tf/jI///nPYzfIinRUneUMzgYaf71RGpbH81MgN7IMaruuiYiIyN6UTLI/ffp08vPz90JvRNpWZ0n2ywnmyo/KCcubcPed7l5avwC7//vJIiIisk/6n//5H3r37s3tt9/+mdrJy8trnQ6JtKHO8gu6K4BcM8t3981h2QTg9+3YJxEREWkHw4cPZ+vWrZ95iM7jjz/Oueeeyz333NNKPRNpfZ1pNp5ngULgm8CpBNNwajYeERERaTIVZrL50Y4dO8jKymqLLok0K9nZeDrLlX2A64E/ANsJbri9OJlEX0RERCQRJfqyr+s0yb67bwXObO9+iIiIyL4tPT29vbsg0mo6yw26IiIiIklZsWJFe3dBpNUo2RcRERGJGD58eHt3QaTVKNkXEREREUlRSvZFRERERFKUkn0RERERkRSlZF9EREREJEUp2RcRERERSVFK9kVEREREUpSSfRERERGRFKVkX0RERDo9/WqupCol+yIiItLppaUpJZLUpDNbREREOj1d2ZdUpWRfREREOj0l+5KqlOyLiIhIp6dkX1KVkn0RERHp9JTsS6pSsi8iIiKdnpJ9SVV7Pdk3s8vMbLmZlZrZAjM7JhLrbmYzzKzMzD41s8sarTvFzNaH6z5mZl0jsZFmNtvMKszsPTObEImlmdn9ZlZsZlvM7Ka9s7ciIiLSESjZl1S1V5N9MxsAPAZcA+QCDwPPRqpMA3oD+wOXAr8xs1HhuuOB+4DzgMHAMODWyLpPAzPD9R8F/mpmXcLYtcDngFHhv981s5NbfQdFRESkQ1KyL6lqb1/ZHwhscfc33N2BGcAgM+sZxicD09y91N3nAH8nSPoBLgf+6O4L3L0E+BFwJYCZjQZGAz919yp3/xWQDhwXafdudy9w92XA7+rXjcfMMs0sp34BeiaqKyIiIh2fkn1JVXs72V8IrDGzU8wsHZgCzHP3MjPrBeQDiyP1FwEHhY/HxYkNN7PuYWyZu++KxD9oYd2DSOwWoCSyrE92B0VERKTjUbIvqapLy1Vaj7vXmtkfCa7YdyVIpE8Jw9lArbtXRFYpDcvr46WNYvXljWPJrJtNYj8lGDJUrydK+EVERFKWkn1JVa2a7JvZTIIx8fH8GHiHYJz9kcDHBOPv/zcchlMOpJtZj0jCnxOWE/6bE2kvJ1LeOJbMuuUk4O47gZ2R/UpUVURERFKAkn1JVa06jMfdT3P3bgmWHwOHAK+6+xJ3r3X3vwAOjHX3ImAzMD7S5ARgSfh4aZzYanevDGOjzSwjEj+khXWXICIiIoKSfUlde3vM/gLgZDMbbYEvAnnAyjA+A7jNzHqGU3KeC/wxjD0FXGxmh5tZLvCDsD7hTbfLgKnhzbXXA7XAnEi7N5tZ33B2n/8CnmzrnRUREZGOQcm+pKq9PWb/dTO7D/gn0AdYA1wWXtUHuB34PbAJKAKuDxN53H2xmX0beIFgGM5fgLsizV8O/AH4PsEQoQvcvSaM/QY4EFgB7AL+291fbav9FBERkY4lLU2/MyqpyYIZMKU54fSbJSUlJeTkNL41QERERDq6ww8/nPfffx8A5UbSEZSWlpKbmwuQ6+6NJ6qJ0cdYERER6fQ0jEdSlZJ9ERER6fSU7EuqUrIvIiIinZ6SfUlVSvZFRESk01OyL6lKyb6IiIh0ekr2JVUp2RcREZFOT8m+pCol+yIiItLpKdmXVKVkX0RERDo9JfuSqpTsi4iISKenZF9SlZJ9ERER6fSU7EuqUrIvIiIinZ6SfUlVSvZFRESk01OyL6lKyb6IiIh0ekr2JVUp2RcREZFOLy1NKZGkJp3ZIiIi0unpyr6kKiX7IiIi0ukp2ZdU1arJvpl1MbO/mNkGM3Mzy49T5yozW2FmO8zsIzMbFYlNMbP1ZlZqZo+ZWddIbKSZzTazCjN7z8wmRGJpZna/mRWb2RYzu6nRNr9gZivDbT5vZr1ac79FRESkY1OyL6mqLa7svwlcGC9gZucA3wbOA7KBs4FtYWw8cF8YGwwMA26NrP40MBPoDTwK/NXMuoSxa4HPAaPCf79rZieH7fYDngK+AfQDyoAHWmE/RUREJEUo2ZdU1arJvrvXuPsD7j4vQZXbgJvcfYkHVrl7YRi7HPijuy9w9xLgR8CVAGY2GhgN/NTdq9z9V0A6cFy47mTgbncvcPdlwO/q1wXOB+a5+8vuvgO4A7jIzDJbc99FRESk41KyL6lqr43ZN7N04DBgfDhUZ7WZ3W5mFlYZByyOrLIIGG5m3cPYMnffFYl/ABzUzLpxY+6+CqgBRjTT10wzy6lfgJ67ubsiIiLSgSjZl1S1N2/Q7Q90AT4PHAz8B3AJMCWMZwOlkfqlkfLGsfp4djPrJoo1jsdzC1ASWdY3U1dEREQ6OCX7kqp2K9k3s5lmVpVgubWF1SvDf+9292J3/xR4EDgzLC8HciL1cyLljWP18fJm1k0UaxyP56dAbmQZ1ExdERER6eCU7Euq6tJylX9z99P2dEPuXmRmGxsVW+TxUmB85PkEYLW7V5rZUmC0mWW4e3UYPwT4n0brLo2suyQSOy+2QbMRBPv9STN93QnsjKzT4v6JiIhIx6VkX1JVqw/jCce7dwufRh8DTCeYKaenmQ0kmEXnxTD2FHCxmR1uZrnAD4AZAOFNt8uAqWH71wO1wJxw3RnAzWbWN5zK87+AJ8PYX4FjzOx0M+sBTAOeDRN6ERERESX7krLaYsz+Mv49ZGdN5DEEifYmgjHw7wDPAX8AcPfFBNNyvhDG1wF3Rda9HDgDKAa+Clzg7jVh7DfALGBF+O+97v5q2G4BcAXwa4JpPnOBG1tpX0VERCQFKNmXVLVbw3iS4e7DmontIkjUv5ogPp3g6n+82Erg+ASxOuCb4RIv/hLwUqJ+iYiISOemZF9S1d6cjUdERERkn6RkX1KVkn0RERHp9NLSlBJJatKZLSIiIp2eruxLqlKyLyIiIp2ekn1JVUr2RUREpNNTsi+pSsm+iIiIdHpK9iVVKdkXERGRTm/ChAnt3QWRNtHq8+yLiIiIdDQnnXQSTz/9NGPGjGnvroi0KnP39u7DPs/McoCSkpIScnJy2rs7IiIiItLJlZaWkpubC5Dr7qWJ6mkYj4iIiIhIilKyLyIiIiKSopTsi4iIiIikKN2guxtKSxMOhxIRERER2WuSzUt1g24SzGx/YH1790NEREREpJFB7r4hUVDJfhLMzICBQFl790WkGT0JPpQOQueq7Lt0nkpHoPNUOoqewEZvJqHXMJ4khAcw4ScmkX1B8JkUgLLmpuASaU86T6Uj0HkqHUiL56du0BURERERSVFK9kVEREREUpSSfZHUsROYFv4rsq/SeSodgc5TSRm6QVdEREREJEXpyr6IiIiISIpSsi8iIiIikqKU7IuIiIiIpCgl+yIiIiIiKUrJvoiIiIhIilKyL52WmWWa2WNmtt7MSszsDTMbH4lPNbOtZlZoZvdY+JOKZjbazP7XzLaF8Rlm1iuy3vVmttDMasxsahL9ONLMFplZhZn9n5kNbRS/ysxWmNkOM/vIzEbtSVtmNt3MdppZebgsCcu/HynbaWbVkee/Det8wcxWhn14Prq/kWO1zszKwn3Pa6aPcY9rMsdC9g1mNs3MlppZnZldGimfEr7+ZWb2iZld20I7zZ2v3cP3VpmZfWpml7XlPknqaeY8zTSzR8K/Q9vM7Akzy2qmHZ2n0qEp2ZfOrAvwCXAM0Bv4O/A3ADM7E7gOOBo4CDgbuCpcLxf4EzASGAZ0Be6NtLsRuDVsr1lmlgk8BzwQ9mEe8EQkfg7wbeA8IDvsx7Y9aSt0h7tnh8tBAO7+k/oy4A7gyUida82sH/AU8A2gH1AWbqN+u98AvgCcAOQAVwJVCfqY8Lgm2X/ZN6wAbgTeblSeCVwL9ALOAaaZ2efiNZDE6z0tLN8fuBT4TXMfdEXiSHSefgMYD4wGhgP9gbgXZnSeSkpwdy1atLhDkLTXAX2Ap4HvRWJXA68lWO9UYHGc8unA1Ba2eTrwUeR5FlABDA2fvw2cmmT/W2ormf5MBaY3Kvsa8HLk+UigkiCxSwc2AQcm2ceEx7Wl/mvZ9xbgDeDSZuJPAt9OEGvpfN0EHB2JPw7c3t77rKXjLY3PU+CXwLTI8+uBvyRYV+eplg6/6Mq+yL8dC2xx9+3AOGBxJLaI4Ep0PMcBS/Zwmw224+47gFXAODNLBw4Dxlsw1Gi1md0eHfaSbFuROjeb2XYzm5PoimsS7a4CaoARwCCgO3CRmW0xs2XRoRtmdoKZFSdqi4bHNZn+SwcRnr9HEXlvmNkHZnZ5+LS5c78XkE/y70GR3fEH4HNm1jcccngh8Ep9UOeppJou7d0BkX2BmeUCDwE/CIuygdJIldKwrPF6hwI3AMkmzo013k50W/0J3qOfBw4mGCLzMrAOeGw324Lga+ibgB3ARcALZnawu69Loo9bE7Tbi2BYU/2QphHAv8xsmbu/7u6zgLxm+hjtX0v9l47lx8AG4J/1Be5+SCTe3OudDdS6e0WcmMhntQooBrYADrwK/L4+qPNUUo2u7EunZ2bdCMbqv+juj4bF5QTJdb2csCy63nDgBeAr7p7UlX0zWxK5+XVInO1Et1UZPr/b3Yvd/VPgQeDMPWgLd3/f3YvcfZe7PwnMJRiC1JJk+jjN3SvD4/BEfR+TaCt6XJvtv3Qc4bc7FwBfcndPUK2517scSDezHnFiIp/VrwkueuQSXLDYBvw8QV2dp9LhKdmXTs3MugDPENxU+51IaCnBDVz1JtBwOEI+wde+P3L3vyW7PXc/yP998+unjbcTzggxEljq7kVhvxp0eU/aStCduiS73bjdEfz75ublwK5EfWypLRoe193tv+yDzOwSgm/ITnf3uDeTh1o69zfTzHtQ5DM4BHjM3cvcvRR4FDg5QV2dp9LhKdmXzu53BGPOpzS6AjkDuM7MhoeJ/bfCsvohP/8EHnf3hxs3aGZdwm8L0oEuZtYtHL8czxtAdjhlYSbBLD4L3H1tGJ8OfNfMeprZQIKZTl7ck7bM7EIzywr7dwnB7DmvtXiE4K/AMWZ2engFaxrwrLvvDMev/hm4NZzObjRwBfBSgrYSHtckjoXsI8wsIzzH04CM8BxPM7PTCG5+PMvd17TQzBs0/3rPAG4Lz/1jgHOBP7bF/khqSnSeAguA/zSzHmHyPoWG4+6j3kDnqXR07X2HsBYt7bUAQwnGa1by769ky4FJYfwWgq93i4B7AAvLvxyuF12nPNLunWE8ukxpph9HAh+E/XiTyOwzBDME/Q4oIRj/PK2+H3vQ1qywnWJgPvD5OOs3mY0nLD+TYJxrBcGUor0isTzgLwRTcq4BrovEJkWPTXPHtaX+a9l3FoIPoY3P8ROB14HqRu+N30bWWwJckeT52p1gNp9ygvtULm/v/dbSsZZmztN+BBcptgHbgeeBgZH1dJ5qSamlPnkREREREZEUo2E8IiIiIiIpSsm+iIiIiEiKUrIvIiIiIpKilOyLiIiIiKQoJfsiIiIiIilKyb6IiIiISIpSsi8iIiIikqKU7IuIiIiIpCgl+yIiIiIiKUrJvoiIiIhIilKyLyIiIiKSov4/ybs17wil02MAAAAASUVORK5CYII=\n",
      "text/plain": [
       "<Figure size 800x250 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "for i in waveform_IRIS_files:\n",
    "    if i != '.DS_Store':\n",
    "        wf = read(\"waveforms_IRIS/\" + i)\n",
    "        wf.plot()\n",
    "        break"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 14,
   "id": "5618ca5c",
   "metadata": {},
   "outputs": [],
   "source": [
    "os.mkdir(\"waveforms_SCEDC\")"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 15,
   "id": "e356b867",
   "metadata": {},
   "outputs": [],
   "source": [
    "out_station_json_file='station_list_SCEDC.json'\n",
    "\n",
    "station_list = {}\n",
    "\n",
    "for ev in SCEDC_inventory:\n",
    "    for st in ev:\n",
    "        station_list[st.code] = {\"latitude\": st.latitude, \"longitude\": st.longitude}\n",
    "        \n",
    "with open(out_station_json_file, 'w') as fp:\n",
    "    json.dump(station_list, fp)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 16,
   "id": "ab373c53",
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "No data available for AZ.TONN\n",
      "No data available for BC.CEA0\n",
      "No data available for CE.01711\n",
      "No data available for CE.01794\n",
      "No data available for CE.11023\n",
      "No data available for CE.11217\n",
      "No data available for CE.11369\n",
      "No data available for CE.11591\n",
      "No data available for CE.11627\n",
      "No data available for CE.11628\n",
      "No data available for CE.11684\n",
      "No data available for CE.12626\n",
      "No data available for CI.ERR\n",
      "No data available for CI.FRK\n",
      "No data available for CI.IMP\n",
      "No data available for CI.IRG5\n",
      "No data available for CI.NSS2\n",
      "No data available for CI.Q0027\n",
      "No data available for CI.Q0082\n",
      "No data available for CI.RXH\n",
      "No data available for CI.SAL\n",
      "No data available for CI.SGL\n",
      "No data available for CI.SLB\n",
      "No data available for CI.SLV\n",
      "No data available for CI.SNR\n",
      "No data available for CI.SWS\n",
      "No data available for CI.WMD\n",
      "No data available for CI.WWF\n",
      "No data available for EN.ELM\n",
      "No data available for EN.ENG\n",
      "No data available for EN.HAT\n",
      "No data available for EN.LIN\n",
      "No data available for EN.OBS\n",
      "No data available for EN.RED\n",
      "No data available for EN.SIM\n",
      "No data available for EN.YOU\n",
      "No data available for NP.117\n",
      "No data available for NP.412\n",
      "No data available for NP.464\n",
      "No data available for NP.5028\n",
      "No data available for NP.5051\n",
      "No data available for NP.5053\n",
      "No data available for NP.5055\n",
      "No data available for NP.5057\n",
      "No data available for NP.5060\n",
      "No data available for NP.5061\n",
      "No data available for NP.5158\n",
      "No data available for NP.5165\n",
      "No data available for NP.5274\n",
      "No data available for NP.5413\n",
      "No data available for NP.5490\n",
      "No data available for NP.931\n",
      "No data available for NP.952\n",
      "No data available for NP.955\n",
      "No data available for PB.DHL2\n",
      "No data available for PB.SCS1\n",
      "No data available for PB.SCS2\n",
      "No data available for SB.WLA01\n",
      "No data available for SB.WLA02\n",
      "No data available for SB.WLA03\n",
      "No data available for SB.WLA04\n",
      "No data available for SB.WLA05\n",
      "No data available for SB.WLA06\n",
      "No data available for SB.WLA09\n",
      "No data available for SB.WLA10\n"
     ]
    }
   ],
   "source": [
    "st_list = []\n",
    "# wave_form = ''\n",
    "\n",
    "for i in SCEDC_contents[\"stations\"]:\n",
    "    curr = i.split()\n",
    "\n",
    "    network = curr[0].split('.')[0]\n",
    "    station = curr[0].split('.')[1]\n",
    "    \n",
    "    try:\n",
    "        wave_form = client2.get_waveforms(network, station, \"*\", \"*\", start_time, end_time)\n",
    "\n",
    "    except:\n",
    "        print(f\"No data available for {network}.{station}\")\n",
    "\n",
    "\n",
    "    for s in wave_form:\n",
    "\n",
    "        if s.stats.channel in chan_priority_list:\n",
    "            s.write(f\"waveforms_SCEDC/{s.stats.network}.{s.stats.station}.{s.stats.channel}__{start_time_file}__{end_time_file}.mseed\", format=\"MSEED\")\n",
    "\n",
    "            if s.stats.station not in st_list:\n",
    "                st_list.append(s.stats.station)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 17,
   "id": "1259dbb2",
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAuEAAADtCAYAAAARHavAAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjYuMCwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy89olMNAAAACXBIWXMAAA9hAAAPYQGoP6dpAAA5HElEQVR4nO3de3ycZZ3//9cn5zTnpG2a5tCkpUdoKS20tJTS0hMgAuqiKK6tpxXxpysiK6CCVYHF9evZFXUVRGBB1l0UEAQFFpSji0BtS4v0mLY0Sdsc27RNev3+uO8Md6YzkzSdSSbT9/PxuB7N3Nd1X/fnvuea9JNrrrnHnHOIiIiIiMjgSRvqAERERERETjRKwkVEREREBpmScBERERGRQaYkXERERERkkCkJFxEREREZZErCRUREREQGmZJwEREREZFBpiRcRERERGSQKQkXERERERlkSsJFRERERAaZknCRCMws28xuN7N6M2sxs6fMbHpYm2vNrNHM9prZN8zM/O2TzewhM2vy6+8ys5LAflea2Stm1mVm1/YjljPM7FUz229m/2tm48LqP2xmb5hZh5mtN7NJA+nLzO4ws4Nm1u6Xtf726wPbDprZ4cDj2/w255vZ3/0YfhN2vqPM7GH/mBvMbEkf5xvxuvbnWoT1k+tf+zYz22Zm7w+rX+U/v63+c50Vo68BnZ+ZXWxmr/tjaJeZfcvM0v26DDP7tZntMDNnZmNiXZeBMrPRZvYrM9vtX9MHzawmUJ9mZt8xs2a/zVWBunlm9qSZ7fPj/0HwOpnZajNbZ2ZHzOyyfsQS9Tr69dea2Xb/OXvFzIoH0pd5r9fOwDh9xN9+W2DbYTM7FHh8vZldaGbP+c9XvZndGHbMh82swa9/wczmxYjvK2Gvlfaw+n82s61+X3eYWXaMviKej183wsx+7Me128w+H6OfNDO72R9ze83sm2a9Xl/nmvf6ajfv9VUbqLvCzDb5r5dtZnZdoK7cvN95jWbWGe34IhLGOaeiohJWgDzgy0AVkA58DngzUH8BsBUYD1QA64CP+HVzgA8BRX4/vwJ+Ftj3EuBC4L+Ba/uIIxvYDnwEyAFuBZ4O1L8T+BtwMmDABKB0gH3d0Y94rgXuCNs2GtgHnO+f713AnYH6XwE/BUb4574XKInSf6zrGjP+CH19A/gdUAjMB5qBSX7ddD+O0/3n6Ungq1H6GfD5AZXAKP/nEuAPwCf9xxnAPwNnAg4Yk6CxXAd8BhgJZAHfBv4YqL8SeNk/z8nALuBcv+48/5zygDLgaeDLgX0/CCwDngcu6yOOvq7jp4H/Bcb5Y/kUIGeAfT3Vj3huA74Stu39wBJ/rNUA64F/DNSfDGQEXntvARal/68At0WpWwbUAyf54+a/gVtixBr1fICbgMf8cV4NvAGsiNL248CrQDlQCjwLfMKvG+mP3QvwfuddAzwX2LcWKPZ/rsD7vXO+/3gUcIV/TToTMY5VVFKxDHkAKirDofjJyxGgzH/8n8AXAvUfAZ6Isu8yYE2E7XfQd9K7AlgfeJwH7AfG+Y9fBJb18xz66qs/8URKwj8BPBJ4PAE44Ccy+cBBoCJQ/zTwoSj9R72ufcUfoa9dwNzA4zuBG/yfbwF+FKg7F9gUpZ+4nB9vJ+HfiFCXsCQ8wrEmAm2Bx88B7ws8/irw8yj7fhx4MML2p+g76Y11HdP952tiP88hal/HEM9RSXiENjcB34+w3YB3+M9bXpR9v0L0JPz/AbcGHs8DdsSII+r5AP+Hnwz7j68D7o7S9tf4fwT6j98P/Nn/+Z3AC4G6NP+aHvWc4CXha4Arw7bXoiRcRaXfRctRRPpnHrDbObfHfzwN7z+hHq/izZJFMh9YO8Dj9jqOc64DeBOY5i9rOA2Y7r91vtnMbgi+vdzfvgJtrjGzPWb2rJktHGCMbwJdeLPZE4EW59yuQPvQtTKzBWbWHK0vel/XmPH7Sxke8n8uAcb0ty+/rs7McuN5foFzbMGbZZwO/DzCMQZT+HgcqrEcvI5VQC5wqb+kYoOZXTHAvnp8318e8biZzRhgzEedrz/GOoGHgO/54xAz+4CZvRa2/2X+so+/mtm7w+os7OexPctvzOzfzezfw9rHOp/wvqI9f321jfS7IziWP2BmbcBOvHej/ivGcUSkD0rCRfpgZkXAj4EvBjbnA62Bx63+tvB9Z+ItBfjaAA8ffpzgscrxljQswXvr/hzgfcCqAfQF8F28t8crgB8CD5pZ9XHGGPOYzrk/OeeKY/QVjK+vvv7VOXdhoG23c25/P/tqDWwPN+Dz8+P6k3OuCC9BvA1vGcWQ8J/PfwW+FNjc37F8Pt7yj28P8PCxrlUl3rKgCXizqe8GbjSzxQPoC+Bf8Jbh1ACPA78zs0jPbVRm9gm818Ivgtv9MVYA/APeMp6e7fc454LJ8a+AKXhLZ74A3G5mp/t1jwEfMLNJflxX+9tH+H1d6Zy7MtBXrPN5DPismRX7a7hX9fQTwWPAJ8yswsxGAp8KtH0OmGRm7zCzTLzlKNnBvvxzLMD7Y/IuvHeiRGSAlISLxGBmOcADwMPOueAMZjveGswehf624L51wIPAR51z/Zo9NLO1gQ9f1UQ4TvBYB/zHtzrnmp1z2/CS5wsG0BfOub865/Y55w455+7G+095WT/CjtVvzGP2o69g22Ppqx1IN7MRUdpGOg5E7+u4z885txlv9vY7EY5xTKJ92K+PfUqBR4GbnXN/CFT1ZyyfAdwOXOKc2z3AGPszllc75w74r5df8vZYPpa+cM696Jxr9/v6hr99Tn/i9o93IXADcKFz7kB4vf8a+TVwtZlNjdSHc26dc+4t51yXc+4xvKVWF/l1vwe+BTwCbMRbm90FRLy2fZzP14EteOvXHwHuB3ZEObX/8Nu8CLyE95mJHf4xmoBL/f7ewnt3Ym2kvpxzfwM66D0xISLHSEm4SBRmlgHci/fWa/gdB9bhzQb1OJXA29bm3eniceBrzrkH+ntM59zJzrl8v2wLP46Z5eHNFq5zzu3zY+sV9kD6ihLOkX6GHd7veLwZ+k14HxIrst53/uh1rWL1Fda23/H71+at/vbl122OlHDF+fzS/JiPS+B57dfsrj9r+jvgN86574dV9zWWpwK/xfuA7PPHEWOs67gROBQe9gD7iqS/Yxl/GdbPgIucc3/vo3kG3gx1f/SKwTn3LefcBOfcWLwlQK8457qPtS/nXIdz7p+ccxXOual469T/Emkn51y3c+5a51y1c64OL+n/S6D+cefcac65Mrx3S3o+gBlJXMayyAltqBajq6gke8Gb+fs9kBmh7h14s091eGuP/8bbd/EowvtP9YYo/Wbgraf8Jd5/dDlAepS2PXcEWeX/fAu972hyE97a1AJgLN5M66oB9vUevA87ZuAta2kFasP6iHV3lBV4b13/kt53qrgfbzlPLt5MYKy7o8S6rjHjj9DXvwWuzZl4d0eZ7NdNB/YAs/zn64/0fXeUYz4/vJnFGv/nk/CWLwQ/kJftP/8O764gEe8GcpzjOAtvGcLPotR/Ci8RGwVMwvvDbolfVw1sA1ZG2TfTj/9p4B/9n9MGeB3vBn7iX5PJfhyLj7UvoBjvHZxs/9yvwvuDrDCsj0h3R5kJNABLIxxzHN5djXL8vv8/vNdIeZQYL/LHVhreB39bgDP9ulxgKt4fGpPwfl9cEqWfmOeDN2M9Bu/DrcvwZq7HRumrDG+5j+HdGWgzMDvs/DPw7pxyD/CdQN0/+tfd/HZb6f0hzxz/eevsuUbxHssqKqlWhjwAFZVkLP5/uA7vbfL2QDk70OY6oMlPBr6Bf6syYKW/b3C/9sB+X/Hrg2VVjFjOAF7zY3mawN1A/P+Uf+r/B78DWN0TxwD6+pPfTzPwAn4iFrb/UUm4v/0CvA9J7sebNS0J1I3Cm4XdjzfjuTRQd3bw2sS6rv2I/3p63zEjFy+xa8dL3j8QdpxV/jVrw7szTHagbi1weRzO73q8W9F14CWz/49Aoo33B0evsZCAsXyO33dH2Jjs+eMgDW+JTDNeAvq5wL434s26BvdbG6i/Izx+YFGMWGJdx2K8u3e0+dflk32cV8S+/OfjL36se4EngFkR9o+UhN8OdIed7yOB3wnP4SXe+4BngHMC+14edm3u9Y/fhjdmg3egKfXHWId/rldEiO22/pwPsBhvfO/H+yPv7LC+Qr+38BL/v/ttXwfeHdb2v/zz2wN8j96vie/jzZy3+318kd6vzfBxsCXeY1lFJdVKT9IgIiIiIiKDRGvCRUREREQGmZJwEREREZFBpiRcRERERGSQKQkXERERERlkSsJFRERERAaZknARERERkUGmJFxEREREZJApCRcRERERGWRKwkVEREREBpmScBERERGRQaYkXERERERkkCkJFxEREREZZErCRUREREQGmZJwEREREZFBpiRcRERERGSQKQkXERERERlkSsJFRERERAaZknARERERkUGmJFxEREREZJApCRcRERERGWQpk4Sb2UIze9DMdpqZM7NLBtCHmdnnzWyjmR00s+1mdn0CwhURERGRE1jGUAcQR3nAq8DtwK8H2Md3geXA54E1QBEwMi7RiYiIiIj4zDk31DHEnZk54F3OuQcC27KArwOXA8XA34AvOOee8uunAq8BpzjnNgxyyCIiIiJyAkmZ5Sj9cDtwFnAZMAO4H3jUzCb69e8ENgEXmtlmM9tiZv9hZqVDE66IiIiIpKoTIgk3swnA+4FLnXPPOOfedM59E/gT8GG/2XhgHHAp8CFgFTAb+K/Bj1hEREREUlkqrQmPZRZgwEYzC27PBvb4P6f5jz/knNsIYGYfBf7PzCZriYqIiIiIxMuJkoSnAd14M9vdYXXt/r+7gK6eBNy33v+3BlASLiIiIiJxcaIk4X8F0oHRzrlnorT5M5BhZhOcc2/62yb5/25NdIAiIiIicuJImbujmFk+cJL/8K/A54Angb3OuW1mdhfeBzOv9utHAucCa5xzvzOzNOAlvJnxz+LNnv8QaHXOLR/McxERERGR1JZKSfgivKQ73C+cc6vMLBP4Et6HLivx1oI/B9zonFvj9zEW+D7evcI7gEeAq51zexN+AiIiIiJywkiZJFxEREREZLg4IW5RKCIiIiKSTIb1BzPNu9/gWKBtqGMREREREfEVADtdjCUnwzoJx0vA64c6CBERERGRMFXAjmiVwz0JbwPYvn07hYWFQx2LiIiIiJzgWltbqa6uhj5Wagz3JByAwsJCJeEiIiIiMmykRBLeX4cPH+bll19m9+7d6K4wsZkZ5eXlzJo1i8zMzKEOR0RERCSlJCQJN7P2sE15wD84534doe0dwPuBw/6mrc65k+Md0wsvvMB1111Ha2trvLtOaQUFBdx0003Mnz9/qEMRERERSRkJv0+4mc0AngfKnXNHrY3xk/DXnXP/OoC+C4GWlpaWmMtRmpqauOiii5g9ezaf+cxnmDBhAmlpujtjLEeOHGHTpk384Ac/4MUXX+SBBx5g9OjRQx2WiIiISFJrbW2lqKgIoMg5F3X2dzAy0cuB30RKwI+VmWWbWWFPwbv9S5+eeOIJjhw5wk033cTEiROVgPdDWloaJ510EjfddBPOOf7whz8MdUgiInKCcs5xwQUXsHLlyqEORSRuEpqN+vfxfj9wdx9NrzGzPWb2rJktjNHuOqAlUPp1e8ItW7ZQW1urD28OQF5eHhMmTGDLli1DHYqIiJygXn/9dR555BHuvPNOfaZLUkaip4QXAiOA38do813gJKAC+CHwoJlVR2l7C1AUKFX9CaKrq0sfLjwOmZmZdHV1DXUYIiIiSsIlZSQ6Cb8c+JVz7nC0Bs65vzrn9jnnDjnn7gaeA5ZFaXvQOdfaU4jDN2X+9Kc/Zfr06eTl5VFTU8PKlSvZsmULixYt4t577424z6pVq8jOzqagoIDCwkJOP/10vvWtb/VKVBctWkROTg75+fmMHj2alStX0t7eflQf+fn5obJq1ap+xVxbW8vzzz/fa9sVV1zBV77yFQCeeuoppkyZctR+OTk5oRntVatWkZ6ezvr160P19957L4sWLepXDCIiIoMluIy0u7t7CCMRiZ+EJeFmlgX8A30vRQl3JAHhRPT1r3+dG264gVtvvZU9e/awbt06zjrrLJ544ok+9129ejVtbW3s3LmTm2++mR//+Md88IMf7NXmjjvuoL29nTVr1vDKK69wyy23HNVHe3t7qNxxxx3xPL0+FRUV8bWvfW1QjykiInKs0tPTQz8fOTJoaYJIQiVyJvwCvHXbz8ZqZGbvMbM8M8sws/cBC4C+s+Dj1NzczM0338yPfvQjLrjggtCs9T/90z/xkY98pN/95Ofns3z5cu677z5+9atf8eqrrx7Vpry8nBUrVvDKK6/E8QyO38c+9jEeeeQRXn/99aEORUREJCrNhEsqSmQSfjlwjwtbvGVml5vZ2sCmq4CdQBPwOeBdzrktCYwLgOeee45Dhw5x4YUXxqW/mTNnUl1dzZ///Oej6nbt2sWjjz7K+PHj43KseCktLeXKK6/UbLiIiCQ1zYRLKkpYEu6cu9Q598UI2+8OfhmPc26Bc67IOVfsnJvrnPtjomIK2rNnDyNHjiQjI37fV1RRUcG+fftCjz/60Y9SWFjI2LFjKSwsZPXq1b3ar169muLi4lD55je/2e9jLVu2rNe+t99+e6/6N954o1d9cXExBw8ePKqfz33uczz88MOaDRcRkaQVnAlXEi6p4oS9YXZZWRlNTU1xvevHrl27KCkpCT3+2c9+RmtrK3/+85/ZsmULTU1NvdrfeOONNDc3h8rnP//5fh/r8ccf77Xvhz/84V71EydO7FXf3NxMdnb2Uf2UlZVx5ZVX8vWvf/0Yz1ZERGRwaDmKpKITNgmfN28emZmZPPzww3Hp79VXX2X79u2cddZZR9XNnz+fT3ziE1xzzTVxOVa8XX311Tz00ENs2LBhqEMRERE5ipajSCqK31qMYaa4uJgvfvGLXHnllWRnZ7N48WK6u7u55557jukbNTs6Onjuuef49Kc/zXvf+15OPfXUiO0+9alPUVNTw5o1a5g+fXq8TiMuysrK+OQnP8n3vve9pItNRETE++4/j2bCJVWcsDPhAF/60pe48cYbueaaaygpKWHy5Mk888wzLFmypFe7Z555hvz8/F7bbrzxRgoKCigvL+cLX/gCH//4x7nrrruiHqu0tJSVK1dy66239uojeJ/w5cuXRzzezTffzPnnnx+PU47q6quv5tChQwk9hoiIyEAEk3DNhEuqsOH8zVNmVgi0tLS0xPxK+ptvvpn169fzy1/+cvCCSyEf/vCHqaur44YbbhjqUERE5ATU0NBAeXk5APX19VRWVg5xRCLRtba2UlRUBFDkf7lkRCf0TLiIiIgML5oJl1RxwiThw3nGf6jp2omISLLQmnBJFSdEEp6fn8/evXuVTA6Ac469e/cetSZeRERkKGgmXFLFCZGEn3HGGTQ0NPDSSy8NdSjDzssvv8zOnTuZM2fOUIciIiInqOAkmmbCJVWcELcoPOOMM5g9ezZXXXUVS5cuZcKECb3uOSpH6+7uZvPmzTz++OPMnDmTuXPnDnVIIiJygrrwwgtDP2smXFLFCXF3FIADBw5w55138sQTT9DQ0KAXcR/S0tIYNWoUixcvZuXKlYwYMWKoQxIRkRNU8BaFa9euZdq0aUMYjUhs/b07ygmThIuIiMjwFEzC16xZwymnnDKE0YjEplsUioiISMrRO9mSKhKWhJvZU2bWaWbtfnkkSrtcM7vLzNrMbJuZvT9RMYmIiMjwpg9mSqpI9AczVznn7u2jzWqgFKgETgF+Z2b/55zbmODYREREZJjRTLikimRYjvKPwGrnXKtz7lngt8BlkRqaWbaZFfYUoGAwAxUREZGhpZlwSRWJTsK/b2aNZva4mc0IrzSzEmAMsCaw+VXg5Cj9XQe0BEp9nOMVERGRJKaZcEkViUzC/wWoA2qAx/GWmYR/7WI+0O2c2x/Y1upvj+QWoChQquIasYiIiCQ1zYRLqkhYEu6ce9E51+6cO+Cc+wbQDoR/7WI7kG5mwZtQF/rbI/V50F+20urf8qUtIcGLiIhIUtJMuKSKwVwTftSrxjm3D3gLmB7YfCqwdrCCEhERkeFDM+GSKhKShJtZsZkt8z9ImWVmV+HdAeUvEZrfBXzZzArM7EzgIuC+RMQlIiIiw5tmwiVVJGomPBNv/fYevJnudwIXOOdazexyMwvOdN+A9yHLXcD9wJXOuQ0JiktERESGMc2ES6pIyH3CnXONwOlR6u4G7g48PgBcnog4REREJLVoJlxSRTLcJ1xERESkX5SES6pQEi4iIiLDhpajSKpQEi4iIiLDhmbCJVUoCRcREZFhQzPhkiqUhIuIiMiwoZlwSRVKwkVERGTY0Ey4pAol4SIiIjJsaCZcUoWScBERERk2NBMuqUJJuIiIiAwbmgmXVKEkXERERIYNzYRLqlASLiIiIsOGZsIlVZhzbqhjGDAzKwRaWlpaKCwsHOpwREREJAHMLPRzRkYGhw8fHsJoRGJrbW2lqKgIoMg51xqtXUJmws0s28xuN7N6M2sxs6fMbHqUtneY2UEza/fL2kTEJCIiIsNfV1fXUIcgEheJWo6SAWwCzgRKgd8CD8Rof6NzLt8vJycoJhERERGRpJCQJNw51+Gc+5pzrt451w38AKgzs7JEHE9EREREZDgZrA9mzgN2O+f2RKm/xsz2mNmzZrYwWif+MpfCngIUJCRaEREREZEESngSbmZFwI+BL0Zp8l3gJKAC+CHwoJlVR2l7HdASKPXxjVZEREREJPESmoSbWQ7eWvCHnXM/j9TGOfdX59w+59wh59zdwHPAsihd3gIUBUpV/KMWEREREUmsjER1bGYZwL3ATuDzx7Br1BuAOucOAgcDxxhwfCIiIjI8dXd3k56ePtRhiByXRM6E/xTIBVa5GDcjN7P3mFmemWWY2fuABcATCYxLREREhrFDhw4NdQgixy0hM+FmNg5YBXQC+wIz1ucDNcD1gVsRXgX8HHDABuBdzrktiYhLREREhr+DBw+Sm5s71GGIHJeEJOHOua1ArLUidwfaLkhEDCIiIpKaNBMuqWCwblEoIiIiEhcHDx7su5FIklMSLiIiIsOKZsIlFSgJFxERkWFFM+GSCpSEi4iISFJ5/vnn2bZtW9R6zYRLKkjYfcJFREREjtW6deuYN28eANHucKwkXFKBZsJFREQkabz22mt9ttFyFEkFSsJFREQkaRQVFYV+PnIk8pdoayZcUoGScBEREUkahYWFoZ87OjoitulrJvzIkSOsXLmS73znO/EMTSSulISLiIhI0gh+E2Zra2vENn3NhD/99NPceeedXHXVVXGNTSSelISLiIhI0jB7+wu3W1paIrZ517vexVtvvRW1j/T09NDPWj8uyUpJuIiIiCSlaDPhAOeff37UuuLi4tDPTU1N8QxJJG6UhIuIiEhSipWEv/LKK1Hr0tLeTm8aGxvjGZJI3CgJFxERkaTx6KOPhn6OlYQD1NfX99lfQ0PDccckyampqYmurq6hDmPAEpKEm9koM3vYzPab2QYzWxKlXa6Z3WVmbWa2zczen4h4REREJLkdOXKE7du3c/3114e2XXrppZSXl0fd593vfnef/fYnUZfhZ+3atVRUVHDJJZcMdSgDlqiZ8B8CO4GRwBeA+82sJEK71UApUAlcBvzIzCYlKCYREZEh45yju7t7qMNIWp/+9Kepqak5anusmeyXXnqJVatWcf/999Pe3s7LL7/MDTfcQHt7e6jNRz/60bjE95vf/Ib3vOc9xzWzXl9fz7333jusZ2+Pxa5duygrK8PMWLduXVz7/uMf/0hXVxcPP/xwr+d7OLFoXwk74A7N8oE9QK1zbpe/7WngP5xzd4a13QVc4px7wX98J/B359xX+3msQqDlsssuY9q0aTz77LO0tbVRWVnJGWecwfr168nLyyM7O5uRI0eye/dutm7dSmlpKRUVFRw5coR9+/bR1tbGjBkz+Nvf/kZZWRnd3d1UVVXx5ptvsmfPHkaNGkVeXh47duxg1qxZrF27lvT0dCoqKnjjjTdobm4mPz+f/fv3U11dTVlZGZs3b2b69OnU19dTVVXFG2+8wb59+2hubiYtLY2xY8dy8skn88orr7B161bmzJnDoUOHqK+vp729naysLHJzcyksLOTUU09l27ZtjBw5km3btlFdXU1zczPjxo1jy5YtFBUVkZGRwb59+6iqqmLbtm2MGTOGNWvWUFNTw8GDB2lra2PSpEls2LCBzMxMMjIy2LNnD7Nnz2bt2rW0tLRQV1fHkSNHGDduHJs2baKiooLOzk7a2tr4y1/+wnnnnRc699raWvLy8ti2bRt79+5l5syZ7N27l/T0dLZv305jYyNtbW2Ul5czc+ZMzIzW1lb27NlDd3c3kyZNoqOjg66uLkaPHs327dvZvn07nZ2dPc8tCxYsoLGxkY6ODioqKjAz2tvb2bdvH1OmTGHXrl0UFhaG+m5sbGTKlCns2bOHvLw8cnNzaWhooKamhu3bt1NaWkp3dzfNzc1UVFSwc+dOxo4dS319PRkZGTQ2NnLyySezbds2Kisr2bx5M+3t7UybNo2mpiaqq6tpamoiKyuLgoICGhoaQte7uLiY7u5u9u/fT11dHQ899BANDQ2MGTOGcePGkZ6ezvjx40lLS2P37t0ArFmzhlmzZtHS0kJlZSUbN25k9+7dTJ8+ndbWVsaNG0d9fX3oA0atra3U1NTwt7/9DeccY8eOJSsri+bmZmpqati6dSujRo0KPWdFRUXU19fT2trKyJEjKSgooKSkhIaGBqqrq9m+fTv5+fmkp6fT0tLCnj17mDJlCg0NDeTn55Obmxs67+3bt1NSUsLhw4dDz0d9fT1jxoxh7969pKWlUVJSQmNjIxs3bqSoqIja2trQuOx5DkaOHMmhQ4c4cOAAY8aMYfv27YwZM4Z9+/aRnp5OUVERu3btorKykh07doTG9t69e6mqqmLr1q1MmDCBPXv2kJmZSWNjI+PGjevVd0dHB2PHjmXr1q1UVFSwe/duzIycnBxaW1sZP348O3fuJC8vj8zMTJqbm9m1axeTJ09m79695OXlkZ+fT2NjI1VVVaGx03PuPWOmvLycLVu2UFpaSmdnJ9XV1WzdupW6ujr27t1LZmYme/bsobq6mm3btlFaWhqKr7y8nB07dlBZWcm+fftIS0ujuLiYhoYGKisr2bRpE4WFheTm5lJaWsqWLVsoLy+nsbGRzMxMtm7dSkNDA/Pnz6ehoYFRo0ZhZmzevDn0Gr344otpa2vDOUdJSQm7d+9m79691NXVsXnzZioqKmhqamLLli0sXbqUN954g6ysLIDQ2H7yySepqKjglFNOYffu3Tz66KPs3LkTgClTpnDRRRfR1dVFXV0dGzdupKmpidGjR5OdnU1jYyNNTU1kZ2dTXl5OUVERW7ZsoaWlhe7ubtLS0igqKqK8vJw33ngD5xxHjhzpVcC7y8WIESMAb6Y0KyuL9PT0XgltWloaR44c4cCBAzQ3N5OXl8eBAwdIT08nOzub9PT0UJ8FBQVkZWXR1NTE5s2baW5upra2lsrKSg4dOsShQ4fIyckhLS2NjIwMuru76ejoYMSIERQXF9PW1kZraysdHR3k5+eTlpbG4cOHSU9PJyMjAzMjPT2drKws9u/fT05ODllZWRw6dIgHH3ww9P/XyJEjQx8YPOmkk5g4cSIZGRns37+frq4ucnNzycnJYdeuXYwYMSL0HDc3N4fG0MGDB+ns7OTw4cN0dnaSlpZGTk4OTU1NbNu2jY6ODurq6qisrCQtLY19+/aFxlVOTg4tLS0cOHCA0tJS2tvbycnJISMjg46OjlBfBQUF5OTkcMopp/DKK6+EYh4xYgTp6ekcOHCAzMxMurq6cM4xYsQI2tra6OrqIj09ne7ubvLy8jAzGhoaaGlpYePGjYwaNYrJkyezdetWNm/e3L8E4zjl5OSQmZlJW1sbY8eOZc6cOWzcuJHm5mYKCgrYsGFDqG1JSQkzZ84kNzeX/fv3s23bNjZt2hS171mzZpGVlcXBgwcZN24cM2bMYP369XR0dJCTk0NnZycHDhzgySefDO1jZjjnOPnkk5k/fz65ubk0NjbinKOlpQXnHHv37mX//v0UFhaG/r/LyMggPT2dadOmhV5np512Gq+//jqdnZ1kZGTQ2dlJc3MzJSUldHV1sWPHDkpLSykrK6O1tRUzIy0tjRkzZnDo0CG2b98OeGNx06ZNTJs2LZRT7Nixg4yMDF5//XXMjOzsbNra2mhra+Oll14KvY4nT57MnDlzeOKJJzh06BArVqxg9erVHDhwIHTOc+fOpbi4mJycnNBruLy8nLS0NH7961+TmZnJ0qVLmTp1Ki0tLezfv59NmzaFznvz5s1kZ2czYsQInnnmmV7XcsmSJWRlZbFkyRKefvppsrKyqK2tZevWrcycOZP169eH8petW7fS2NjI4cOHKSgo4NRTT2Xt2rXU1tbS2dnJjh07yM3NZcqUKbz66qtMnTqVF154gT/+8Y/MmDGDd73rXTzzzDO0tLTwjne8g927d1NRUcGmTZvIz8/n6aef5uWXXwYocs5FXVOViCT8NOD3zrnRgW3fB/Y7574Q2FYC7AXynHP7/W1XA3Occ++L0nc2kB3YVADofSYRERERSTYxk/BELEfJB8IP2OpvD2/X3ZOAx2gXdB3QEiihBLygoKBXw9zcXGbMmAHAtGnTGDFiBOeee26ovqqqirS0NM444wzA+6u3trYWgNNPP52cnBzOPvvsUPtJk7xVMjk5OaF9xo8fT0ZGxlFB1tXVAZCRkcGSJUsoLS1lwoQJR7ULj3nevHkRT7qoqIhFixZRXV1NRUUFAOeeey5FRUXMmzePcePGMX78eObOnUtBQQFLly5l5MiRlJWVAbBo0SIASktLOemkkwCYPn064N3GaebMmaFjLViwgLKyMpYuXUp5eTm1tbWhrxDOyspi4sSJgPfXckVFBWeeeSYAY8aMYe7cuUyePPmo+Gtraxk/fjxz5swJbauoqGDGjBnMmDGDoqIiFi5ceNR+I0aMYMkS7+MElZWVTJo0KdRHUVERCxYsYNy4cdTU1ISuXUFBAfPnz2f8+PFUVlaycOFCSkpKWLRoEVVVVUyePJnTTjuNESNGsGzZMoqKipgwYUIo7uzsbJYsWRI69x7nnXceo0aNYt68eUyePJnq6moWL15MWVkZixYtoqamhokTJzJ37tzQtQo6++yzGTVqFBdddBGzZ89m9OjRoXYrVqxg9OjRvcbIkiVLKC4uZuHChdTU1DBhwgTmzp1Lfn5+6JqPHz+e2tpazj77bEpLS1m8eHFo5r3neZ4/f36oz8mTJzNmzBiWLl1KaWkp8+bNo6amhrq6ulC7nus3YcIExo0bxznnnENpaSmLFi2isrKSqVOncvrpp1NaWsqKFSsoLy9n1qxZTJs2jVGjRnHOOeeEjheMb+TIkSxdupSqqiqmT5/OaaedRllZGeeddx4VFRXMmTOHKVOmUF5ezrnnnktZWRkLFy6krq6Ok046iXnz5lFSUsKyZcuorq7mwgsvpKamhkWLFlFWVsbixYsZO3YsU6dOZfbs2aG2ZWVlTJo0idraWnJycjj//PMpLy9nwYIF1NXVUVtbGxo7+fn5nHXWWUycOJGxY8eyYMECioqKWLx4MdXV1UycOJHZs2dTXFzM8uXLGTNmTGjczJs3j9LSUpYtW0ZNTQ3vfOc7qamp6fXcVFVVMWXKFE477TTy8/NZvnw5o0eP5vTTT2fKlCmMHj2apUuXUlJSwpw5c6itrWXKlCmcc845nH/++YwePZp58+YxceLE0PgqKChgyZIl1NTUUFNTw5gxY0LXv7KyklmzZjF16tTQNQJvBnb8+PGh12FPP1VVVYwcOTL0mg2aMGFCr/suBy1cuJDS0lIWLFgQ2nbaaacd1a7nmOGqqqoibpfkkpeXp+fqGARfi/11yimnHPM+p556KuD9n5ifHyt9ii6Y6wTjXrFiBaWlpUyfPj30u6Iv0X5P9CUvL6/X49GjR5OTk8OUKVOOua/s7GzGjh0LvJ27paWlcfrppwNejhaeu40aNSr0rnPwepSURFpJDZmZmb3aLF++nLFjx4Z+d86ePbtfsR6dQR6/dqAwbFuhvz28XbqZjQgk4pHaBd0CfCvwuACob2lp6fU1tyIiIsmmvr6eP/zhD6E/9uRtzjlmz57NX//61wHtP2PGDD72sY/x1ltv0d7ezve+971QXc+Skr40Njby1ltvMWXKlF5JVo/rr7+eW265hccee4xly5YB8NZbb1FeXt7rC4Zi2blzJ8899xwXX3xxxEm84cg5F/X8X3rppdDk2X333cd73/veuB03eMx4r+o4Xq2traEJzFgSuSZ8nHPuLX9bQteEKwkXEREZ3rq6unjsscd4xzvecUz7tba2HvXO8pQpU0JrvDdv3nzUuzsDFSvhlMgOHjxIdnZ23w2P0amnnsprr70GJHUSPrjLUZxz7cBvgdX+LQgvAk4BHozQ/C7gy2ZWYGZnAhcB98U7JhEREUluGRkZrFixote2H/zgB6EllZEcOHDgqAQcvOVYPeKVgANKwAcgEQk4wK233gqgWxRGcCVQjTcj/k3gvc65fWZ2uZmtDbS7AW9t9y7gfuBK59yGo3oTERGRlJeent4rqfrUpz7V6yvog8aMGUNOTk7EumASLqnpvPPOY926ddxzzz1DHcqAJWRBknOuEbggwva7gbsDjw8AlyciBhERERl+em4q0CPS+myARx55JGof+hDpiWHq1KlDHcJx0dfWi4iISNIoLS3t9ThaEh68s1c4zYTLcKAkXERERJJGf5LwBx54IGYfPUtaBnKLO5HBkhr3xxEREZGUEL7EINKt/C6++OKYfVRWVoa+/VckWSkJFxERkaRx0UUX8eUvfzm03CTacpS+jBo1Ko5RicSfknARERFJGmbGV7/69teFDDQJF0l2WhMuIiIiSUtJuKQqJeEiIiKStJSES6pSEi4iIiJJS0m4pCol4SIiIpK0lIRLqlISLiIiIklLSbikKiXhIiIikrSUhEuqUhIuIiIiSUtJuKQqJeEiIiKStJSES6qKaxJuZpPN7CEzazKzRjO7y8xKYrTfYmb7zazdL7fFMx4REREZ3pSES6qK90x4EfArYAJQC2QB3+xjn3Odc/l+uSLO8YiIiMgwFp6Ef/vb3x6iSETiK65fW++cexF4seexmf0U+FY8jyEiIiInjvAk/LOf/ezQBCISZ4leEz4fWNtHmwfMbLeZ/Y+ZjYvV0MyyzaywpwAFcYtUREREko6Wo0iqSlgSbmYzgc8AX4vR7AN4y1YmAtvwEnKL0f46oCVQ6uMRq4iIiCQnJeGSqo4pCTezx8ysM0r5UqBdHfAg8FHnXNSZcOfcs865TudcK/A5vGS8LkYIt+CtO+8pVccSv4iIiAwvSsIlVR3TmnDn3PK+2pjZGOBx4GvOuQeOpXu/xDr+QeBg4FjH0L2IiIgMNxkZcf34mkjSiPctCouA3wN3Oud+0kfbGjObZ2aZZpYH/BuwFdgSz5hERERk+NJMuKSqeK8JvwSYAfxL4N7f7T2VZnZb4F7gBcBPgGa8xPsk4GLn3JE4xyQiIiLDlJJwSVXxvkXhL4BfxKi/IvDzWmB6PI8vIiIiqUVJuKQqfW29iIiIJC0l4ZKqlISLiIhI0lISLqlKSbiIiIgkLSXhkqqUhIuIiEjSUhIuqUpJuIiIiCQtJeGSqpSEi4iISNJSEi6pSkm4iIiIJC0l4ZKqlISLiIhI0lISLqlKSbiIiIgkLSXhkqqUhIuIiEjSUhIuqUpJuIiIiCQtJeGSqpSEi4iISNJKS1OqIqlJI1tEREREZJDFPQk3M2dmHWbW7pfrY7SdYGZ/NrP9ZvaymZ0a73hERERERJJNRoL6neCce6sf7f4TeBhYAnwM+B8zm+Sc60pQXCIiIiIiQ27IlqOY2WRgMnCLc67TOfcDIB2YH2OfbDMr7ClAwSCFKyIiIkPkvvvuG+oQROIuUUn4y2a2w8zuMLOyKG2mARucc4cC214DTo7R73VAS6DUxyVaERERSVpmNtQhiMRdIpLwhcA4YCYwAvh5lHb5QGvYtlZ/ezS3AEWBUnU8gYqIiEjy0x1SJBUd06g2s8fMrDNK+RKAc+4Z59xh51wj8BngAjPLitBdO1AYtq3Q3x6Rc+6gc661pwBtxxK/iIiIDD9KwiUVHdMHM51zy4+x/yP+v5HeR1oHTDazTOfcYX/bDODfjvEYIiIiksKUhEsqiuuoNrOTzexUM0s3sxLgO8DjzrmD4W2dcxuADcC1/gcurwS6gWfjGZOIiIgMb0rCJRXFe1SXA/fjre1ej5dUr+qpNLPbzOy2QPsPAOcBzcDHgXfr9oQiIiISpA9mSiqK633CnXNPAJNi1F8R9vjvwFnxjEFERERSi2bCJRVpVIuIiEhSUxIuqUijWkRERJKaknBJRRrVIiIiktSUhEsq0qgWERGRpKYkXFKRRrWIiIgkNd0dRVKRknARERFJapoJl1SkUS0iIiJJTUm4pCKNahEREUlqSsIlFWlUi4iISFJTEi6pSKNaREREkpqScElFGtUiIiKS1HR3FElFSsJFREQkqWkmXFKRRrWIiIgkNSXhkoriOqrN7Gwzaw+U/WZ2xMxGRWm/xW/T0/62eMYjIiIiw5+ScElFGfHszDn3DJDf89jMPgNc6pxrjLHbuc655+MZh4iIiKQOJeGSiuKahEdwOfDzeHVmZtlAdmBTQbz6FhERkeRUXl4+1CGIxF3CknAzOwmYCdzfR9MHzPvY87PAZ51zW2O0vQ64MT4RioiIyHBQXl7Oww8/TF5e3lCHIhI35pxLTMdmNwKznHMXx2gzH3gZyAK+Biz094kYVJSZ8PqWlhYKCwvjFruIiIiIyEC0trZSVFQEUOSca43W7pgWWZnZY2bWGaV8Kaz5B4C7Y/XnnHvWOdfpB/g5YCJQF6P9Qedca08B2o4lfhERERGRZHBMy1Gcc8v7087MzgAqgAePpXu/iIiIiIiktER93Phy4L+dcweiNTCzGjObZ2aZZpYH/BuwFdiSoJhERERERJJC3D+YaWbpwPuAD0Wouw3AOXcF3nrunwDjgf14H8y82Dl35FiP2doadbmNiIiIiMig6W9emrAPZg4GM6sE6oc6DhERERGRMFXOuR3RKod7Em7AWPQBTUluBXh/LFahsSrJS+NUhgONUxkuCoCd0e74B4n/sp6E8k8s6l8YIsnA+1sRgLZYtyoSGUoapzIcaJzKMNLn+NT3wIqIiIiIDDIl4SIiIiIig0xJuEjiHQRW+/+KJCuNUxkONE4lZQzrD2aKiIiIiAxHmgkXERERERlkSsJFRERERAaZknARERERkUGmJFxEREREZJApCRcRERERGWRKwiXpmFm2md1uZvVm1mJmT5nZ9ED9tWbWaGZ7zewb5n+FmplNNrOHzKzJr7/LzEoC+11pZq+YWZeZXduPOM4ws1fNbL+Z/a+ZjQur/7CZvWFmHWa23swmDaQvM7vDzA6aWbtf1vrbrw9sO2hmhwOPb/PbnG9mf/dj+E3wfAPXaruZtfnnXhwjxojXtT/XQpKDma02s3VmdsTMLgtsX+U//21mtsnMruijn1jjNdd/bbWZ2TYze38iz0lST4xxmm1mP/N/DzWZ2S/NLC9GPxqnMqwpCZdklAFsAs4ESoHfAg8AmNkFwCeBucDJwIXAh/39ioBfAROAWiAL+Gag353Al/z+YjKzbOC/ge/6MTwP/DJQ/07gauASIN+Po2kgffludM7l++VkAOfczT3bgBuBuwNtrjCz0cA9wKeB0UCbf4ye434aOB9YABQCHwQ6o8QY9br2M35JDm8A/wy8GLY9G7gCKAHeCaw2s4WROujH873a314JXAb8KNYfoCIRRBunnwamA5OBOqAciDhhonEqKcE5p6KS1AUvmT4ClAH/CXwhUPcR4Iko+y0D1kTYfgdwbR/HXAGsDzzOA/YD4/zHLwLL+hl/X331J55rgTvCtn0CeCTweAJwAC/hSgd2ARP7GWPU69pX/CrJV4CngMti1N8NXB2lrq/xuguYG6i/E7hhqM9ZZfiV8HEKfB9YHXh8JfDrKPtqnKoM+6KZcBkO5gG7nXN7gGnAmkDdq3gzt5HMB9YO8Ji9juOc6wDeBKaZWTpwGjDdvCUzm83shuDyjf72FWhzjZntMbNno81Q9qPfN4EuYDxQBeQCl5rZbjPbEFyCYGYLzKw5Wl/0vq79iV+GCX/8ziHw2jCz18zsA/7DWGO/BBhD/1+DIsfiF8BCMxvlL517D/B4T6XGqaSajKEOQCQWMysCfgx80d+UD7QGmrT628L3mwl8BuhvQhsu/DjBY5XjvXaWAKfgLfV4BNgO3H6MfYH3dupVQAdwKfCgmZ3inNvejxgbo/Rbgrc8p2dpznjgD2a2wTn3pHPuT0BxjBiD8fUVvwwvXwd2AL/v2eCcmxGoj/V85wPdzrn9EepEjtebQDOwG3DAH4H/6KnUOJVUo5lwSVpmloO3Fvxh59zP/c3teElvj0J/W3C/OuBB4KPOuX7NhJvZ2sCHHmsiHCd4rAP+41udc83OuW3AD4ELBtAXzrm/Ouf2OecOOefuBp7DW0rTl/7EuNo5d8C/Dr/sibEffQWva8z4Zfjw3w15N/APzjkXpVms57sdSDezERHqRI7Xv+NNRhThTSQ0Ad+O0lbjVIY9JeGSlMwsA7gX78OUnw9UrcP74E6PU+n9tvoYvLcvv+ace6C/x3POneze/tDjtvDj+J/QnwCsc87t8+PqFfJA+ooSzpF+hh3e73je/lDrRuBQtBj76ove1/VY45ckZGbvw3tHaYVzLuKHiH19jf23iPEaFDkOM4DbnXNtzrlW4OfAuVHaapzKsKckXJLVT/HWNK8Km7G7C/ikmdX5Cffn/G09S1d+D9zpnPtJeIdmluHPrqcDGWaW46+PjeQpIN+/tVs23l1V/uKc2+rX3wH8i5kVmNlYvDtPPDyQvszsPWaW58f3Pry7mTzR5xWC/wHONLMV/ozPauB+59xBf33kfwFf8m/7NRm4HPhdlL6iXtd+XAtJEmaW6Y/xNCDTH+NpZrYc70Nv73DObemjm6eI/XzfBXzZH/tnAhcB9yXifCQ1RRunwF+AD5nZCD+pXkXvdd1BT6FxKsPdUH8yVEUlvADj8NYDHuDttxbbgbP9+uvw3qbcB3wDMH/7Sn+/4D7tgX6/4tcHy6oYcZwBvObH8TSBu4Hg3bHlp0AL3vra1T1xDKCvP/n9NAMvAEsi7H/U3VH87RfgraPcj3frxZJAXTHwa7xbF24BPhmoOzt4bWJd177iV0megvfHYfgYXwQ8CRwOe23cFthvLXB5P8drLt7dVdrxPgfxgaE+b5XhVWKM09F4kwdNwB7gN8DYwH4apyopVXqSFxERERERGSRajiIiIiIiMsiUhIuIiIiIDDIl4SIiIiIig0xJuIiIiIjIIFMSLiIiIiIyyJSEi4iIiIgMMiXhIiIiIiKDTEm4iIiIiMggUxIuIiIiIjLIlISLiIiIiAwyJeEiIiIiIoPs/wfZMpjavfrFsAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 800x250 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "waveform_SCEDC_files = os.listdir(\"waveforms_SCEDC\")\n",
    "\n",
    "for i in waveform_SCEDC_files:\n",
    "    if i != '.DS_Store':\n",
    "        wf = read(\"waveforms_SCEDC/\" + i)\n",
    "        wf.plot()\n",
    "        break"
   ]
  }
 ],
 "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.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 5
}