Source code :: main

[Return]
[Download]#!/usr/bin/python2 import matplotlib matplotlib.rcParams['backend'] = 'Agg' matplotlib.rc('font', family='arial', size='10') matplotlib.rc('text', usetex=True) import sys sys.path.append('/home/michal/Desktop/Prace/Golem/web_framework') from numpy import * from pygolem.config import * from pygolem.modules import * from basic_diagn import * def plot_data(file_type): shot = 1000 plasma = 1 t0 = time.time() #GPfPlasma => home page data = [ get_data('Uloop_dp', 'Loop voltage', 'U [V]', ylim = [0,None]), get_data('Btoroidal_dp', 'Toroidal mag. field', 'B$_t$ [T]' , ylim = [0,None]), get_data('Iplasma_dp', 'Plasma current', 'I$_{p}$ [kA]' , ylim = [0,None], data_rescale = 1e-3) if plasma else \ get_data('Irogowski_dp', 'Plasma current', 'I$_{p}$+I$_{ch}$ [kA]', data_rescale = 1e-3, ylim = [0,None]) , [get_data('Photod_dp', 'Visible', 'Photod [a.u.]' , ylim = [0,None]), get_data('PhotodHalpha_dp', 'H$_\\alpha', 'Photod [a.u.]' , ylim = [0,None])] , ] paralel_multiplot(data, 'Golem shot No:' + str(shot), 'graphpres', (9,7), 100, 'vertical', file_type) # GPf data = [ get_data('Uloop_dp', 'Loop voltage', 'U [V]'), get_data('Btoroidal_dp', 'Toroidal mag. field', 'B$_t$ [T]'), get_data('Irogowski_dp', 'Total current', 'I$_{ch}$ [kA]' , data_rescale = 1e-3), [get_data('Photod_dp', 'Visible', 'Photod [a.u.]' ), get_data('PhotodHalpha_dp', 'H$_\\alpha', 'Photod [a.u.]' )] , ] paralel_multiplot(data, 'Golem shot No:' + str(shot), 'graphprint', (9,7), 100, 'vertical', file_type) data = [ get_data('Uloop_dp', 'Loop voltage', 'U [V]', xlim = [0,25]), get_data('Btoroidal_dp', 'Toroidal mag. field', 'B$_t$ [T]' , xlim = [0,25]), get_data('Iplasma_dp', 'Plasma current', 'I$_{p}$ [kA]' , xlim = [0,25], data_rescale = 1e-3) if plasma else \ get_data('Irogowski_dp', 'Plasma current', 'I$_{p}$+I$_{ch}$ [kA]', data_rescale = 1e-3) , [get_data('Photod_dp', 'Visible', 'Photod [a.u.]', xlim = [0,25]), get_data('PhotodHalpha_dp', 'H$_\\alpha', 'Photod [a.u.]' , xlim = [0,25])] , ] paralel_multiplot(data, 'Golem shot No:' + str(shot) + ' final data ', 'graphpresfull', (9,7), 100, 'vertical', file_type) ##GPitegrated data = [ [get_data('Uloop_dp', 'Loop voltage', 'U [V]', xlim = [0,25])], get_data('Btoroidal_dp', 'Toroidal mag. field', 'B$_t$ [T]', xlim = [0,25]), [get_data('Irogowski_dp', 'Total current', 'I$_p$+I$_{ch}$ [kA]' , xlim = [0,25], data_rescale = 1e-3)] + \ [get_data('Iplasma_dp', 'Plasma current', 'I$_{p}$ [kA]' , xlim = [0,25], data_rescale = 1e-3)] if plasma else [], [get_data('Photod_dp', 'Visible', 'Photod [a.u.]', xlim = [0,25]), get_data('PhotodHalpha_dp', 'H$_\\alpha', 'Photod [a.u.]' , xlim = [0,25])] , ] paralel_multiplot(data, 'Golem shot No:' + str(shot) + " - integrated data", 'graphpresi', (9,7), 100, 'vertical', file_type) ###raw data = [ get_data('Uloop_dp', 'Loop voltage', 'U [V]', xlim = [0,25] ), get_data('dBdt_toroidal_dp', 'Derivative of mag. field', 'dB$_t$/dt [DAS V]', xlim = [0,25] ), get_data('dIdt_rogowski_dp', 'Derivative of current', 'dI$_{p+ch}$/dt [DAS V]', xlim = [0,25] ), [get_data('Photod_dp', 'Visible', 'Photod [a.u.]' , xlim = [0,25] ), get_data('PhotodHalpha_dp', 'H$_\\alpha', '[DAS V]', xlim = [0,25] )] , ] paralel_multiplot(data, 'Golem shot No:' + str(shot) + " - raw data", 'graphpresb', (9,7), 100, 'vertical', file_type) # icon data = [ get_data('Uloop_dp', '', '', xlabel = ""), [get_data('Irogowski_dp', '', '', xlabel = "" )] +\ [get_data('Iplasma_dp','', '', xlabel = "" )] if plasma else [] ] paralel_multiplot(data, '', 'icon', (4,3), 40) # GPic data = [ get_data('Uloop_dp', '', 'U [V]', ylim = [0,None], xlim = [0,25], xlabel=""), [get_data('Irogowski_dp', '', 'I$_p$+I$_{ch}$ [kA]' , ylim = [0,None], xlim = [0,25], xlabel = "", data_rescale = 1e-3)] + \ [get_data('Iplasma_dp', '', 'I$_{p}$ [kA]' , ylim = [0,None], xlim = [0,25], xlabel = "", data_rescale = 1e-3)] if plasma else [] ] paralel_multiplot(data, "", 'graphic', (9,1.5), 100, 'horizontal', file_type) print ' Time: %g' % ( time.time() - t0) def prepare_data(): t = time.time() Aktual_PfeifferMerkaVakua = loadconst("Aktual_PfeifferMerkaVakua") [Btor, BtMax, BtMean] = getBtoroidal() [Uloop, UloopMax, UloopMean] = getUloop() [Irog, dIdt_rogMax, IrogMax] = getIrogowski() [Ipla, Ich] = getIplasma(Uloop,Irog) Plasma, PlasmaStart, PlasmaEnd, PlasmaTimeLength = PlasmaDetect(Ipla) Photod = getPhotod(PlasmaStart) PhotodHalpha = getPhotodHalpha(PlasmaStart) MeanBt = getMeanBt(Btor,PlasmaStart, PlasmaEnd ) TotalCharge = getTotalCharge(Ipla, PlasmaStart,PlasmaEnd ) MeanUloop = getMeanUloop(Uloop, PlasmaStart,PlasmaEnd ) MeanIpla = getMeanCurrent(Ipla, PlasmaStart,PlasmaEnd ) OhmicHeatingPower = getOhmicHeatingPower(MeanUloop,MeanIpla) Qedge = getQedge(MeanBt,MeanIpla) ElectronTemperature = getElectronTemperature(MeanUloop, MeanIpla) Umax, Btime, Bbreak, Bipla = getBreakDownVoltage(Uloop, Btor, Ipla, PlasmaStart, PlasmaEnd) StateEqElectronDensity = getStateEqElectronDensity(Aktual_PfeifferMerkaVakua) ChamberResistance = getChamberResistance(Plasma) Plasma = Failures(Plasma, UloopMax, dIdt_rogMax, MeanUloop, BtMax, MeanBt, PlasmaStart, PlasmaEnd ) GreenwaldDensity = getGreenwaldDensity(Ipla) OhmicHeatingPower = getOhmicHeatingPowerTime(Ipla, Uloop) Qedge = getQedgeTime(Btor,Ipla) MagneticFlux = getMagneticFlux(Uloop) print "time" , time.time() - t def main(): prepare_data() plot_data('png') plot_data('svg') if __name__ == "__main__": main()[Return]

Navigation