The file user.f.timedep contains an example of solving a (parabolic) time dependent problem. The Crank Nicolson method is used for discretizing the first order time derivative. Each iteration advances the solution one time step. This also illustrates the use of subroutines save and ssolut to save the state and evaluate the saved solution. Additionally, it shows how to call gnuplot directly to create graphics.