Rewrite plot functionality
The way to plot from the FORTRAN code through an input file routed to gnuplot is cumbersome and it is difficult to make changes without breaking anything. This functionality should be handled differently, maybe in the python part using matplotlib.