Add ipynb to create four-digit year of entry and plot