%matplotlib inline #If interactive mode

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt


shot_no = 33993

data_URL = "http://golem.fjfi.cvut.cz/shots/{shot_no}/DASs/StandardDAS/{identifier}.csv"  
ds = np.DataSource(destpath='/tmp')

def open_remote(shot_no, identifier, url_template):
    return ds.open(url_template.format(shot_no=shot_no, identifier=identifier))

def read_signal(shot_no, identifier): 
    file = open_remote(shot_no, identifier, data_URL)
    return pd.read_csv(file, names=['Time',identifier],
                     index_col='Time', squeeze=True)  # squeeze makes simple 1-column signals a Series

loop_voltage = read_signal(shot_no, 'LoopVoltageCoil_raw')
ax = loop_voltage.plot(grid=True)
ax.set(xlabel="Time [s]", ylabel="$U_l$ [V]", title="Loop voltage $U_l$ #{}".format(shot_no));

plt.show()
plt.savefig('graph.jpg')


#Run it: save it as script.py and run "python script.py" or execute in a ceel in a Jupyter Notebook

ShotNo = '33993'
DAS='DASs/StandardDAS/'
identifier = 'LoopVoltageCoil_raw.csv' ;
baseURL='http://golem.fjfi.cvut.cz/shots/'
DataURL= baseURL.ShotNo.'/'.DAS.identifier
set datafile separator ',';
set title "Uloop@ShotNo";
! wget -q @DataURL ;
plot identifier u 1:2 w l t 'Uloop'


#Run it interactively or  save it as script.gp and run "gnuplot --persist script.gp"