fix['value'] = np.append(fix['value'],0)
fix['BC'] = np.append(fix['BC'],'Bz')
'''
del sf.coil['Coil9']
#del sf.coil['Coil10']
#del sf.coil['Coil11']
#del sf.coil['Coil12']
#del sf.coil['Coil13']
eq.coils(delta=0.25) # multi-filiment coils
conf.TF(sf)
rb.TFopp(False,objF=conf.TFopp) # L==length, V==volume
rb.TFfill()
eq.coil_currents(fix)
sf.plot_coils(Color,coils=sf.coil,label=False,plasma=False,current=True)
sf.plot_coils(Color,coils=eq.coil,label=False,plasma=False)
eq.update_psi() # levels=sf.cs.levels
sf.contour(Nstd=1.5)
#sf.sol(plot=True)
eq = EQ([3,14],[-10,10],8e3,sf)
#eq = EQ([2,16],[-10,10],2e3,sf)
sf.plot_coils(Color)
conf.TF(sf)
rb = RB(conf, sf, Np=200)
rb.divertor_outline(True)
rb.fill(dt=conf.tfw, alpha=0.7, color=next(Color))
rb.fill(dt=conf.BB, alpha=0.7, color=next(Color))
rb.fill(dt=conf.tBBsupport, alpha=0.7, color=next(Color))
rb.BBsheild_fill(dt=conf.sheild,
ref_o=2 / 8 * np.pi,
alpha=0.7,
color=next(Color))
rb.fill(dt=conf.VV, alpha=0.7, color=next(Color), loop=True)
rb.set_TFbound() # TF boundary conditions
#rb.TFbound['ro_min'] -= 0.25
#rb.plot_TFbounds()
rb.TFopp(False, objF='L') # L==length, V==volume
rb.TFfill()
pl.text(16, 9.5, config, fontsize=36, ha='right', va='top')
'''
Color = cycle(sns.color_palette('Set2'))
for leg in conf.targets.keys():
R,Z = rb.sol.legs(leg)
pl.plot(R,Z,color=next(Color))
pl.tight_layout()
col_labels=['volume','length']
row_labels=['TF','LCFS','ratio']
table_vals=[[r'{:1.0f}m$^3$'.format(rb.TFvol),r'{:1.1f}m'.format(rb.TFlength)],
[r'{:1.0f}m$^3$'.format(rb.Pvol),r'{:1.1f}m'.format(rb.Plength)],
['{:1.2f}'.format(rb.Rvol),r'{:1.2f}'.format(rb.Rlength)]]
P6support['bottom'] = {}
P5support = {}
P5support['top'] = {}
P5support['bottom'] = {}
if config is 'DD3':
rb.VVfill(connect=[0.36, 0.802],
r=[13.56, 13.5, 10, 6.8, 5.2],
z=[-1.5, -5, -10, -9.5, -3.5],
w=[0.9, 2, 2, 3.5, 2],
dt=VV,
ref_o=3 / 8 * np.pi,
alpha=0.3,
color='b',
label='vacuum vessel')
rb.TFopp(TFopp)
rb.rzGet() # VV ports
rb.fill(trim=[0.795, 0.795 + 0.04],
dt=VV,
ref_o=3 / 8 * np.pi,
alpha=1,
color='r')
rb.rzGet()
rb.fill(trim=[0.745, 0.745 + 0.04],
dt=VV,
ref_o=3 / 8 * np.pi,
alpha=1,
color='r',
label='ports')
r95,z95 = sf.get_boundary(alpha=0.95)
ru = r95[np.argmax(z95)] # triangularity
rl = r95[np.argmin(z95)]
del_u = (R-ru)/a
del_l = (R-rl)/a
kappa = (np.max(z95)-np.min(z95))/(2*a)
print('R',R,'a',a,'AR',AR,'del_u',del_u,'del_l',del_l,'kappa',kappa)
Pvol = rb.loop_vol(rb.Rp,rb.Zp,plot=False)
print(Pvol)
'''
rb.TFopp(False,config='SXex') # load TF outline
rb.TFfill() # construct TF geometory
'''
conf.TFopp = 'L'
rb.set_TFbound() # TF boundary conditions
#rb.TFbound['ro_min'] -= 0.5
rb.plot_TFbounds()
rb.TFopp(True,objF=conf.TFopp) # L==length, V==volume
rb.TFfill() # construct TF geometory
'''
fig.tight_layout(rect=[-0.3,-0.54,1.1,1.5])
pl.savefig('../Figs/SXcoils_const.png',dpi=300)