checkpoint is somewhat improper because checkpointing functionality was
implemented in vorticity_equation instead of
direct_numerical_simulation, so that would be the next step.

preliminary test works with no obvious problems.