eq.run(update=False)
sf.contour()
pf.plot(subcoil=False, current=True, label=True)
pf.plot(subcoil=True)
inv.ff.plot()
mc = main_chamber('dtt', date='2017_03_08')
mc.load_data(plot=True) # load from file
mc.shp.plot_bounds()
rb = RB(sf, Setup(name))
rb.generate(mc, debug=False)
#profile = Profile(config['TF'],family='S',part='TF',nTF=nTF,obj='L')
shp = Shape(tf.profile, eqconf=config['eq_base'], ny=3)
shp.add_vessel(rb.segment['vessel_outer'])
#shp.minimise(ripple=ripple,verbose=True)
shp.tf.fill()
sf.eqwrite(sf, config='SXex')
'''
inv.Lnorm = Lnorm
sw = SWING(inv,sf)
sw.flat_top()
sw.output()
for ends,name in zip([0,-1],['SOF','EOF']):
inv.solve_slsqp(inv.swing['flux'][ends])
eq.run(update=False)
sf.eqwrite(pf,config='SXex_{}'.format(name))
'''
import matplotlib.animation as manimation
import numpy as np
nTF = 18
family = 'S'
config = {'TF': 'dtt', 'eq': 'SN'}
config['TF'] = '{}{}{:d}'.format(config['eq'], config['TF'], nTF)
demo = DEMO()
demo.fill_part('Blanket')
demo.fill_part('Vessel')
#demo.fill_part('TF_Coil')
profile = Profile(config['TF'], family=family, part='tmp', load=False)
shp = Shape(profile, nTF=nTF, obj='L', eqconf=config['eq'], ny=3)
shp.add_vessel(demo.parts['Vessel']['out'])
shp.minimise(ripple=False, verbose=True)
#shp.update()
shp.loop.set_input(x=shp.xo[0])
shp.tf.fill()
shp.plot_bounds()
'''
filename = 'movie'
moviename = '../Movies/{}'.format(filename)
moviename += '.mp4'
print(moviename)
with writer.saving(fig,moviename,100):