'lines.linewidth': 0.75 } sns.set(context='paper', style='white', font='sans-serif', palette='Set2', font_scale=7 / 8, rc=rc) Color = cycle(sns.color_palette('Set2')) pl.figure() pl.axis('equal') pl.axis('off') conf = Config('SN') sf = SF(conf) conf.TF(sf) rb = RB(conf, sf, Np=500) a = (sf.LFPr - sf.HFPr) / 2 R = (sf.LFPr + sf.HFPr) / 2 AR = R / a 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) pl.plot(sf.Mpoint[0], sf.Mpoint[1], 'o') pl.plot(sf.LFPr, sf.LFPz, 'o')
'savefig.jpeg_quality':100,'savefig.pad_inches':0.1, 'lines.linewidth':0.75} sns.set(context='paper',style='white',font='sans-serif',palette='Set2', font_scale=7/8,rc=rc) Color = cycle(sns.color_palette('Set2')) pl.figure() pl.axis('equal') pl.axis('off') conf = Config('SN') sf = SF(conf) conf = Config('X') sf.conf = conf conf.TF(sf) #rb = RB(conf,sf,Np=200) #rb.divertor_outline(False,plot=False,debug=False) eq = EQ(sf,dCoil=0.5,limit=[5,12,-8,5.5],n=1e4) eq.set_sf_psi() # set psi eq.gen() inv = INV(sf,eq,configTF='X',config='X') Lpf = inv.grid_PF(nPF=5) Lcs = inv.grid_CS(nCS=5) Lo = np.append(Lpf,Lcs) inv.eq.coils() # re-grid inv.update_coils() inv.fix_boundary_psi(N=31,alpha=1-1e-4,factor=1) # add boundary points #inv.fix_boundary_feild(N=31,alpha=1-1e-4,factor=1) # add boundary points