Other results ...
Procedure (This notebook to download)
In [1]:
# -*- coding: utf-8 -*-
"""
Created on Wed May 6 17:58:19 2020
@author: martina
"""
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import rc,rcParams
from urllib.request import urlopen
rc('font', weight='bold')
plt.rcParams['axes.labelweight'] = 'bold'
In [2]:
shot_no = 33105
shot = str(shot_no)
BPP = np.loadtxt(urlopen('http://golem.fjfi.cvut.cz/shots/' + (shot) +'/DASs/PetiSonda/3-BPP.csv'),delimiter = ',')
BPP[:,0] *= 1000
BPP[:,1] *= 100
TP_Lim = np.loadtxt(urlopen('http://golem.fjfi.cvut.cz/shots/' + (shot) +'/DASs/PetiSonda/1-Lim.csv'),delimiter = ',')
TP_Lim[:,0] *= 1000
TP_Lim[:,1] *= 100
TP_OLim = np.loadtxt(urlopen('http://golem.fjfi.cvut.cz/shots/' + (shot) +'/DASs/PetiSonda/2-OLim.csv'),delimiter = ',')
TP_OLim[:,0] *= 1000
TP_OLim[:,1] *= 100
LP = np.loadtxt(urlopen('http://golem.fjfi.cvut.cz/shots/' + (shot) +'/DASs/PetiSonda/4-LP.csv'),delimiter = ',')
LP[:,0] *= 1000
LP[:,1] *= 100
In [3]:
fig,ax = plt.subplots(1)
ax.axhline(y = 0, linestyle = '--', color = 'black')
ax.plot(BPP[:,0],BPP[:,1], label = 'BPP')
ax.plot(TP_Lim[:,0],TP_Lim[:,1], label = 'TP_Lim')
ax.plot(TP_OLim[:,0],TP_OLim[:,1], label = 'TP_OLim')
ax.plot(LP[:,0],LP[:,1], label = 'LP')
ax.set_xlabel('t [ms]')
ax.set_ylabel('U [V]')
ax.set_xlim(4,17)
ax.set_ylim(-40,20)
ax.set_title('#'+str(shot),fontweight = 'bold')
ax.legend()
plt.savefig('icon-fig')
In [4]:
Te = (BPP[:,1] - TP_OLim[:,1])/2.5
In [5]:
fig,ax = plt.subplots(1)
ax.plot(BPP[:,0],Te[:], label = r'$T_\mathrm{e}$')
ax.set_xlabel('t [ms]')
ax.set_ylabel(r'$T_\mathrm{e} \ \mathrm{[eV]}$')
ax.axhline(y = 0, linestyle = '--', color = 'black')
ax.set_xlim(4,17)
ax.set_ylim(-5,20)
ax.set_title('#'+str(shot),fontweight = 'bold')
ax.legend()
Out[5]:
In [6]:
fig,ax = plt.subplots(5,figsize=(8, 8))
ax[0].axhline(y = 0, linestyle = '--', color = 'black',linewidth = 1)
ax[0].plot(TP_Lim[:,0],TP_Lim[:,1], label = 'TP_Lim')
ax[0].set_ylabel('U [V]')
ax[0].set_xlim(4,17)
ax[0].label_outer()
ax[0].set_ylim(-40,20)
ax[0].set_title('#'+str(shot),fontweight = 'bold')
ax[0].legend()
ax[1].axhline(y = 0, linestyle = '--', color = 'black',linewidth = 1)
ax[1].plot(TP_OLim[:,0],TP_OLim[:,1], label = 'TP_OLim')
ax[1].set_xlabel('t [ms]')
ax[1].set_ylabel('U [V]')
ax[1].set_xlim(4,17)
ax[1].label_outer()
ax[1].set_ylim(-40,20)
ax[1].legend()
ax[2].axhline(y = 0, linestyle = '--', color = 'black',linewidth = 1)
ax[2].plot(BPP[:,0],BPP[:,1], label = 'BPP')
ax[2].set_xlabel('t [ms]')
ax[2].set_ylabel('U [V]')
ax[2].set_xlim(4,17)
ax[2].label_outer()
ax[2].set_ylim(-40,20)
ax[2].legend()
ax[3].axhline(y = 0, linestyle = '--', color = 'black',linewidth = 1)
ax[3].plot(LP[:,0],LP[:,1], label = 'LP')
ax[3].set_xlabel('t [ms]')
ax[3].set_ylabel('U [V]')
ax[3].set_xlim(4,17)
ax[3].label_outer()
ax[3].set_ylim(-40,20)
ax[3].legend()
ax[4].axhline(y=0,color = 'black',linestyle = '--',linewidth = 1)
ax[4].plot(BPP[:,0],Te[:], label = r'$T_\mathrm{e}$')
ax[4].set_xlabel('t [ms]')
ax[4].set_ylabel(r'$T_\mathrm{e} \ \mathrm{[eV]}$')
ax[4].set_xlim(4,17)
ax[4].set_ylim(-5,20)
ax[4].legend()
Out[6]:
