print(f"k24 = {k24}, Lambda = {Lambda}, q = {q}")
psistuff = PsiData(scan=scan,
loadsuf=loadsuf,
savesuf=savesuf,
name=f"psivsr",
sfile_format="pdf")
rs = psistuff.r() / q
psis = psistuff.psi()
ax1.plot(rs, psis, markertypes[i], label=f'{type} cornea', lw=2)
print("Radius = ", rs[-1], " nm.")
print("surface twist = ", psis[-1], " rad.")
ax1.set_xlabel(r'$\tilde{r}$' + ' (' + r'$\si{\nano\meter}$' + ')',
fontsize=10)
ax1.set_ylabel(r'$\psi(\tilde{r})$' + ' (' + r'$\si{\radian}$' + ')',
fontsize=10)
ax1.set_xlim(left=0)
ax1.set_ylim(bottom=0)
ax1.legend(frameon=False, fontsize=10)
fig.subplots_adjust(left=0.2, bottom=0.2)
psistuff.name = f"twocorneas"
fig.savefig(psistuff.psivsr_sname())
for i, type in enumerate(types):
psistuff = PsiData(scan=scan,
loadsuf=loadsuf,
savesuf=savesuf,
name=f"psivsr_{type}",
sfile_format="pdf")
rs = psistuff.r() / q
psis = psistuff.psi()
ax1.plot(rs, psis, markertypes[i], label=f'{type} tendon', lw=2)
print("Radius = ", rs[-1], " nm.")
print("surface twist = ", psis[-1], " rad.")
ax1.set_xlabel(r'$\tilde{r}$' + ' (' + r'$\si{\nano\meter}$' + ')',
fontsize=10)
ax1.set_ylabel(r'$\psi(\tilde{r})$' + ' (' + r'$\si{\radian}$' + ')',
fontsize=10)
ax1.set_xlim(left=0)
ax1.set_ylim(bottom=0)
ax1.legend(frameon=False, fontsize=10)
fig.subplots_adjust(left=0.2, bottom=0.2)
psistuff.name = f"coexistingtendons"
fig.savefig(psistuff.psivsr_sname())