A-614 opened this issue 3 years ago · 0 comments
Plotting zundel.out:
axs[0].plot(time, (KE - KE[0]), linewidth = 1.5, label = r'$\Delta PE$') axs[0].plot(time, (PE - PE[0]), linewidth = 1.5, label = r'$\Delta KE$')