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):