NK_4.ipynb

{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": [
    "import gps_particle_data\n",
    "import wget\n",
    "import os\n",
    "import numpy as np \n",
    "import time\n",
    "import json\n",
    "import sys\n",
    "import math\n",
    "from datetime import datetime, timedelta, date\n",
    "from itertools import compress\n",
    "from random import randint\n",
    "import urllib2\n",
    "import aacgmv2\n",
    "import scipy\n",
    "import matplotlib\n",
    "from matplotlib import gridspec\n",
    "import numpy as np\n",
    "import matplotlib.pyplot as plt\n",
    "%matplotlib inline"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "====================================\n",
      "Loading data for satellite  53\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_061029_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  54\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_061029_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  56\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_061029_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  59\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_061029_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  60\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_061029_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  61\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_061008_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_061015_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_061022_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_061029_v1.03.ascii\n"
     ]
    }
   ],
   "source": [
    "#Script simple plots for NK nuclear tests on the following dates:\n",
    "#9 October 2006, 25 May 2009, 12 February 2013, 6 January 2016, 9 September 2016\n",
    "\n",
    "start_date = datetime(2006,10,7,0,0,0);\n",
    "end_date = datetime(2006,10,30,0,0,0);\n",
    "ms = gps_particle_data.meta_search();\n",
    "ms.load_data(start_date,end_date);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "\"for ns in output_data.keys():\\n    if len(output_data[ns]) !=0 and ns != 53:\\n        for i in range(11):\\n            cor = np.correlate(a0, np.asarray(output_data[ns][signal])[:,i], 'full')\\n            c_dif = np.argmax(cor)-np.argmax(ctrl)\\n            #cor_dict.update({ns:c_dif})\\n            print ns,'channel:',i,'--->', c_dif\""
      ]
     },
     "execution_count": 3,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "#We could look at not normalised cross correlation by following the shift of the peak.\n",
    "\n",
    "output_data = ms.get_all_data_by_satellite()\n",
    "signal = 'rate_electron_measured'\n",
    "a0 = np.asarray(output_data[53][signal])[:,2]\n",
    "#cor_dict ={}\n",
    "ctrl = np.correlate(a0,a0,'full')\n",
    "'''for ns in output_data.keys():\n",
    "    if len(output_data[ns]) !=0 and ns != 53:\n",
    "        for i in range(11):\n",
    "            cor = np.correlate(a0, np.asarray(output_data[ns][signal])[:,i], 'full')\n",
    "            c_dif = np.argmax(cor)-np.argmax(ctrl)\n",
    "            #cor_dict.update({ns:c_dif})\n",
    "            print ns,'channel:',i,'--->', c_dif'''\n",
    "        "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 6,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SATELLITE 54 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.13602126699\n",
      "channel: 1 X-cor  --->    0.117981837073\n",
      "channel: 2 X-cor  --->    0.0911004588186\n",
      "channel: 3 X-cor  --->    0.0655564934594\n",
      "channel: 4 X-cor  --->    0.073383917756\n",
      "channel: 5 X-cor  --->    0.0331515803175\n",
      "channel: 6 X-cor  --->    0.0298903936414\n",
      "channel: 7 X-cor  --->    0.0425039337645\n",
      "channel: 8 X-cor  --->    0.128672071677\n",
      "channel: 9 X-cor  --->    0.367403219871\n",
      "channel: 10 X-cor  --->    0.798417339149\n",
      "\n",
      "\n",
      "SATELLITE 56 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.7259604756\n",
      "channel: 1 X-cor  --->    0.704419554565\n",
      "channel: 2 X-cor  --->    0.665912151527\n",
      "channel: 3 X-cor  --->    0.658587360571\n",
      "channel: 4 X-cor  --->    0.686517720377\n",
      "channel: 5 X-cor  --->    0.661703205998\n",
      "channel: 6 X-cor  --->    0.662673335235\n",
      "channel: 7 X-cor  --->    0.604801346331\n",
      "channel: 8 X-cor  --->    0.530716454962\n",
      "channel: 9 X-cor  --->    0.6741933951\n",
      "channel: 10 X-cor  --->    0.92270558598\n",
      "\n",
      "\n",
      "SATELLITE 59 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.128083688345\n",
      "channel: 1 X-cor  --->    0.113253440273\n",
      "channel: 2 X-cor  --->    0.0944193886717\n",
      "channel: 3 X-cor  --->    0.0745819038011\n",
      "channel: 4 X-cor  --->    0.0836780250456\n",
      "channel: 5 X-cor  --->    0.0425661291043\n",
      "channel: 6 X-cor  --->    0.0377796392822\n",
      "channel: 7 X-cor  --->    0.0454269038983\n",
      "channel: 8 X-cor  --->    0.127433792955\n",
      "channel: 9 X-cor  --->    0.359896911394\n",
      "channel: 10 X-cor  --->    0.8006011118\n",
      "\n",
      "\n",
      "SATELLITE 60 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.23280812681\n",
      "channel: 1 X-cor  --->    0.220812810444\n",
      "channel: 2 X-cor  --->    0.218454411861\n",
      "channel: 3 X-cor  --->    0.201117551785\n",
      "channel: 4 X-cor  --->    0.187169741869\n",
      "channel: 5 X-cor  --->    0.141339519246\n",
      "channel: 6 X-cor  --->    0.118888033349\n",
      "channel: 7 X-cor  --->    0.109221912817\n",
      "channel: 8 X-cor  --->    0.214439182971\n",
      "channel: 9 X-cor  --->    0.475902567441\n",
      "channel: 10 X-cor  --->    0.863289011235\n",
      "\n",
      "\n",
      "SATELLITE 61 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.152921264895\n",
      "channel: 1 X-cor  --->    0.142494643973\n",
      "channel: 2 X-cor  --->    0.130535534026\n",
      "channel: 3 X-cor  --->    0.114356286692\n",
      "channel: 4 X-cor  --->    0.118340993333\n",
      "channel: 5 X-cor  --->    0.0723524117532\n",
      "channel: 6 X-cor  --->    0.0594830766859\n",
      "channel: 7 X-cor  --->    0.0613802274961\n",
      "channel: 8 X-cor  --->    0.163216518828\n",
      "channel: 9 X-cor  --->    0.437114223341\n",
      "channel: 10 X-cor  --->    0.854712373058\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "#Normalised cross correlation values as compared to satellite 53.\n",
    "\n",
    "for ns in output_data.keys():\n",
    "    if len(output_data[ns]) !=0 and ns != 53:\n",
    "        print 'SATELLITE', ns, '\\n', '===============\\n'\n",
    "        for i in range(11):\n",
    "            a0 = np.asarray(output_data[53][signal])[:,i]\n",
    "            ar2 = np.asarray(output_data[ns][signal])[:,i]\n",
    "            cor = np.correlate(a0, ar2)/(math.sqrt(np.sum(a0**2)*np.sum(ar2**2)))\n",
    "            print 'channel:', i, 'X-cor  --->   ', cor[0]\n",
    "        print'\\n'"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 7,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "\n",
      "====================================\n",
      "Loading data for satellite  53\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns53/ns53_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  54\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns54/ns54_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  55\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns55/ns55_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns55/ns55_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns55/ns55_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns55/ns55_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  56\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns56/ns56_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  57\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns57/ns57_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns57/ns57_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns57/ns57_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns57/ns57_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  58\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns58/ns58_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns58/ns58_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns58/ns58_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns58/ns58_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  59\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns59/ns59_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  60\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns60/ns60_090614_v1.03.ascii\n",
      "\n",
      "====================================\n",
      "Loading data for satellite  61\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_090524_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_090531_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_090607_v1.03.ascii\n",
      "File : https://www.ngdc.noaa.gov/stp/space-weather/satellite-data/satellite-systems/gps/data/ns61/ns61_090614_v1.03.ascii\n"
     ]
    }
   ],
   "source": [
    "#25 May 2009\n",
    "start_date = datetime(2009,5,23,0,0,0);\n",
    "end_date = datetime(2009,6,14,0,0,0);\n",
    "ms = gps_particle_data.meta_search();\n",
    "ms.load_data(start_date,end_date);"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 8,
   "metadata": {},
   "outputs": [
    {
     "name": "stdout",
     "output_type": "stream",
     "text": [
      "SATELLITE 54 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.13602126699\n",
      "channel: 1 X-cor  --->    0.117981837073\n",
      "channel: 2 X-cor  --->    0.0911004588186\n",
      "channel: 3 X-cor  --->    0.0655564934594\n",
      "channel: 4 X-cor  --->    0.073383917756\n",
      "channel: 5 X-cor  --->    0.0331515803175\n",
      "channel: 6 X-cor  --->    0.0298903936414\n",
      "channel: 7 X-cor  --->    0.0425039337645\n",
      "channel: 8 X-cor  --->    0.128672071677\n",
      "channel: 9 X-cor  --->    0.367403219871\n",
      "channel: 10 X-cor  --->    0.798417339149\n",
      "\n",
      "\n",
      "SATELLITE 56 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.7259604756\n",
      "channel: 1 X-cor  --->    0.704419554565\n",
      "channel: 2 X-cor  --->    0.665912151527\n",
      "channel: 3 X-cor  --->    0.658587360571\n",
      "channel: 4 X-cor  --->    0.686517720377\n",
      "channel: 5 X-cor  --->    0.661703205998\n",
      "channel: 6 X-cor  --->    0.662673335235\n",
      "channel: 7 X-cor  --->    0.604801346331\n",
      "channel: 8 X-cor  --->    0.530716454962\n",
      "channel: 9 X-cor  --->    0.6741933951\n",
      "channel: 10 X-cor  --->    0.92270558598\n",
      "\n",
      "\n",
      "SATELLITE 59 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.128083688345\n",
      "channel: 1 X-cor  --->    0.113253440273\n",
      "channel: 2 X-cor  --->    0.0944193886717\n",
      "channel: 3 X-cor  --->    0.0745819038011\n",
      "channel: 4 X-cor  --->    0.0836780250456\n",
      "channel: 5 X-cor  --->    0.0425661291043\n",
      "channel: 6 X-cor  --->    0.0377796392822\n",
      "channel: 7 X-cor  --->    0.0454269038983\n",
      "channel: 8 X-cor  --->    0.127433792955\n",
      "channel: 9 X-cor  --->    0.359896911394\n",
      "channel: 10 X-cor  --->    0.8006011118\n",
      "\n",
      "\n",
      "SATELLITE 60 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.23280812681\n",
      "channel: 1 X-cor  --->    0.220812810444\n",
      "channel: 2 X-cor  --->    0.218454411861\n",
      "channel: 3 X-cor  --->    0.201117551785\n",
      "channel: 4 X-cor  --->    0.187169741869\n",
      "channel: 5 X-cor  --->    0.141339519246\n",
      "channel: 6 X-cor  --->    0.118888033349\n",
      "channel: 7 X-cor  --->    0.109221912817\n",
      "channel: 8 X-cor  --->    0.214439182971\n",
      "channel: 9 X-cor  --->    0.475902567441\n",
      "channel: 10 X-cor  --->    0.863289011235\n",
      "\n",
      "\n",
      "SATELLITE 61 \n",
      "===============\n",
      "\n",
      "channel: 0 X-cor  --->    0.152921264895\n",
      "channel: 1 X-cor  --->    0.142494643973\n",
      "channel: 2 X-cor  --->    0.130535534026\n",
      "channel: 3 X-cor  --->    0.114356286692\n",
      "channel: 4 X-cor  --->    0.118340993333\n",
      "channel: 5 X-cor  --->    0.0723524117532\n",
      "channel: 6 X-cor  --->    0.0594830766859\n",
      "channel: 7 X-cor  --->    0.0613802274961\n",
      "channel: 8 X-cor  --->    0.163216518828\n",
      "channel: 9 X-cor  --->    0.437114223341\n",
      "channel: 10 X-cor  --->    0.854712373058\n",
      "\n",
      "\n"
     ]
    }
   ],
   "source": [
    "for ns in output_data.keys():\n",
    "    if len(output_data[ns]) !=0 and ns != 53:\n",
    "        print 'SATELLITE', ns, '\\n', '===============\\n'\n",
    "        for i in range(11):\n",
    "            a0 = np.asarray(output_data[53][signal])[:,i]\n",
    "            ar2 = np.asarray(output_data[ns][signal])[:,i]\n",
    "            cor = np.correlate(a0, ar2)/(math.sqrt(np.sum(a0**2)*np.sum(ar2**2)))\n",
    "            print 'channel:', i, 'X-cor  --->   ', cor[0]\n",
    "        print'\\n'"
   ]
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 2",
   "language": "python",
   "name": "python2"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 2
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython2",
   "version": "2.7.13"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 2
}