# coding: utf-8

# In[1]:

import matplotlib.pyplot as plt
import numpy as np
import matplotlib.image as mpimg
import matplotlib.ticker as tic
from pylab import *

#%matplotlib inline
data=np.loadtxt("FinalniData.shots", delimiter=' ')


fig=plt.figure(1)
plt.subplots_adjust(hspace=0.00)
plt.figure(figsize=(13, 10), dpi=100)
plt.rc("font", size=15)

plt.subplot(511)
plt.ylabel('Delka\nvyboje\n$\\tau_{pl}$ [ms]');plt.ylim(0,25)
plt.plot(data[:,1], data[:,2], 'o')
plt.plot(data[:,1], data[:,2]*0+14,'r') #pardon
#plt.axvline(x=22000,c='r')
plt.tick_params(axis='x',labelbottom='off') 

plt.subplot(512)
plt.ylabel('Elektronova\nteplota\n$T_{e}$ [eV]');plt.ylim(0,80)
plt.plot(data[:,1], data[:,3], 'bo')
plt.plot(data[:,1], data[:,3]*0+40,'r') #pardon
plt.yticks(np.arange(10, 80, 20))
plt.tick_params(axis='x',labelbottom='off') 

plt.subplot(513)
plt.ylabel('Elektronova\nhustota\n$n_{e}$ [$\mathrm{10^{18}/m^3}$]');plt.ylim(0,13)
plt.plot(data[:,1], data[:,4], 'bo')
plt.plot(data[:,1], data[:,4]*0+10,'r') #pardon
plt.yticks(np.arange(0, 14, 3))
plt.tick_params(axis='x',labelbottom='off') 

plt.subplot(514)
plt.ylabel('Proud\nplazmatem\n$I_{pl}$ [kA]');plt.ylim(0,6)
plt.plot(data[:,1], data[:,5], 'bo')
plt.plot(data[:,1], data[:,5]*0+4,'r') #pardon
plt.yticks(np.arange(0, 7, 2))
plt.tick_params(axis='x',labelbottom='off') 

plt.subplot(515)
plt.ylabel('Davka\n$D$ [$10^{-2} \mathrm{\mu Sv}$]');plt.ylim(-0.5,12)
plt.yticks(np.arange(0, 20, 4))
plt.plot(data[:,1], data[:,6], 'bo')


plt.xlabel('Cislo vyboje');
plt.xticks(rotation=90)

plt.savefig('DlouhodobaStabilita@tokamakGOLEM.pdf')
plt.savefig('DlouhodobaStabilita@tokamakGOLEM.jpg')
plt.show()


# In[ ]:




# In[ ]: