{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 24,
   "metadata": {
    "collapsed": false
   },
   "outputs": [],
   "source": [
    "import numpy as np\n",
    "from matplotlib.pyplot import *\n",
    "from pylab import *\n",
    "import matplotlib.image as mpimg\n",
    "from urllib import urlopen\n",
    "\n",
    "baseURL = \"http://golem.fjfi.cvut.cz/utils/data/\";ShotNo = 22471\n",
    "fig=figure(1);subplots_adjust(hspace=0.001)\n",
    "\n",
    "sbp1=subplot(511)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/loop_voltage')\n",
    "obj1=np.loadtxt(dataURL, delimiter='\\t')\n",
    "ylim(0,26);yticks(arange(0, 30, 5))\n",
    "title('#'+str(shotno));ylabel('$U_l$ [V]')\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='Loop voltage $U_l$' );legend(loc=0)\n",
    "\n",
    "sbp1=subplot(512, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/toroidal_field')\n",
    "obj1=np.loadtxt(dataURL, delimiter='\\t')\n",
    "yticks(arange(0, 0.5 , 0.1));ylim(0,0.35)\n",
    "ylabel('$B_t$ [T]')\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='Toroidal mag. field $B_t$');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(513, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/plasma_current')\n",
    "obj1=np.loadtxt(dataURL, delimiter='\\t')\n",
    "yticks(arange(0, 4.5, 1));ylim(0,4.5)\n",
    "ylabel('$I_p$ [kA]')\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1]/1000, 'k-',label='Plasma current $I_p$');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(514, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/photodiode_alpha')\n",
    "obj1=np.loadtxt(dataURL, delimiter='\\t')\n",
    "yticks(arange(0, 0.09 , 0.02));ylim(0,0.09)\n",
    "ylabel('Intensity [a.u.]')\n",
    "plt.plot(obj1[:,0]*1000, obj1[:,1], 'k-',label='$H_\\\\alpha$ radiation');legend(loc=0)\n",
    "\n",
    "sbp1=subplot(515, sharex=sbp1)\n",
    "dataURL = urlopen(baseURL+ str(ShotNo) + '/electron_density')\n",
    "obj2=np.loadtxt(dataURL, delimiter=' ')\n",
    "yticks(arange(0, 0.8 , 0.2));ylim(0,0.8)\n",
    "ylabel('$n_e$')\n",
    "xticks(arange(8, 30, 5));xlim(5,25)\n",
    "xlabel('Time [ms]')\n",
    "plt.plot(obj2[:,0]*1000, obj2[:,1]*1e-19, 'k-',label='electron density $n_e$');legend(loc=0)\n",
    "\n",
    "xticklabels = sbp1.get_xticklabels()\n",
    "setp(xticklabels, visible=False)\n",
    "\n",
    "savefig('basicgraph.pdf')\n",
    "savefig('basicgraph.jpg')\n",
    "\n",
    "show()\n"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {
    "collapsed": true
   },
   "outputs": [],
   "source": []
  }
 ],
 "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.6"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 0
}