checkpoints were overwritten by zero functions.
Also PML boundaries and viscoelastic modelling now work with checkpointing.
I didnt't implement higher fd-orders in time because checkpoint files would become very large.