In [9]:
import MDSplus as mds
%pylab inline
In [37]:
# mds.ipp.ac.cn corresponds to public IP 202.127.204.12
# Also working : 202.127.204.42, 202.127.205.8
cnx = mds.Connection('mds.ipp.ac.cn')
In [55]:
shot_nb = 55672
cnx.openTree('east', shot_nb)
In [32]:
Ip = cnx.get('\IPM')
t_Ip = cnx.get('dim_of(\IPM)')
In [33]:
plot(t_Ip, Ip, lw=2)
grid()
xlim(0,8)
Out[33]:
A more pretty (or Pythonesque) approach from Laurent :
In [53]:
def get_signal(signal_name, shot_nb, mdsplus_connection):
"""
Retrieves a signal and its associated time vector
for a given shot number on EAST
"""
mdsplus_connection.openTree('EAST', shot_nb)
return mdsplus_connection.get(signal_name), \
mdsplus_connection.get('dim_of('+signal_name+')')
In [59]:
Ip, t_Ip = get_signal('\IPM', shot_nb, cnx)
In [61]:
Ne, t_Ne = get_signal('TOP:\dfsdev', shot_nb, cnx)
sig_name | description |
---|---|
\IPM |
plasma current |
\dfsdev |
|
\plhi1 |
LHCD 1 Incident power |
\plhi2 |
LHCD 2 Incident power |
\plhr1 |
LHCD 1 Reflected power |
\plhr2 |
LHCD 2 Reflected power |
\nbi1lhv |
BOX1 - LEFT - NBI |
\nbi2lhv |
BOX2 |
\nbi1rhv |
BOX1 - RIGHT - NBI |
\nbi2rhv |
BOX2 |
\icrf9 |
antenna B - Line 1 |
\icrf10 |
|
\icrf11 |
antenna B - Line 2 |
\icrf12 |
|
\icrf13 |
antenna B - Line 3 |
\icrf14 |
|
\icrf15 |
antenna B - Line 4 |
\icrf16 |
|
\te0_hrs |
Te0 |
In [ ]: