#%%
if 1: # plot mesh
fig = plt.figure(8); fig.clf()
tets = bba.get_convergences_max_tets()
varsz = bba.get_variable_vs(swpvar)
Y = {}
for key in tets.keys():
Y[varsz[key]] = tets[key]
y = pd.Series(Y ) #.values(), index = varsz.values())
y.plot(marker = '*', ms = 20)
ax7t = ax7.twinx()
ax7t.plot(y, marker = '*', ms = 10, c = 'g')
#%%
if 1:
fig = plt.figure(21); fig.clf()
tts = bba.get_convergences_Tets_vs_pass()
for key, x in tts.iteritems():
#np.log10(x).plot(label = varsz[key])
x.plot(label = varsz[key])
plt.legend(loc = 0)
#%%
if 0:
variation = '0'; pJ_method = 'J_surf_mag';
#pJ_mj_series = bbp.calc_Pjs_from_I_for_mode(variation, bbp.U_H,bbp.U_E, bbp.LJs, junc_rects, junc_lens, method = pJ_method) # to be implemented
res = bbp.calc_avg_current_J_surf_mag(variation,junc_rects[0], junc_lens[0])
if 0: # for debug
variation = '0'; junc_rect = 'juncV';
print_color(' Setup: ' + bbp.setup.name)
lv = bbp.get_lv(variation)
